I want the total lengths of regions that are exons, splicing, UTRs, and ncRNA. I can't distinguish them from the provided bed file. If there is a way using the bed file, please let me know. Thank you.
•
written
by
anp375
1 answer
You can use a GTF with annotated intervals for the categories you want. For example this one from Ensembl. Use something like unix grep to create files for the categories you are interested. Then use bedtools intersect to get the regions of your BED file overlapping with each category GTF file. Finally use the procedure described here to find total lengths.
Also, if you are only interested on the ratio for each category, you could simplify this by using bedtools jaccard. In this case you would only need your BED file and the categories GTF files.
Log in to answer this question.
-
Drosophila melanogaster dm3 intron percentage too high
written by Christian •
-
determine softmask regions from FASTA
written by Matteo Ungaro •
-
BED files
written by adarsh •
-
How to identify exon sequences
written by Ak •
-
Find count of reads in BAM file that fall within BED region
written by KH •
-
Bam file metrics
written by Sabeen •
-
How to extract exon, intron and UTR regions from gene?
written by ganeshram •
-
Extracting 5'UTR and 3'UTR bed files from gtf file
written by c_u
-
GTF separate coding and non-coding exons
written by igor
-
How to get distingushable CDS and UTR sequences of gene from biomart?
written by sohra •
