Hi Nicola,

I'd recommend you look at our suite of tools, called the GATK, for working with NGS data. It includes a very useful test of tools for comparing and evaluating SNP call quality:

http://www.broadinstitute.org/gsa/wiki/index.php/VariantEval

Best,

On Jul 9, 2010, at 4:45 AM, Nicola Nadeau wrote:

Hi,

I am trying to find SNPs and/or indel variants that differ between two groups of samples. The data are 454 sequence capture results spanning a region of interest for a non-model organism without a complete genome. I have mapped these to my reference sequence of the region and generated BAM and pileup files for them. Does anyone know of a tool or method (either in galaxy or elsewhere) that will allow me to compare datasets and pick out variation between them (as opposed to things that differ from the reference sequence).

Many thanks,

Nicola

---

galaxy-user mailing list galaxy-user@lists.bx.psu.edu http://lists.bx.psu.edu/listinfo/galaxy-user

Mark A. DePristo, Ph.D. Manager, Medical and Population Genetics Analysis Broad Institute of Harvard and MIT depristo@broadinstitute.org mark@depristo.com