Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,594,484 A→C 7.8% R313R (CGT→CGG kgtP ← alpha‑ketoglutarate permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,594,4840AC 7.8% 156.3 / 3.8 51R313R (CGT→CGGkgtPalpha‑ketoglutarate permease
Reads supporting (aligned to +/- strand):  ref base A (22/25);  new base C (2/2);  total (24/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTCTGCATGTAGGTGGTGTAGGTG  >  NC_002947/1594341‑1594631
                                                                                                                                               |                                                                                                                                                   
ttGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGAcggcg                                                                                                                                                >  1:225262/1‑149 (MQ=255)
      gAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGcc                                                                                                                                             <  8:101320/146‑1 (MQ=255)
      gAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGcc                                                                                                                                             >  7:101320/1‑146 (MQ=255)
       aGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGCCGCCGTCCGATc                                                                                                                                         <  1:36093/149‑1 (MQ=255)
       aGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATc                                                                                                                                         <  4:152913/149‑1 (MQ=255)
         gTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGACCtt                                                                                                                                       <  2:225262/149‑1 (MQ=255)
               ccGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCgg                                                                                                                                  >  5:49642/1‑148 (MQ=255)
                 gCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACa                                                                                                                               <  2:56621/149‑1 (MQ=255)
                         tgatCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGTTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCTATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCg                                                                                                                       <  8:207513/149‑1 (MQ=255)
                            tCAGTGCGGCCATGATCAGTAAGAACGCACCACAACAGGTCTGGATGGTGTGCAGGGTTCTGATGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGatc                                                                                                                    <  3:81952/149‑1 (MQ=255)
                               gTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGATGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGatcatc                                                                                                                 >  3:91235/1‑149 (MQ=255)
                                     ccATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGc                                                                                                           <  5:174723/149‑1 (MQ=255)
                                      cATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGTCCACCGATCATCGGCTg                                                                                                            >  7:311238/1‑147 (MQ=255)
                                       aTGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAg                                                                                                         <  6:380412/149‑1 (MQ=255)
                                               gaagaaCGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATg                                                                                                  >  2:226508/1‑148 (MQ=255)
                                                 agaaCGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCTAAGGCGATCAGGATCGGACGGCGGCCg                                                                                                                                            <  4:316488/104‑1 (MQ=255)
                                                  gaaCGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATgaaca                                                                                              <  1:226508/149‑1 (MQ=255)
                                                    aCGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATgaacagg                                                                                            <  6:315394/149‑1 (MQ=255)
                                                    aCGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCAACGGCTTCAGGCACATgaacag                                                                                             <  6:49642/148‑1 (MQ=255)
                                                                 ggTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTggcgg                                                                               >  5:18069/1‑149 (MQ=255)
                                                                       gATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGcc                                                                             >  5:381090/1‑145 (MQ=255)
                                                                       gATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGcc                                                                             >  1:351954/1‑145 (MQ=255)
                                                                          ggTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGAtgg                                                                      >  3:307996/1‑149 (MQ=255)
                                                                              tgCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATggtggt                                                                  <  6:342772/149‑1 (MQ=255)
                                                                                       cTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCAc                                                          <  4:307996/148‑1 (MQ=255)
                                                                                         gAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGAtgg                                                                      <  4:24942/134‑1 (MQ=255)
                                                                                         gAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGAtgg                                                                      >  3:24942/1‑134 (MQ=255)
                                                                                          aGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGAt                                                      <  3:261812/149‑1 (MQ=255)
                                                                                           ggATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATg                                                     <  6:381090/149‑1 (MQ=255)
                                                                                            gATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGc                                                    <  2:351954/149‑1 (MQ=255)
                                                                                              tGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCt                                                   >  5:274813/1‑148 (MQ=255)
                                                                                                 ggACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGAt                                                      <  1:155146/142‑1 (MQ=255)
                                                                                                 ggACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACACCGTGGCGGCCGAGATCCTGGTCGAGTCACTGAt                                                      >  2:155146/1‑142 (MQ=255)
                                                                                                  gACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGCTCATCGGCTGCCGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGc                                              >  5:292015/1‑149 (MQ=255)
                                                                                                     ggTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCa                                           >  4:381058/1‑149 (MQ=255)
                                                                                                       agAACAGGGTGCCGAGCACAACGAAGGCGATCAGGATCGGACGGAGTCCGATCTTGACGGACAGACCAACGATCATCGGGTGCAGGCACATGAACATGAACAGCGTGGCGGCAGAGAAGGGGCTCGAGTCACAGATCCGCCTGccaaag                                         >  4:87498/2‑145 (MQ=255)
                                                                                                          aCAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGc                                              <  3:179082/141‑1 (MQ=255)
                                                                                                          aCAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGc                                              >  4:179082/1‑141 (MQ=255)
                                                                                                              ggTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGa                                  <  2:254111/149‑1 (MQ=255)
                                                                                                                    gAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACCGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGAGTTGACCAgg                             >  7:76546/1‑148 (MQ=255)
                                                                                                                          acCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTCt                      >  2:72247/1‑149 (MQ=255)
                                                                                                                            cGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTCTg                     >  5:1982/1‑148 (MQ=255)
                                                                                                                            cGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTCTGc                    <  4:91235/149‑1 (MQ=255)
                                                                                                                             gAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGTCCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTCTGCa                   <  8:311238/149‑1 (MQ=255)
                                                                                                                              aaGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTCTGCAt                  <  1:72247/149‑1 (MQ=255)
                                                                                                                                       aagATCGGCCGCCGTCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTCTGCATGTAggtggt         >  3:16580/3‑149 (MQ=255)
                                                                                                                                       aGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTCTGCATGTAggtggt         <  6:274813/149‑1 (MQ=255)
                                                                                                                                         gATCGGCCGCCGTCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTCTGCATGTAGgtggtg        >  1:268474/1‑148 (MQ=255)
                                                                                                                                         gATCGGCCGCCGTCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTCTGCATGTAGGTggtgt       <  4:16580/149‑1 (MQ=255)
                                                                                                                                              gACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTCTGCATGTAGGTggtgtaggtg  >  7:332668/1‑149 (MQ=255)
                                                                                                                                              gACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACAGATGCTCATGCCCACGGTGTTGACCAGGTATTTCTGCATGTAGGTggtgtaggtg  >  3:331187/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
TTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTCTGCATGTAGGTGGTGTAGGTG  >  NC_002947/1594341‑1594631

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: