Skip to content

Issue filtering host with >2 read files #11

Description

@LanderDC

Hi,

When I try to filter my reads (3 files) against my host genome without the quality filter (has already been performed beforehand), I get an error from minimap2:

[ERROR] incorrect input: in the sr mode, please specify no more than two query files.

I believe the issue is on these lines:

minimap2 -t $CPUS -ax sr ${FILTER_LOC}/filter_seqs.fna ${READS} | \
samtools fastq -n -f 4 - > ${TEMP_DIR}/${SAMPLE}.EV_input.fastq

and can be fixed by:

cat ${READS} | minimap2 -t $CPUS -ax sr ${FILTER_LOC}/filter_seqs.fna - | \
samtools fastq -n -f 4 - > ${TEMP_DIR}/${SAMPLE}.EV_input.fastq

Thanks!

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Fields

    No fields configured for issues without a type.

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions