This is a test version of Biostars. For the public version, visit https://www.biostars.org.
how to change bam file chromosome names using sed or awk?

I have several bam files with illegal chromosome names, eg: ENSMUST00000193812.1_4933401J01Rik

I wanna visualize these bam files in IGV and those illegal names are causing problems.

I know that I can edit the reference .fasta file header and do mapping again, But I really wanna know how to directly edit the bam file using awk or sed command. I see similar post here: Bam File: Change Chromosome Notation However it is quite difficult for me to understand and apply it to my own situation.

What I wanna do is just replace all the dot in the bam file chromosome name using lower dash.

Thanks in advance for any help!

bam awk
• 2,469 views • link
updated by swbarnes2 swbarnes2 • written by qianxu0517 • qianxu0517

ENSMUST00000193812.1_4933401J01Rik

how do you know it's 'illegal' ? why is it causing any problem ?

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

2 answers

You should use samtools reheader. Do not do this with sed or awk - you could compromise the file.

• 0 views • link
written by Ram Ram

Thanks! I think I need to rename the fast file and then do mapping agian.

• 0 views • link
written by qianxu0517 • qianxu0517

bam files are compressed. You'd have to feed in the uncomprassed version into sed or awk. Easiest way to do that with with samtools, in which case, you might as well use samtools to do it by itself.

• 0 views • link
written by swbarnes2 swbarnes2

Log in to answer this question.

More posts like this