I am interested in Reference based transcriptome assembly for my species of interest for which the reference sequence is available from a different accession. For reference based approach i used tophat/cufflinks/gffread to get the final transcript file. Now i am wondering how can i edit the transcript file based on RNAseq information. Basically i don't want the sequence from the reference and all i want is the transcript sequence information based on RNASeq reads.
Any ideas of how to do that?
Thanks
•
written
by
upendrakumar.devisetty
1 answer
If you need a fully de-novo transcriptome assembly that does not incorporate other information then you should create it as such. You cannot later remove information used in the assembly, nor can you reverse engineer to undo the effects of including the reference.
Log in to answer this question.
-
Problem trying to convert from .gtf to .fasta file
written by nattzy94 •
-
Studying variances in assembled transcripts
written by rmihael •
-
Different annotation gtf file in Tophat2 and Cufflink
written by aditi.qamra
-
how to compare the Ka/Ks ratio for chromosomes
written by upendrakumar.devisetty
-
Cuffcompare: Some Cufflinks Transcripts Are Not Contained By Cuffcompare Output Tmap File
written by zju.whw •
-
How To Extract Cds And Protein Sequences From Cufflinks Transcripts.Gtf File?
written by Rahul Sharma
-
Tuxedo Suite: When To Skip Cufflinks/Cuffmerge?
written by jobinv
-
How To Merge Gff Files
written by upendrakumar.devisetty
-
How To Distinguish Chimeric Transcript From Fusion Gene
written by upendrakumar.devisetty
-
Fused Or Not Fused Gene
written by upendrakumar.devisetty
You want to extract transcript sequences from the reference genome based on assembly annotations generated by tophat/cufflinks?
No i have already done that. Now i have final transcript file but the problem is the sequence of the transcript file is from the reference and not from my RNAseq reads which i don't want because it is from a different accession.