示例代码: bamtools filter -tag 'NM:i:<3' -in input.bam -out output.bam
此命令将读取序列中的NM标记小于3的读取序列写入到output.bam文件中。
示例代码: samtools view -b -s 0.1 -f 4 -o out.bam in.bam chr:1-1000
此命令将输入文件in.bam中的区域chr:1-1000截取出来,并将其写入到out.bam文件中。使用"-b"选项可以将BAM格式输出到文件,使用"-s"选项可以对文件进行子抽样,使用"-f"选项可以控制要过滤掉的读取标志。
示例代码: java -jar picard.jar CompressSam I=input.bam O=output.bam
此命令将输入文件input.bam中的所有数据重新排列并写入到output.bam文件中。压缩后的BAM文件可以大幅减小文件大小。
注意:如果需要确保保留的是完整的信息,请在操作前备份原始数据。