美文网首页Whole Genome Resequencing
【bwa mem比对问题汇总】

【bwa mem比对问题汇总】

作者: 云养江停 | 来源:发表于2021-09-08 08:42 被阅读0次

bwa mem比对出现:skip orientation RF as there are not enough pairs

This is normal. The Illumina reads orientation is FR. 
You got majority FR and skipped other three.

mapping生成sam文件时出现[mem_sam_pe] paired reads have different names错误

方法一(这种方法我在使用后依旧报错):

bbrename.sh in1=read1.fq in2=read2.fq out1=renamed1.fq out2=renamed2.fq
#To install this package with conda run one of the following: conda install -c bioconda bbmap

对于多个fq文件,可以用以下命令:
while read line
do
      nohup bbrename.sh in=${line}_R1.fq in2=${line}_R2.fq out=${line}_R1.sh.fq out2=${line}_R2.sh.fq &    
done < name.txt #####其中,name.txt指的是包含sample名字文件

方法二 使用repair.sh进行修复:(使用这种方法后fastq文件变小,可能是合并部分文件造成的。
此外, fastq.gz也可以,另外就是记得输出文件要和源文件命名不同,以防覆盖。

while read line
do
      nohup repair.sh in=${line}_R1.fastq in2=${line}_R2.fastq out=${line}_R1.sh.fastq out2=${line}_R2.sh.fastq &    
done < name.txt

相关文章

网友评论

    本文标题:【bwa mem比对问题汇总】

    本文链接:https://www.haomeiwen.com/subject/ugprwltx.html