Question: IGV pileup not showing any reads at any position in bam
0
gravatar for bioguy24
24 months ago by
bioguy24190
Chicago
bioguy24190 wrote:

I am not able to figure out why when I load an Ion Torrent bam file into IGV the result is a blank/empty pileup. I routinely run these bam files and have no problems, but am now and not sure why. This is a bam uses a hotspot file, but I am not sure if that is necessary or how to load it properly.

 samtools view -h torrent_316.2017_A.bam | head -20

 @HD    VN:1.4  GO:none SO:coordinate
 @SQ    SN:chr1 LN:249250621
 @SQ    SN:chr2 LN:243199373
 @SQ    SN:chr3 LN:198022430
 @SQ    SN:chr4 LN:191154276
 @SQ    SN:chr5 LN:180915260
 @SQ    SN:chr6 LN:171115067
 @SQ    SN:chr7 LN:159138663
 @SQ    SN:chr8 LN:146364022
 @SQ    SN:chr9 LN:141213431
 @SQ    SN:chr10    LN:135534747
 @SQ    SN:chr11    LN:135006516
 @SQ    SN:chr12    LN:133851895
 @SQ    SN:chr13    LN:115169878
 @SQ    SN:chr14    LN:107349540
 @SQ    SN:chr15    LN:102531392
 @SQ    SN:chr16    LN:90354753
 @SQ    SN:chr17    LN:81195210
 @SQ    SN:chr18    LN:78077248
 @SQ    SN:chr19    LN:59128983
 @SQ    SN:chr20    LN:63025520
 @SQ    SN:chr21    LN:48129895
 @SQ    SN:chr22    LN:51304566
 @SQ    SN:chrX LN:155270560
 @SQ    SN:chrY LN:59373566
 @SQ    SN:chrM LN:16569
 @RG    ID:PP5M3.IonXpress_002  PL:IONTORRENT   PU:PGM/316v2/IonXpress_002  LB:hg19/IonXpress_002   FO:TACGTACGTCTGAGCATCGATCGATGTACAGCTACGTACGTCTGAGCATCGATCGATGTACAGCTACGTACGTCTGAGCATCGATCGATGTACAGCTACGTACGTCTGAGCATCGATCGATGTACAGCTACGTACGTCTGAGCATCGATCGATGTACAGCTACGTACGTCTGAGCATCGATCGATGTACAGCTACGTACGTCTGAGCATCGATCGATGTACAGCTACGTACGTCTGAGCATCGATCGATGTACAGCTACGTACGTCTGAGCATCGATCGATGTACAGCTACGTACGTCTGAGCATCGATCGATGTACAGCTACGTACGTCTGAGCATCGATCGATGTACAGCTACGTACGTCTGAGCATCGATCGATGTACAGCTACGTACGTCTGAGCA DT:2014-10-10T17:17:22-0500 SM:NGS_02_CAP_316   PG:tmap KS:TCAGTAAGGAGAACGAT    CN:TorrentServer/sn274670099

Also, the result of this command shows zero reads at this position:

 samtools view -c torrent_316.2017_A.bam "chr1:115258747"
 0

However, the variant caller in Ion Reporter called a SNV with 1935 reads.

Thank you :).

ngs • 896 views
ADD COMMENTlink modified 24 months ago • written 24 months ago by bioguy24190

Run a flagstat on this not just an intersection of a single coordinate.

ADD REPLYlink written 24 months ago by Istvan Albert ♦♦ 80k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1688 users visited in the last hour