2.3 years ago by
Duarte, CA
If they are closely related genomes, it might be worth considering an RNA-Seq assembly from a related genomic alignment (using an algorithm like cufflinks)? If so, it would probably be good to have multiple samples (where I think the cuffmerge assembly is probably better than any of the individual sample assemblies).
If you have genomic sequence, you could also try using a program like MAKER (where you can provide RNA-Seq data in the annotation process). However, that may take a while for a whole genome sequence (particularly if it is a vertebrate genome).
Otherwise, I sometimes use Oases (or maybe even Velvet contigs) for RNA-Seq de novo assembly, but I would usually be assembling partial transcript sequences (most likely, for a fraction of the total transcripts in any sample). In other words, I would have some concerns about using that assembly for quantification. However, if there is a set of RefSeq sequences for your organism (I would guess from ESTs, and other sources), that may be the best option for transcript quantification.
Hi, Similar question here
Thank you. I see the post, in your case, you have one genome that is of low quality. This is different from mine, the two genome are both of high quality, and they are similar. If I align the two sample reads to one reference genome, the reference genome itself is some different from the two genome. Is is possible to do something more precisely?
I am not an expert, I don't want to tell you a wrong answer, but in your case maybe I would check the genomes similarity by mapping reads from species 2 to species 1 (or/and species 1 to species 2) and see the overall alignment rate. And if it is good enough for what you want to do, use this genome as a reference (this is still a non expert opinion).
OK, thank you for your advice.
Depending on what species you're working on, you could find ortholog information in ensembl. I think that would be a much better approach than mapping the RNA-seq of the two species to the genome of one of them.
All right, thank you~