按照重测序课程中的命令
for i in CB2 CB3; do
echo "RUN CMD: bwa mem Prunus_dulcis.T2T.assembly.fa ${i}_1.fq.gz \
${i}_2.fq.gz -t 10 -M \
-R '@RG\tID:${i}\tLB:${i}\tPL:ILLUMINA\tSM:${i}' \
|samtools view -bS -h - > ${i}.bam"
nohup bwa mem Prunus_dulcis.T2T.assembly.fa ${i}_1.fq.gz \
${i}_2.fq.gz -t 10 -M \
-R '@RG\tID:${i}\tLB:${i}\tPL:ILLUMINA\tSM:${i}' \
|samtools view -bS -h - > ${i}.bam
done &