Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,788,127 A→G 7.1% Q136R (CAG→CGG)  ahpF → alkyl hydroperoxide reductase subunit F

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,788,1270AG 7.1% 160.7 / 2.8 56Q136R (CAG→CGG) ahpFalkyl hydroperoxide reductase subunit F
Reads supporting (aligned to +/- strand):  ref base A (30/22);  new base G (2/2);  total (32/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 = 4.95e-01

CGTTGGTGCTGGCACTGCTGCAAGTGGGTGGCCACCCGTCAAAGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTACCTCAACGGTGAGGTGTTTGGTCAGGGCCGTATGG  >  NC_002947/2787983‑2788276
                                                                                                                                                |                                                                                                                                                     
cGTTGGTGCTGGCACTGCTGCAAGTGGGTGGCCACCCGTCAAAGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCAc                                                                                                                                                  <  1:1200815/150‑1 (MQ=255)
cGTTGGTGCTGGCACTGCTGCAAGTGGGTGGCCACCCGTCAAAGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCAc                                                                                                                                                  <  1:30404/150‑1 (MQ=255)
                          ggtggCCACCCGTCAAAGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAAcccc                                                                                                                        >  1:999721/1‑150 (MQ=255)
                               ccACCCGTCAAAGGCCAGTGCCGAAGTGATCGAGCAGAGCCAGGCCCTGGAAGGTGAATTCAACTTCTGATGCTCCTTCACCCGTACGGGCCAGGACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTGCGGAACCCCAACGt                                                                                                                   >  1:629721/1‑150 (MQ=255)
                                 aCCCGTCAAAGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTgc                                                                                                                 <  2:1525963/150‑1 (MQ=255)
                                   ccGTCAAAGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTgcgc                                                                                                               >  1:1649067/1‑150 (MQ=255)
                                    cGTCAAAGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACgaaaggc                                                                                                              <  1:906234/150‑7 (MQ=255)
                                         aaGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCAAGACACTGATCCGGATGGCGGGAGGCAACCCCAACGTGCGCCATGtt                                                                                                         >  2:1287499/1‑150 (MQ=255)
                                            gCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGcc                                                                                                      >  1:1363882/1‑150 (MQ=255)
                                                    cGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCTGGACGACGTCCGGGCagtt                                                                                                                                                <  1:421649/100‑4 (MQ=255)
                                                    cGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCTGGACGACGTCCGGGCagtt                                                                                                                                                >  2:421649/1‑97 (MQ=255)
                                                     gAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGt                                                                                             <  2:1146084/150‑1 (MQ=255)
                                                           aTCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTg                                                                                       >  2:135282/1‑150 (MQ=255)
                                                           aTCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTg                                                                                       >  1:1500798/1‑150 (MQ=255)
                                                             cGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGtt                                                                                     <  2:1500798/150‑1 (MQ=255)
                                                              gAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAAc                                                                                                                           >  1:468485/1‑111 (MQ=255)
                                                              gAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAAc                                                                                                                           <  2:468485/111‑1 (MQ=255)
                                                                       cAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCAAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGaa                                                                           >  2:276263/1‑150 (MQ=255)
                                                                           cccTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTgg                                                                       >  1:1929041/1‑150 (MQ=255)
                                                                             cTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGaa                                                                     >  2:606715/1‑150 (MQ=255)
                                                                             cTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGaa                                                                     <  2:1363882/150‑1 (MQ=255)
                                                                                   ggTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTAGTACAGGCACTGAACCTGGTGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCgcgc                                                               >  1:725490/1‑150 (MQ=255)
                                                                                       aaTTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGa                                                           >  2:1292222/1‑150 (MQ=255)
                                                                                          tCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGa                                                                                                                                         >  1:531111/1‑69 (MQ=255)
                                                                                          tCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGa                                                                                                                                         <  2:531111/69‑1 (MQ=255)
                                                                                          tCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCa                                                        <  1:606715/150‑1 (MQ=255)
                                                                                            aaCTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATg                                                      >  2:276104/1‑150 (MQ=255)
                                                                                             aCTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATgg                                                     <  2:1742466/150‑1 (MQ=255)
                                                                                                  gAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGCCGAAGGATATTCGCGCAAGATCATGGCGGTa                                                >  2:933465/1‑150 (MQ=255)
                                                                                                          cTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTAt                                        <  1:759834/150‑1 (MQ=255)
                                                                                                           ttCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATc                                       >  2:1758906/1‑150 (MQ=255)
                                                                                                             ctcGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTa                                     >  2:304262/1‑150 (MQ=255)
                                                                                                             ctcGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTa                                     <  2:629721/150‑1 (MQ=255)
                                                                                                              tcGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTAc                                    <  1:276104/150‑1 (MQ=255)
                                                                                                              tcGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTAc                                    <  1:1292222/150‑1 (MQ=255)
                                                                                                              tcGCTGGCGGGCCAGAACTTCCCGGACGTCCTCCAGGCACTGAACCTGATGGCGGTACTCAACACCAACGTGCGCCATGTTGCCACAGACGGTCCGCTGTTCAACGATGATTTGGAATCGCGCAAGACCATGGCGGTACCGCGTATCTAc                                    <  2:559386/150‑1 (MQ=255)
                                                                                                                ctgatCGTGCCAGAACTGCCCCGACGTCGTCCAGGCACTGACCCTGATGGCGGCACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTGCAGGACGAAGTGGACTCGCGCAAGATCGTGGAGGTACCGAGTATCTAcca                                  >  1:559386/5‑149 (MQ=255)
                                                                                                                   gTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTACCTCaa                               <  1:1287499/150‑1 (MQ=255)
                                                                                                                    tCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCGGGACGACGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTACCTCAAc                              >  1:1564947/1‑150 (MQ=255)
                                                                                                                      gTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCt                                                                                        >  2:1429948/1‑90 (MQ=255)
                                                                                                                      gTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCt                                                                                        <  1:1429948/90‑1 (MQ=255)
                                                                                                                      gTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTACCTCAACgg                            >  2:1002634/1‑150 (MQ=255)
                                                                                                                        gCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGc                                                    >  2:499332/1‑124 (MQ=255)
                                                                                                                        gCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGc                                                    <  1:499332/124‑1 (MQ=255)
                                                                                                                           aGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTACCTCAACGGTGAgg                       >  2:1932604/1‑150 (MQ=255)
                                                                                                                             aaCTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTACCTCAACGGTGAGgtg                     >  2:1454131/1‑150 (MQ=255)
                                                                                                                                tGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTACCTCAACGGTGAGGTGttt                  <  1:500763/150‑1 (MQ=255)
                                                                                                                                 gcctggACGACGTCCGGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGa                                            <  1:1764203/119‑1 (MQ=255)
                                                                                                                                 gCCTGGACGACGTCCGGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGa                                            >  2:1764203/1‑123 (MQ=255)
                                                                                                                                      gACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTACCTCAACGGTGAGGTGTTTGGTCAg            <  2:1649067/150‑1 (MQ=255)
                                                                                                                                         gtcgtcCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTACCTCAACGGTGAGGTGTTTGGTCAGGGc         >  1:986346/1‑150 (MQ=255)
                                                                                                                                             tcCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTACCTCAACGGTGAGGTGTTTGGTCAGGGCCGTa     >  1:1248052/1‑150 (MQ=255)
                                                                                                                                               cAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTAc                                               >  1:328674/1‑106 (MQ=255)
                                                                                                                                               cAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTAc                                               <  2:328674/106‑1 (MQ=255)
                                                                                                                                               cAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTACCTCAACGGTGAGGTGTTTGGTCAGGGCCGTATg   >  2:429979/1‑150 (MQ=255)
                                                                                                                                               cAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCACAAGATCATGGCGGTACCGAGTATCTACCTCAACGGTGAGGTGTTTGGTCAGGGCCGTATg   >  1:1577102/1‑150 (MQ=255)
                                                                                                                                                aGGCACTGAACCTGATGGCGGTACTCACCCCCGACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTACCTCAACGGTGAGGTGTTTGGTCAGGGCCGTATgg  >  1:1808390/1‑150 (MQ=255)
                                                                                                                                                aGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTACCTCAACGGTGAGGTGTTTGGTCAGGGCCGTATgg  >  2:487814/1‑150 (MQ=255)
                                                                                                                                                aGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTACCTCAACGGTGAGGTGTTTGGTCAGGGCCGTATgg  >  2:1550231/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
CGTTGGTGCTGGCACTGCTGCAAGTGGGTGGCCACCCGTCAAAGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTACCTCAACGGTGAGGTGTTTGGTCAGGGCCGTATGG  >  NC_002947/2787983‑2788276

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