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