Genome Information From Bam

Question

[ChIP]

Hi!

Is their a way to know the genome build(ie., hg18 or hg19) to which the BAM file is aligned, other than manual inspection.

Thank you

Answer 1

[Pierre Lindenbaum]

get the dictionary in the sam header and test one chromosome

$ samtools view -H file.bam |\
 grep -E 'SN\:(chr)?22' |\
 awk '{if(index($0,"LN:51304566")!=0) printf("hg19\n"); else if(index($0,"LN:49691432")!=0) printf("hg18\n"); else printf("??\n");}'
hg19