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

TCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACC  >  NC_002947/4740677‑4740953                                                                                                                                            |                                                                                                                                          tcCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGcc                                                                                                                                  >  1:44431/1‑148 (MQ=255)  ccAGCCCGTACAAGATCGGCCGTCCGAGCCAGCTGTACACCGGCGAGCAGAAGAGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGACCTCGGGCAG‑CGTTTTTTATTTGTCCTGTGCCg                                                                                                                                 >  8:102762/1‑148 (MQ=255)   cAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCgg                                                                                                                                >  2:128631/1‑148 (MQ=255)          tACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTc                                                                                                                         <  3:53427/148‑1 (MQ=255)                           gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                                        <  8:26375/148‑1 (MQ=255)                           gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                                        <  7:155620/148‑1 (MQ=255)                             gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                      <  4:28902/148‑1 (MQ=38)                             gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                      <  6:124702/148‑1 (MQ=255)                              ccAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAagag                                                                                                     <  1:98834/148‑1 (MQ=255)                              ccAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAagag                                                                                                     <  5:98551/148‑1 (MQ=255)                                  cTGTACACGGGCGATCATAAGCGCGACATCGTTGCCCTTAAACACCGCAAGTATGCTGCTGTAGCCGAACGCAATATGGCTGCCGTCGGGCTG‑CTTTTTTTTTTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCTAATAGGCCg                                                                                                 <  2:176225/148‑1 (MQ=255)                                        acCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCATATCCCGCGAAGAGGCCGGCACAg                                                                                           >  8:73927/1‑148 (MQ=255)                                                                gTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGCAGCGGCCGGCACAGCCACCACAAAACCTCAGTTTTtct                                                                   >  5:98961/1‑148 (MQ=255)                                                                               cGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTTCCCGCGAAGCGGCCGGCCCAGACACCACAAAACCACAGTTTTTCTACGCCCGCTCACGCa                                                    >  8:39352/1‑148 (MQ=255)                                                                                     tAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCATGTGCAGGCCTCTTAGCGGCCTTGCCCGCGAAGAGGAAGGCACAGACACCACAAAACCTAAGTTTTTATCCGCCCGCTCACGCAGcccat                                              >  7:58356/1‑146 (MQ=255)                                                                                        gCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCTTGGCCGCGAAGCGGCCGGGCCAGACCCCACAAAACCGCAGTTTTTATCCCCCCGCTCCCGCCCCCCGTTCa                                           >  1:61857/1‑148 (MQ=255)                                                                                          tGCAGTAGCCTAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                         <  6:47016/148‑1 (MQ=255)                                                                                          tGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGCCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                         <  5:47655/148‑1 (MQ=255)                                                                                                          aagaagGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCCCAGACCCCACAAAACCTCAGGTTTTCTACGCCCGCTCACGCAGCCCTTTCAGGGTATTTAACGgcgcct                         >  3:58447/1‑146 (MQ=255)                                                                                                           atacggCTGCCCTCGGGCAG‑CTTTTTTTATTTGCCCTGTCCCGGCCGCTACGCGCCGTGCCCCGCTAACAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATc                        <  7:73927/144‑1 (MQ=255)                                                                                                             aagGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATcca                      <  2:126404/148‑1 (MQ=255)                                                                                                                 cTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAAACACCACAAAACCACAGTTTTTATCCGCCCGCTAACGCAGCCCGTTAACGGTATTGAACGGCGCAAAcaccac                  >  8:25152/1‑148 (MQ=255)                                                                                                                  tGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCATGCCCGCGAAGACGCCGGCACACACACCACAAAACCTCCCTTTTTATCCCCCCGCGAACGCACCCCAGAGATCGTATTTAACGGCGCATCCACCACg                 >  4:104490/1‑148 (MQ=255)                                                                                                                           ccAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCTGCCTGGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgttg        <  2:62671/147‑1 (MQ=255)                                                                                                                                  ttttttATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGCGGCCGGCACAGACACCACAAAACCGCAGTTTTTCTCCGCCCGCTCACGCAGCCCGTTCCGGGTATTGCACGGCGCATCCCCCACGAACTTGTTGGccacc  >  6:183861/1‑148 (MQ=255)                                                                                                                                            |                                                                                                                                          TCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACC  >  NC_002947/4740677‑4740953

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑