Update: adaptation of a default run bam->pileup on my small test dataset gives the same result in my case using mpileup instead of pileup in the sam_pileup.py script (line 97). Not tested any additional parameters though so there might be subtle differences in the outcome.
#prepare basic pileup command cmd = 'samtools mpileup %s -f %s %s > %s' try:
Van: galaxy-dev-bounces@lists.bx.psu.edu [mailto:galaxy-dev-bounces@lists.bx.psu.edu] Namens Bossers, Alex Verzonden: vrijdag 20 juli 2012 12:20 Aan: 'galaxy-dev@lists.bx.psu.edu' Onderwerp: [galaxy-dev] samtools pileup deprecated
Hi All,
Maybe I missed it on the list here, but it seems that for samtools the pileup command has been deprecated without backward compatibility from the more recent samtools packages... http://massgenomics.org/2012/03/5-things-to-know-about-samtools-mpileup.html
Did anyone rewrite some of the samtools tool wrappers for this....? I am currently testing whether the Bam->pileup command can be easily changed in the sam_pileup.py script. But some functionality has gone and some additional is there. It seems I am up-to-date with the latest wrappers in the galaxy_dist repository.
Any pointers to recent wrappers would be welcomed.
Alex