Issue with bam-readcount and fpfilter.pl: Saying "failed to get readcounts for variant allele"
0
0
Entering edit mode
2 days ago
1phamcha2 • 0

Hi everyone,

I'm using Somatic Sniper and working on somatic mutation indel filtering using bam-readcount and fpfilter.pl, but I keep running into an issue where fpfilter.pl fails to retrieve read counts for all variants.

Below is an outline of how I processed and filtered SomaticSniper variant calls using standard filtering, bam-readcount, and fpfilter.pl:

> perl snpfilter.pl –snp-file my_sniper_file –indel-file my_indel_pileup
> perl prepare_for_readcount.pl –snp-file my_sniper_file.SNPfilter
> bam-readcount -b 15 -f my_ref.fasta -l my_sniper_file.SNPfilter.pos my_tumor.bam > my_readcounts.rc
> perl fpfilter.pl –snp-file my_sniper_file.SNPfilter –readcount-file my_readcounts.rc

However, I keep getting:
514 variants
514 failed to get readcounts for variant allele
0 had read position < 0.1
0 had strandedness < 0.005
0 had var_count < 2
0 had var_freq < 0.02
0 had mismatch qualsum difference > 50
0 had mapping quality difference > 40
0 had read length difference > 25
0 had var_distance_to_3' < 0.1
0 passed the strand filter

What I've checked so far:
-Checked that positions match in both files (my_readcounts.rc and my_sniper_file.SNPfilter)
ex. of output for readcount
chr1    759567  T   3   =:0:0.00:0.00:0.00:0:0:0.00:0.00:0.00:0:0.00:0.00:0.00  A:0:0.00:0.00:0.00:0:0:0.00:0.00:0.00:0:0.00:0.00:0.00      C:0:0.00:0.00:0.00:0:0:0.00:0.00:0.00:0:0.00:0.00:0.00   G:3:55.33:30.33:55.33:0:3:0.53:0.02:455.00:0:0.00:76.00:0.74        T:0:0.00:0.00:0.00:0:0:0.00:0.00:0.00:0:0.00:0.00:0.00 :0.00:0.00:0.00:0:0:0.00:0.00:0.00:0:0.00:0.00:0.00
chr1    915173  C   5   =:0:0.00:0.00:0.00:0:0:0.00:0.00:0.00:0:0.00:0.00:0.00  A:0:0.00:0.00:0.00:0:0:0.00:0.00:0.00:0:0.00:0.00:0.00      C:5:45.60:29.00:45.60:5:0:0.24:0.02:425.80:5:0.74:76.00:0.74     G:0:0.00:0.00:0.00:0:0:0.00:0.00:0.00:0:0.00:0.00:0.00      T:0:0.00:0.00:0.00:0:0:0.00:0.00:0.00:0:0.00:0.00:0.00 :0.00:0.00:0.00:0:0:0.00:0.00:0.00:0:0.00:0.00:0.00

ex. of output for SNPfilter
chr1    759567  .   A   G   .   .   .   GT:IGT:DP:DP4:BCOUNT:GQ:JGQ:VAQ:BQ:MQ:AMQ:SS:SSC    /0:4:0,4,0,0:4,0,0,0:39:.:0:30:60:60:0:.    1/1:3:0,0,0,3:0,0,3,0:36:.:36:30:55:55:2:44
chr1    915173  .   A   C   .   .   .   GT:IGT:DP:DP4:BCOUNT:GQ:JGQ:VAQ:BQ:MQ:AMQ:SS:SSC    /0:2:2,0,0,0:2,0,0,0:33:.:0:30:57:57:0:.    1/1:5:0,0,5,0:0,5,0,0:42:.:42:29:46:46:2:4

8

What I've tried: -Removing the normal column from the SNP file, didn't work (same results) -Adjusting the fpfilter.pl parameters so its less selective, didn't work (same results)

Any help would be greatly appreciated, thank you!

bam-readount somaticsniper • 235 views
ADD COMMENT

Login before adding your answer.

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