Genome Information From Bam

Latest

Open

RNA-Seq

ChIP-Seq

SNP

Assembly

Tutorials

Tools

Jobs

Forum

Planet

All »

Home

Welcome to Biostar! about • faq • rss

Community

Log In

Sign Up

Add New Post

Question: Genome Information From Bam

5.2 years ago by

ChIP • 500

Netherlands

ChIP • 500 wrote:

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

bam chip-seq genomics • 1.1k views

ADD COMMENT • link •

modified 5.2 years ago by Pierre Lindenbaum ♦ 120k • written 5.2 years ago by ChIP • 500

5.2 years ago by

Pierre Lindenbaum ♦ 120k

France/Nantes/Institut du Thorax - INSERM UMR1087

Pierre Lindenbaum ♦ 120k wrote:

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

ADD COMMENT • link written 5.2 years ago by Pierre Lindenbaum ♦ 120k

Please log in to add an answer.

Similar posts • Search »