Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AE015451	4,740,819	2 bp→CT	intergenic (+56/+696)	gltA → / ← yeiW	citrate synthase/putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	4,740,819	0	T	C	100.0%	61.9 / NA	18	intergenic (+56/+697)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (8/10);  total (8/10)
*	AE015451	4,740,820	0	G	T	100.0%	60.0 / NA	18	intergenic (+57/+696)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (8/10);  total (8/10)

CGCTGTACCCGAACGTCGACTTCTACTCGGGGATCATCCTCAAGGCGATCGGCATTCCGACCAGCATGTTCACCGTGATCTTCGCCCTGGCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAGGTCACTTCAGTGCTGTCACCCTTGGGCACCAGCTTCCAGAAGCCGTCGACTTTTGCCACCCGAACAAAACCCTTCTCTTCCGGCAGGTAAGTCGGCTCT  >  AE015451/4740546‑4741094                                                                                                                                                                                                                                                                                  ||                                                                                                                                                                                                                                                                                   cgctGTACCCGAACGTCGACTTCTACTCGGGGATCATCCTCAAGGCGATCGGCATTCCGACCAGCATGTTCACCGTGATCTTCGCCCTGGCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGTCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCcg                                                                                                                                                                                                                                                         <  2:865343/301‑1 (MQ=255)        cccGAACGTCGACTTCTACTCGGGGATCATCCTCAAGGCGATCGGCATTCCGACCAGCATGTTCACCGTGATCTTCGCCCTGGCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCTGCCTCTTCGCCGCCTTGCCCGCGAaggg                                                                                                                                                                                                                                                  >  1:167006/1‑299 (MQ=255)               gTCGACTTCTACTCGGGGATCATCCTCAAGGCGATCGGCATTCCGACAAGCATGTTCACCGTGATCTTCGCCCTGGCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGca                                                                                                                                                                                                                                           <  1:97778/301‑1 (MQ=255)                                 atcatcCTCAAGGCGATCGCCCTTCCAACCAGCATGTGCACCGTGATCTTCGCCCTGGCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTc                                                                                                                                                                                                                         <  1:854760/301‑1 (MQ=255)                                       cTCAAGGCGATCGGCATTCCGACCAGCATGGTCCCCGTGATCTTCGCCCTGGCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGtttt                                                                                                                                                                                                                   <  2:276972/301‑1 (MQ=255)                                                                           gggATCTTCGCCCTGGCACGTACCATGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGa                                                                                                                                                                               <  1:357511/299‑1 (MQ=255)                                                                           gTGATCTTCGCCCTGGCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGCCCGGCACAGCCACCACAAAACCTCAGTTTTTCTCCGCCCGCCCACGCAGCCCTCTCAGGGTATTGa                                                                                                                                                                               >  1:575086/1‑301 (MQ=255)                                                                                                                  cACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTTGCCACCATCGCCg                                                                                                                                        >  1:43527/1‑301 (MQ=255)                                                                                                                                                                gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGTCCACCATCGCCGGCACGCTACCACCCGTCTCGGTGTGCCCCTGGTAGTGCACTTCatt                                                                                          >  2:94441/1‑299 (MQ=255)                                                                                                                                                                                   gCAGAAGCGCGCCATCGATGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATT‑AACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAGGTCACTTCAGTGCTGTCACCCTTGGGCACCa                                                                      <  2:166940/301‑1 (MQ=255)                                                                                                                                                                                                       cccTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAGGTCACTTCAGTGCTGTCACCCTTGGGCACCAGCTTCCAGAAGCCGTCGAc                                                   <  1:709543/301‑1 (MQ=255)                                                                                                                                                                                                                 ccGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAGGTCACTTCAGTGCTGTCACCCTTGGGCACCAGCTTCCAGAAGCCGTCGACTTTTGCCAcc                                         <  1:730606/301‑1 (MQ=255)                                                                                                                                                                                                                      aGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAGGTCACTTCAGTGCTGTCACCCTTGGGCACCAGCTTCCAGAAGCCGTCGACTTTTGCCACCCGAAc                                    <  2:43521/301‑1 (MQ=255)                                                                                                                                                                                                                                    cccGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGCGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAGGTCACTTCAGTGCTGTCACCCTTGGGCACCAGCTTCCAGAAGCCGTCGACTTTTGCCACCCGAACAAAACCCTTCTCtt                      <  1:234325/300‑1 (MQ=255)                                                                                                                                                                                                                                            cAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAGGTCACTTCAGTGCTGTCACCCTTGGGCACCAGCTTCCAGAAGCCGTCGACTTTTGCCACCCGAACAAAACCCTTCTCTTCCGGCAgg              >  1:665911/1‑301 (MQ=255)                                                                                                                                                                                                                                                    cTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAGGTCACTTCAGTGCTGTCACCCTTGGGCACCAGCTTCCAGAAGCCGTCGACTTTTGCCACCCGAACAAAACCCTTCTCTTCCGGCAGGTAAGTCgg      >  1:210035/1‑301 (MQ=255)                                                                                                                                                                                                                                                     tGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAGGTCACTTCAGTGCTGTCACCCTTGGGCACCAGCTTCCAGAAGCCGTCGACTTTTGCCACCCGAACAAAACCCTTCTCTTCCGGCAGGTAAGTCGGc     >  1:230018/1‑301 (MQ=255)                                                                                                                                                                                                                                                        ccTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAGGTCACTTCAGTGCTGTCACCCTTGGGCACCAGCTTCCAGAAGCCGTCGACTTTTGCCACCCGAACAAAACCCTTCTCTTCCGGCAGGTAATTCGGctct  >  1:429230/1‑301 (MQ=255)                                                                                                                                                                                                                                                                                  ||                                                                                                                                                                                                                                                                                   CGCTGTACCCGAACGTCGACTTCTACTCGGGGATCATCCTCAAGGCGATCGGCATTCCGACCAGCATGTTCACCGTGATCTTCGCCCTGGCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAGGTCACTTCAGTGCTGTCACCCTTGGGCACCAGCTTCCAGAAGCCGTCGACTTTTGCCACCCGAACAAAACCCTTCTCTTCCGGCAGGTAAGTCGGCTCT  >  AE015451/4740546‑4741094

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑