Extract assembled transcripts from Cufflinks gtf without to use reference genome
Entering edit mode
3.9 years ago

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

RNA-Seq cufflinks gtf bcftools bam • 1.0k views

Login before adding your answer.

Traffic: 1738 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6