Can we use the SAM tools pileup tool in Galaxy to get an accurate count of the coverage across the genome? I ask because I have run the pileup tool in Galaxy with bam files generated from tophat and noticed that the coverage reported in the pileup data was consistently lower then what I see in IGV. For example, I look at chr_1, coordinate 3515 and I see 10 reads mapped to that coordinate. I look in the pileup and it reads 1x coverage. I have not applied any filter to the pileup because all I want is a readout of the coverage at each NT including the number of variants that occur. My settings for generate pileup were: do not print mapping quality, print all lines, cap mapping quality to- 0, call consensus according to MAQ model- NO. Any suggestions? Is there a tool in Galaxy other than pileup that can do this? Thanks, Patrick
___________________________________________________________ The Galaxy User list should be used for the discussion of Galaxy analysis and other features on the public server at usegalaxy.org. Please keep all replies on the list by using "reply all" in your mail client. For discussion of local Galaxy instances and the Galaxy source code, please use the Galaxy Development list: http://lists.bx.psu.edu/listinfo/galaxy-dev To manage your subscriptions to this and other Galaxy lists, please use the interface at: http://lists.bx.psu.edu/