Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,271,600 (C)5→6 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.C94.7% 62.7 / ‑2.8 19coding (176/1146 nt)garKglycerate kinase I
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (10/8);  total (10/9)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

TTTTCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCA  >  NC_000913/3271463‑3271714
                                                                                                                                     |                                                                                                                       
ttttCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGcccccc                                                                                                                   >  2:809945/1‑138 (MQ=255)
ttttCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGcccccc                                                                                                                   >  2:378251/1‑138 (MQ=255)
   tCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGcccccctgg                                                                                                                <  2:881444/139‑5 (MQ=255)
    ccGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCttggg                                                                                                               <  1:357295/139‑5 (MQ=255)
          gTACCAGCTCCTGCTCACTTGCCCCCCCCATTTCAATAAACCCGCTTTTCCCATCGCCCGAGATCCCCCAACTCGCATTCACTTTCTCGCCCAGCGCCTCTCTAACCCATGCGTGACGTTCAGCCCCCCTGGGTGGCTg                                                                                                          <  1:330283/139‑1 (MQ=255)
            aCCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCCTGGGTGGCTGCa                                                                                                        >  2:536662/1‑139 (MQ=255)
                                 gccgccATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCt                                                                                   >  1:113597/1‑139 (MQ=255)
                                   cgccATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTc                                                                                 >  2:910039/1‑139 (MQ=255)
                                                  cgcgGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAAc                                                                  <  1:356632/139‑1 (MQ=255)
                                                   gcgGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAAcc                                                                 <  1:378251/139‑1 (MQ=255)
                                                      gTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGa                                                              >  2:57797/1‑139 (MQ=255)
                                                                             ccAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAgg                                       >  1:1040189/1‑139 (MQ=255)
                                                                               aaCTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTAc                                                  >  1:610/1‑126 (MQ=255)
                                                                               aaCTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTAc                                                  <  2:610/126‑1 (MQ=255)
                                                                                           cTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGaa                         <  2:30455/139‑1 (MQ=255)
                                                                                               ctcGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCAt                                           <  2:590877/117‑1 (MQ=255)
                                                                                               ctcGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCAt                                           >  1:590877/1‑117 (MQ=255)
                                                                                                        cGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGaa                                                             <  2:536617/90‑1 (MQ=255)
                                                                                                        cGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGaa                                                             >  1:536617/1‑90 (MQ=255)
                                                                                                                  aCCCAGGCGTGACGTTCAGCCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCa  >  1:478895/1‑139 (MQ=255)
                                                                                                                     cAGGCGTGACGTTCAGCCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCa  >  1:912315/1‑136 (MQ=255)
                                                                                                                     cAGGCGTGACGTTCAGCCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCa  <  2:912315/136‑1 (MQ=255)
                                                                                                                                     |                                                                                                                       
TTTTCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCA  >  NC_000913/3271463‑3271714

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