Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,032,086 G→A 8.3% A53A (GCC→GCT)  kdsA‑II ← 2‑dehydro‑3‑deoxyphosphooctonate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,032,0860GA 8.3% 133.7 / 2.4 48A53A (GCC→GCT) kdsA‑II2‑dehydro‑3‑deoxyphosphooctonate aldolase
Reads supporting (aligned to +/- strand):  ref base G (23/21);  new base A (2/2);  total (25/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.32e-01

GGGCCGTTTGGTAGATTTCGTGAACGTCGGTGATTATCGGCACGCCAAATTCTGCCTTCACCTTTTCAAAGATCCGCAGCCCTTCTTCCATGCCGGGGCCACGGTAGGAATGGATCGATGAGCGATTGGCTTTA‑TCAAAGCTGGCTTTGAACACATACGGGATGCCCAGCTTCTGCGTGACGCGCACATATTCTGCACAGGCAGTAAGGGCCAGGTCTTCAGACTCGAGCACGTTGATGCCGCCGAATAGCACAAAAGGAGCCGAGTTGGAGCATTCGATTGTGGA  >  NC_002947/2031944‑2032229
                                                                                                                                               |                                                                                                                                               
gggCCGTTTGGTAGATTTCGTGAACGTCGGTCATAATCGGCACCCCAAATTATCCATACACCTTTTCAAAAAACCGCAGCCATACTTCCATGCCGGGGCAACGGTAGGAATGGAACGAAGAGCGATTGGAGTTA‑TCAAACATGGGttt                                                                                                                                            >  6:248349/1‑148 (MQ=255)
           tAGATTTCGTGAACGTCGGTGATTATCGGCACGCCAAATTCTGCCTTCACCTTTTCAAAGATCCGCAGCCCTTCTTCCATGCCGGGGCCACGGTAGGAATGGATCGATGAGCGATTGGCTTTA‑TCAAAGCTGGCTTTGAACACATACgg                                                                                                                                >  2:10611/1‑149 (MQ=255)
                   gTGAACGTCGGTGATTATCGGCACGCCAAATTCTGCCTTCACCTTTTCAAAGATCCGCAGCCCTTCTTCCATGCCGGGGCCACGGTAGGAATGGATCGATGAGCGATTGGCTTTA‑TCAAAGCTGGCTTTGAACACATACGGGATGCCCa                                                                                                                        <  7:98731/149‑1 (MQ=255)
                   gTGAACGTCGGTGATTATCGGCACGCCAAATTCTGCCTTCACCTTTTCAAAGATCCGCAGCCCTTCTTCCATGCCGGGGCCACGGTAGGAATGGATCGATGAGCGATTAGCTTTA‑TCAAAGCTGGCTTTAAACACATACGGGATGcc                                                                                                                          >  3:294447/1‑147 (MQ=255)
                      aaCGTCGGTGATTATCGGCACGCCAAATTCTGCCTTCACCTATTCAAAGATCCGCAGCCCTTCTTCCATGCCGGGGCCACGGTAGGAATGGATCGATGAGCGATTGGCTTTA‑TCAAAGCTGGCTTTGAACACATACGGGATGCCCAGCt                                                                                                                     <  5:306177/149‑1 (MQ=255)
                          tCGGTGATTATCGGCACGGCAAATTCTGCCTTCACCTTTTCAAAGATCCGCAGCCCTTCTTCCATGCCGGGGCCACGGTAGGAAGGGATCGATGAGCGATTGGCTTTA‑TCTAAGCTGGCTTTGAACACATACGGGATGCCCAGCTTCTg                                                                                                                 <  5:248349/149‑1 (MQ=255)
                                            ccAAATTCTGCCTTCACCTTTTCAAAGATCCGCAGCCCTTCTTCCATGCCGGGGCCACGGTAGGAATGGATCGATGAGCGATTGGCTTTA‑TCAAAGCTGGCTTTGAACACATACGGGATGCCCAGCTTCTGCGTGACGCGCACATATTc                                                                                               >  6:47138/1‑149 (MQ=255)
                                             cAAATTCTGCCTTCACCTTTTCAAAGATCCGCAGCCCTTCTTCCATGCCGGGGCCACGGTAGGAATGGATCGATGAGCGATTGGCTTTA‑TCAAAGCTGGCTTTGAACACATACGGGATGCCCAGCTTCTGCGTGACGCGCACATATTCt                                                                                              >  1:240976/1‑149 (MQ=255)
                                              aaaTTCTGCCTTCACCTTTTCAAAGATCCGCAGCCCTTCTTCCATGCCGGGGCCACGGTAGGAATGGATCGATGAGCGATTGGCTTTA‑TCAAAGCTGGCttt                                                                                                                                            <  1:265199/102‑1 (MQ=255)
                                              aaaTTCTGCCTTCACCTTTTCAAAGATCCGCAGCCCTTCTTCCATGCCGGGGCCACGGTAGGAATGGATCGATGAGCGATTGGCTTTA‑TCAAAGCTGGCttt                                                                                                                                            >  2:265199/1‑102 (MQ=255)
                                                   cTGCCTTCACCTTTTCAAAGATCCGCAGCCCTTCTTCCATGCCGGGGCCACGGTAGGAATGGATCGATGAGCGATTGGCTTTA‑TCAAAGCTGGCTTTGAACACATACGGGATGCCCAGCTTCTGCGTGACGCGCACATATTCTGcaca                                                                                         <  5:47138/148‑1 (MQ=255)
                                                       cTTCACCTTTTCAAAGATCCGCAGCCCTTCTTCCATGCCGGGGCCACGGTAGGAATGGATCGATGAGCGATTAGCTTTA‑TCAAAGCTGGCTTTAAACACATACGGGATGCCCAGCTTCTGCGTGACGCGCACATATTCTGCACAGGCAg                                                                                    <  4:294447/149‑1 (MQ=255)
                                                        ttCACCTTTTCAAAGATCCGCAGCCCTTCTTCCATGCCGGGGCCACGGTAGGAATGGATCGATGAGCGATTGGCTTTA‑TCAAAGCCAGCtt                                                                                                                                             >  7:238952/1‑91 (MQ=255)
                                                        ttCACCTTTTCAAAGATCCGCAGCCCTTCTTCCATGCCGGGGCCACGGTAGGAATGGATCGATGAGCGATTGGCTTTA‑TCAAAGCCAGCtt                                                                                                                                             <  8:238952/91‑1 (MQ=255)
                                                                  cAAAGATCCGCAGCCCTTCTTCCATGCCGGGGCCACGGTAGGAATGGATCGATGAGCGATTGGCTTTA‑TCAAAGCTGGCTTTGAACACATACGGGATGCCCAGCTTCTGCGTGACGCGCACATATTCTGCACAGGCAGt                                                                                   >  8:31012/1‑139 (MQ=255)
                                                                  cAAAGATCCGCAGCCCTTCTTCCATGCCGGGGCCACGGTAGGAATGGATCGATGAGCGATTGGCTTTA‑TCAAAGCTGGCTTTGAACACATACGGGATGCCCAGCTTCTGCGTGACGCGCACATATTCTGCACAGGCAGt                                                                                   <  7:31012/139‑1 (MQ=255)
                                                                        tCCGCAGCCCTTCTTCCATGCCGGGGCCACGGTAGGAATGGATCGATGAGCGATTGGCTTTA‑TCAAAGCTGGCTTTGAACACATACGGGATGCCCAGCTTCTGCGTGACGCGCACATATTCTGCACAGGCAGTAAGGGCCAGGTCTTCa                                                                   <  6:143711/149‑1 (MQ=255)
                                                                                ccTTCTTCCATGCCGGGGCCACGGTAGGAATGGATCGATGAGCGATTGGCTTTA‑TCAAAGCTGGCTTTGAACACATACGGGATGCCCAGCTTCTGCGTGACGCGCACATATTCTGCACAGGCAGTAAGGGCCAGGTCTTCAGACTCGAg                                                           <  8:202345/149‑1 (MQ=255)
                                                                                ccTTCTTCCATGCCGGGGCCACGGTAGGAATGGATCGATGAGCGATTGGCTTTA‑TCAAAGCTGGCTTTGAACACATACGGGATGCCCAGCTTCTGCGTGACGCGCACATATTCTGCACAGGCAGTAAGGGCCAGGTCTTCAGACTCGAg                                                           <  4:55365/149‑1 (MQ=255)
                                                                                 cttcttCCATGCCGGGGCCACGGTAGGAATGGATCGATGAGCGATTGGCTTTA‑TCAAAGCTGGCTTTGAACACATACGGGATGCCCAGCTTCTGCGTGACGCGCACATATTCTGCACAGGCAGTAAGGGCCAGGTCTTCAGACTCGAg                                                           >  5:161380/1‑148 (MQ=255)
                                                                                    cttcCATGCCGGGGCCACGGTAGGAATGGATCGATGAGCGATTGGCTTTA‑TCAAAGCTGGCTTTGAACACATACGGGATGCCCAGCTTCTGCGTGACGCGCACATATTCTGCACAGGCAGTAAGGGCCAGGTCTTc                                                                    >  6:283172/1‑136 (MQ=255)
                                                                                    cttcCATGCCGGGGCCACGGTAGGAATGGATCGATGAGCGATTGGCTTTA‑TCAAAGCTGGCTTTGAACACATACGGGATGCCCAGCTTCTGCGTGACGCGCACATATTCTGCACAGGCAGTAAGGGCCAGGTCTTc                                                                    <  5:283172/136‑1 (MQ=255)
                                                                                    cttcCATGCCGGGGCCACGGTAGGAATGGATCGATGAGCGATTGGCTTTA‑TCAAAGCTGGCTTTGAACACATACGGGATGCCCAGCTTCTGCGTGACGCGCACATATTCTGCACAGGCAGTAAGGGCCAGGTCTTCAGACTCGAGCAc                                                        <  1:10611/148‑1 (MQ=255)
                                                                                         aTGCCGGGGCCACGGTAGGAATGGATCGATGAGCGATTGGCTTTA‑TCAAAGCTGGCTTTGAACACATACGGGATGCCCAGCTTCTGCGTGACGCGCACATATTCTGCACTGGCAGTAAGGGCCAGGTCTTCAGACTCGAGCACGTTGAt                                                  <  2:240976/149‑1 (MQ=255)
                                                                                           gCCGGGGCCACGGTAGGAATGGATCGATGAGCGATTGGCTTTAGTCAAAGCTGGCTTTGAACACATACGGCATGCCCAGCTTCTGCATGACGCGCACATATTCTGCACAGGCAGTAAGGGCCAGGTCTTCAGACTCGACCACGTTGATg                                                 >  1:68529/1‑149 (MQ=255)
                                                                                           gCCGGGGCCACGGTAGGAATGGATCGATGAGCGATTGGCTTTA‑TCAAAGCTGGCTTTGAACACATACGGGATGCCCAGCTTCTGCGTGACGCGCACATATTCTGCACAGGCAGTAAGGGCCAGGTCTTCAGACTCGAGCACGTTGATgc                                                >  1:54162/1‑149 (MQ=255)
                                                                                                           gAATGGATCGATGAGCGATTGGCTTTA‑TCAAAGCTGGCTTTGAACACATACGGGATGCCCAGCTTCTGCGTGACGCGCACATATTCTGCACAGGCAGTAAGGGCCAGGTCTTCAGACTCGAGCACGTTGATGCCGCCGAATAGCACaa                                 >  5:14590/1‑148 (MQ=255)
                                                                                                                 aTCGATGAGCGATTGGCTTTA‑TCAAAGCTGGCTTTGAACACATACGGGATGCCCAGCTTCTGCGTGACGCGCACATATTCTGCACAGGCAGTAAGGGCCAGGTCTTCAGACTCGAGCACGTTGATGCCGCCGAATAGCACAAAAGGAGc                          >  4:15812/1‑149 (MQ=255)
                                                                                                                  tCGATGAGCGATTGGCTTTA‑TCAAAGCTGGCTTTGAACACATACGGGATGCCCAGCTTCTGCGTGACGCGCACATATTCTGCACAGGCAGTAAGGGCCAGGTCTTCAGACTCGAGCACGTTGATGCCGCCGAATAGCACAAAAGGAGcc                         >  7:86682/1‑149 (MQ=255)
                                                                                                                            aTTGGCTTTA‑TCAAAGCTGGCTTTGAACACATACGGGATGCCCAGCTTCTGCGTGACGCGCACATATTCTGCACAGGCAGTAAGGGCCAGGTCTTCAGACTCGAGCACGTTGATGCCGCCGAATAGCACAAAAGGAGCCGAGTTGGAGc               <  2:54162/149‑1 (MQ=255)
                                                                                                                                gCTTTA‑TCAAAGCTGGCTTTGAACACATACGGGATGCCCAGCTTCTGCGTGACGCGCACATATTCTGCACAGGCAGTAAGGGCCAGGTCTTCAGACTCGAGCACGTTGATGCCGCCGAATAGCACAAAAGGAGCCGAGTTGGAGCAtt            <  6:14590/148‑1 (MQ=255)
                                                                                                                                 cTTTA‑TCAAAGCTGGCTTTGAACACATACGGGATGCCCAGCTTCTGCGTGACGCGCACATATTCTGCACAGGCAGTAAGGGCCAGGTCTTCAGACTCGAGCACGTTGATGCCGCCGAAt                                        >  4:168604/1‑119 (MQ=255)
                                                                                                                                  tttA‑TCAAAGCTGGCTTTGAACACATACGGGATGCCCAGCTTCTGCGTGACGCGCACATATTCTGCACAGGCAGTAAGGGCCAGGTCTTCAGACTCGAGCACGTTGATGCCGCCGAATAGCACAAAAGGAGCCGAGTTGGAGCATTCGa         >  1:19341/1‑149 (MQ=255)
                                                                                                                                    tgttCAAAGCCAGCTTTGAACACATACGGGATGCCCAGCTTCTGCGTGACGCGCACATATTCTGCACAGGCAGTAAGGGCCAGGTCTTCa                                                                   <  8:307803/87‑1 (MQ=255)
                                                                                                                                    tgttCAAAGCCAGCTTTGAACACATACGGGATGCCCAGCTTCTGCGTGACGCGCACATATTCTGCACAGGCAGTAAGGGCCAGGTCTTCa                                                                   >  7:307803/4‑90 (MQ=255)
                                                                                                                                         aaaGCTGGCTTTGAACACATACGGGATGCCCAGCTTCTGCGTGACGCGCACATATTCTGCACAGGCAGTAAGGGCCAGGTCTTCAGACTCGAGCACGTTGATGCCGCCGAATAGCACaaaa                               >  2:208683/1‑121 (MQ=255)
                                                                                                                                         aaaGCTGGCTTTGAACACATACGGGATGCCCAGCTTCTGCGTGACGCGCACATATTCTGCACAGGCAGTAAGGGCCAGGTCTTCAGACTCGAGCACGTTGATGCCGCCGAATAGCACaaaa                               <  1:208683/121‑1 (MQ=255)
                                                                                                                                         aaaGCTGGCTTTGAACACATACGGGATGCCCAGCTTCTGCGTGACGCGCACATATTCTGCACAGGCAGTAAGGGCCAGGTCTTCAGACTCGAGCACGTTGATGCCGCCGAATAGCACAAAAGGAGCCGAGTTGGAGCATTCGATtgtg    >  7:78754/1‑148 (MQ=255)
                                                                                                                                         aaaGCTGGCTTTGAACACATACGGGAGGCCCAGCTTCTGCGTGACGCGCACATATTCTGCACAGGCAGTAAGGGCCAGGTCTTCAGACTCGAGCACGTTGATGCCGCCGAATAGCACAAAAGGAGCCGAGTTGGAGCATTCGATtgtg    >  1:186209/1‑148 (MQ=255)
                                                                                                                                          aaGCTGGCTTTGAACACATACGGGATGCCCAGCTTCTGCGTGACGCGCACATATTCTGCACAGGCAGTAAGGGCCAGGTCTTCAGACTCGAGCACGTTGATGCCGCCGAATAGCACAAAAGGAGCCGAGTTGGAGCATTCGATTGTGGa  <  3:15812/149‑1 (MQ=255)
                                                                                                                                            gCTGGCTTTGAACACATACGGGATGCCCAGCTTCTGCGTGACGCGCACATATTCTGCACAGGCAGTAAGGGCCAGGTCTTCAGACTCGAGCACgg                                                      <  7:229558/95‑2 (MQ=255)
                                                                                                                                            gCTGGCTTTGAACACATACGGGATGCCCAGCTTCTGCGTGACGCGCACATATTCTGCACAGGCAGTAAGGGCCAGGTCTTCAGACTCGAGCACGt                                                      >  8:229558/1‑95 (MQ=255)
                                                                                                                                             cTGGCTTTGAACACATACGGGATGCCCAGCTTCTGCGTGACGCGCACATATTCTGCACAGGCAGTAAGGGCCAGGTCTTCAGACTCGAGCACGTTGATGCCGCCGAATAGCACAAAAGGAGCCGAGTTGGAGCATTCGATTGTGGa  <  7:323443/146‑1 (MQ=255)
                                                                                                                                             cTGGCTTTGAACACATACGGGATGCCCAGCTTCTGCGTGACGCGCACATATTCTGCACAGGCAGTAAGGGCCAGGTCTTCAGACTCGAGCACGTTGATGCCGCCGAATAGCACAAAAGGAGCCGAGTTGGAGCATTCGATTGTGGa  >  8:323443/1‑146 (MQ=255)
                                                                                                                                              tGGCTTTGAACACATACGGGATGCCCAGCTTCTGCGTGACGCGCACATATTCTGCACAGGCAGTAAGGGCCAGGTCTTCa                                                                   <  4:168624/80‑1 (MQ=255)
                                                                                                                                              tGGCTTTGAACACATACGGGATGCCCAGCTTCTGCGTGACGCGCACATATTCTGCACAGGCAGTAAGGGCCAGGTCTTCa                                                                   <  4:163379/80‑1 (MQ=255)
                                                                                                                                              tGGCTTTGAACACATACGGGATGCCCAGCTTCTGCGTGACGCGCACATATTCTGCACAGGCAGTAAGGGCCAGGTCTTCa                                                                   >  3:168624/1‑80 (MQ=255)
                                                                                                                                              tGGCTTTGAACACATACGGGATGCCCAGCTTCTGCGTGACGCGCACATATTCTGCACAGGCAGTAAGGGCCAGGTCTTCa                                                                   >  3:163379/1‑80 (MQ=255)
                                                                                                                                               |                                                                                                                                               
GGGCCGTTTGGTAGATTTCGTGAACGTCGGTGATTATCGGCACGCCAAATTCTGCCTTCACCTTTTCAAAGATCCGCAGCCCTTCTTCCATGCCGGGGCCACGGTAGGAATGGATCGATGAGCGATTGGCTTTA‑TCAAAGCTGGCTTTGAACACATACGGGATGCCCAGCTTCTGCGTGACGCGCACATATTCTGCACAGGCAGTAAGGGCCAGGTCTTCAGACTCGAGCACGTTGATGCCGCCGAATAGCACAAAAGGAGCCGAGTTGGAGCATTCGATTGTGGA  >  NC_002947/2031944‑2032229

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