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
Running preseq from the latest release preseq-3.1.2.tar.gz
lc_extrap give the following error:
However, converting the .bam to .bed first and then running lc_extrap on the .bed file completes as expected:
Any idea what's going wrong?
Example .bam and .bed files attached.
bams.zip