將SAM文件轉(zhuǎn)換為BAM文件：首先要拂，將SAM文件轉(zhuǎn)換為BAM文件，這樣更容易處理÷塘埽可以使用samtools進行轉(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）闯传。
