Is it possible using the out of the box samtools to list what reads (read IDs) were used in determining a specific SNP call? Poking around it seems like it might be possible using the BioPerl or pysam libraries but was wondering if the samtools binary could maybe do this when using pileup or mpileup.
Thanks!
Thanks!