$ grep "^>" my_fasta.fna | cut -c 2- > fasta.txt
$ samtools view -H my_bam.bam | grep "^@SQ" | samtools view -H bam10k/10.bam | grep "^@SQ" | cut -f 2 | cut -c 4- > bam.txt
$ diff fasta.txt bam.txt
用来查看 BAM 和 FASTA文件header的差异
最新推荐文章于 2024-06-04 23:11:16 发布