preseq icon indicating copy to clipboard operation
preseq copied to clipboard

ERROR: max count before zero is less than min required count (4) duplicates removed

Open didillysquat opened this issue 3 years ago • 2 comments

Running preseq from the latest release preseq-3.1.2.tar.gz

lc_extrap give the following error:

Singularity> preseq lc_extrap -P scaf1.1m.bam
ERROR:	max count before zero is less than min required count (4) duplicates removed

However, converting the .bam to .bed first and then running lc_extrap on the .bed file completes as expected:

Singularity> bedtools bamtobed -i scaf1.1m.bam > scaf1.1m.bed
Singularity> preseq lc_extrap -P scaf1.1m.bed > lc_extrap.out

Any idea what's going wrong?

Example .bam and .bed files attached. bams.zip

didillysquat avatar Feb 15 '22 08:02 didillysquat