Question: Filter vcf file using coordinate of bed , the result is SnpSift had filtered some wrong items, why ?
gravatar for duolin1
3 months ago by
duolin110 wrote:

when I using snpsift software filter vcf by coordinate of bed, I found the result of vcf have wrong item that non_exist in bed ;

this is bed's coordinate

chr8    29608637    29609316
chr8    29609651    29609846

this is the "wrong" item that I found in result that is filtered by bed

< chr7  123666723   .   AT  A   207.73  PASS    AC=1;AF=0.500;AN=2;BaseQRankSum=-0.749;ClippingRankSum=0.860;DP=28;FS=0.000;MLEAC=1;MLEAF=0.500;MQ=60.39;MQ0=0;MQRankSum=-0.416;QD=7.42;ReadPosRankSum=-0.250;SOR=0.495 GT:AD:DP:GQ:PL  0/1:13,15:28:99:245,0,194

< chr8  29609647    .   TAAAAAGA    T   179.73  PASS    AC=1;AF=0.500;AN=2;BaseQRankSum=0.167;ClippingRankSum=-1.303;DP=23;FS=2.066;MLEAC=1;MLEAF=0.500;MQ=60.00;MQ0=0;MQRankSum=0.702;QD=7.81;ReadPosRankSum=0.634;SOR=1.402   GT:AD:DP:GQ:PL  0/1:16,7:23:99:217,0,976

as you can see ,actually this is "wrong" item because this "chr8...." is not accord with bed's coordinate; by the way, I found the wrong item often happened in the first item of the new "chr"( it can also see in my provided item)

this is my code like this

 cat variants.vcf | java -jar SnpSift.jar intervals my_intervals.bed > variants_intersecting_intervals.vcf

so,my question is why , is it a bug in snpsift???

snp snpeff vcf snpsift • 198 views
ADD COMMENTlink modified 3 months ago by zx87548.2k • written 3 months ago by duolin110

which version of snpsift ?

how about just using bcftools view?

ADD REPLYlink written 3 months ago by Pierre Lindenbaum122k

I'd report it to the SnpSift author.

It is a 7bp deletion which overlaps your specified coordinate so it might make sense depending on how you look at it. I would prefer to have such hits in my results.

ADD REPLYlink written 3 months ago by colindaven1.7k

thank you very much. I understand..

ADD REPLYlink written 3 months ago by duolin110
Please log in to add an answer.


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