Skip to content

Don't use random FASTQ access during hcid check - #115

Merged
BioWilko merged 2 commits into
mainfrom
hcid_bugfixc
May 18, 2026
Merged

Don't use random FASTQ access during hcid check#115
BioWilko merged 2 commits into
mainfrom
hcid_bugfixc

Conversation

@BioWilko

Copy link
Copy Markdown
Member

No description provided.

@BioWilko
BioWilko requested a review from Sam-Sims May 18, 2026 11:04

@Sam-Sims Sam-Sims left a comment

Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

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

minor, non-blocking comment - otherwise LGTM

Comment thread bin/check_hcid.py
for name in hcid_dict[taxid]["mapped_read_ids"]
}
read_records = {}
for name, seq, qual in pyfastx.Fastq(read_file, full_name=False, build_index=False):

Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

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

very minor but could consider an early return here once all needed reads are found e.g all hcid reads were in first 10% of a 5gb fastq file then no need to continue parsing

as I type this though probably not worth it in our use case - so not blocking

Copy link
Copy Markdown
Member Author

Choose a reason for hiding this comment

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

Yeah that's a good idea, let me have a look

@BioWilko
BioWilko merged commit 6a63bec into main May 18, 2026
1 check passed
@BioWilko
BioWilko deleted the hcid_bugfixc branch May 18, 2026 12:14
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

Successfully merging this pull request may close these issues.

2 participants