nf-core / eager

A fully reproducible and state-of-the-art ancient DNA analysis pipeline
https://nf-co.re/eager
MIT License
142 stars 80 forks source link

Documentation for qualimap doesn't account for large genomes #578

Closed TCLamnidis closed 3 years ago

TCLamnidis commented 3 years ago

https://github.com/nf-core/eager/blob/dev/docs/output.md#qualimap

The expected coverage histograms and cumulative coverage histograms for larger genomes (e.g. human) would include a majority of 0x coverage positions, with some positions having 1x and fewer having 2x etc. (See reference image below)

This will be the case for the majority of human aDNA data, since we can rarely produce 30x whole genome sequencing data. The docs should be updated to mention/acknowledge this possibility.

Screenshot 2020-10-16 at 15 59 40

The above is results from human aDNA.

jfy133 commented 3 years ago