Skip to content

How do I skip demultiplexing while running curator.py!! #48

Description

@mrqasimasif

Since on long read data, demultiplexing takes a long time to process, and if the process fails during SAM or BAM generation for one FASTQ file, it has to run on the entire processed.fastq.gz to demultiplex and perform minimap alignment.
This is highly inefficient since the cell barcode demultiplexing is a time-intensive task.

Also, in curator_io.py the curation_master() function fails to handle the SAM file generation if --inc_contig is not used, since the code does not account for the header lines.

if options.inc_contig:
SAM.write(out_msg.decode("utf-8"))
#---filter out non-autosome---
else: ---------------------------------------> # changed it with the code below
lines = out_msg.decode("utf-8").split("\n")
for line in lines:
if line == "":
break
if line.split("\t")[2].startswith("chr"):
SAM.write(line + "\n")

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