Question: Extract 3'utr from annotation file
0
gravatar for najibveto
4.1 years ago by
najibveto40
Korea, Republic Of
najibveto40 wrote:

hello

i am looking to extract the 3'utr sequences of fathead minnow fish using the annotation file and the genome. http://www.ncbi.nlm.nih.gov/genome/?term=fathead%20minnow

can someone help me or advice me on the method to do it?

genome • 2.1k views
ADD COMMENTlink modified 4.1 years ago by Devon Ryan94k • written 4.1 years ago by najibveto40
1

If UTR features are annotated in your annotation file, you can use bedtools getfasta utility to extract the sequences. You must provide the genome file and a BED/GFF/VCF file of ranges to extract (3'UTR in your case).

If UTR are not annotated yet, you'd need to annotate them using some ORF prediction tool ...

ADD REPLYlink modified 4.1 years ago • written 4.1 years ago by iraun3.7k
1
gravatar for Devon Ryan
4.1 years ago by
Devon Ryan94k
Freiburg, Germany
Devon Ryan94k wrote:

If the UTRs are unannotated, you might have luck with the GenomicFeatures package in R. UTRs would be the difference between exon and CDS entries per transcript (there are some functions that enable doing that, read through the vignette). If you're really luck then the txdb object that you create will already do this, but I've never tried.

The you can either write those to a file and use bedtools getfasta or use an equivalent method in R.

ADD COMMENTlink modified 4.1 years ago • written 4.1 years ago by Devon Ryan94k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1014 users visited in the last hour