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

AVITI manifest improvements #374

Merged
merged 9 commits into from
Nov 4, 2024

Conversation

kedhammar
Copy link

For AVITI manifest generation: make PhiX manifest variant, fix udf typo, remove unused func, clarify var names

@chuan-wang
Copy link
Member

chuan-wang commented Oct 28, 2024

Hi @kedhammar, there is one more thing for us to fix. For the following cases index 2 need to be converted to the reverse complementary form:

(1) All inhouse-prepared libraries
(2) User submitted libraries with special indexes that match with TENX_DUAL_PAT or SMARTSEQ_PAT.

Please refer to the codes for MiSeq as example.

Copy link
Member

@chuan-wang chuan-wang left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Changes required.

@kedhammar
Copy link
Author

kedhammar commented Nov 4, 2024

Changes required.

How's this @chuan-wang? Not tested.

2bc45ad

@chuan-wang
Copy link
Member

``> > Changes required.

How's this @chuan-wang? Not tested.

2bc45ad

Hi, the logic should be the opposite:

# Special cases to reverse-complement index2
if not user_library or (
    user_library
    and (
        TENX_DUAL_PAT.findall(lims_label)
        or SMARTSEQ_PAT.findall(lims_label)
    )
):
    logging.info(f"Reverse-complementing index2 of {sample.name}.")
    row["Index2"] = revcomp(row["Index2"])

@kedhammar
Copy link
Author

``> > Changes required.

How's this @chuan-wang? Not tested.
2bc45ad

Hi, the logic should be the opposite:

# Special cases to reverse-complement index2
if not user_library or (
    user_library
    and (
        TENX_DUAL_PAT.findall(lims_label)
        or SMARTSEQ_PAT.findall(lims_label)
    )
):
    logging.info(f"Reverse-complementing index2 of {sample.name}.")
    row["Index2"] = revcomp(row["Index2"])

Oop, that's right 👀

ab77d5e

@chuan-wang chuan-wang merged commit 9497e52 into SciLifeLab:master Nov 4, 2024
7 checks passed
@kedhammar kedhammar deleted the aviti-manifest-dev branch November 4, 2024 15:05
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants