Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,817,653 A→G 7.2% P258P (CCA→CCG nrfA → nitrite reductase, formate‑dependent, cytochrome

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,817,6530AG 7.2% 233.2 / 8.8 83P258P (CCA→CCGnrfAnitrite reductase, formate‑dependent, cytochrome
Reads supporting (aligned to +/- strand):  ref base A (39/38);  new base G (4/2);  total (43/40)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGACCT  >  minE/2817585‑2817723
                                                                    |                                                                      
gAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCaa                                                                       <  1:3209401/70‑1 (MQ=255)
gAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAAt                                                                      >  1:3286527/1‑71 (MQ=255)
gAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAAt                                                                      >  1:2905041/1‑71 (MQ=255)
  aaaTATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATg                                                                     <  1:4047776/70‑1 (MQ=255)
   aaTATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCt                                                                   <  1:1346247/71‑1 (MQ=255)
   aaTATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCt                                                                   <  1:2375288/71‑1 (MQ=255)
     tatGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGa                                                                 <  1:3953330/71‑1 (MQ=255)
     tatGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGa                                                                 <  1:3600725/71‑1 (MQ=255)
      atGGAGCAGTATTAAGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCt                                                                   <  1:1974422/68‑1 (MQ=255)
       tGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGaaa                                                               <  1:73764/71‑1 (MQ=255)
          aGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAgcg                                                            <  1:4042071/71‑1 (MQ=255)
          aGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAgcg                                                            <  1:2239684/71‑1 (MQ=255)
          aGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAgcg                                                            <  1:2975046/71‑1 (MQ=255)
           gCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAgcg                                                            <  1:3702113/70‑1 (MQ=255)
               tATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGa                                                                 <  1:3663386/61‑1 (MQ=255)
               tATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGa                                                                 <  1:218424/61‑1 (MQ=255)
               tATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGa                                                                 <  1:2290929/61‑1 (MQ=255)
               tATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGa                                                                 <  1:114285/61‑1 (MQ=255)
                 ttACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCa                                                                        <  1:3527908/52‑1 (MQ=255)
                 ttACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCGATGCTGAAAGCgca                                                          >  1:3793347/1‑66 (MQ=255)
                 ttACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCGATGCTGAAAGCgca                                                          >  1:858926/1‑66 (MQ=255)
                 ttACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGc                                                                    >  1:1243392/1‑56 (MQ=255)
                   aCGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCGATGCTGAAAGCGCAGCACCCg                                                   >  1:2814798/1‑71 (MQ=255)
                   aCGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCg                                                   >  1:1406366/1‑71 (MQ=255)
                   aCGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCg                                                   >  1:1005533/1‑71 (MQ=255)
                   aCGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCg                                                   >  1:936040/1‑71 (MQ=255)
                   aCGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCg                                                   >  1:1163269/1‑71 (MQ=255)
                    cGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCg                                                   <  1:1330173/70‑1 (MQ=255)
                     gACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCgg                                                  >  1:1372048/1‑70 (MQ=255)
                     gACAAAATTGCCT‑‑TCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAat                                               <  1:1652844/71‑1 (MQ=255)
                      aCAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGaa                                                <  1:1280558/71‑1 (MQ=255)
                      aCAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGaa                                                <  1:2544144/71‑1 (MQ=255)
                       cAAAATTGCCTTCTCTTACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGaa                                                >  1:3289683/1‑70 (MQ=255)
                       cAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCg                                                   <  1:3338866/67‑1 (MQ=255)
                       cAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGaa                                                <  1:1032216/70‑1 (MQ=255)
                       cAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGaa                                                >  1:80981/1‑70 (MQ=255)
                       cAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGaa                                                >  1:1294872/1‑70 (MQ=255)
                        aaaaTTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAata                                              >  1:1899087/1‑71 (MQ=255)
                        aaaaTTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAata                                              >  1:2498032/1‑71 (MQ=255)
                         aaaTTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAatat                                             <  1:3661853/71‑1 (MQ=255)
                          aaTTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATg                                            <  1:3169788/71‑1 (MQ=255)
                          aaTTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATg                                            <  1:2671822/71‑1 (MQ=255)
                            ttGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGaa                                          <  1:943314/71‑1 (MQ=255)
                                cTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCGATGCTGAAAgcg                                                            <  1:388214/49‑1 (MQ=255)
                                 ttCTCTGACTGGACTAACTCCCTGTCGAAAACGCCGATGCTGAAAGCGCAGCACCCGGAATATGAAACCTg                                     >  1:1491079/1‑71 (MQ=255)
                                  tctcTGACTGGACTAACTCCCTGTCGAAAACGCCAATg                                                                     >  1:3379992/1‑38 (MQ=255)
                                   ctctGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTg                                     <  1:396463/69‑1 (MQ=255)
                                   ctctGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTg                                     <  1:2945246/69‑1 (MQ=255)
                                   ctctGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGa                                   >  1:3155985/1‑71 (MQ=255)
                                   ctctGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGa                                   >  1:2106484/1‑71 (MQ=255)
                                    tctGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCgcag                                                         >  1:1338583/1‑48 (MQ=255)
                                    tctGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCgcag                                                         >  1:1860034/1‑48 (MQ=255)
                                    tctGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCgcag                                                         >  1:800823/1‑48 (MQ=255)
                                    tctGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCgcag                                                         >  1:255607/1‑48 (MQ=255)
                                    tctGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGa                                   >  1:1591137/1‑70 (MQ=255)
                                      tgacttgactAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAg                                >  1:2101561/1‑71 (MQ=255)
                                       gACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCgca                                                          >  1:1984732/1‑44 (MQ=255)
                                        aCTGGACTAACTCCCTGTCGAAAACGCCGATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAg                                <  1:743604/69‑1 (MQ=255)
                                            gACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCa                          <  1:2740383/71‑1 (MQ=255)
                                              cTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGc                           >  1:1427986/1‑68 (MQ=255)
                                              cTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCa                          <  1:3671617/69‑1 (MQ=255)
                                              cTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCa                          <  1:1681001/69‑1 (MQ=255)
                                              cTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCAtt                        <  1:775087/71‑1 (MQ=255)
                                              cTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCAtt                        >  1:2269735/1‑71 (MQ=255)
                                              cTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCAtt                        >  1:2624053/1‑71 (MQ=255)
                                              cTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCAtt                        >  1:230485/1‑71 (MQ=255)
                                                 aCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTc                       <  1:3092498/69‑1 (MQ=255)
                                                 aCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCa                      >  1:3910039/1‑70 (MQ=255)
                                                 aCTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCAc                     >  1:3070418/1‑71 (MQ=255)
                                                   tCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGa                                   >  1:1102035/1‑55 (MQ=255)
                                                      cTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGaa                                          <  1:2232144/45‑1 (MQ=255)
                                                      cTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTa                 >  1:2311299/1‑70 (MQ=255)
                                                        gTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAAc                                        >  1:971150/1‑45 (MQ=255)
                                                         tCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTaaaaa             >  1:96925/1‑71 (MQ=255)
                                                         tCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTaaaaa             <  1:3295755/71‑1 (MQ=255)
                                                           gAAAACGCCAATGCTGAAAGCGCAGCACCCGGAata                                              >  1:3158234/1‑36 (MQ=255)
                                                           gAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCa                      >  1:412672/1‑60 (MQ=255)
                                                           gAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAaaca           <  1:1925816/71‑1 (MQ=255)
                                                           gAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAaaca           <  1:1938158/71‑1 (MQ=255)
                                                           gAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAaac            >  1:1606798/1‑70 (MQ=255)
                                                            aaaaCGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTc                       <  1:1306865/58‑1 (MQ=255)
                                                              aaCGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCAc                     >  1:3426427/1‑58 (MQ=255)
                                                               aCGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGt       <  1:2053077/71‑1 (MQ=255)
                                                                  ccAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGAc    >  1:1905037/1‑71 (MQ=255)
                                                                   cAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGAcc   >  1:768898/1‑71 (MQ=255)
                                                                    aaTGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGACCt  <  1:2607568/71‑1 (MQ=255)
                                                                    |                                                                      
GAAAATATGGAGCAGTATTACGACAAAATTGCCTTCTCTGACTGGACTAACTCCCTGTCGAAAACGCCAATGCTGAAAGCGCAGCACCCGGAATATGAAACCTGGACAGCGGGCATTCACGGTAAAAACAACGTGACCT  >  minE/2817585‑2817723

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