pub fn parse_sequences<P, F>(path: P, callback: F) -> Result<()>Expand description
Parse a FASTA/FASTQ file and call a function for each valid DNA sequence
§Arguments
path- Path to input file (may be gzipped)callback- Function called for each sequence, receives (name, sequence)
§Errors
Returns error if:
- File cannot be opened
- File format is invalid
- Sequence contains non-DNA characters