Read coverage shows a difference in IGV
Hello,
This is my bam file in IGV. There are way more reads mapped to the region 2 in the panel displaying reads. On the other hand, the coverage panel shows the opposite - there are more reads mapped to the region 1.
I also ran samtools depth and it said there are more reads mapped to the region 1. I am confused by this. Did IGV not show me all the reads mapped to the region 1 for some reason? I didn't change any setting for IGV.
1 answer
Log in to answer this question.
More posts like this
-
How to calculate the depth of nanopore reads in IGV
written by kirankumareripogu •
-
How can I find reads that have multi-mapped with multiple genes overlapping in the same region
written by Berk •
-
samtools depth -r
written by lz1903 •
-
Remove reads from bam file that partially covers a region
written by arbarharw •
-
No coverage in macs2 peak summit!
written by blur
-
Extracting only reads mapped to specific genomic region
written by agata •
-
Detect Transgene Insertion Site using SAM format
written by kspata
-
BAM File Export Coverage Values from IGV or using PySam/Samtools
written by EOX123 •
-
Pysam - Different amount of mapped and aligned reads
written by AlchemistMoz •
-
view of a circular viral genome and its annotation
written by tonja.r
may be the IGV algorithm collapsing reads generates a visual bias.
I see. I will post this on the IGV help forum.
perhaps it has something to do with the mapping quality? Image is a bit small but it seems there are much more gaps/indels/mismatches in region2 in stead of 1 ...
perhaps samtools only counts the well mapped reads or such?
Well, I filtered reads out in my bam by mapping quality. Yeah, I will look into
samtools depthmore.I think I found why. If I unclick
Downsamplingoption in alignment in preferences, it looks like IGV shows every read...