美文网首页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