Feature motifs
This feature adds a motif finder from a fasta file. It does not take into account ambiguity codes and overlapping matches.
It adds a requirement for numpy. It uses array.array('I')
in stead of a simple list of tuples to keep the memory footprint small (< 1GB
for hg19
). It outputs in bed format. To call:
fastools famotif2bed <INPUT_FASTA> <OUTPUT_BED> --motif <sequence_to_be_found>