Question: Filter vcf file using coordinate of bed , the result is SnpSift had filtered some wrong items, why ?
gravatar for duolin1
26 days 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 • 125 views
ADD COMMENTlink modified 25 days ago by zx87547.5k • written 26 days ago by duolin110

which version of snpsift ?

how about just using bcftools view?

ADD REPLYlink written 25 days ago by Pierre Lindenbaum120k

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 25 days ago by colindaven1.3k

thank you very much. I understand..

ADD REPLYlink written 24 days 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: 1135 users visited in the last hour