454 - Transcriptome SNPs Detection

[Giorgio C]

Hy everybody,
I have assembled a Transcriptome with Newbler and than blasted all the isotig.fna. I would estimate the number of SNPs and their frequency but manually is impossible. Is there any automate tool to do this? Can anyone help me?
Thank you very much

[westerman]

I suggest that you use the Newbler 'mapper' instead of the assembler. This will automatically provide variants.

[Giorgio C]

Thank you, but i need to map with "EST database" or with for exemple "Ensemble trasncriptome 37"? Where can i find the frequency in the output files or the statistic number of the work (total number of SNP..etc) ???

[Kaurh5]

Hi all!

I am trying to assemble a low coverage 454 data of a plant using Newbler. I have two raw sff files from two different genotypes of my experimental plant. newbler completes the assembly step without a considerable error for the individual sffs. But when I try to assemble the sff files of both genotypes together(using incremental denovo assembly) it just adds up the total contigs and the singletons for that matter neglecting the possible common contigs between the two genotypes. To my understanding newbler is treating every read in both the sff files as unique which is very unlikely to happen. My basic aim is to find the SNPs and repeats in the genome and if newbler is assembling every read into a unique contig then this could be a matter of concern to me. Please provide the necessary explaination for this behaviour.