can pass in custom reference seq if needed
use Bio::SeqIO module