Skip to main content

write_gff_streaming

Function write_gff_streaming 

Source
pub fn write_gff_streaming<P>(
    path: &Path,
    feature: &str,
    produce: P,
) -> Result<(), Box<dyn Error>>
where P: FnOnce(&mut dyn FnMut(&str, &[u8], &[NucleosomeCall], usize) -> Result<()>) -> Result<()>,
Expand description

Write nucleosome calls as GFF, pulling records from produce as they are needed.

One feature per call. The attribute column carries the called sequence, as the reference implementation does.

Coordinates follow the reference rather than the GFF specification. It prints dyad - half_width directly, which is a zero-based index into the record, where GFF expects one-based inclusive coordinates. Everything it emits is therefore one base to the left of where a browser will read it. Reproduced so that positions can be compared against the reference; see the Algorithm Issues page.