This is a test version of Biostars. For the public version, visit https://www.biostars.org.
fasta generation for functional elements

Hi,

I have got several suggestion how to extract fasta sequence from assembled cufflinks like to use gff read program or extract genomic DNA. Both of these programs help in generating fasta sequence with respect to the reference annotation file but I want to extract fasta sequence from mapped reads so that I have a better picture to find the functional elements present in my samples.This will also help to identify novel genes and transcripts.

Is there a better way to extract fasta format?

Thanks

NISHA

mapped-reads fasta-generation
• 2,496 views • link
updated by Ram Ram • written by ajmaninisha • ajmaninisha

2 answers

I answered a question similar to this, please check my comment in this C: Tool to retrieve regions with reads?

• 0 views • link
written by James Ashmore James Ashmore

Thanks james and is it possible to get fasta from assembled cufflinks transcript file too.

I am using galaxy so the tools suggested by you are available in the tool shed.

Thanks

• 0 views • link
updated by Ram Ram • written by ajmaninisha • ajmaninisha

The bedtools getfasta command requires your query file to be in .bed format, so you could convert the transcript .gtf file into a .bed file and use that, again check the awk command I listed on the other thread to see how to do that, hope that helps

• 0 views • link
written by James Ashmore James Ashmore

Log in to answer this question.

More posts like this