多样性测序单端barcode分

1.trimmotic质控

trimmomatic-0.30.jar PE -phred33 all.R1.fastq all.R2.fastq trim.1.fastq unpair.1.fastq trim.2.fastq unpair.2.fastq LEADING:0 TRAILING:20 SLIDINGWINDOW:50:20 MINLEN:50 2> trim.log

2.flash拼接

flash -m 10 -M 100 -x 0.2 -p 33 -t 6 -d merge -o  all.trim  trim.200.1.fastq trim.200.2.fastq >merge.log 

3.创建一个barcode.list文件,样品名+ barcode; 然后直接进行下列命令

barcode list格式

4、拆分

less 3.trim.extendedFrags.fastq | fastx_barcode_splitter.pl --bcfile barcode.list --bol --mismatches 0 --prefix data --suffix ".fq"

使用软件:

http://hannonlab.cshl.edu/fastx_toolkit/commandline.html

https://ccb.jhu.edu/software/FLASH/

亲亲们支持一下!这里是小编编的画图店,【不会画图请找我】https://shop113546122.taobao.com/?spm=2013.1.1000126.d21.7c4ae3e8LU3Kmg,有生信问题和R问题欢迎来咨询~请小编编吃外卖就可以啦~

©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

推荐阅读更多精彩内容