Skip to main content

parse_sequences

Function parse_sequences 

Source
pub fn parse_sequences<P, F>(path: P, callback: F) -> Result<()>
where P: AsRef<Path>, F: FnMut(&[u8], &[u8]) -> 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