trimmomatic的介绍

trimmomatic的官网trimmomatic
trimmomatic可以对测序数据进行过滤
java -jar trimmomatic-0.35.jar PE -phred33 input_forward.fq.gz input_reverse.fq.gz output_forward_paired.fq.gz output_forward_unpaired.fq.gz output_reverse_paired.fq.gz output_reverse_unpaired.fq.gz ILLUMINACLIP:$file_path/TruSeq3-PE.fa:2:30:10 LEADING:3 TRAILING:3 SLIDINGWINDOW:4:15 MINLEN:36

运行上面的命令可以完成以下任务
Remove adapters (ILLUMINACLIP:TruSeq3-PE.fa:2:30:10) #去掉接头
Remove leading low quality or N bases (below quality 3) (LEADING:3) #去掉开头质量低于3或N碱基
Remove trailing low quality or N bases (below quality 3) (TRAILING:3) #去掉末尾质量低于3或N碱基
Scan the read with a 4-base wide sliding window, cutting when the average quality per base drops below 15 (SLIDINGWINDOW:4:15) #以4个碱基为滑动窗口对reads进行扫描,当平均质量值低于15时进行剪切
Drop reads below the 36 bases long (MINLEN:36) #去掉长度小于36的reads

©著作权归作者所有,转载或内容合作请联系作者
【社区内容提示】社区部分内容疑似由AI辅助生成,浏览时请结合常识与多方信息审慎甄别。
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

相关阅读更多精彩内容

友情链接更多精彩内容