interpreting inversions on IGV
Entering edit mode
4 months ago
Shaghayegh ▴ 20

I just started working with PacBio and Nanopore long reads, after some preliminary QCing I mapped the reads using minimap2. enter image description here

I am looking at the inversions for now and I have igv images for one of my regions. I want to know is it how I should expect to see the inversions on IGV? I do not understand whether the gaps are oK or should I expect the reads below the gaps or should I expect to see them filling the gaps? Thanks for any help interpreting it.

inversion long-reads igv • 424 views
Entering edit mode
4 months ago
cmdcolin ★ 3.2k

you may be able to turn on a setting "Link supplementary alignments" to show multiple segments of the alignment as a single, continuous row, with the inversion in the middle. in the SAM format, the inversion is likely split into e.g. 3 lines for a long read that goes through the whole inversion, one line is e.g. forward, one like is reverse, and then the other is forward, and one of the lines will be marked as primary (no flag 2048) and the other two will be marked as supplementary (flag 2048) and they will all share the same read name and have an SA tag linking them


Entering edit mode

Thanks a lot, very helpful. I activated "linked supplementary alignments" and can see gaps are filled.


Login before adding your answer.

Traffic: 1071 users visited in the last hour
Help About
Access RSS

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

Powered by the version 2.3.6