Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,788,121 T→A 8.0% V134D (GTC→GAC)  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,1210TA 8.0% 142.3 / 5.5 50V134D (GTC→GAC) ahpFalkyl hydroperoxide reductase subunit F
Reads supporting (aligned to +/- strand):  ref base T (25/21);  new base A (2/2);  total (27/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 = 3.24e-01

GAGTTCACATCGTTGGTGCTGGCACTGCTGCAAGTGGGTGGCCACCCGTCAAAGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTACCTCAACGGTGAGGTGTTTGGTCAGGGC  >  NC_002947/2787973‑2788269
                                                                                                                                                    |                                                                                                                                                    
gAGTTCACATCGTTGGTGCTGGCACTGCTGCAAGTGGGTGGCCACCCGTCAAAGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGAcgtc                                                                                                                                                     <  1:525288/150‑1 (MQ=255)
          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                                                                                                                                         >  2:421649/1‑97 (MQ=255)
                                                              cGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCTGGACGACGTCCGGGCagtt                                                                                                                                         <  1:421649/100‑4 (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:1363882/150‑1 (MQ=255)
                                                                                       cTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGaa                                                              >  2:606715/1‑150 (MQ=255)
                                                                                             ggTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTAGTACAGGCACTGAACCTGGTGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCgcgc                                                        >  1:725490/1‑150 (MQ=255)
                                                                                                 aaTTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGa                                                    >  2:1292222/1‑150 (MQ=255)
                                                                                                    tCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGa                                                                                                                                  <  2:531111/69‑1 (MQ=255)
                                                                                                    tCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGa                                                                                                                                  >  1:531111/1‑69 (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)
                                                                                                                                                    |                                                                                                                                                    
GAGTTCACATCGTTGGTGCTGGCACTGCTGCAAGTGGGTGGCCACCCGTCAAAGGCCAGTGCCGAAGTGATCGAGCAGATCCAGGCCCTGGAAGGTGAATTCAACTTCGAAACCTACTTCTCGCTGTCGTGCCAGAACTGCCCGGACGTCGTCCAGGCACTGAACCTGATGGCGGTACTCAACCCCAACGTGCGCCATGTTGCCATCGACGGTGCGCTGTTCCAGGACGAAGTGGAATCGCGCAAGATCATGGCGGTACCGAGTATCTACCTCAACGGTGAGGTGTTTGGTCAGGGC  >  NC_002947/2787973‑2788269

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: