This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Bedtools intersect error

Hi, I am trying to use bedtools intersect to intersect between a vcf and a bed file. My command is this

bedtools intersect -header -u -a ABC.vcf -b XYZ.bed > check.vcf

My bed file looks like this

chr6    29000001        33000000

But everytime I try this. It gives me the following error:

Error: Invalid record in file ABC.vcf. Record is  
chr6        60002        .        A        <*>        0        .        DP=5;I16=5,0,0,0,191,7301,0,0,208,8656,0,0,47,495,0,0;QS=1,0;MQ0F=0     PL  0,15,104

I tried to remove this particular position and redo. But it gives error for the next position. What could be the issue here ? Any help would be appreciated.

software error sequence vcf

Can you first check if all the file formats you are using are OK for input to bedtools? Perhaps it's the <*> that is causing problems?

Yes, I have rechecked all the file formats from the original sam file to the final input vcf file. Those seem to be fine

I'm far from an expert here but I don't see the <*> as a valid option in the vcf file format . Did you try removing/altering this

0 answers

No answers yet.

Log in to answer this question.