Bedtools intersect behavior under different coordinate systems?
0
0
Entering edit mode
6.3 years ago
edsouza ▴ 30

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.

gene genome • 1.2k views
ADD COMMENT
0
Entering edit mode

Apologies but your question appears strange... Why do you not just execute these commands to see what they produce? Then, you will be enlightened.

ADD REPLY
0
Entering edit mode

I figured I'd ask in case anybody else had the question and wanted to quickly find an answer by googling it.

ADD REPLY
0
Entering edit mode

Okay, and is everything okay now?

ADD REPLY

Login before adding your answer.

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