將SAM文件轉(zhuǎn)換為BAM文件：首先，將SAM文件轉(zhuǎn)換為BAM文件，這樣更容易處理?？梢允褂胹amtools進(jìn)行轉(zhuǎn)換：

samtools view -bS input.sam > input.bam

從BAM文件中提取比對上的reads：使用samtools，你可以從BAM文件中提取比對上的reads。在此示例中，我們假設(shè)比對上的reads具有flag值不為4：

samtools view -b -F 4 input.bam > mapped.bam

將比對上的reads轉(zhuǎn)換為雙端FASTQ文件：使用samtools的fastq命令，你可以將比對上的reads從BAM文件中提取并輸出為雙端FASTQ文件：

samtools fastq -1 mapped_1.fq -2 mapped_2.fq -s mapped_singles.fq mapped.bam

這將從mapped.bam中提取比對上的reads，并將它們保存為兩個雙端FASTQ文件（mapped_1.fq和mapped_2.fq），以及一個單端FASTQ文件（mapped_singles.fq）。