similar bams with rescale and non rescale about paleomix HOT 2 CLOSED

Hi,

mapDamage performs rescaling by decreasing the base quality scores of C>T/G>A substitutions based on the estimated probability that they are the result of DNA damage. No other changes are made and the coverage should therefore be identical before and after rescaling. However, if you don't have much of a damage signal, then rescaling probably won't make much of a difference to your base qualities either.

And take care if you are running with/without rescaling using the same YAML file and in the same directory: The current version of paleomix does not always detect changes in parameters like that and BAM file(s) might not be updated if you just change the YAML file and re-run paleomix in the same directory. In this particular case, it should be enough to (re)move the final BAM files generated by the pipeline after you've tweaked the mapDamage feature.

Best,
Mikkel

from paleomix.

Hi Mikkel,

Thank you very much for the quick reply. I have also check intersection and complement in SNPs for both vcf files with bcftools stats - virtually all SNPs are shared between rescaled and non rescale indicating that nothing was filtered out while generating and filtering the vcfs.

Also run it in separate folders to be sure.

With all best wishes,

Flávia

from paleomix.