Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 1,092,469 A→C 21.9% V50G (GTC→GGC)  kdsD ← Arabinose 5‑phosphate isomerase KdsD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1221,092,4690AC21.9% 80.7 / 4.8 53V50G (GTC→GGC) kdsDArabinose 5‑phosphate isomerase KdsD
Reads supporting (aligned to +/- strand):  ref base A (17/22);  new base C (10/2);  total (29/24)
Fisher's exact test for biased strand distribution p-value = 2.15e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATCTTGTTGCCGATGTGCCCGGACTTGCCCATGCCGACCACGACGACCCGGCCCTTGCTGGCCAGAATCAGTTCGCAGGCTTTGACGAAATTGTCGTCGATGCGCGCCAAAAGGGCCTCTACGGCCTCGAGTTCCAGGCGCAGAGTGCGTTGGGCGGATTGGATCAGCTC  >  ACB122/1092340‑1092602
                                                                                                                                  |                                                                                                                                     
tcatcaCGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGGGCTCTTGTTGCCGATGTGCCCGGCCTTGGCCCTGCCGGCCccg                                                                                                                                  >  1:302565/1‑135 (MQ=255)
tcatcaCGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGGGCTTGTTGCCGATGTGCCCGGCCTGGCCCCTGCCCGCCacc                                                                                                                                  >  6:247993/1‑134 (MQ=255)
tcatcaCGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGCTCTTGTTGCCGCTGTGGCCGGGCTTGCCCCTGCAGAccccc                                                                                                                                  >  1:37612/1‑132 (MQ=255)
tcatcaCGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATCTTGTTGCCGATGTGCCCGGCCGTGCCCATGCCGCCCccg                                                                                                                                  >  6:346277/1‑135 (MQ=255)
     aCGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATCTTGTTGCCGATGTGGCCGGACTTGCCCATGCCGACCACcacgac                                                                                                                             >  2:27904/1‑135 (MQ=255)
     aCGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATCTTGTTGCCGATGTGCCCGGACTTGCCCATGCCGACCacgacgcc                                                                                                                             >  2:356373/1‑133 (MQ=255)
     aCGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATCTTGTTGCCGATGTGCCCGGACTTGCCCATGCCGACCACgacgac                                                                                                                             >  8:82379/1‑135 (MQ=255)
     aCGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATCTTGTTGCCGATGTGCCCGGACTTGCCCATGCCGACCACgacgac                                                                                                                             >  4:437766/1‑135 (MQ=255)
      cGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGCTCTTTTTGCCGATGTGCCCGGACGTGCCCCTGCCGCCCCCGACGccc                                                                                                                            >  5:371398/1‑135 (MQ=255)
      cGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATCTTGTTGCCCATGTGCCCGGGCTTGCCCATGCCGACCACGaccacc                                                                                                                            >  6:424869/1‑135 (MQ=255)
      cGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGCCGATCTGGTTGCCGATGTGCCCGGACTTGCCCATGCCGCCCCCGACGAcc                                                                                                                            >  8:120181/1‑135 (MQ=255)
            aTCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATCTTGTTGCCGATGTGCCCGGACTTGCCCATGCCGCCCACGACGACCCCGccc                                                                                                                      >  5:141635/1‑135 (MQ=255)
               aTGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGCTCTTGTTGCCGATGTGCCCGGACTTGCCCAAGCCGCCCACGCCCACCCGGCCCGTg                                                                                                                   >  2:369007/1‑135 (MQ=255)
                 gCCCAGGCCGCCATGGCGGGCTGCGGCCGGATGCACCAAACACGCCGGGGTGCCGGTGCTGGCCAGGGTCGCGGCGATCTTGTTGCCGATGTGCCCGGACTTGCCCATGCCGACCACGACGACCCGGCCCTTGCt                                                                                                                 <  3:180409/135‑1 (MQ=255)
                    cATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATCTTGTTGCCGATGTGCCCGGCCTTGGCCATGCCGACCACGACGACCCGGCCCTTGGTGGc                                                                                                              >  8:91605/1‑135 (MQ=255)
                    cATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGAGCTTGTTGCCGCTGTGCCCCGACTTGCCCATGCCGACCACGCCGACCCGCCCCTTGCGGGc                                                                                                              >  3:178274/1‑135 (MQ=255)
                    cATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGACGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATCTTGTTGCCGATGTGCCCGGACTTTCCCCAGCCGACCACGACCCCCCCGCCCGTGCAgg                                                                                                               >  8:441675/1‑134 (MQ=255)
                                        cGGCCGGATGCACGAAAAACGCCGGTGTCCCGGTGCTCGCCAGGGTCGCGGCGATCTTGTTGCCGATGTGCCCGGACATGCCCATGCCCCCCCCCACCCCCCCGCCCATGCTGGCCACGATCAGTTCGCCGGcgt                                                                                          >  6:188142/1‑133 (MQ=255)
                                              ccgccACTAAACACGCCGGGGTGCCGGTGCTGGCCTGGGGCGCGGCGATCTTGTTGCCGATGTGCCCGGACTGGCCCATGCCGACCACGACGACCCGGCCCTTGCTGGCCAGAATCAGTTCGCAGGCTTTGACGa                                                                                    <  1:397420/131‑1 (MQ=255)
                                                  cACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGAGCTTGTTGCCGATGTGGCCGGGCTTGCCCCTCCCGACCACGACGCCCCGCGCCCTGGGGGCCAGAATCAGTTCGCAGGCTTTGACGATAtt                                                                                >  2:197712/1‑135 (MQ=255)
                                                   aCGAAAAACGCCGGTATGCCGGTGCTGGCCAGGGTCGCGGCGCTCGTGTTTCCGCTGTGCCCGGCCTTGCCCCTGCCGCCCACCACGCCCCGGCCCAGGCGGGGCAGGAGCAGGTCGGAGGCTTTTGCGAAATTg                                                                               >  1:184465/1‑135 (MQ=255)
                                                       aaaaCGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATATTGTTGCCGATGTTGCCGGCCTTGCCCATGCCGCCCCCCACGCCCCCGCCCCTGCTGGGCAGAAAGAGGTGGCCGGGGTTGACGAAAttttagg                                                                           >  5:363023/1‑130 (MQ=255)
                                                         aaCGCCGGTGTGCCGGTGCTGGCCAGGGGCGCGGCGATCTTGTTGCCGATGTGCCCGGACTTGACCATGCCGACCACGACGACCCGGCCCTTGCTGGCCAGAATCAGTTCGCAGGCTTTGACGAAATTGtcgtcg                                                                         <  4:65440/135‑1 (MQ=255)
                                                                  gtgCCGGTGCTGGCCAGGGTCGCGGCGTTCTTGTTGCGGTTGGGCCCGGAGTTGCCCATGCCTACCACGACGACCCGGCCCTTGCTGGCCAGAATCAGTTCGCAGGCTTTGACGAAATTGTCGTCGATGCGCGcc                                                                <  3:202078/135‑1 (MQ=255)
                                                                    gCCGGTGCTGGCCAGGGTCGCGGCGATCTTGTTGCCGATGTGCCCGGACTTGCCCATGCCGACCACGACGACCCGGCCCTTGCTGGCCAGAATCAGTTCGCAGGCTTTGACGAAATTGTCGTCGATGCGCGCCaa                                                              <  7:82379/135‑1 (MQ=255)
                                                                    gCCGGTGCTGGCCAGGGTCGCGGCGATCTTGTTGCCGATGTGCCCGGACTTGCCCATGCCGACCACGACGACCCGGCCCTTGCTGGCCAGAATCAGTTCGCAGGCTTTGACGAAATTGTCGTCGATGCGCGCCaa                                                              <  3:437766/135‑1 (MQ=255)
                                                                     ccGGTGCTGGCCAGGGTCGCGGCGATCTTGTTGCCGATGTGCCCGGACTTGCCCATGCCGACCACGACGACCCGGCCCTTGCTGGCCAGAATCAGTTCGCAGGCTTTGACGAAATTGTCGTCGATGCGCGCCaaa                                                             <  7:91605/135‑1 (MQ=255)
                                                                       ggTGCTGGCCAGGGTCGCGGCGATCTTGTTGCCGATGTGCCCGGACTTGCCCATGCCGACCACGACGACCCGGCCCTTGCTGGCCAGAATCAGTTCGCAGGCTTTGACGAAATTGTCGTCGATGCGCGCCAAAAg                                                           <  7:20766/135‑1 (MQ=255)
                                                                       ggTGCTGGCCAGGGTCGCGGCGATCTTGTTGCCGATGTGCCCGGACTTGCCCATGCCGACCACGACGACCCGGCCCTTGCTGGCCAGAATCAGTTCGCAGGCTTTGACGAAATTGTCGTCGATGCGCGCCAAAAg                                                           <  8:316790/135‑1 (MQ=255)
                                                                       ggTGCGTGCCAGAGTCGCGGCGAGCTGGTTGCGGATGGGCCCGGACTTGCCGATGCCGACCACGACGACCCGGCCCTTGCTGGCCAGAATCAGTTCGCAGGCTTTGACGAAATTGTCGTCGATGCGCGCCAAAAg                                                           <  5:308880/135‑1 (MQ=255)
                                                                            tGGCCCGGGGCGCGGCGATCTTGTGGCCTTGGGGCCCGTCCTGGCCCATGCCCCCCCCGACGACCCGGCCCTTGCTGGCCAGAATCAGTGCGCAGGCTTTGACGAAATTGTCGTCGATGCGCGCCAAAAGGGCct                                                      <  6:22928/135‑1 (MQ=255)
                                                                            tGGCCAGGGGCGCGGCGATCTTGTTGCCGATGTGCCCGGACTTGCCCATGCCGACCACGACGACCCGGCCCTTGCTGGCCAGAATCAGTTCGCAGGCTTTGACGAAATTGTCGTCGATGCGCGCCAAAAGGGCct                                                      <  1:369007/135‑1 (MQ=255)
                                                                              gCCTGGGTCGGGGGGATTTTGTGCCCTTTGTGCCGGCACTTGCCCATCACGACCCCGACGACCCGGCCCTTGCTGGCCAGAATCAGTTCGCAGGCTTTGACGAAATTGTCGTCGATGCGCGCCAAAAGGGCctct                                                    <  5:424869/135‑1 (MQ=255)
                                                                              gCCCGGGGCGCGGCTAGCTTGTTGCCGATGTGCCCGGACTTGCCCATGCCGACCACGACGACCCGGCCCTTGCTGGCCAGAATCAGTTCGCAGGCTTTGACGAAATTGTCGTCGATGCGCGCCAAAAGGGCctct                                                    <  1:356373/135‑1 (MQ=255)
                                                                                cAGGGTCGCGGGGATCGTGGTGGCGAAGTGGCCGGACTTGCCCATGCCGCCCCCGACGCCCCCGCCCCTGCGGGCCAGAATCAGTTCGCAGGCGTGGACGAAATTGTCGTCGATGCGCGCCAAAACGGCCTCTAc                                                  >  2:307129/1‑135 (MQ=255)
                                                                                 aGGGTCGCGGCGATCTTGTTGCCGCTGTGCCCGGACTTGCCCATGCCGCCCCCGACGACCCGGCCCGTGCTGGGCAGAATCAGTTCGCAGGCTTTGACGAAATTGTCGTCGCTGCGCGCCAAAAGGGCCTCTTCg                                                 >  6:234613/1‑135 (MQ=255)
                                                                                     tCGCGGCGATCTTGTTGCCGAGGTGCCCGGACTTGCCCATGCCGACCACGACGACCCGGCCCTTGCTGGCCAGAATCAGTTCGCAGGCTTTGACGAAATTGTCGTCGATGCGCGCCAAAAGGGCCTCTACGGCCt                                             <  7:441675/135‑1 (MQ=255)
                                                                                      cgcgGCGATCTTGTTGCCGCTGTGCCCGGACTTGCCCATGCCGACCACGACGCCCCGGCCCGTGCTGGCCAGCATCAGTTCGCAGGCTTTGACGAAATTGTCGTCGATGCGCGCCAAAAGGGCCTCTACGGccgc                                            >  7:313218/1‑133 (MQ=255)
                                                                                                      ccGATGGGCCCGGACTTGCCCATGCCGACCACGACGACCCGGCCCTTGCTGGCCAGAATCAGTTCGCAGGCTTTGACGAAATTGTCGTCGATGCGCGCCAAAAGGGCCTCTACGGCCTCGAGTTCCAGGCGCag                             <  2:78529/134‑1 (MQ=255)
                                                                                                             gcgcGGACTTGCCATGCCGGAGCACGACGAGCGGGCCCTTGCTGGCAAGAATCAGTTCGCTGGGTGTGCCGAAATTGTCGTCGATGCGCGCCAAAAGGGCTTCTACGGCCTCGAGTTCCAGGCGCAGAGTGCGtt                     <  4:128539/132‑1 (MQ=255)
                                                                                                             gCCCGGCCTGGCCCATGCCGACCACCACGACCCGGCCCTTGCTGGCCAGAATCAGTTCGCAGGCTTTGACGAAATTGTCGTCGATGCGCGCCAAAAGGGCCTCTACGGCCTCGAGTTCCAGGCGCAGAGTGCGtt                     <  2:184465/135‑1 (MQ=255)
                                                                                                             gCCCGGACTTGCCCATGCCGACCACGACGACCCGGCCCTTGCTGGCCAGAATCAGTTCGCAGGCTTTGACGAAATTGTCGTCGATGCGCGCCAAAAGGGCCTCTACGGCCTCGAGTTCCAGGCGCAGAGTGCGtt                     <  3:183870/135‑1 (MQ=255)
                                                                                                             gCCCGCACTCCCCAATGCCGACCACGACGACCCGGCCCTTGCTGGCCAGAATCAGTTCGCAGGCTTTGACGAAATTGTCGGCGATGCGCGCCAAAAGGGCCTCTACGGCCGCGAGTTCCAGGCGCAGAGTGCGtt                     <  1:451951/135‑1 (MQ=255)
                                                                                                                  gACTTGCCCTTGCCGCCCACGACGACCCGGCCCTTGCTGGCCAGAATCAGTTCGCAGGCTTTGACGAAATTGTCGTCGATGCGCGCCAAAAGGGCCTCTACGGCCTCGAGTTCCAGGCGCAGAGTGCGTTGGGCg                <  3:327355/135‑1 (MQ=255)
                                                                                                                         ccATGCCGACCACGACGACCCGGCCCTTGCTGGCCAGAATCAGTTCGCAGGCTTTGACGAAATTGTCGTCGATGCGCGCCAAAAGGGCCTCTACGGCCTCGAGTTCCAGGCGCAGAGTGCGTTGGGCGGATTGGa         >  4:109725/1‑135 (MQ=255)
                                                                                                                         ccATGCCGACCACGACGACCCGGCCCTTGCTGGCCAGAATCAGTTCGCAGGCTTTGACGAAATTGTCGTCGATGCGCGCCAAAAGGGCCTCTACGGCCTCGAGTTCCAGGCGCAGAGTGCGTTGGGCGGATTGGa         >  1:274725/1‑135 (MQ=255)
                                                                                                                          cATGCCGACCACGACGACCCGGCCCTTGCTGGCCAGAATCAGTTCGCAGGCTTTGACGAAATTGTCGTCGATGCGCGCCAAAAGGGCCTCTACGGCCTCGAGTTCCAGGCGCAGAGTGCGTTGGGCGGATTGGAt        >  3:276675/1‑135 (MQ=255)
                                                                                                                               cccaCCCCGACGACCCGGCCCTTGCTGGCCAGAATCAGTTCGCAGGCTTTGACGAAATTGTCGTCGATGCGCGCCAAAAGGGCCTCTACGGCCTCGAGTTCCAGGCGCAGAGTGCGTTGGGCGGATTGGATCAGc    <  5:254686/132‑1 (MQ=255)
                                                                                                                               cGACCACGACGACCCGGCCCTTGCTGGCCAGAATCAGTTCGCAGGCTTTGACGAAATTGTCGTCGATGCGCGCCAAAAGGGCCTCTACGGCCTCGAGTTCCAGGCGCAGAGTGCGTTGGGCGGATTGGATCAGc    <  8:240648/134‑1 (MQ=255)
                                                                                                                               cGACCACGACGACCCGGCCCTTGCTGGCCAGAATCAGTTCGCAGGCTTTGACGAAATTGTCGTCGATGCGCGCCAAAAGGGCCTCTACGGCCTCGAGTTCCAGGCGCAGAGTGCGTTGGGCGGATTGGATCAGCt   <  8:117524/135‑1 (MQ=255)
                                                                                                                                 gACCACGACGACCCGGCCCTTGCTGGCCAGAATCAGTTCGCAGGCTTTGACGAAATTGTCGTCGATGCGCGCCAAAAGGGCCTCTACGGCCt                                             >  2:280249/1‑92 (MQ=255)
                                                                                                                                 gACCACGACGACCCGGCCCTTGCTGGCCAGAATCAGTTCGCAGGCTTTGACGAAATTGTCGTCGATGCGCGCCAAAAGGGCCTCTACGGCCTCGAGTTCCAGGCGCAGAGTGCGTTGGGCGGATTGGATCAGCTc  >  2:102406/1‑135 (MQ=255)
                                                                                                                                 gACCACGACGACCCGGCCCTTGCTGGCCAGAATCAGTTCGCAGGCTTTGACGAAATTGTCGTCGATGCGCGCCAAAAGGGCCTCTACGGCCTCGAGTTCCAGGCGCAGAGTGCGTTGGGCGGATTGGATCAGCTc  >  5:435417/1‑135 (MQ=255)
                                                                                                                                 gACCACGACGACCCGGCCCTTCCTGGCCAGAATCAGTTCGCAGGCTTTGACGAAATTGTCGTCGATGCGCGCCAAAAGGGCCTCTACGGCCt                                             <  1:280249/92‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                     
TCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATCTTGTTGCCGATGTGCCCGGACTTGCCCATGCCGACCACGACGACCCGGCCCTTGCTGGCCAGAATCAGTTCGCAGGCTTTGACGAAATTGTCGTCGATGCGCGCCAAAAGGGCCTCTACGGCCTCGAGTTCCAGGCGCAGAGTGCGTTGGGCGGATTGGATCAGCTC  >  ACB122/1092340‑1092602

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