Hi. I have a folder full of .fa files, and a .gff. The gff file contains information about which loci look like they code for RNA sequences. The .fa contain the DNA sequences for a set of human chromosomes. I want to get all the sequences which code for RNA, as defined by the gff file, out of the DNA in the fasta files. I also have a file telling me which RNA types have higher priority (lincRNA is higher priority than miRNA for example), this tells me which are more important and how I should decided between RNAs for overlapping reads in the gff.
I have been trying to code my own little program in F# that will read these files and give me each RNA read defined in the gff, and its corresponding DNA. However I am a bit confused about how it works. Do the start and end of each feature in the gff file define a character in the corresponding .fa file? Are they 1 or 0 indexed? Does it matter what strand they are ('+' or '-') for my purposes?
Ultimately my goal is to get a bunch of RNAs with their corresponding types (miRNA, lincRNA, snRNA... etc) to do some computations on.
My question is this: what is the easiest way to get it out of the data I have?
The data I am using is freely available here: http://wanglab.pcbi.upenn.edu/coral/ under the heading "Annotation packages" if anyone is interested or needs specifics.
Thank you!
•
written
by
MarkR •
1 answer
Perhaps my answer to to A: Extract cds fastas from a gff annotation + reference sequence will serve your purposes....???
hmmm - your Q said nothing about "reads".... anyway, good luck with the f# ... you might want to reword you q to indicate this is required.
Log in to answer this question.
-
Rows of read count files
written by Abhishek •
-
Full GTF file vs Subset GTF file
written by Arindam Ghosh
-
GTF file in Ensembl database for RNA-Seq analysis
written by wangdp123
-
Retrive orthologs for coral housekeeping genes
written by chetana •
-
non-coding RNAs target predictions
written by Leite
-
Structure prediction on DNA folding form or RNA folding form
written by The Last Word
-
Linc and miRNA in GTF file
written by Sandeep
-
Extracting regions from .bam file using a .gff or .gtf file
written by Tom A •
-
Mutual Information
written by ftp
-
Exonerate Gff In Artemis
written by Daniel Fernández •
