seismic cleanfa
Purpose
seismic cleanfa fixes FASTA files that are incompatible with SEISMIC-RNA.
Common issues from databases like NCBI and Gencode include metadata in header
lines and non-standard characters in sequences.
It trims headers to valid path characters, uppercases sequence bases, replaces
ambiguous IUPAC codes with N, and converts U``↔``T as needed.
Inputs
- FASTA files
One or more FASTA files to clean.
Outputs
Cleaned FASTA files with the same names as the inputs, written to the output
directory (default ./out).
Quick example
Clean a reference FASTA:
seismic cleanfa refs.fa
Output goes to out/refs.fa.
Specify a different directory with -o:
seismic cleanfa -o clean refs.fa
Options
--inplace/--no-inplaceOverwrite the original files instead of writing to
--out-dir(default off).Warning
--inplaceis irreversible — back up your files first.--out-dir DIR(-o)Write cleaned files here (default
./out). Ignored when--inplaceis set.
--num-cpus N — multiprocessing; see Parallelize Tasks.
--force — overwrite existing outputs.
The auto-generated Command Line Reference lists every option with its current default.
Common unexpected results
- Line contains sequence characters other than whitespace or IUPAC codes
The sequence has characters that cannot be interpreted automatically. Remove them manually and rerun. The error message names the affected line(s).
See also
FASTA: Reference sequences — the FASTA format required by SEISMIC-RNA
seismic align — uses the cleaned FASTA as input