I have de novo assembled a consensus ("reference") sequence for a non-model organism using RNAseq from Illumina and I want to annotate the contigs with CDS information. I have BLASTx'd everything against NCBI's protein data base (translating everything frist, of course) and now I would like to annotate the consensus contigs so that I can do a proper SNP analysis (i.e. I want to know which SNPs are synonymous/not, but I can't do this unless CLC knows where the coding regions are - and presumably this information comes from the BLAST results, but how do I put these together??)
Does anyone know how to do this? Or am I missing something?
Does anyone know how to do this? Or am I missing something?
Comment