visualize my .BAM file in IGV
2
0
Entering edit mode
4.5 years ago
chrisgr ▴ 20

Hi,

I got a little confused trying to visualize my bowtie2 mapped reads. I mapped some short (synthetic) RNA reads to a self-made (short) reference sequence with bowtie2, and I have aquired the .bam file.

Now I want to visualize this .bam file. I wanted to do this with IGV, however they want me to also provide a 'genome sequence' (which I understood is the reference sequence). Isn't the reference sequence already provided in the .bam file? Why do I need to provide a reference sequence again? I just want to visualize the mapping bowtie2 did.

I'd have to read more of the manuals of IGV but I can't seem to find an answer to this question. If there is an informaive page that could help me with this specific question I'd also be very happy!

Thanks in advance!

Chris

RNA bowtie2 IGV alignment BAM • 3.2k views
ADD COMMENT
0
Entering edit mode

what is the reference organism? @ christiangriffioen

ADD REPLY
0
Entering edit mode

Not an organism, the reads mapped are synthetic, so I used a reference genome made up of these sequences pasted after each other.

ADD REPLY
0
Entering edit mode

IGV has few inbuilt reference genomes and allows to load external sequences as well. In your case, you may have to follow custom genome load instructions for IGV.

ADD REPLY
2
Entering edit mode
4.5 years ago
GenoMax 146k

however they want me to also provide a 'genome sequence' (which I understood is the reference sequence)

No. The reference can be any sequence. In your case it will be self-made (short) reference sequence if that is what you used to create the bowtie index and for alignments. Create a custom genome (.genome file) in IGV with that.

ADD COMMENT
2
Entering edit mode
4.5 years ago

Isn't the reference sequence already provided in the .bam

no, the bam only stores the positions and sequences of the reads.

The reference sequence is the fasta sequence you used to build the bowtie index. Most of the time, it must be indexed with samtools faidx.

ADD COMMENT
0
Entering edit mode

Many thanks for the reply, I thought/read about this and don't understand yet, I really want to grasp this.

Like you said, the .bam file stores the positions and sequences of the reads on the reference sequence I provided right? I provided the reference sequence when I mapped them to the reference sequence using bowtie2, why is the reference sequence not included in the .bam file? If in IGV I select a reference genome and it would be slighty different, won't that mess with the .bam 'positions 'of the reads?

ADD REPLY
0
Entering edit mode

why is the reference sequence not included in the .bam

because the BAM format just don't do this.

If in IGV I select a reference genome and it would be slighty different, won't that mess with the .bam 'positions 'of the reads?

if slighty means a few subsitutions you'll see some mismatches in the reads.

if slighty means a few insertions/deletions : you'll get a mess.

ADD REPLY
0
Entering edit mode

why is the reference sequence not included in the .bam

imagine your mapping only one read to the human genome. The small bam would contains the 3G from the human genome ? no.

ADD REPLY

Login before adding your answer.

Traffic: 819 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6