Reverse-Complement FASTA in Seqqio: Orientation, Ambiguous Bases, and Batch Exports
Transform DNA orientation without losing record identity. Worked FASTA examples explain the three operations, ambiguity codes, U-to-T normalization and export checks.
Changing orientation is simple arithmetic but an easy place to lose experimental meaning. Reversing the characters is not the same operation as finding the opposite strand. The examples below make the operation explicit and show how ambiguous bases and an invalid record behave in a batch.
Distinguish the three operations with four bases
Take the input 5′-AGTC-3′. Complementing it replaces each base while preserving its position; reversing it changes order without replacing bases. Doing both gives the opposite strand in 5′-to-3′ orientation.
| Operation | Output string | How to read it |
|---|---|---|
| Reverse | CTGA | Reversed character order; no base substitution |
| Complement | TCAG | Complement aligned to the input: 3′-TCAG-5′ |
| Reverse complement | GACT | Complementary strand written 5′-GACT-3′ |
The plain sequence string does not carry visible direction labels by itself. Keep the operation in the saved run or surrounding metadata. Do not call every transformed string a 5′-to-3′ opposite strand: that description applies to the reverse complement of the stated input orientation.
Run a batch that includes ambiguity and an error
Open Reverse Complement, choose batch FASTA input and select the reverse-complement operation. Use these four short control records. The third deliberately contains U, and the fourth contains an unsupported question mark.
>plain
AGTC
>ambiguous
ACGTRYMKWSBVDHN
>rna_input
AUGC
>invalid
ACG?
| Record | Result | Status |
|---|---|---|
| plain | GACT | Success |
| ambiguous | NDHBVSWMKRYACGT | Success |
| rna_input | GCAT | Success; U normalized to T |
| invalid | No output sequence | unsupported_symbol |
The records are transformed separately. Their boundaries must not be removed by concatenating every sequence line into one molecule before the operation. In a larger dataset, check that successful records plus failed records account for the original batch, and retain the failure report alongside the exported valid FASTA.
Understand what ambiguous-base complements preserve
An ambiguity code represents a set of possible bases. Complementing it transforms that set rather than selecting one base. R, meaning A or G, complements to Y, meaning T or C. M and K exchange; B and V exchange; D and H exchange. W, S and N remain their own complements.
For the second control, reversing and complementing every position gives NDHBVSWMKRYACGT. The ambiguity remains unresolved. A transformed N is still unknown; it is not evidence that the underlying position has been determined.
Biopython’s DNA complement documentation describes the same base transformations and the distinction from RNA-specific operations. We used Biopython 1.87 as an independent check of the valid outputs. Its handling of arbitrary non-base characters is not a substitute for Seqqio’s own input validation.
Treat U-to-T normalization as part of the method
In this Seqqio tool, AUGC is normalized to the DNA string ATGC before transformation. Its reverse complement is therefore GCAT. This is DNA output, not an RNA reverse complement containing U.
Record the normalization when the original source used RNA letters. Do not assume that every Seqqio analysis accepts or normalizes U identically: input rules belong to the selected tool. If the next step requires RNA, choose an explicitly appropriate conversion or RNA workflow and verify its alphabet.
Check the transformation before a downstream analysis
For a valid normalized DNA record, applying reverse complement twice returns the normalized original sequence. Length is preserved. These are useful sanity checks, but they do not determine whether the original record was supplied in the biologically relevant orientation.
After transformation, any positional annotations also need review. For a sequence of length L, an original one-based inclusive interval a–b maps to L−b+1 through L−a+1 in the reverse-complement coordinates. The coordinate arithmetic does not automatically update a separate annotation file. Keep that step explicit if you carry features into another tool.
If your goal is to find where a primer binds, use a mapping workflow rather than assuming that reverse-complementing a primer locates a binding site. The Seqqio primer-mapping guide covers that separate task.
Export sequences and the decisions that produced them
Save the output FASTA together with the original input, operation, normalization note, Seqqio version and complete report. Stable record IDs make it possible to join transformed sequences back to sample or construct metadata. Check for repeated identifiers before using an ID alone as a database key.
Seqqio fits a Windows user who wants local transformations, record-level results and retained exports in the same workspace as other DNA analyses. A short script may be preferable inside an automated pipeline; the desktop workflow is useful when you want to inspect and preserve individual analysis runs.
References
- Biopython contributors. Bio.Seq module Biopython 1.87 documentation Sequence complement and reverse-complement operations.