Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Paired read detection not working for ReadQC reports from megsap pipeline #1

Open
caspargross opened this issue Jan 17, 2023 · 0 comments

Comments

@caspargross
Copy link
Member

Description of bug

Field for pair-end is determined based on source file names in .qcML file:

def check_paired_end(qcml_contents):
"""Check if both R1 and R2 are present in input files."""
found_r1 = False
found_r2 = False
root = xml.etree.cElementTree.fromstring(qcml_contents)
for el in root.findall(".//{http://www.prime-xs.eu/ms/qcml}metaDataParameter"):
if el.attrib["name"] == "source file":
if "R1" in el.attrib["value"]:
found_r1 = True
if "R2" in el.attrib["value"]:
found_r2 = True
return found_r1 and found_r2

ReadQC Files from RNA-Seq pipeline dont have correct source file names:

    <metaDataParameter ID="md0004" name="source file" value="analyze_rna_nMIt6T_bc1.fastq.gz" cvRef="QC" accession="QC:1000005"/>
    <metaDataParameter ID="md0005" name="source file" value="analyze_rna_b9JIQq_bc2.fastq.gz" cvRef="QC" accession="QC:1000005"/>

File that triggers the error

No response

MultiQC Error log

No response

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant