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.9 / NA	26	intergenic (+53/+700)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (13/13);  total (13/13)

CGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCT‑CTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACC  >  NC_002947/4740684‑4740953                                                                                                                                     |                                                                                                                                           cGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCct‑ct                                                                                                                           >  6:93630/1‑148 (MQ=255)      aaGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCT‑CTTcgcg                                                                                                                      >  4:154004/1‑148 (MQ=255)        gATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCT‑CTTCGCGGc                                                                                                                    <  5:42406/148‑1 (MQ=255)          tCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCT‑CTTCGCGGCCt                                                                                                                  >  5:12390/1‑148 (MQ=255)              ccGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCT‑CTTCGCGGCCTTGcc                                                                                                              <  6:149831/148‑1 (MQ=255)                      gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCT‑CTTCGCGGGCAAGGCCGCGAaga                                                                                                      <  8:88914/148‑1 (MQ=38)                      gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCT‑CTTCGCGGCCTTGCCCGCGAaga                                                                                                      <  6:12390/148‑1 (MQ=255)                       ccAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCT‑CTTCGCGGCCTTGCCCGCGAagag                                                                                                     <  3:6893/148‑1 (MQ=255)                             gTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCT‑CTTCGCGGCCTTGCCCGCGAAGAGGCCGGc                                                                                               >  6:173223/1‑148 (MQ=255)                                   cGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCT‑CATAGCCGCCATTCCCCCGAATAGGCCCGCACAGac                                                                                         >  6:123639/1‑148 (MQ=255)                                          cAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCA‑CATAGCCGCCGTGCCCGCGAAGACGGCGGCACACAccccacc                                                                                   >  2:186798/1‑146 (MQ=255)                                               gcgcGACATCGTTGCCCTGAAAGACCGCAAGTATGCTGCTGTAGCCGACCGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGGGCCGGCCT‑CTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCt                                                                             <  6:93903/148‑1 (MQ=255)                                                        cGTTGCCCTTAAAGACCGCAAGTAAGCTGCAGTAGCCCAACGCAAGAAGGCTGCCCTCGGGCAG‑CATTTTTTATTTGTCGTGTGCCGGCCT‑CTTCGCGGCCTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTtc                                                                    <  6:191646/148‑1 (MQ=255)                                                            gCCCTGAAAGACCGCAAGTACGCTGCAGTAGCCTAACGCAATACGGCTGCCCGGGGGCAG‑CGTTTTTTATTGGTCCTGTGCCGGCCT‑CTTCGCGGCCTTGCCCGCTAATAGGCCGGCACAGACACCACAAAACCTCAGTTTTtct                                                                   <  2:129129/145‑1 (MQ=255)                                                                tGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCT‑CTTAGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACCAAACCTCAGTTTTTCTccgccc                                                             >  4:31828/1‑147 (MQ=255)                                                                      aCCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGGCTAATTGCCGGCCTTGCCC‑CGAAGTGGGCGGCACCAACAACACAAAAACCAAGGTTTTCTaagcccccgacac                                                      >  4:2614/1‑141 (MQ=255)                                                                           aaGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑ATTTTTATTTGTCCTGTGCCGGACT‑CTTCGCGGACTTTCACGCGAAGAGGCCGGCACAGCCACCACAAAACCTCAGTTTTTATCCGCCCCCTCACGCAGcc                                                 >  4:146588/1‑148 (MQ=255)                                                                                  cTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCT‑CTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAg                                          >  2:4307/1‑148 (MQ=255)                                                                                   ggCAGTACCCCAACGCAAGACGGCGGCCCTCGGGCAG‑CGTTTTTTATTTGCCCGGTGCCGGCCG‑CTGCGCGGCCTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                         <  5:169512/147‑1 (MQ=255)                                                                                             gcaCGCATGACGGGTGCCGTGGGGCTGCTTTTTTTTTTTGGCCTGG‑GCCGGCGG‑CTGCGCGGCCTGGCCCGGGATGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACg                               <  2:122042/146‑1 (MQ=255)                                                                                                 gCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCT‑ATTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGCCACCCCAAAACCTCAGTTTTTCTCCGCCCGCAAACGCAGCCCATTCAGGGTATTGAACGgcgc                           >  6:56293/1‑148 (MQ=38)                                                                                                    agaagGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCT‑CTTCGCGGCCTTGCCCGCGCAGAGGCCGGCCCAGACACCACAAAACCTCAGTGTTTCTCCGCCCGCTCACGCCGCCCGTTCAGGGTATTGAACGGCGCATc                        >  8:170902/1‑148 (MQ=255)                                                                                                           tGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGACGGCCT‑ATTAGCGGCCTTGCCCGCGAATAGGACGGCACACACAAAAAAAAAACAAAGTTTTTATCACCCACCCCACGCAACCCATTCAGGGTATTGAACGGCGCATCcaccac                  >  8:23773/1‑147 (MQ=255)                                                                                                                    ggaagcTTTTTTTATTGGGCCTGGGCCGGCCT‑CTTCGGGGCCTGGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgtt         <  3:8519/143‑1 (MQ=255)                                                                                                                         ccttttttATTTGTCCTGTGCCGGCCT‑CTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGcca    <  5:93630/146‑1 (MQ=255)                                                                                                                           ttttttATTTGTCCTGTGCCGGCCT‑CTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGccacc  <  3:45694/148‑1 (MQ=255)                                                                                                                                     |                                                                                                                                           CGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCT‑CTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACC  >  NC_002947/4740684‑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: ‑