Sorry to bother you with this, but I am trying to process the KB1 Bushman data, through a pipeline that starts with samtools mpileup -C50 -uf ~/Bushman/hg18/chromFa/chr1.fa ../KB1illumChr1.bam ../KB1454chr1.bam | bcftools view -c - >Chrm1_illum_454.vcf
Processing the data one cromosome at a time, aligning the data against hg18, using mpileup to combine the illumina and 454 data. All goes very well, everything looks good all the way through, results are stable, hi-quality, and conceivably even meaningful, except for KB1454chr5.bam, the 454 data for Chrm5 When combined with the illumina data for Chrm5 it results in garbage, even though if I run the illumina data alone, it looks like all the other Chromosomes. I confess that I am worse than a Newbie, I'm an old guy working on my own, so may in fact just be in way over my head. But if there is any chance you could take a look at KB1454chr5.bam to see if there were a lot of bad reads, or an error was made in processing, or anything else is wrong, I'd be most grateful.
___________________________________________________________ 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/