Generate reverse-complemented FASTQ reads while preserving quality scores
Reverse Complement
Generates the reverse complement of each read in a FASTQ file. Sequences are reversed and complemented (A↔T, C↔G), while quality score strings are reversed to maintain base-quality correspondence.
.fastq,
.fq)
.fastq.gz,
.fq.gz)
Original
@read1
ATCG
+
FFFF
↓
@read1
CGAT
+
FFFF
seqkit seq -r -p \
input.fastq.gz \
-o revcomp.fastq.gz
fastx_reverse_complement \
-i input.fastq \
-o revcomp.fastq
from Bio import SeqIO
with open("revcomp.fastq", "w") as out:
for record in SeqIO.parse(
"input.fastq",
"fastq"):
SeqIO.write(
record.reverse_complement(
id=record.id,
description=""
),
out,
"fastq"
)