lmdu/pyfastx

How to read two fastq files simultaneously?

Opened this issue · 1 comments

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?

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