批量bwa比对和samtools排序

最新推荐文章于 2024-12-31 09:43:05 发布

爱笑的小牙

最新推荐文章于 2024-12-31 09:43:05 发布

阅读量6.1k

点赞数 1

CC 4.0 BY-SA版权

分类专栏：生物信息分析

本文链接：https://blog.csdn.net/Cassiel60/article/details/89025882

生物信息分析专栏收录该内容

45 篇文章

订阅专栏

本文介绍如何处理FASTQ测序数据文件，包括为参考序列建立索引、使用bwa工具进行批量比对及利用samtools进行排序和索引等步骤。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

拿到fastq文件后，要进行比对和排序

第一步是对reference建立索引

bwa index -a bwtsw $reference         #对于大基因组建立FM-Index

#bwa index -a is ref.fasta             #对小基因组建立index，速度快，内存消耗大

第二步，批量比对和排序

# 2.使用bwa men 比对
mkdir -p ../fina1_bam
cat SRR_Acc_List.txt|while read line 
do 
$bwa mem -t 16 -M -Y  -R "@RG\tID:$line\tPL:ILLUMINA\tLB:$line\tSM:$line" $reference ../bam/${line}.fastq.gz | $samtools view -Sb -t 16 -f bam > ../fina1_bam/${line}.bam &&
echo "** ${line} BWA MEM done **"
# 3.排序
$samtools sort -@ 16 -m 4G  -O bam ../fina1_bam/${line}.bam -o ../fina1_bam/${line}.sorted.bam && \
$samtools index ../fina1_bam/${line}.sorted.bam
echo "** ${line} sorted raw bam file done **"
done