去核糖体后的文件用hisat2比对的问题

2021-03-05  本文已影响0人  宗肃書

当用hisat2比对会报错

(base) [jychu@localhost cleandata]$ hisat2 -p 20 -x /public/jychu/zhengxt/refs/Oar_v3.1/hisat2_index/Oar_v3.1_index -1 Born-3_nrRNA_1.fq -2 Born-3_nrRNA_2.fq | samtools view -S -b -> ../hisat2/Born-3.bam
Error: Read AFFFJJJJFJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJFJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJFJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJFFFF has more read characters than quality values.
terminate called after throwing an instance of 'int'
(ERR): hisat2-align died with signal 6 (ABRT) (core dumped)
(base) [jychu@localhost cleandata]$ hisat2 -p 20 -x /public/jychu/zhengxt/refs/Oar_v3.1/hisat2_index/Oar_v3.1_index -1 E120-1_nrRNA_1.fq -2 E120-1_nrRNA_2.fq -S E120-1.sam
Error: Read E120-1.14.2 14 length=150 has more read characters than quality values.
terminate called after throwing an instance of 'int'
(ERR): hisat2-align died with signal 6 (ABRT) (core dumped)

在网上搜索是因为文件里的reads双端不一样长,检查文件也没有问题

解决方法:用bwa比对

上一篇 下一篇

猜你喜欢

热点阅读