-
Notifications
You must be signed in to change notification settings - Fork 17
Description
Dear author
There are some problems in calculating library complexity using preseq lc_extrap. Versions over 3.0 do not have -B parameters. Using this version(preseq lc_extrap -e 10000 -P -D -v -o A1.test.dat SRR13810470-f2.q30.bam 2>A1.test.log) will encounter an error (ERROR: max count before zero is less than min required count (4) duplicates removed) in A1.test.log file and without generating A1.test.dat file. However, Versions below 2.0 have -B parameters. Using this version(preseq lc_extrap -e 10000 -P -B -D -v -o A1.test.dat SRR13810470-f2.q30.bam 2>A1.test.log) will generate two files without any errors.
I want to know why versions over 3.0 have this error.
Error as following
PAIRED_END_BED_INPUT
TOTAL READS = 3213007
DISTINCT READS = 3.21293e+06
DISTINCT COUNTS = 3
MAX COUNT = 3
COUNTS OF 1 = 3.21286e+06
MAX TERMS = 2
OBSERVED COUNTS (4)
1 3212858
2 73
3 1
ERROR: max count before zero is less than min required count (4) duplicates removed