Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,661,664 A→G 8.0% K161E (AAG→GAG)  PP_2333 → GntR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,661,6640AG 8.0% 139.0 / 3.1 50K161E (AAG→GAG) PP_2333GntR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (24/22);  new base G (2/2);  total (26/24)
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.46e-01

ACATCGACGAACTGCGCCGGGTGCTCGACACCCATGAGCGCGATGCCGCGTTCCAGGCCGGCCTGGGCTACTACCAGCAGGAAGGCGACTACGACTTCCATTACCGGATTATTCAGGGTAGCGGCAACCAGACCCTGGTCAAGATGCTCTGCGGCGAGCTGTACCAGCTGGTGCGCATGTACCGTATCCAGTTCTCTGCCACGCCTAACCGGCCACGCCAGGCCTTTGCCGAACACCATCGCATTCTTGATGCCATCGCCGACCGTGACGGCGAACTGGCCGAACTCCT  >  NC_002947/2661524‑2661812
                                                                                                                                            |                                                                                                                                                    
aCATCGACGAACTGCGCCGGGTGCTCGACACCCATGAGCGCGATGCCGCGTTCCAGGCCGGCCTGGGCTACTACCAGCAGGAAGGCGACTACGACTTCCATTACCGGATTATTCAGGGTAGCGGCAACCAGACCCTGGTCAAGATGctc                                                                                                                                              <  5:221547/149‑1 (MQ=255)
     gacgaACTGCGCCGGGTGCTCGACACCCATGAGCGCGATGCCGCGTTCCAGGCCGGCCTGGGCTACTACCAGCAGGAAGGCGACTACGACTTCCATTACCGGATTATTCAGGGTAGCGGCAACCAGACCCTGGTCAAGATGCTCTgcgg                                                                                                                                         >  7:64138/1‑149 (MQ=255)
           cTGCGCCGGGTGCTCGACACCCATGAGCGCGATGCCGCGTTCCAGGCCGGCCTGGGCTACTACCAGCAGGAAGGCGACTACGACTTCCATTACCGGATTATTCAGGGTAGCGGCAACCAGACCCTGGTCAAGATGCTCTGCGGCGAGCt                                                                                                                                   <  4:268102/149‑1 (MQ=255)
              cgcCGGGTGCTCGACACCCATGAGCGCGATGCCGCGTTCCAGGCCGGCCAGGGCTACTACCAGCAGGAAGGCGACAACGACATCCAGTACCGGATTATTCAGGGTAGCGGCAACCAGACCCATCGAAAGATGCTATGCGgagatccga                                                                                                                                 >  8:283337/1‑143 (MQ=255)
               tccGGGTGCTCGACACCCATGAGCGCGATGCCGCGTTCCAGGCCGGCCTGGGCTACTACCAGCAGGAAGGCGACTACGACTTCCATTACCGGATTATTCAGGGTAGCGGCAACCAGACCCTGGTCAAGATGCTCTGCGGCGAGCTGTAc                                                                                                                               <  5:186504/148‑1 (MQ=255)
                     tGCTCGACACCCATGATCTCGATGCCGCGTTCCAGGCCGGCCTGTTCTACTACCAGCAGGAAGGCGACTACGACTTCCATTACCGGATTATTCAGGTTAGCGCCAACCAGACCCTGGTCAAGATTCTCTGCGGCGAGCTGTACCAGCTg                                                                                                                         <  8:143672/149‑1 (MQ=255)
                        tCGACACCCATGAGCGCGATGCCGCGTTCCAGGCCGGCCTGGGCTACTACCAGCAGGAAGGCGACTACGACTTCCATTACCGGATTATTCAGGGTAGCGGCAACCAGACCCTGGTCAAGATGCTCTGCGGCGAGCTGTACCAGCTGGTg                                                                                                                      >  8:259168/1‑149 (MQ=255)
                               ccATGAGCGCGATGCCGCGTTCCAGGCCGGCCTGGGCTACTACCAGCAGGAAGGCGACTACGACTTCCATTACCGGATTATTCAGGGTAGCGGCAACCAGACCCTGGTCAAGATGCTCTGCGGCGAGCTGTACCAGCTGGTGCGCATg                                                                                                                >  1:137944/1‑148 (MQ=255)
                                cATGAGCGCGATGCCGCGTTCCAGGCCGGCCTGGGCTACTACCAGCAGGAAGGCGACTACGACTTCCATTACCGGATTATTCAGGGTAGCGGCAACCAGACCCTGGTCAAGATGCTCTGCGGCGAGCTGTACCAGCTGGTGCGCATGTa                                                                                                              <  1:186030/149‑1 (MQ=255)
                                       gcgATGCCGCGTTCCAGGCCGGCCTGGGCTACTACCAGCAGGAAGGCGACTACGACTTCCATTACCGGATTATTCAGGGTAGCGGCAACCAGACCCTGGTCAAGATGCTCTGCGGCGAGCTGTAcc                                                                                                                              <  6:144900/126‑1 (MQ=255)
                                       gcgATGCCGCGTTCCAGGCCGGCCTGGGCTACTACCAGCAGGAAGGCGACTACGACTTCCATTACCGGATTATTCAGGGTAGCGGCAACCAGACCCTGGTCAAGATGCTCTGCGGCGAGCTGTAcc                                                                                                                              >  5:144900/1‑126 (MQ=255)
                                        cgATGCCGCGTTCCAGGCCGGCCTGGGCTACTACCAGCAGGAAGGCGACTACGACTTCCATTACCGGATTATTCAGGGTAGCGGCAACCAGACCCTGGTCAAGATGCTCTGCGGCGAGCTGTACCAGCTGGTGCGCATGTACCGTATcc                                                                                                      <  3:141766/149‑1 (MQ=255)
                                            gCCGCGTTCCAGGCCGGCCTGGGCTACTACCAGCAGGAAGGCGACTACGACTTCCATTACCGGATTATTCAGGGTAGCGGCAACCAGACCCTGGTCAAGATGCTCTGCGGCGAGCTGTACCAGCTGGTGCGCATGTACCGTATCCAGt                                                                                                   >  3:190059/1‑148 (MQ=255)
                                             ccGCGTTCCAGGCCGGCCTGGGCTACTACCAGCAGGAAGGCGACTACGACTTCCATTACCGGATTATTCAGGGTAGCGGCAACCAGACCCTGGTCAAGATGCTCTGCGGCGAGCTGTACCAGCTGGTGCGCATGTACCGTATCCAGtt                                                                                                  <  1:129994/148‑1 (MQ=255)
                                             ccGCGTTCCAGGCCGGCCTGGGCTACTACCAGCAGGAAGGCGACTACGACTTCCATTACCGGATTATTCAGGGTAGCGGCAACCAGACCCTGGTCAAGATGCTCTGCGGCGAGCTGTACCAGCTGGTGCGCATGTACCGTATCCAGTtc                                                                                                 >  8:218234/1‑149 (MQ=255)
                                                cgTTCCAGGCCGGCCTGGGCTACTACCAGCAGGAAGGCGACTACGACTTCCATTACCGGATTATTCAGGGTAGCGGCAACCAGACCCTGGTCAAGATGCTCTGCGGCGAGCTGTACCAGCTGGTGCGCATGTACCGTATCCAGTTctct                                                                                              <  2:83716/149‑1 (MQ=255)
                                                 gTTCCAGGCCGGCCTGGGCTACTACCAGCAGGAAGGCGACTACGACTTCCATTACCGGATTATTCAGGGTAGCGGCAACCAGACCCTGGTCAAGATGCTCTGCGGCGAGCTGTACCAGCTGGTGCGCATGTACCGTATCCAGTTCTCTg                                                                                             <  8:72907/149‑1 (MQ=255)
                                                  ttCCAGGCCGGCCTGGGCTACTACCAGCAGGAAGGCGACTACGACTTCCATTACCGGATTATTCAGGGTAGCGGCAACCAGACCCTGGTCAAGATGCTCTGCGGCGAGCTGTACCAGCTGGTGCGCATGTACCGTATCCAGTTCTCTGc                                                                                            >  4:86334/1‑149 (MQ=255)
                                                                    tactacCAGCAGGAAGGCGACTACGACTTCCATTACCGGATTATTCAGGGTAGCGGCAACCAGACCCTGGTCAAGATGCTCTGCGGCGAGCTGTACCAGCTGGTGCGCATGTACCGTATCCAGTTCTCTGCCACGCCTAAc                                                                                  >  6:255624/1‑141 (MQ=255)
                                                                    tactacCAGCAGGAAGGCGACTACGACTTCCATTACCGGATTATTCAGGGTAGCGGCAACCAGACCCTGGTCAAGATGCTCTGCGGCGAGCTGTACCAGCTGGTGCGCATGTACCGTATCCAGTTCTCTGCCACGCCTAAc                                                                                  <  5:255624/141‑1 (MQ=255)
                                                                    tactacCAGCAGGAAGGCGACTACGACTTCCATTACCGGATTATTCAGGGTAGCGGCAACCAGACCCTGGTCAAGATGCTCTGCGGCGAGCTGTACCAGCTGGTGCGCATGTACCGTATCCAGTTCTCTGCCACGCCTAACCGGCCACg                                                                          >  8:179927/1‑149 (MQ=255)
                                                                         ccAGCAGGAAGGCGACTACGACTTCCATTACCGGATTATTCAGGGTAGCGGCAACCAGACCCTGGTCGAGATGCTCTGCGGCGAGCTGTACCAGCTGGTGCGCATGTACCGTATCCAGTTCTCTGCCACGCCTAACCGGCCACGCCAgg                                                                     >  6:32042/1‑149 (MQ=255)
                                                                         ccAGCAGGAAGGCGACTACGACTTCCATTACCGGATTATTCAGGGTAGCGGCAACCAGACCCTGGTCGAGATGCTCTGCGGCGAGCTGTACCAGCTGGTGCGCATGTACCGTATCCAGTTCTCTGCCACGCCTAACCGGCCACGCCAgg                                                                     >  6:27189/1‑149 (MQ=255)
                                                                          cagcagGAAGGCGACTACGACTTCCATTACCGGATTATTCAGGGTAGCGGCAACCAGACCCTGGTCGAGATGCTCTGCGGCGAGCTGTACCAGCTGGTGCGCATGTACCGTATCCAGTTCTCTGCCACGCCTAACCGGCCACGCCAGGc                                                                    <  5:27189/149‑1 (MQ=255)
                                                                          cagcagGAAGGCGACTACGACTTCCATTACCGGATTATTCAGGGTAGCGGCAACCAGACCCTGGTCGAGATGCTCTGCGGCGAGCTGTACCAGCTGGTGCGCATGTACCGTATCCAGTTCTCTGCCACGCCTAACCGGCCACGCCAGGc                                                                    <  5:32042/149‑1 (MQ=255)
                                                                                   ggCGACTACGACTTCCATTACCGGATTATTCAGGGTAGCGGCAACCAGACCCTGGTCAAGATGCTCTGCGGCGAGCTGTACCAGCTGGTGCGCATGTACCGTATCCAGTTCTCTGCCACGCCTAACCGGCCACGCCAGGCCTTTGCCg                                                            >  2:23885/1‑148 (MQ=255)
                                                                                     cGACTACGACTTCCATTACCGGATTATTCAGGGTAGCGGCAACCAGACCCTGGTCAAGATGCTCTGCGGCGAGCTGTACCAGCTGGTGCGCATGTACCGTATCCAGTTCTCTGCCACGCCTAACCGGCCACGCCAGGCCTTTGCCGAac                                                         <  2:137944/149‑1 (MQ=255)
                                                                                         tACGACTTCCATTACCGGATTATTCAGGGTAGCGGAAACCAGACCCTGGTCAAGATGCTCTGCGGCGAGCTGTACCAGCTGGTGCCCATGTACCGTATCCAGTTCTATGCCACGACTAACAGGCCACGCCAGGCCTTTGCCGAACACCa                                                     >  5:221366/1‑149 (MQ=255)
                                                                                           cGACTTCCATTACCGGATTATTCAGGGTAGCGGCAACCAGACCCTGGTCAAGATGCTCTGCGGCGAGCTGTACCAGCTGGTGCGCATGTACCGTATCCAGTTCTCTGCCACGCCTAACCGGCCACGCCAGGCCTTTGCCGAACACCa                                                     <  2:94297/147‑1 (MQ=255)
                                                                                            gACTTCCATTACCGGATTATTCAGGGTAGCGGCAACCAGACCCTGGTCAAGATGCTCTGCGGCGAGCTGTACCAGCTGGTGCGCATGTACCGTATCCAGTTCTCTGCCACGCCTAACCGGCCACGCCAGGCCTTTGCCGAACACCATCg                                                  <  7:259168/149‑1 (MQ=255)
                                                                                            gACTTCCATTACCGGATTATTCAGGGTAGCGGCAACCAGACCCTGGTCAAGATGCTCTGCGGCGAGCTGTACCAGCTGGTGCGCATGTACCGTATCCAGTTCTCTGCCACGCCTAACCGGCCACGCCAGGCCTTTGCCGAACACCATCg                                                  <  7:179927/149‑1 (MQ=255)
                                                                                               ttCCATTACCGGATTATTCAGGGTAGCGGCAACCAGACCCTGGTCAAGATGCTCTGCGGCGAGCTGTACCAGCTGGTGCGCATGTACCGTATCCAGTTCTCTGCCACGCCTAACCGGCCACGCCAGGCCTTTGCCGAACACCATCGCAt                                               >  6:94661/1‑149 (MQ=255)
                                                                                                   aTTACCGGATTATTCAGGGTAGCGGCAACCAGACCCTGGTCAAGATGCTCTGCGGCGAGCTGTACCAGCTGGTGCGCATGTACCGTATCCAGTTCTCTGCCACGCCTAACCGGCCACGCCAGGCCTTTGCCGAACACCATCGCattatt                                           >  7:316869/1‑149 (MQ=255)
                                                                                                   aTTACCGGATTATTCAGGGTAGCGGCAACCAGACCCTGGTCAAGATGCTCTGCGGCGAGCTGTACCAGCTGGTGCGCATGTACCGTATCCAGTTCTCTGCCACGCCTAACCGGCCACGCCAGGCCTTTGCCGAACACCATCGCATTCtt                                           <  3:86334/149‑1 (MQ=255)
                                                                                                       ccGGATTATTCAGGGTAGCGGCAACCAGACCCTGGTCAAGATGCTCTGCGGCGAGCTGTACCAGCTGTTGCGCATGTACCGTATCCAGTTCTCTGCCACGCCTAACCGGCCACGCCAGGCCTTTGCCGAACACCATCGCATTCTTGAt                                        <  7:218234/148‑1 (MQ=255)
                                                                                                       ccGGATTATTCAGGGTAGCGGCAACCAGACCCTGGTCAAGATGCTCTGCGGCGAGCTGTACCAGCTGGTGCGCATGTACCGTATCCAGTTCTCTGCCACGCCTAACCGGCCACGCCAGGCCTTTGCCGAACACCATCGCATTCTTGAt                                        <  4:190059/148‑1 (MQ=255)
                                                                                                       ccGGATTATTCAGGGTAGCGGCAACCAGACCCTGGTCAAGATGCTCTGCGGCGAGCTGTACCAGCTGGTGCGCATGTACCGTATCCAGTTCTCTGCCACGCCTAACCGGCCACGCCAGGCCTTTGCCGAACACCATCGCATTCTTGAt                                        <  4:97215/148‑1 (MQ=255)
                                                                                                          gATTATTCAGGGTAGCGGCAACCAGACCCTGGTCAAGATGCTCTGCGGCGAGCTGTACCAGCTGGTgcgc                                                                                                                   <  3:262501/70‑1 (MQ=255)
                                                                                                          gATTATTCAGGGTAGCGGCAACCAGACCCTGGTCAAGATGCTCTGCGGCGAGCTGTACCAGCTGGTgcgc                                                                                                                   >  4:262501/1‑70 (MQ=255)
                                                                                                            ttattCAGGGTAGCGGCAACCAGACCCTGGTCAAGATGCTCTGCGGCGAGCTGTACCAGCTGGTGCGCATGTACCGTATCCAGTTCTCTGCCACGCCTAACCGGCCACGCCAGGCCTTTGCCGAACACCATCGCATTCTTGATGCCATc                                  >  5:318015/1‑149 (MQ=255)
                                                                                                             tattCAGGGTCGCGGCAACCAGACCCTGGTCAAGATGCTCTGCGGCGAGCTGTACCAGCTGGTGCGCATGTACCGTATCCAGTTCTCTGCCACGCCTAACCg                                                                                <  5:180044/102‑1 (MQ=255)
                                                                                                             tattCAGGGTAGCGGCAACCAGACCCTGGTCAAGATGCTCTGCGGCGAGCTGTACCAGCTGGTGCGCATGTACCGTATCCAGTTCTCTGCCACGCCTAACCg                                                                                >  6:180044/1‑102 (MQ=255)
                                                                                                                         cGGCAACCAGACCCTGGTCAAGATGCTCTGCGGCGAGCTGTACCAGCTGGTGCGCATGTACCGTATCCAGTTCTCTGCCACGCCTAACCGGCCACGCCAGGCCTTTGCCGAACACCATCGCATTCTTGATGCCATCGCCGACCGTGACg                     >  1:306047/1‑149 (MQ=255)
                                                                                                                          ggCAACCAGACCCTGGTCAAGATGCTCTGCGGCGAGCTGTACCAGCTGGTGCGCATGTACCGTATCCAGTTCTCTGCCACGCCTAACCGGCCACGCCAGGCCTTTGCCGAACACCATCGCATTCTTGATGc                                      >  6:284401/1‑131 (MQ=255)
                                                                                                                          ggCAACCAGACCCTGGTCAAGATGCTCTGCGGCGAGCTGTACCAGCTGGTGCGCATGTACCGTATCCAGTTCTCTGCCACGCCTAACCGGCCACGCCAGGCCTTTGCCGAACACCATCGCATTCTTGATGc                                      <  5:284401/131‑1 (MQ=255)
                                                                                                                          ggCAACCAGACCCTGGTCAAGATGCTCTGCGGCGAGCTGTACCAGCTGGTGCGCATGTACCGTATCCAGTTCTCTGCCACGCCTAACCGGCCACGCCAGGCCTTTGCCGAACACCATCGCATTCTTGATGCCATCGCCGACCGTGACgg                    >  4:161476/1‑149 (MQ=255)
                                                                                                                           gCAACCAGACCCTGGTCAAGATGCTCTGCGGCGAGCTGTACCAGCTGGTGCGCATGTACCGTATCCAGTTCTCTGCCACGCCTAACCGGCCACGCCAGGCCTTTGCCGAACACCATCGCATTCTTGATGCCATCGCCGACCGTGACGGc                   >  8:174342/1‑149 (MQ=255)
                                                                                                                              aCCAGACCCTGGTCAAGATGCTCTGCGGCGAGCTGTACCAGCTGGTGCGCATGTACCGTATCCAGTTCTCTGCCACGCCTAACCGGCCACGCCAGGCCTTTGCCGAACACCATCGCATTCTTGATGCCATCGCCGACCGTGACGGCGaa                >  1:332685/1‑149 (MQ=255)
                                                                                                                                         gTCAAGATGCTCTGCGGCGAGCTGTACCAGCTGGTGCGCATGTACCGTATCCAGTTCTCTGCCACGCCTAACCGGCCACGCCAGGCCTTTGCCGAACACCATCGCATTCTTGATGCCATCGCCGACCGTGACGGCGAACTGGCCGAACt     >  6:271139/1‑149 (MQ=255)
                                                                                                                                          tCAAGATGCTCTGCGGCGAGCTGTACCAGCTGGTGCGCATGTACCGTATCCAGTTCTCTGCCACGCCTAACCGGCCACGCCAGGCCTTTGCCGAACACCATCGCATTCTTGATGCCATCGCCGACCGTGACGGCGAACTGGCCGAACTc    >  1:188583/1‑149 (MQ=255)
                                                                                                                                            aaGATGCTCTGCGGCGAGCTGTACCAGCTGGTGCGCATGTACCGTATCCAGTTCTCTGCCACGCCTAACCGGCCACGCCAGGCCTTTGCCGAACACCATCGCATTCTTGATGCCATCGCCGACCGTGACGGCGAACTGGCCGAACTCCt  >  3:267308/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                    
ACATCGACGAACTGCGCCGGGTGCTCGACACCCATGAGCGCGATGCCGCGTTCCAGGCCGGCCTGGGCTACTACCAGCAGGAAGGCGACTACGACTTCCATTACCGGATTATTCAGGGTAGCGGCAACCAGACCCTGGTCAAGATGCTCTGCGGCGAGCTGTACCAGCTGGTGCGCATGTACCGTATCCAGTTCTCTGCCACGCCTAACCGGCCACGCCAGGCCTTTGCCGAACACCATCGCATTCTTGATGCCATCGCCGACCGTGACGGCGAACTGGCCGAACTCCT  >  NC_002947/2661524‑2661812

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: ‑