Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 263,957 A→T 14.0% E236D (GAA→GAT)  gabD‑I → succinate‑semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947263,9570AT14.0% 129.8 / 13.4 50E236D (GAA→GAT) gabD‑Isuccinate‑semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (24/19);  new base T (3/4);  total (27/23)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.27e-01

CCGCTCTGGCCCTGGTCGAGCTGGCCCACCGTGCCGGCATCCCGGCTGGCGTGCTGAGTGTGGTTACCGGCAGCGCCGGCGAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACG  >  NC_002947/263812‑264105
                                                                                                                                                 |                                                                                                                                                    
ccGCTCTGGCCCTGGTCGAGCTGGCCCACCGTGCCGGCATCCCGGCTGGCGTGCTGAGTGTGGTTACCGGCAGCGCCGGCGAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCg                                                                                                                                                  <  2:647653/150‑1 (MQ=255)
 cGCTCTGGCCCTGGTCGAGCTGGCCCACCGTGCCGGCATCCCGGCTGGCGTGCTGAGTGTGGTTACCGGCAGCGCCGGCGAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCgg                                                                                                                                                 >  1:603794/1‑150 (MQ=255)
   ctctGGCCCTGGTCGAGCTGGCCCACCGTGCCGGCATCCCGGCTGGCGTGCTGAGTGTGGTTACCGGCAGCGCCGGCGAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTc                                                                                                                                               <  2:603794/150‑1 (MQ=255)
     ctGGCCCTGGTCGAGCTGGCCCACCGTGCCGGCATCCCGGCTGGCGTGCTGAGTGTGGTTACCGGCAGCGCCGGCGAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGATTTCGGTCga                                                                                                                                             <  2:543637/150‑2 (MQ=255)
      tGGCCCTGGTCGAGCTGGCCCACCGTGCCGGCATCCCGGCTGGCGTGCTGAGTGTGGTTACCGGCAGCGCCGGCGAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGcc                                                                                                                                            >  2:509612/1‑150 (MQ=255)
               tCGAGCTGGCCAACCGTGCCGGCAACCTCGCTGGGGGGGTGAGGGTGGCTGCCCGCAGCGCCGGCGCAGATGGCGGCGGAGTGACCGTCAACTCCCTGCGACGCAAGCTGTCCTTCACCGGCTCCAGCGCAATGGATAGACAGCtcgtgg                                                                                                                                   >  2:282820/1‑145 (MQ=255)
                 gAGCTGGCCCACCGTGCCGGCATCCCGGCTGGCGTGCTGAGTGTGGTTACCGGCAGCGCCGGCGAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGgaa                                                                                                                                 <  1:509612/150‑1 (MQ=255)
                 gAGCTGGCCCACCGTGCCGGCATCCCGGCTGGCGTGCTGAGTGTGGTTACCGGCAGCGCCGGCGAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGgaa                                                                                                                                 >  1:541197/1‑150 (MQ=255)
                         ccACCGTGCCGGCATCCCGGCTGGCGTGCTGAGTGTGGTTACCGGCAGCGCCGGCGAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATgcgc                                                                                                                         >  2:774910/1‑150 (MQ=255)
                         ccACCGTGCCGGCATCCCGGCTGGCGTGCTGAGTGTGGTTACCGGCAGCGCCGGCGAAGTTGGCGGCGAACTAACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGATTTCGGTCgagc                                                                                                                                           <  1:623800/132‑4 (MQ=255)
                         ccACCGTGCCGGCATCCCGGCTGGCGTGCTGAGTGTGGTTACCGGCAGCGCCGGCGAAGTTGGCGGCGAACTAACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGATTTCGGTCgagc                                                                                                                                           >  2:623800/1‑129 (MQ=255)
                                  cGGCATCCCGGCTGGCGTGCTGAGTGTGGTTACCGGCAGCGCCGGCGAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACAt                                                                                                                >  2:532870/1‑150 (MQ=255)
                                      aTCCCGGCTGGCGTGCTGAGTGTGGTTACCGGCAGCGCCGGCGAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCaag                                                                                                            >  2:242310/1‑150 (MQ=255)
                                       tCCCGGCTGGCGTGCTGAGTGTGGTTACCGGCAGCGCCGGCGAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCaaga                                                                                                           >  1:684272/1‑150 (MQ=255)
                                                gCGTGCTGAGTGTGGTTACCGGCAGCGCCGTCGAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAgg                                                                                                        >  1:544032/1‑144 (MQ=255)
                                                gCGTGCTGAGTGTGGTTACCGGCAGCGCCGGCGAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAgg                                                                                                        <  2:544032/144‑1 (MQ=255)
                                                 cGTGCTGAGTGTGGTTACCGGCAGCGCCGGCGAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCt                                                                                                 <  1:774910/150‑1 (MQ=255)
                                                   tGCTGAGTGTGGTTACCGGCAGCGCCGGCGAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGttt                                                                                                     <  2:98766/144‑1 (MQ=255)
                                                   tGCTGAGTGTGGTTACCGGCAGCGCCGGCGAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGttt                                                                                                     >  1:98766/1‑144 (MQ=255)
                                                         gtgtGGTTACCGGCAGCGCCGGCGAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTgg                                                                                         >  2:694291/1‑150 (MQ=255)
                                                         gtgtGGTTACCGGCAGCGCCGGCGAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGCAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTgg                                                                                         >  1:649827/1‑150 (MQ=255)
                                                               ttACCGGCAGCGCCGGCGAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGATTTCGGTCgagc                                                                                                                                           <  1:60496/94‑4 (MQ=255)
                                                               ttACCGGCAGCGCCGGCGAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGATTTCGGTCgagc                                                                                                                                           >  2:60496/1‑91 (MQ=255)
                                                                 aCCGGCAGCGCCGGCGCATTGGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAAc                                                                                 <  1:242310/150‑1 (MQ=255)
                                                                     gCAGCGCCGGCGAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACTGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAg                                                                                             >  1:438972/1‑134 (MQ=255)
                                                                     gCAGCGCCGGCGAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAg                                                                                             <  2:438972/134‑1 (MQ=255)
                                                                     gCAGCGCCGGCGAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGccc                                                                             >  1:642331/1‑150 (MQ=255)
                                                                          gcCGGCGAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCg                                                                                                                                              >  1:391302/1‑80 (MQ=255)
                                                                          gcCGGCGAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCg                                                                                                                                              <  2:391302/80‑1 (MQ=255)
                                                                                gAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCaagaa                                                                                                          >  1:345451/1‑110 (MQ=255)
                                                                                gAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCaagaa                                                                                                          <  2:345451/110‑1 (MQ=255)
                                                                                   gTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTCCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTc                                                               <  1:282820/150‑1 (MQ=255)
                                                                                    tgttcggcgAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTcg                                                              >  2:596248/5‑150 (MQ=255)
                                                                                         ggcgAACTGACCGGCAACTCCCTGGTACGCAAGCTGTTCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCGGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCgacgac                                                         <  2:278634/150‑1 (MQ=255)
                                                                                                 gACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACAACGCCGACCt                                                 <  2:684272/150‑1 (MQ=255)
                                                                                                  aCCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGAGTTCGACGACGCCGACCTg                                                >  2:526142/1‑150 (MQ=255)
                                                                                                    cGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGa                                              <  1:694291/150‑1 (MQ=255)
                                                                                                    cGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGa                                              >  1:1529/1‑150 (MQ=255)
                                                                                                           tCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCg                                       <  2:541197/150‑1 (MQ=255)
                                                                                                           tCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCg                                       >  1:5215/1‑150 (MQ=255)
                                                                                                                  tACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAgg                                <  2:1529/150‑1 (MQ=255)
                                                                                                                       aaGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGcgcg                           >  2:752142/1‑150 (MQ=255)
                                                                                                                            gTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGatcat                      <  1:532870/150‑1 (MQ=255)
                                                                                                                            gTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGatc                        >  2:557669/1‑148 (MQ=255)
                                                                                                                            gTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGatc                        <  1:557669/148‑1 (MQ=255)
                                                                                                                                   aCCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAg               >  1:726256/1‑150 (MQ=255)
                                                                                                                                     cGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTa             >  2:163272/1‑150 (MQ=255)
                                                                                                                                       gCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTAcc           <  1:752142/150‑1 (MQ=255)
                                                                                                                                        ggcGACCGATTTCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTcga                                                             <  1:720194/97‑1 (MQ=255)
                                                                                                                                        ggcGACCGATTTCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTcga                                                             >  2:720194/3‑99 (MQ=255)
                                                                                                                                                aaaTCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACg  <  1:526142/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CCGCTCTGGCCCTGGTCGAGCTGGCCCACCGTGCCGGCATCCCGGCTGGCGTGCTGAGTGTGGTTACCGGCAGCGCCGGCGAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACG  >  NC_002947/263812‑264105

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