Question: How to load BAM files that do not have chromosome name in the header?
0
gravatar for JC
10 months ago by
JC0
JC0 wrote:

I cannot view a BAM file with IGV

It complains that "File RNAseq_no1.bam does not contain any sequence names which match the current genome. File:      hg19_refGene_NM_032291, hg19_refGene_NM_032785, hg19_refGene_NM_018090, hg19_refGene_NM_052998, ... Genome: chr1, chr2, chr3, chr4, ..."

I have take look at the head file with samtools:

samtools view -H RNAseq_no1.bam

Here are the first couple lines in the Header:

@HD VN:1.0  GO:none SO:queryname
@SQ SN:hg19_refGene_NM_032291   LN:4894
@SQ SN:hg19_refGene_NM_032785   LN:3188
@SQ SN:hg19_refGene_NM_018090   LN:2281
@SQ SN:hg19_refGene_NM_052998   LN:2382
@SQ SN:hg19_refGene_NM_001145278    LN:2203

It doesn't contain any chromosome name. I think this is what IGV complains about. Does anyone know how I can possibly load the BAM files into IGV?

Thank you.

igv rna-seq • 447 views
ADD COMMENTlink modified 10 months ago by WouterDeCoster40k • written 10 months ago by JC0

I added (code) markup to your post for increased readability. You can do this by selecting the text and clicking the 101010 button. When you compose or edit a post that button is in your toolbar, see image below:

101010 Button

ADD REPLYlink written 10 months ago by WouterDeCoster40k
1
gravatar for Chris Miller
10 months ago by
Chris Miller21k
Washington University in St. Louis, MO
Chris Miller21k wrote:

It looks like you have a sequence data that was aligned to a transcriptome, rather than a genome. If you want to use IGV, you will either need to a) create a new "genome" for IGV with one contig per gene (this is almost certainly not what you want) or b) align your RNAseq data to a genome, using something like HISAT2, STAR, etc.

ADD COMMENTlink written 10 months ago by Chris Miller21k
0
gravatar for swbarnes2
10 months ago by
swbarnes26.5k
United States
swbarnes26.5k wrote:

Those bam files were made by aligning to a list of transcripts. You cannot view them in IGV as aligned to chromsomes, but you can load the fasta of transcripts into IGV, and look at them by transcript.

ADD COMMENTlink written 10 months ago by swbarnes26.5k

Could you tell me step-by-step how to do it?

ADD REPLYlink written 10 months ago by JC0
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 576 users visited in the last hour