trimmomatic 去接头,处理reads

<span style="font-size:10px;"></span><pre name="code" class="plain"><span style="font-size:10px;">nohup java -jar</span>
/home/lixiangyong/software/Trimmomatic/trimmomatic-0.36.jar PE -threads 2 -phred33
 
<span style="font-size:10px;">./Sample_hua_20160524_CGTACG_L001_R1.fastq ./Sample_hua_20160524_CGTACG_L001_R2.fastq </span>
<span style="font-size:10px;">Sample_hua_20160524_CGTACG_L001_R1_paired.fq <span style="font-family: Arial, Helvetica, sans-serif;">Sample_hua_20160524_CGTACG_L001_R1_unpaired.fq </span></span>
<span style="font-family: Arial, Helvetica, sans-serif;"><span style="font-size:10px;">Sample_hua_20160524_CGTACG_L001_R2_paired.fq Sample_hua_20160524_CGTACG_L001_R2_unpaired_fq  </span></span>
<span style="font-size:10px;">ILLUMINACLIP:/data/public/wanglei/Lj/hua_from_wanglei_analysis/illuminaTruSeq3-PE-2.fa:2:30:10:1 LEADING:5 TRAILING:5 SLIDINGWINDOW:4:15 MINLEN:36 &</span>



nohup java -jar trimmomatic-0.36.jar   PE -threads 2 -phred33 
input_forward.fq input_reverse.fq output_forward_paired.fq output_forward_unpaired.fq output_reverse_paired.fq output_reverse_unpaired.fq
<span style="color: rgb(55, 96, 146); font-family: simsun; font-size: 14px; line-height: 21px; background-color: rgb(240, 231, 180);"></span>ILLUMINACLIP:1.adapter.list:2:30:10 LEADING:3 TRAILING:3 SLIDINGWINDOW:4:15 MINLEN:36<span style="color: rgb(55, 96, 146); font-family: simsun; font-size: 14px; line-height: 21px; background-color: rgb(240, 231, 180);"></span>

同一个指令下的不同参数可以用':'来界定

phred33    设置碱基的质量格式,如果不设置,默认的是-phred64

trimlog file就是产生日志,包括如下部分内容:read的名字,留下来的序列的长度,第一个碱基的起始位置,从开始trimmed的长度,最后的一个碱基位于初始read的位置,最后trimmed的数量,其他步骤产生的日志

LEADING:3   切除首端碱基质量小于3的碱基或者N

TRAILING:3   切除末端碱基质量小于3的碱基或者

ILLUMINACLIP: 1.adapter.lis:2:30:10 1.adapter.listadapter文件,允许的最大mismatch 数,palindrome模式下匹配碱基数阈值:simple模式下的匹配碱基数阈值

SLIDINGWINDOW:4:15  Windowssize4个碱基,其平均碱基质量小于15,则切除

MINLEN:36  最低reads长度为36

CROP:  保留的reads长度

HEADCROP: reads的首端切除指定的长度


原文来自于铁汉 http://blog.sina.com.cn/s/blog_670445240101kq1r.html

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
提供的源码资源涵盖了Java应用等多个领域,每个领域都包含了丰富的实例和项目。这些源码都是基于各自平台的最新技术和标准编写,确保了在对应环境下能够无缝运行。同时,源码中配备了详细的注释和文档,帮助用户快速理解代码结构和实现逻辑。 适用人群: 适合毕业设计、课程设计作业。这些源码资源特别适合大学生群体。无论你是计算机相关专业的学生,还是对其他领域编程感兴趣的学生,这些资源都能为你提供宝贵的学习和实践机会。通过学习和运行这些源码,你可以掌握各平台开发的基础知识,提升编程能力和项目实战经验。 使用场景及目标: 在学习阶段,你可以利用这些源码资源进行课程实践、课外项目或毕业设计。通过分析和运行源码,你将深入了解各平台开发的技术细节和最佳实践,逐步培养起自己的项目开发和问题解决能力。此外,在求职或创业过程中,具备跨平台开发能力的大学生将更具竞争力。 其他说明: 为了确保源码资源的可运行性和易用性,特别注意了以下几点:首先,每份源码都提供了详细的运行环境和依赖说明,确保用户能够轻松搭建起开发环境;其次,源码中的注释和文档都非常完善,方便用户快速上手和理解代码;最后,我会定期更新这些源码资源,以适应各平台技术的最新发展和市场需求。 所有源码均经过严格测试,可以直接运行,可以放心下载使用。有任何使用问题欢迎随时与博主沟通,第一时间进行解答!

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值