Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AE015451	4,740,816	T→G	intergenic (+53/+700)	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,816	0	T	G	100.0%	56.4 / NA	18	intergenic (+53/+700)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (8/10);  total (8/10)

GCCCTGACCGATCCTTACTTCATCGAGCGCTCGCTGTACCCGAACGTCGACTTCTACTCGGGGATCATCCTCAAGGCGATCGGCATTCCGACCAGCATGTTCACCGTGATCTTCGCCCTGGCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAGGTCACTTCAGTGCTGTCACCCTTGGGCACCAGCTTCCAGAAGCCGTCGACTTTTGCCACCCGAACAAAACCCTTCTCTTCCGGCAGGTAAGTCGGCTCT  >  AE015451/4740515‑4741094                                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                       gCCCTGACCGATCCTTACTTCATCGAGCGCTCGCTGTACCCGAACGTCGACTTCTACTCGGGGATCATCCTCAAGGCGATCGGCATTCCGACCAGCATGTTCACCGTGATCTTCGCCCTGGCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTAtttt                                                                                                                                                                                                                                                                                        >  1:133293/1‑301 (MQ=255)                                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)                                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                       GCCCTGACCGATCCTTACTTCATCGAGCGCTCGCTGTACCCGAACGTCGACTTCTACTCGGGGATCATCCTCAAGGCGATCGGCATTCCGACCAGCATGTTCACCGTGATCTTCGCCCTGGCACGTACCGTGGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGGTAGGTCACTTCAGTGCTGTCACCCTTGGGCACCAGCTTCCAGAAGCCGTCGACTTTTGCCACCCGAACAAAACCCTTCTCTTCCGGCAGGTAAGTCGGCTCT  >  AE015451/4740515‑4741094

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