In rho_est.nf, subsample process, samtools pileup is used to determine the maximum while this method seems to work well, perhaps samtools coverage or a different tool that is designed for this task might be better? Furthermore by default samtools pileup has a default max depth of 8000.
In rho_est.nf, subsample process, samtools pileup is used to determine the maximum while this method seems to work well, perhaps samtools coverage or a different tool that is designed for this task might be better? Furthermore by default samtools pileup has a default max depth of 8000.