In the package Stacks, the -m parameter specifies the read depth, i.e. the minimal number of reads required to identify and accept a locus. I'm using the denovo_map.pl script and subsequently the populations script to output an input file for Structure.

With respect to denovo, the manual (p. 26) says:

m — specify a minimum number of identical, raw reads required to create a stack.

Now, in the 'populations' script, the manual (p.24) says

m — specify a minimum stack depth required for individuals at a locus.

Now I wonder, if -m in 'denovo' sets the min nr of read per allele, and -m in 'populations' sets it per individual, should not this parameter in populations have at leats 2x its value in denovo? Otherwise, the filtering would favour homozygotes over heterozygotes.

Is my reasoning correct?

Similar questions and discussions