Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	4,743,226	2 bp→GC	intergenic (+495/+256)	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,743,226	0	C	G	100.0%	129.1 / NA	36	intergenic (+495/+257)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (20/16);  total (20/16)
*	NC_002947	4,743,227	0	A	C	100.0%	131.8 / NA	36	intergenic (+496/+256)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (20/16);  total (20/16)

GTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AACCTCCGCCAAC‑TGCATTCGTGAACCAAGCAACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTG  >  NC_002947/4743104‑4743345                                                                                                                             ||                                                                                                                       gTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGCATTCGTGAATCCAGGCACAGGCCCCgg                                                                                                             <  1:135637/136‑1 (MQ=255)   cGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGCATTCGTGAATCCAGGCACAGGCCCCGGac                                                                                                           <  1:660599/136‑1 (MQ=255)   cGTAATGAGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGCATTCGTGAATCCAGGCACAGGCCCCGGac                                                                                                           <  2:1294829/136‑1 (MQ=255)    gTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGCATTCGTGAATCCAGGCACAGGCCCCGGaca                                                                                                          >  1:921730/1‑136 (MQ=255)    gTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGCATTCGTGAATCCAGGCACAGGCCCCGGaca                                                                                                          >  1:554983/1‑136 (MQ=255)       aTGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGCATTCGTGAATCCAGGCACAGGCCCCGGACACGt                                                                                                       >  1:1262084/1‑136 (MQ=255)                   gAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAgcg                                                                                           >  2:679239/1‑136 (MQ=255)                   gAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAgcg                                                                                           <  1:726329/136‑1 (MQ=255)                       cACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTc                                                                                       >  1:905258/1‑136 (MQ=255)                                 cACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGcc                                                                             >  1:1247698/1‑136 (MQ=255)                                  aCCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCa                                                                            <  2:1208203/136‑1 (MQ=255)                                                 ttGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCa                                                             >  1:1152405/1‑136 (MQ=255)                                                      cTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTc                                                        <  1:1094299/136‑1 (MQ=255)                                                      cTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTc                                                        <  1:999061/136‑1 (MQ=255)                                                      cTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTc                                                        <  1:1084607/136‑1 (MQ=255)                                                       tGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTca                                                       <  2:1262084/136‑1 (MQ=255)                                                         ggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTcaca                                                     >  2:1326037/1‑136 (MQ=255)                                                            ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgcggc                                                                       >  1:1098152/1‑115 (MQ=255)                                                            ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgcggc                                                                       <  2:1098152/115‑1 (MQ=255)                                                            ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcg                                                  >  1:757731/1‑136 (MQ=255)                                                            ataGGTCCAGCTCTGGTCACCATCGTGCTTTAGCAGCTTCTGCGACT‑TTCATTCGTGAATCCAGGCACAGTCACCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGACTTGATCACGGTCACAccc                                                  >  1:847562/1‑135 (MQ=255)                                                             taGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc                                                 <  1:729924/136‑1 (MQ=255)                                                             taGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc                                                 <  2:921730/136‑1 (MQ=255)                                                             taGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc                                                 <  1:639019/136‑1 (MQ=255)                                                                      ctctGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGcca                                        >  1:1236861/1‑135 (MQ=255)                                                                      ctctGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGcca                                        >  1:176733/1‑135 (MQ=255)                                                                      ctctGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCt                                        >  2:1077631/1‑136 (MQ=255)                                                                      ctctGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCt                                        >  1:1120055/1‑136 (MQ=255)                                                                      ctctGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCt                                        >  1:983995/1‑136 (MQ=255)                                                                      ctctGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCt                                        >  2:944430/1‑136 (MQ=255)                                                                             cACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACtt                                 >  2:278525/1‑136 (MQ=255)                                                                              aCCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTg                                <  2:267030/136‑1 (MQ=255)                                                                              aCCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTg                                <  2:512928/136‑1 (MQ=255)                                                                              aCCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTg                                >  1:179978/1‑136 (MQ=255)                                                                                   cGTGCTTGAG‑CAGCTTCTGCGACTTGCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGAccc                           >  2:1075219/1‑136 (MQ=255)                                                                                                             gCATTCGTGAATCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTg  <  1:944430/136‑1 (MQ=255)                                                                                                                             ||                                                                                                                       GTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AACCTCCGCCAAC‑TGCATTCGTGAACCAAGCAACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTG  >  NC_002947/4743104‑4743345

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