How to include strand information in the result when running bedtools complement?
Hi everyone,
I am using bedtools complement to find the regions not covered in a GFF file.
Genome file is as follows.
chr1 12222
GFF file likes below.
chr1 db exon 1 45 . + . ID=gene1
chr1 db exon 90 450 . + . ID=gene2
chr1 db exon 1 30 . - . ID=gene99
I run complementBed -i <gff_file> -g <genome_file> and got the following result.
chr1 45 89
chr1 450 12222
This is not what I expected. I hope to get the following information, which includes strand symbol. Could anyone help? Thank you very much!
chr1 45 89 +
chr1 450 12222 +
chr1 30 12222 -
1 answer
Log in to answer this question.
More posts like this
-
Remove multiple gene id from bed file
written by anii •
-
Tumor mutation burden calculation based on bed file?
written by lincaijin1994 •
-
Average length of each exon and intron
written by kashiff007
-
How can I number all exons and then subsequently subdivide all exons into bins in a GFF file?
written by lasra •
-
converting dexseq_count.py 's output of Gene-exon ID to Transcript Id
written by Ice P •
-
Taking genotypes out of a vcf file
written by Jautis
-
use bedtools to find overlaps of bed and gff files
written by biolab
-
RSeQC (infer_experiment) problem
written by biolab
-
how to produce mapping graph based on blat result?
written by biolab
-
Identify Overlapping And Non Overlapping Regions For Paired-End Data
written by ancient_learner
