How to read two fastq files simultaneously?
Opened this issue · 1 comments
yc-pro commented
A single fastq file can be read as follows:
example:for name,seq,qual in pyfastx.Fastq('XXX.fq.gz')
But,If the sequencing data consists of a pair of PE files, how can they be read simultaneously?
mr-eyes commented
x = pyfastx.Fastq('XXX_R1.fq.gz')
y = pyfastx.Fastq('XXX_R2.fq.gz')
# Iterate over both sequences simultaneously
for (name1, seq1), (name2, seq2) in zip(x, y):
# ...
pass