If there are too many reads, the PMDtools rule gives this error:
Traceback (most recent call last):
File "/proj/nobackup/metagenomics/pochonz/scripts/pmdtools.0.60.py", line 944, in <m
odule>
L_D = L_D * L_match(z,ancient_model_deam,qualsrev,options.polymorphismancient)
File "/proj/nobackup/metagenomics/pochonz/scripts/pmdtools.0.60.py", line 148, in L
match
P_error= phred2prob( (ord(fquals[fposition])-33))/3.0
KeyboardInterrupt
But this can be fixed by limiting the number of reads used to e.g. 100000 using the option --number for PMDtools.
If there are too many reads, the PMDtools rule gives this error: Traceback (most recent call last): File "/proj/nobackup/metagenomics/pochonz/scripts/pmdtools.0.60.py", line 944, in <m odule> L_D = L_D * L_match(z,ancient_model_deam,qualsrev,options.polymorphismancient) File "/proj/nobackup/metagenomics/pochonz/scripts/pmdtools.0.60.py", line 148, in L match P_error= phred2prob( (ord(fquals[fposition])-33))/3.0 KeyboardInterrupt
But this can be fixed by limiting the number of reads used to e.g. 100000 using the option --number for PMDtools.