Question: Extract assembled transcripts from Cufflinks gtf without to use reference genome
0
gravatar for marija.durdevic
2.2 years ago by
marija.durdevic0 wrote:

Hi all,

I know there were similar questions, but unfortunately I couldn't find proper answer for my problem. I am doing de novo transcriptome assembling of unknown fish. I used Trinity for that purpose and I acceptable assemblies. I also did genome-guided assembling with STAR->Cufflinks using close related organism O.niloticus. My idea is to do augmentation of trinity assemblies with results from Cufflinks (guided by this article) and annotate it afterwords.

The problem is to extract generated Cufflinks transcripts into fasta format. I read about several options, but all of them use reference genome. I cannot use it in my case since it's close related species and I have a lot of SNPs.

Is there any method to extract fasta file from alignment only (from bam files) or from reference but to change all inappropriate nucleotides?

Thanks in advance, Marija

bam rna-seq bcftools cufflinks gtf • 702 views
ADD COMMENTlink modified 14 months ago by Biostar ♦♦ 20 • written 2.2 years ago by marija.durdevic0
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: 1771 users visited in the last hour