Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,967,855	T→C	D216G (GAC→GGC)	cheB ←	protein‑glutamate methylesterase/protein glutamine deamidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,967,855	0	T	C	100.0%	144.0 / NA	40	D216G (GAC→GGC)	cheB	protein‑glutamate methylesterase/protein glutamine deamidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (17/23);  total (17/23)

CATCTACCGAAGGCCGATGACGGTTAACCGCCGGGCCATCGTGAATTTTGATTTGGTAATTTGCGCCACTACGCGACAGCTCCATATGCCGATCGCCCGGCGCAATATAGGCATGCCCCGGCAAGACACGTTCTCCGTCTTCGGCTTCTTTAACCCCGATCTGGCAAAGCTTATTAAGTCTGTCGGCAAAAGAGCGGGTGAAACCGGGCGGCATATGCTGGGTAATTAACAGTGCCGGGCTGGAAAGCGGCAACGGTT  >  NC_000913/1967718‑1967975                                                                                                                                          |                                                                                                                         caTCTACCGAAGGCCGATGACGGTTAACCGCCGGGCCATCGTGAATTTTGATTTGGTAATTTGCGCCACTACGCGACAGCTCCATATGCCGATCGCCCGGCGCAATATAGGCATGCCCCGGCAAGACACGTTCTCCGCCt                                                                                                                        >  2:136959/1‑140 (MQ=255)    cTACCGAAGGCCGATGACGGTTAACCGCCGGGCCATCGTGAATTTTGATTTGGTAATTTGCGCCACTACGCGACAGCTCCATATGCCGATCGCCCGGCGCAATATAGGCATGCCCCGGCAAGACACGTTCTCCGCCTTCg                                                                                                                     >  2:253652/1‑140 (MQ=255)      aCCGAAGGCCGATGACGGTTAACCGCCGGGCCATCGTGAATTTTGATTTGGTAATTTGCGCCACTACGCGACAGCTCCATATGCCGATCGCCCGGCGCAATATAGGCATGCCCCGGCAAGACACGTTCTCCGCCTTCGGc                                                                                                                   <  1:997968/140‑1 (MQ=255)      aCCGAAGGCCGATGACGGTTAACCGCCGGGCCATCGTGAATTTTGATTTGGTAATTTGCGCCACTACGCGACAGCTCCATATGCCGATCGCCCGGCGCAATATAGGCATGCCCCGGCAAGACACGTTCTCCGCCTTCGGc                                                                                                                   <  1:811306/140‑1 (MQ=255)      aCCGAAGGCCGATGACGGTTAACCGCCGGGCCATCGTGAATTTTGATTTGGTAATTTGCGCCACTACGCGACAGCTCCATATGCCGATCGCCCGGCGCAATATAGGCATGCCCCGGCAAGACACGTTCTCCGCCTTCGGc                                                                                                                   <  1:804165/140‑1 (MQ=255)        cGAAGGCCGATGACGGTTAACCGCCGGGCCATCGTGAATTTTGATTTGGTAATTTGCGCCACTACGCGACAGCTCCATATGCCGATCGCCCGGCGCAATATAGGCATGCCCCGGCAAGACACGTTCTCCGCCTTCGGctt                                                                                                                 >  1:43968/1‑140 (MQ=255)         gAAGGCCGATGACGGTTAACCGCCGGGCCATCGTGAATTTTGATTTGGTAATTTGCGCCACTACGCGACAGCTCCATATGCCGATCGCCCGGCGCAATATAGGCATGCCCCGGCAAGACACGTTCTCCGCCTTCGGcttc                                                                                                                >  2:610846/1‑140 (MQ=255)          aaGGCCGATGACGGTTAACCGCCGGGCCATCGTGAATTTTGATTTGGTAATTTGCGCCACTACGCGACAGCTCCATATGCCGATCGCCCGGCGCAATATAGGCATGCCCCGGCAAGACACGTTCTCCGCCTTCGGcttct                                                                                                               <  1:119340/140‑1 (MQ=255)           aGGCCGATGACGGTTAACCGCCGGGCCATCGTGAATTTTGATTTGGTAATTTGCGCCACTACGCGACAGCTCCATATGCCGATCGCCCGGCGCAATATAGGCATGCCCCGGCAAGACACGTTCTCCGCCTTCGGcttctt                                                                                                              <  2:605754/140‑1 (MQ=255)              ccGATGACGGTTAACCGCCGGGCCATCGTGAATTTTGATTTGGTAATTTGCGCCACTACGCGACAGCTCCATATGCCGATCGCCCGGCGCAATATAGGCATGCCCCGGCAAGACACGTTCTCCGCCTTCGGCTTCTTTaa                                                                                                           >  2:824508/1‑140 (MQ=255)                       gTTAACCGCCGGGCCATCGTGAATTTTGATTTGGTAATTTGCGCCACTACGCGACAGCTCCATATGCCGATCGCCCGGCGCAATATAGGCAGGCCCCGGCGAGAGACGGTCTCCGCCCTCGGGTTCTTTAACCCCGATCt                                                                                                  >  2:308049/1‑140 (MQ=255)                            ccgccgGGCCATCGTGAATTTTGATTTGGTAATTTGCGCCACTACGCGACAGCTCCATATGCCGATCGCCCGGCGCAATATAGGCATGCCCCGGCAAGACACGTTCTCCGCCTTCGGCTTCTTTAACCCCGATCTGGCaa                                                                                             <  1:864806/140‑1 (MQ=255)                            ccgccgGGCCATCGTGAATTTTGATTTGGTAATTTGCGCCACTACGCGACAGCTCCATATGCCGATCGCCCGGCGCAATATAGGCATGCCCCGGCAAGACACGTTCTCCGCCTTCGGCTTCTTTAACCCCGATCTGGCaa                                                                                             <  2:10315/140‑1 (MQ=255)                              gccgGGCCATCGTGAATTTTGATTTGGTAATTTGCGCCACTACGCGACAGCTCCATATGCCGATCGCCCGGCGCAATATAGGCATGCCCCGGCAAGACACGTTCTCCGCCTTCGGCTTCTTTAACCCCGATCTGGCAAAg                                                                                           <  1:695767/140‑1 (MQ=255)                              gccgGGCCATCGTGAATTTTGATTTGGTAATTTGCGCCACTACGCGACAGCTCCATATGCCGATCGCCCGGCGCAATATAGGCATGCCCCGGCAAGACACGTTCTCCGCCTTCGGCTTCTTTAACCCCGATCTGGCAAAg                                                                                           <  1:610846/140‑1 (MQ=255)                                  ggCCATCGTGAATTTTGATTTGGTAATTTGCGCCACTACGCGACAGCTCCATATGCCGATCGCCCGGCGCAATATAGGCATGCCCCGGCAAGACACGTTCTCCGCCTTCGGCTTCTTTAACCCCGATCTGGCaa                                                                                             >  1:911106/1‑134 (MQ=255)                                  ggCCATCGTGAATTTTGATTTGGTAATTTGCGCCACTACGCGACAGCTCCATATGCCGATCGCCCGGCGCAATATAGGCATGCCCCGGCAAGACACGTTCTCCGCCTTCGGCTTCTTTAACCCCGATCTGGCaa                                                                                             <  2:911106/134‑1 (MQ=255)                                            aaTTTTGATTTGGTAATTTGCGCCACTACGCGACAGCTCCATATGCCGATCGCCCGGCGCAATATAGGCATGCCCCGGCAAGACACGTTCTCCGCCTTCGGCTTCTTTAACCCCGATCTGGCAAAGCTTATTAAGTCTGt                                                                             <  2:43968/140‑1 (MQ=255)                                                      tGGTAATTTGCGCCACTACGCGACAGCTCCATATGCCGATCGCCCGGCGCAATATAGGCATGCCCCGGCAAGACACGTTCTCCGCCTTCGGCTTCTTTAACCCCGATCTGGCAAAGCTTATTAAGTCTGTCGGCAAAaga                                                                   <  1:800302/140‑1 (MQ=255)                                                                   cACTACGCGACAGCTCCATATGCCGATCGCCCGGCGCAATATAGGCATGCCCCGGCAAGACACGTTCTCCGCCTTCGGCTTCTTTAACCCCGATCTGGCAAAGCTTATTAAGTCTGTCGGCAAAAGAGCGGGTGAAACCg                                                      >  1:823261/1‑140 (MQ=255)                                                                    ccTACGCGACAGCTCCATATGCTGATCGCTCGGCGCAATATAGGCATGCCCCGGCAAGACACGTTCTCCGCCTTCGGCTTCTTTAACCCCGATCTGGCAAAGCTTATTAAGTCTGTCGGCAAAAGAGCGGGTGAAACCgg                                                     <  1:385164/139‑1 (MQ=255)                                                                    aCTACGCGACAGCTCCATATGCCGATCGCCCGGCGCAATATAGGCATGCCCCGGCAAGACACGTTCTCCGCCTTCGGCTTCTTTAACCCCGATCTGGCAAAGCTTATTAAGTCTGTCGGCAAAAGAGCGGGTGAAACCgg                                                     <  1:29046/140‑1 (MQ=255)                                                                     cTACGCGACAGCTCCATATGCCGATCGCCCGGCGCAATATAGGCATGCCCCGGCAAGACACGTTCTCCGCCTTCGGCTTCTTTAACCCCGATCTGGCAAAGCTTATTAAGTCTGTCGGCAAAAGAGCGGGTGAAACCggg                                                    <  1:918718/140‑1 (MQ=255)                                                                         gcgACAGCTCCATATGCCGATCGCCCTGCGCAATATAGGCATGCCCCGGCAAGACACGTTCTCCGCCTTCGGCTTCTTTAACCCCGATCTGGCAAAGCTTATTAAGTCTGTCGGCAAAAGAGCGGGTGAAACCGggcggc                                                <  1:92237/140‑1 (MQ=255)                                                                         gcgACAGCTCCATATGCCGATCGCCCGGCGCAATATAGGCATGCCCCGGCAAGACACGTTCTCCGCCTTCGGCTTCTTTAACCCCGATCTGGCAAAGCTTATTAAGTCTGTCGGCAAAAGAGCGGGTGAAACCGggcggc                                                <  1:929534/140‑1 (MQ=255)                                                                         gcgACAGCTCCATATGCCGATCCCCTGGCTCAATATATGCATGCCCCGGCAAGACACGTTCTCCGCCTTCGGCTTCTTTAACCCCGATCTGGCAAAGCTTATTAAGTCTGTCGGCAAAAGAGCGGGTGAAACCGggcggc                                                <  1:308049/140‑1 (MQ=255)                                                                             cAGCTCCATATGCCGATCGCCCGGCGCAATATAGGCATGCCCCGGCAAGACACGTTCTCCGCCTTCGGCTTCTTTAACCCCGATCTGGCAAAGCTTATTAAGTCTGTCGGCAAAAGAGCGGGTGAAACCGGGCGGCatat                                            >  1:578094/1‑140 (MQ=255)                                                                                cTCCATATGCCGATCGCCCGGCGCAATATAGGCATGCCCCGGCAAGACACGTTCTCCGCCTTCGGCTTCTTTAACCCCGATCTGGCAAAGCTTATTAAGTCTGTCGGCAAAAGAGCGGGTGAAACCGGGCGGCATATGCt                                         <  2:823261/140‑1 (MQ=255)                                                                                cTCCATATGCCGATCGCCCGGCGCAATATAGGCATGCCCCGGCAAGACACGTTCTCCGCCTTCGGCTTCTTTAACCCCGATCTGGCAAAGCTTATTAAGTCTGTCGGCAAAAGAGCGGGTGAAACCGGGCGGCATATGCt                                         <  1:136959/140‑1 (MQ=255)                                                                                           gATCGCCCGGCGCAATATAGGCATGCCCCGGCAAGACACGTTCTCCGCCTTCGGCTTCTTTAACCCCGATCTGGCAAAGCTTATTAAGTCTGTCGGCAAAAGAGCGGGTGAAACCGGGCGGCATATGCTGGGTAATTAAc                              >  1:168999/1‑140 (MQ=255)                                                                                                                cATGCCCCGGCAAGACACGTTCTCCGCCTTCGGCTTCTTTAACCCCGATCTGGCAAAGCTTATTAAGTCTGTCGGCAAAAGAGCGGGTGAAACCGGGCGGCATATGCTCGGTAATTAACAGTGCCGGGCTGGAAAGCGGc         <  2:578094/140‑1 (MQ=255)                                                                                                                   gCCCCGGCAAGACACGTTCTCCGCCTTCGGCTTCTTTAACCCCGATCTGGCAAAGCTTATTAAGTCTGTCGGCAAAAGAGCGGGTGAAACCGGGCGGCATATGCTGGGTAATTAACAGTGCCGGGCTGGAAAGCGGCAAc      >  1:179528/1‑140 (MQ=255)                                                                                                                    ccccGGCAAGACACGTTCTCCGCCTTCGGCTTCTTTAACCCCGATCTGGCAAAGCTTATTAAGTCTGTCGGCAAAAGAGCGGGTGAAACCGGGCGGCATATGCTGGGTAATTAACAGTGCCGGGCTGGAAAGCGGCAACg     >  1:884576/1‑140 (MQ=255)                                                                                                                    ccccGGCAAGACACGTTCTCCGCCTTCGGCTTCTTTAACCCCGATCTGGCAAAGCTTATTAAGTCTGTCGGCAAAAGAGCGGGTGAAACCGGGCGGCATATGCTGGGTAATTAACAGTGCCGGGCTGGAAAGCGGCAACg     <  1:824508/140‑1 (MQ=255)                                                                                                                    ccccGGCAAGACACGTTCTCCGCCTTCGGCTTCTTTAACCCCGATCTGGCAAAGCTTATTAAGTCTGTCGGCAAAAGAGCGGGTGAAACCGGGCGGCATATGCTGGGTAATTAACAGTGCCGGGCTGGAAAGCGGCAACg     >  1:707620/1‑140 (MQ=255)                                                                                                                    ccccGGCAAGACACGTTCTCCGCCTTCGGCTTCTTTAACCCCGATCTGGCAAAGCTTATTAAGTCTGTCGGCAAAAGAGCGGGTGAAACCGGGCGGCATATGCTGGGTAATTAACAGTGCCGGGCTGGAAAGCGGCAACg     >  2:904006/1‑140 (MQ=255)                                                                                                                    ccccGGCAAGACACGTTCTCCGCCTTCGGCTTCTTTAACCCCGATCTGGCAAAGCTTATTAAGTCTGTCGGCAAAAGAGCGGGTGAAACCGGGCGGCATATGCTGGGTAATTAACAGTGCCGGGCTGGAAAGCGGCAACg     >  2:930567/1‑140 (MQ=255)                                                                                                                      ccGGCAAGACACGTTCTCCGCCTTCGGCTTCTTTAACCCCGATCTGGCAAAGCTTATTAAGTCTGTCGGCAAAAGAGCGGGTGAAACCGGGCGGCATATGCTGGGTAATTAACAGTGCCGGGCTGGAAAGCGGCAACGGt   <  1:253652/140‑1 (MQ=255)                                                                                                                       cGGCAAGACACGTTCTCCGCCTTCGGCTTCTTTAACCCCGATCTGGCAAAGCTTATTAAGTCTGTCGGCAAAAGAGCGGGTGAAACCGGGCGGCATATGCTGGGTAATTAACAGTGCCGGGCTGGAAAGCGGCAACGGtt  >  2:684913/1‑140 (MQ=255)                                                                                                                       cGGCAAGACACGTTCTCCGCCTTCGGCTTCTTTAACCCCGATCTGGCAAAGCTTATTAAGTCTGTCGGCAAAAGAGCGGGTGAAACCGGGCGGCATATGCTGGGTAATTAACAGTGCCGGGCTGGAAAGCGGCAACGGtt  >  2:956658/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                         CATCTACCGAAGGCCGATGACGGTTAACCGCCGGGCCATCGTGAATTTTGATTTGGTAATTTGCGCCACTACGCGACAGCTCCATATGCCGATCGCCCGGCGCAATATAGGCATGCCCCGGCAAGACACGTTCTCCGTCTTCGGCTTCTTTAACCCCGATCTGGCAAAGCTTATTAAGTCTGTCGGCAAAAGAGCGGGTGAAACCGGGCGGCATATGCTGGGTAATTAACAGTGCCGGGCTGGAAAGCGGCAACGGTT  >  NC_000913/1967718‑1967975

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