I'm doing de novo assembly on RNA-seq data and I do a stand-alone BLAST against a reference transcriptome. I want to know how much of my assembly aligned, preferably as a percentage so I can estimate quality and have simple statistics to present along with the usual average and median contig lengths-
Velvet-Oases 6.3%
Trinity 7.9%
etc. Is there a tool out there that will estimate how much of a reference database I have covered? Essentially I could add up the lengths of each aligned transcript and divide by the length of all the reference transcripts....not sure if this makes sense, but I could use some advice.
Velvet-Oases 6.3%
Trinity 7.9%
etc. Is there a tool out there that will estimate how much of a reference database I have covered? Essentially I could add up the lengths of each aligned transcript and divide by the length of all the reference transcripts....not sure if this makes sense, but I could use some advice.
Comment