see_chrom_read_input.c
This file contains functions for reading DNA sequences out of a .fa format file.
The core of the logic is see_chrom_read_go(), which
reads a .fa file, skips over comments, and passes
each nucleotide one at a time to see_chrom_gen_model.
(Conceptually, this is lower-level code, but callbacks are often slower than
normal calls because they are harder to inline/optimize,
and this is a performance-critical loop. So I'll just have read_input.c
be at a higher level than see_chrom_gen_model and feed it data).
For simplicity, this tool treats lower-case bases
the same as upper-case, and skips over variable bases like n.
For maximum performance, input files must be in .fa format with
only one > comment at the start of the file.
See documentation for more information.
Currently using data from USCS, see the page
http://hgdownload.cse.ucsc.edu/goldenPath/hg19/chromosomes/