Skip to content

secapr clean_reads and secapr assemble_reads stall #41

@faustovrz

Description

@faustovrz

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

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions