-
Notifications
You must be signed in to change notification settings - Fork 12
Open
Description
Hello
I've tried to run secapr in a sequence capture library, and I found that the tutorial instructions aren't working for my data.
both clean_reads and assemble_reads start and stall without giving an error.
I know they stall because 1) I can get no I/O nor significant memory usage from those processes and 2) by using fastp and spades.py I can clean and I can do de novo assembly in minutes on the same data.
It's been 2 years since the last commit so probably there is no longer support for this software, but maybe somebody has some suggestion.
This is an example for clean_reads
secapr clean_reads \
--input ../testin \
--read_min 190000 \
--sample_annotation_file sample_annotation.csv \
--output testout/clean
/share/maize/frodrig4/conda/env/secapr_env/lib/python3.8/site-packages/Bio/Application/__init__.py:40: BiopythonDeprecationWarning: The Bio.Application modules and modules relying on it have been deprecated.
Due to the on going maintenance burden of keeping command line application
wrappers up to date, we have decided to deprecate and eventually remove these
modules.
We instead now recommend building your command line and invoking it directly
with the subprocess module.
warnings.warn(
[Info:] Samples with a total read-count of less than 190000 (forward + reverse reads) are not being processed. If required you can set a different threshold, using the --read_min flag.
Zx0540_P3_P5_P1111_S27: Counting all reads (forward + reverse) belonging to this sample...
221222
Zx0550_P4_P3_P5311_S20: Counting all reads (forward + reverse) belonging to this sample...
196898
Zx0580_P2_P5_P2411_S29: Counting all reads (forward + reverse) belonging to this sample...
236377
##################################################
Processing Zx0540_P3_P5_P1111_S27...
Metadata
Metadata
Assignees
Labels
No labels