I'm currently waiting for a pipeline I'm running to produce some BED files and VCF files, which I will then intersect with a GFF file of genes that I want. The command I'm running will be:
bedtools intersect -wa -wb -a file.bed -b searchSpace.gff > out.bed
bedtools intersect -wa -wb -a snps.vcf -b searchSpace.gff > out.vcf
I know that GFF runs under a different coordinate system than BED or VCF format. I also know that bedtools natively accepts these types of files. My question is, when I use the -wa -wb option to write the coordinate features from both input files, will there be a discrepancy in the output? Or will bedtools automatically take care of that? I need to know whether it's necessary to take the extra step of converting all of my GFF files to BED coordinate system, or if doing so will result in a double-correction by bedtools.
0 answers
No answers yet.
Log in to answer this question.
-
Bedtools merge but only if intersection exists
written by jrose •
-
bedtools interset doesn't return a VCF file?
written by ManuelDB
-
Determine the significance of peak overlaps (bedtools)
written by kstangline •
-
intersection of two GFFs
written by BioDH •
-
Bedtools intersect vcf and bed files
written by User000
-
terminate called after throwing an instance of 'std::bad_alloc?
written by star
-
Error using bedtools "intersect" command
written by hila •
-
bedtools closest - result for only the first query
written by yjyj •
-
bedtools intersect "complete outer join"
written by graeme.thorn
-
Get completely overlapping features
written by iraun
Apologies but your question appears strange... Why do you not just execute these commands to see what they produce? Then, you will be enlightened.
I figured I'd ask in case anybody else had the question and wanted to quickly find an answer by googling it.
Okay, and is everything okay now?