Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 635,912 G→A G235D (GGC→GAC)  ahpF → alkyl hydroperoxide reductase subunit F

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273635,9120GA91.1% 167.5 / 6.7 56G235D (GGC→GAC) ahpFalkyl hydroperoxide reductase subunit F
Reads supporting (aligned to +/- strand):  ref base G (3/2);  new base A (23/28);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.53e-01

AAATCGTTGCCAAAATTGATACTGGCGCGGAAAAACGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATA  >  NZ_CP009273/635777‑636047
                                                                                                                                       |                                                                                                                                       
aaaTCGTTGCCAAAATTGATACTGGCGCGGAAAAACGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACAt                                                                                                                                      <  1:888723/139‑1 (MQ=255)
aaaTCGTTGCCAAAATTGATACTGGCGCGGAAAAACGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACAt                                                                                                                                      <  2:790792/139‑1 (MQ=255)
aaaTCGTTGCCAAAATTGATACTGGCGCGGAAAAACGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACAt                                                                                                                                      <  2:748659/139‑1 (MQ=255)
     gTTGCCAAAATTGATACTGGCGCGGAAAAACGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTa                                                                                                                                 <  2:238971/139‑1 (MQ=255)
      ttGCCAAAATTGATACTGGCGCGGAAAAACGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTAc                                                                                                                                >  2:456384/1‑139 (MQ=255)
        gCCAAAATTGATACTGGCGCGGAAAAACGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTAcc                                                                                                                               <  2:73677/138‑1 (MQ=255)
                    aCTGGCGCGGAAAAACGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTAATCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCg                                                                                                                  >  2:342803/1‑139 (MQ=255)
                        gcgcGGAAAAACGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACg                                                                                                              >  1:299632/1‑139 (MQ=255)
                               aaaaCGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTtggt                                                                                                       <  2:299632/139‑1 (MQ=255)
                                   cGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTc                                                                                                   <  1:879022/139‑1 (MQ=255)
                                     tGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAg                                                                                                 <  2:443669/139‑1 (MQ=255)
                                     tGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAg                                                                                                 <  1:986221/139‑1 (MQ=255)
                                        ggCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGGGCAGCGGCAGCCATTTAATCCGCACGTAAAGACATCCGTACCGGTGTGATGGGCGAAAGTTTTGGGGGGCAGATc                                                                                              >  2:915895/1‑139 (MQ=255)
                                        ggCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGGGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGGGGTCAGATc                                                                                              >  2:915840/1‑139 (MQ=255)
                                          cAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCt                                                                                            <  2:1093621/139‑1 (MQ=255)
                                            gaagaGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCg                                                                                          <  2:155704/139‑1 (MQ=255)
                                                   tGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGt                                                                                   >  1:986970/1‑139 (MQ=255)
                                                            gTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCg                                                                           >  2:1135331/1‑138 (MQ=255)
                                                             tGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGaa                                                                         <  2:986970/139‑1 (MQ=255)
                                                             tGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGaa                                                                         <  1:700328/139‑1 (MQ=255)
                                                               aTGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCCGCAATTTACTCCGCACGTAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGGGGTCAGATCCTCGATACCGTTGATATCGaaaa                                                                       >  1:211439/1‑139 (MQ=255)
                                                               aTGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGaaaa                                                                       >  2:294054/1‑139 (MQ=255)
                                                               aTGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGaaaa                                                                       >  1:297765/1‑139 (MQ=255)
                                                               aTGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGGGCAGCAGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGaaa                                                                        >  2:1198529/1‑138 (MQ=255)
                                                                  cTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTa                                                                    <  1:342803/139‑1 (MQ=255)
                                                                     aTGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACAt                                                                 >  1:568506/1‑139 (MQ=255)
                                                                      tGACGTATTAATCGTCGGTTCCGGCCCGGCGGGGGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACAtt                                                                >  2:1062409/1‑139 (MQ=255)
                                                                         cGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCt                                                             <  2:568506/139‑1 (MQ=255)
                                                                            aTTAATCGTCGGTTCCGGCCCGGCGGGTGCCGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATTGAAAACTACATTTCTGTa                                                          >  1:484079/1‑139 (MQ=255)
                                                                            aTTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCg                                                                                          >  1:689465/1‑107 (MQ=255)
                                                                            aTTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCg                                                                                          <  2:689465/107‑1 (MQ=255)
                                                                            aTTAATCGTCGGTTCCGGCCCGGCGGGGGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTa                                                          >  1:1230871/1‑139 (MQ=255)
                                                                            aTTAATCGTCGGTTCCGGCCCGGCGGGGGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGGTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTa                                                          >  2:918448/1‑139 (MQ=255)
                                                                             ttAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTAc                                                         >  2:342239/1‑139 (MQ=255)
                                                                               aaTCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTAcc                                                        <  1:862130/138‑1 (MQ=255)
                                                                               aaTCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTAcc                                                        <  1:456384/138‑1 (MQ=255)
                                                                                 tcgtcgGTTGCGGCACGGCGGGTGCAGCCGCAGCAAGTGACTCTGCACGTAAAGACATCCGCACCGGTCTGATGGACGAACGTTTTGGTGGTCAGATCCTCGAGACCGTTGATATCGAAAACTTCATTTCTGTACCGaa                                                     <  2:1219159/139‑1 (MQ=255)
                                                                                 tcgtcgGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGaa                                                     <  2:297765/139‑1 (MQ=255)
                                                                                     cgGTTCCGGTCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACt                                                 >  1:882448/1‑139 (MQ=255)
                                                                                     cgGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACt                                                 >  1:1070960/1‑139 (MQ=255)
                                                                                       gTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGa                                               <  1:342239/139‑1 (MQ=255)
                                                                                       gTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGa                                               <  1:1135331/139‑1 (MQ=255)
                                                                                               ccggcGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGa                                       <  1:116698/139‑1 (MQ=255)
                                                                                                  gcggGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGc                                    >  1:1290965/1‑139 (MQ=255)
                                                                                                         cAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAgg                             <  1:294054/139‑1 (MQ=255)
                                                                                                           gCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGgcg                           <  2:665518/139‑1 (MQ=255)
                                                                                                            cGGCAGCAATTTACTCCGCACGTAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGgcg                           >  2:1097537/1‑138 (MQ=255)
                                                                                                                 gcaATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGaaaa                                                                       >  1:1135639/1‑89 (MQ=255)
                                                                                                                 gcaATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGaaaa                                                                       <  2:1135639/89‑1 (MQ=255)
                                                                                                                 gcaATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTg                      >  2:31990/1‑138 (MQ=255)
                                                                                                                         cTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCa              <  1:1062409/138‑1 (MQ=255)
                                                                                                                            cGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACg            <  2:484079/137‑1 (MQ=255)
                                                                                                                            cGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGtt          >  2:205593/1‑139 (MQ=255)
                                                                                                                             gCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGtt          >  2:14093/1‑138 (MQ=255)
                                                                                                                             gCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTtg         <  2:291950/139‑1 (MQ=255)
                                                                                                                                    aaGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATa  <  2:105295/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                       
AAATCGTTGCCAAAATTGATACTGGCGCGGAAAAACGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATA  >  NZ_CP009273/635777‑636047

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