Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,973,169	G→A	G57G (GGC→GGT)	cheW ←	purine‑binding chemotaxis protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,973,169	0	G	A	100.0%	68.7 / NA	21	G57G (GGC→GGT)	cheW	purine‑binding chemotaxis protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (11/10);  total (11/10)

TCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGAC  >  NC_000913/1973058‑1973298                                                                                                                |                                                                                                                                  tCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACACCTTTGATAAACGCTGGCGTGTTCGCa                                                                                                        >  1:474543/1‑139 (MQ=255)  ccGACATTCAGGACGATAACTACCGTGTTGTCGTTACAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGCCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACACCTTTGATAAACGCTGGCGTGTTCGCaa                                                                                                       <  1:321061/139‑1 (MQ=255)    gagaTTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACACCTTTGATAAACGCTGGCGTGTTCGCAATc                                                                                                     <  2:118104/139‑1 (MQ=255)     agaTTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACACCTTTGATAAACGCTGGCGTGTTCGCAATcc                                                                                                    >  1:76012/1‑139 (MQ=255)                                         tCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACACCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACtt                                                                >  2:66170/1‑139 (MQ=255)                                          ccACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACACCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACttt                                                               >  2:491477/1‑139 (MQ=255)                                             caTCCACCTGGCCGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACACCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAg                                                            >  1:265625/1‑139 (MQ=255)                                                  aCCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACACCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCTTGCACTTTCAGGatat                                                       <  2:63407/139‑1 (MQ=255)                                                  aCCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACACCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGatat                                                       <  1:350822/139‑1 (MQ=255)                                                   ccTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACACCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATc                                                      >  1:352223/1‑139 (MQ=255)                                                   ccTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACACCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATc                                                      >  2:259304/1‑139 (MQ=255)                                                          gAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACACCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCg                                               >  2:254526/1‑139 (MQ=255)                                                             cTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACACCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTAc                                            >  1:47645/1‑139 (MQ=255)                                                                 aTTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACACCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATa                                                  >  2:378724/1‑129 (MQ=255)                                                                 aTTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACACCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATa                                                  <  1:378724/129‑1 (MQ=255)                                                                   tCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACACCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTATCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTtca                                      <  1:378277/139‑1 (MQ=255)                                                                        aGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACACCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCAcc                                 <  1:322638/139‑1 (MQ=255)                                                                              caaTCGGCGCAATAACGCCGCGCAGATTCGTGACACCTTTGGTAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTc                                                              <  1:333126/104‑1 (MQ=255)                                                                              caaTCGGCGCAATAACGCCGCGCAGATTCGTGACACCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTc                                                              >  2:333126/1‑104 (MQ=255)                                                                                                cgcgCAGATTCGTGACACCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTg         <  2:265625/139‑1 (MQ=255)                                                                                                       aTTCGTGACACCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGAc  <  2:474543/139‑1 (MQ=255)                                                                                                                |                                                                                                                                  TCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGAC  >  NC_000913/1973058‑1973298

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