Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,271,600 (C)5→7 coding (171/1146 nt) garK ← glycerate kinase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,271,5951.C100.0% 89.0 / NA 24A59G (GCT→GGT) garKglycerate kinase I
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (12/12);  total (12/12)
*NC_0009133,271,5952.C100.0% 88.2 / NA 24A59G (GCT→GGT) garKglycerate kinase I
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (12/12);  total (12/12)

CGCGTTTTTCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAG‑‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGG  >  NC_000913/3271458‑3271724
                                                                                                                                          ||                                                                                                                                 
cgcgTTTTTCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAG‑‑c                                                                                                                                  <  1:382629/139‑1 (MQ=255)
                     gCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCCCTGGGTGGCTGCAatc                                                                                                               <  1:175021/139‑1 (MQ=255)
                       tCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCTGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCCCTGGGTGGCTGCAatcat                                                                                                             <  2:418397/139‑1 (MQ=255)
                           gCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCCCTGGGTGGCTGCAATCATCGCt                                                                                                         >  2:847829/1‑139 (MQ=255)
                                cTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCCCTGGGTGGCTGCAATCATCGCTTCCAc                                                                                                    >  2:609437/1‑139 (MQ=255)
                                     cgccgccATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTc                                                                                               >  2:947984/1‑139 (MQ=255)
                                                           gTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAAc                                                                            >  1:306721/1‑136 (MQ=255)
                                                           gTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAAc                                                                            <  2:306721/136‑1 (MQ=255)
                                                             tttGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGaa                                                                       <  1:947984/139‑1 (MQ=255)
                                                             tttGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGaa                                                                       <  1:1018747/139‑1 (MQ=255)
                                                             tttGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGaa                                                                       <  1:509665/139‑1 (MQ=255)
                                                                 ccATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGaa                                                                       <  1:1080947/135‑1 (MQ=255)
                                                                 ccATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGaa                                                                       >  2:1080947/1‑135 (MQ=255)
                                                                                ccccAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATc                                                    <  2:8867/139‑1 (MQ=255)
                                                                                ccccAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAG‑‑ccccccc                                                                                                                            >  2:611284/1‑63 (MQ=255)
                                                                                ccccAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAG‑‑ccccccc                                                                                                                            <  1:611284/65‑3 (MQ=255)
                                                                                      cTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAAc                                                                            >  1:50956/1‑109 (MQ=255)
                                                                                      cTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAAc                                                                            <  2:50956/109‑1 (MQ=255)
                                                                                                cTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGa                                    >  2:352927/1‑139 (MQ=255)
                                                                                                   tctcGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAAt                                 >  1:323617/1‑139 (MQ=255)
                                                                                                        cccAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCg                                                                                 >  2:699102/1‑86 (MQ=255)
                                                                                                        cccAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCg                                                                                 <  1:699102/86‑1 (MQ=255)
                                                                                                                 ccTGTAACCCAGGCGTGACGTTCAGCCCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCa                                                                                       <  2:349989/71‑1 (MQ=255)
                                                                                                                 ccTGTAACCCAGGCGTGACGTTCAGCCCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCa                                                                                       >  1:349989/1‑71 (MQ=255)
                                                                                                                             gCGTGACGTTCAGCCCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTc       >  1:690268/1‑139 (MQ=255)
                                                                                                                                tGACGTTCAGCCCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCt    <  2:323617/139‑1 (MQ=255)
                                                                                                                                  aCGTTCAGCCCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTgg  >  1:390654/1‑139 (MQ=255)
                                                                                                                                          ||                                                                                                                                 
CGCGTTTTTCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAG‑‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGG  >  NC_000913/3271458‑3271724

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