IGV deletion insert size
1
0
Entering edit mode
14 months ago
robinycfang ▴ 20

Hi,

I am trying to extract variant reads for deletion (large SVs) from DNA-seq bam files. In IGV, variant reads for deletions are labelled in red (insert size larger than expected): IGV link. I am currently using the pysam package to parse bam files. I have a look at those deletion variant reads in red, but couldn't find unique features to extract them out from the bam files...not sure how IGV identifies those reads. Any suggestions? Thanks!

deletion IGV variantion structural DNA • 730 views
ADD COMMENT
1
Entering edit mode
14 months ago

but couldn't find unique features to extract them out from the bam files...

  • reads not-properly-paired with a large TLEN (distance) between the read and it's mate.
    • reads with supplementary alignments mapping far from the original read.
ADD COMMENT
0
Entering edit mode

Thanks! This gave me some ideas! I guess I will start working on measuring the inner distance between two paired reads (segments), and if the inner distance is unexpected large (e.g., >>100bp), then they are potential variant reads. Is that correct?

ADD REPLY

Login before adding your answer.

Traffic: 1304 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