Question: Problem to Understand the visualization of IGV
0
Entering edit mode

I use BWA to create a bam file and visualize the alignment with IGV. I didn't split the window. But I get two alignments exact at the position of my gene of interest. What does it mean? Untitled_2

ADD COMMENTlinkeditmoderate 15 months ago saadleeshehreen • 60 • updated 15 months ago Devon Ryan 90k
Entering edit mode
0

Can you highlight the portion which you are looking for in the image?

ADD REPLYlinkeditmoderate 15 months ago
toralmanvar
• 760
Entering edit mode
0

red one ... IGV_screenshot

ADD REPLYlinkeditmoderate 15 months ago
saadleeshehreen
• 60
Entering edit mode
0

Does your BAM file have read groups?

ADD REPLYlinkeditmoderate 15 months ago
Devon Ryan
90k
Entering edit mode
0

Its my first time I am working with IGV , BWA and samtools. I tried the same procedure with other genome and didn't see such alignment there. That seems quite normal. I don't know what happens with it. What do you mean by 'read groups'?

ADD REPLYlinkeditmoderate 15 months ago
saadleeshehreen
• 60
Entering edit mode
0

In situations where you might want to do variant calling you end up needing to merge multiple biological samples together in 1 BAM file. You can then label each alignment in that file to designate the sample it came from. These labels are call "read groups". My expectation is that this BAM file has read groups in it. You can see that by posting the output of samtools view -H your_file.bam | tail.

ADD REPLYlinkeditmoderate 15 months ago
Devon Ryan
90k

Login before adding your answer.

Powered by the version 2.0