I am looking for a tool that calculates span of GFF file (for example total number of bp that are composed of repeats or genes, given these are described in GFF file). I wrote awk script to do that (sums lengths of all features), but since this way I have no validation of inputs, I would prefer tool similar to bedtools to do this.
•
written
by
Biomonika (Noolean)
1 answer
It seems to me that the tool you're looking for is bedtools genomecov.
You'll probably want to use the option "-max 1" if you're only interested to know how many bases are spanned by some type of features.
For example, you might do something like this:
bedtools genomecov -max 1 -i <(awk '$3 == "start_codon" {print}' file.gff) -g hg19.genome > start_codon.txt
The output file will tell you how many bases on each chromosome are spanned by start codons, and also how many bases in the whole genome are spanned by start codons.
Log in to answer this question.
-
Drosophila melanogaster dm3 intron percentage too high
written by Christian •
-
random list of coordinates of different lengths
written by Chironex •
-
cwlexec and cwltool not interpreting baseCommand the same way
written by drkennetz
-
bp overlap between multiple beds
written by varsha619
-
identify length of repeat arrays with given motif
written by Biomonika (Noolean)
-
Bedtools coverage adding feature scores
written by jaime alvarez •
-
length of combined transcripts
written by bsp017 •
-
How to get gene names from tophat results
written by sebastianzeki0
-
per bp depth intersected with gff file
written by Biomonika (Noolean)
-
Finding Overlapping Variants (I.E. Indels, Snps) Using Annovar Format.
written by youngcsong •
