Hello Shashi, You have some choices:
Using SAM Tools on the command line, samtools pileup -c will call variants from the SAM/BAM alignments with BED output. This could be loaded at a custom track. This may be the quickest method if your dataset is very large. Or, you can use Galaxy and load the results as a custom track. Send your existing BAM custom track data from the Table browser (if not too large) or upload files directly, then after the analysis, view results in the UCSC browser. At Galaxy, see "NGS: SAM Tools -> Generate pileup & Filter pileup". Links from UCSC: http://genome.ucsc.edu/ -> Galaxy (left blue bar) from Table browser: Check "Galaxy" box in output section Or directly: http://galaxy.psu.edu/ Hopefully this helps, Jennifer --------------------------------- Jennifer Jackson UCSC Genome Informatics Group http://genome.ucsc.edu/ On 6/2/10 7:03 AM, Shashikant Pujar wrote: > Hi > > I have loaded a custom track (Illumina NGS reads in BAM format of a 2Mb > region) on the UCSC Dog Genome Browser. Is there a way I can extract all > SNPs and Indels between the custom track and canFam2? > > Thanks > > Shashi Pujar > _______________________________________________ > Genome maillist - [email protected] > https://lists.soe.ucsc.edu/mailman/listinfo/genome _______________________________________________ Genome maillist - [email protected] https://lists.soe.ucsc.edu/mailman/listinfo/genome
