IGV pileup not showing any reads at any position in bam
0
0
Entering edit mode
6.8 years ago
bioguy24 ▴ 230

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 • 2.3k views
ADD COMMENT
0
Entering edit mode

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

ADD REPLY

Login before adding your answer.

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