Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,629,895 C→A 12.5% V633V (GTG→GTT aegA ← fused predicted FeS binding subunit and predicted NAD/FAD‑binding subunit of oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,629,8950CA12.5% 212.3 / 9.3 70V633V (GTG→GTTaegAfused predicted FeS binding subunit and predicted NAD/FAD‑binding subunit of oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (30/31);  new base A (4/5);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACTGATTTTACCCCCAGCCAGTCAATAATCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAACCACCAGATCCGCACCACGCACGGCGTCACCACCAGCGAAGATTTTCGGATTGGTGGTCTGGTAACGGTAC  >  minE/1629826‑1629964
                                                                     |                                                                     
gACTGATTTTACCCCCAGCCAGTCAATAATCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAacca                                                                      >  1:2190950/1‑71 (MQ=255)
  cTGATTTTACCCCCAGCCAGTCAATAATCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAaccacc                                                                    <  1:889878/71‑1 (MQ=255)
  cTGATTTTACCCCCAGCCAGTCAATAATCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAaccacc                                                                    <  1:4151787/71‑1 (MQ=255)
  cTGATTTTACCCCCAGCCAGTCAATAATCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAaccacc                                                                    <  1:5292382/71‑1 (MQ=255)
     aTTTTACCCCCAGCCAGTCAATAATCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAACCACCAg                                                                  <  1:3077788/70‑1 (MQ=255)
     aTTTTACCCCCAGCCAGTCAATAATCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAACCACCAGa                                                                 <  1:424007/71‑1 (MQ=255)
      ttttACCCCCAGCCAGTCAATAATCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAaccac                                                                     >  1:329004/1‑66 (MQ=255)
            ccccAGCCAGTCAATAATCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAACCACCAGa                                                                 <  1:5186335/64‑1 (MQ=255)
            ccccAGCCAGTCAATAATCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAACCACCAGa                                                                 <  1:2443061/64‑1 (MQ=255)
            ccccAGCCAGTCAATAATCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAACCACCAGATCCGcac                                                          <  1:2857216/71‑1 (MQ=255)
              ccagccagTCAATAATCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAACCACCAGa                                                                 >  1:3663022/1‑62 (MQ=255)
              ccagccagTCAATAATCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAACCACCAGATCCGcacca                                                        >  1:1959509/1‑71 (MQ=255)
              ccagccagTCAATAATCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAACCACCAGATCCGcacca                                                        >  1:2145963/1‑71 (MQ=255)
              ccagccagTCAATAATCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAACCACCAGATCCGcacca                                                        >  1:3935993/1‑71 (MQ=255)
                    agTCAATAATCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAAccacca                                                                   >  1:320343/1‑54 (MQ=255)
                        aatgatCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAACCACCAGa                                                                 <  1:2454803/52‑1 (MQ=255)
                        aataatCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAACCACCAGAt                                                                >  1:4407690/1‑53 (MQ=255)
                        aataatCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAACCACCAGATCCGCACCACGCACgg                                                 >  1:3839548/1‑68 (MQ=255)
                        aataatCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAACCACCAGATCCGCACCACGCACgg                                                 >  1:4137203/1‑68 (MQ=255)
                        aataatCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAACCACCAGATCCGCACCACGCACGGCGt                                              <  1:5431478/71‑1 (MQ=255)
                        aataatCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAACCACCAGATCCGCACCACGCACGGCGt                                              >  1:2434561/1‑71 (MQ=255)
                        aataatCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAACCACCAGATCCGCACCACGCACGGCGt                                              <  1:1828145/71‑1 (MQ=255)
                        aataatACCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAACCACCAGATCCGCACCACGCACGGCGt                                              >  1:3448/1‑71 (MQ=255)
                         ataatCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAACCACCAGATCCGCACCACGCACGGCGt                                              <  1:5437777/70‑1 (MQ=255)
                         ataatCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAACCACCAGATCCGCACCACGCACGGCGTc                                             <  1:2688501/71‑1 (MQ=255)
                           aatCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAACCACCAg                                                                  >  1:1523683/1‑48 (MQ=255)
                           aatCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAACCACCAg                                                                  >  1:1565966/1‑48 (MQ=255)
                           aatCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAACCACCAGATCCGCACCACGCACGGCGTcac                                           >  1:1581383/1‑71 (MQ=255)
                             tCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAACCACCAGATCCGCACCACGCACGGCGTcacca                                         >  1:288794/1‑71 (MQ=255)
                               cccTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAACCACCAGATcc                                                              <  1:3379958/48‑1 (MQ=255)
                               cccTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAACCACCAGATCCGCACCACGCACGGCGTcaccac                                        >  1:2983440/1‑70 (MQ=255)
                               cccTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAACCACCAGATCCGCACCACGCACGGCGTcaccac                                        >  1:2276367/1‑70 (MQ=255)
                               cccTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAACCACCAGATCCGCACCACGCACGGCGTCaccacc                                       <  1:4186157/71‑1 (MQ=255)
                                  cgcgccgccTGCCTTCCTGCGGCCATTGCAGTGACAACCAGATCCGCGCCATGAACTGCATCACCACCAGc                                    <  1:3800687/68‑1 (MQ=255)
                                    tgCCGCATGACGTCCTTCTGCCATTGCGGTAACCACCAGATCCGCACCACGCACGGCGTcaccac                                        >  1:2781339/1‑65 (MQ=255)
                                    tgCCGCATGACGTCCTTCTGCCATTGCGGTAACCACCAGATCCGCACCACGCACGGCGTcaccac                                        >  1:885890/1‑65 (MQ=255)
                                    tgCCGCATGACGTCCTTCTGCCATTGCGGTAACCACCAGATCCGCACCACGCACGGCGTCACCACCAGCg                                   <  1:2499490/70‑1 (MQ=255)
                                    tgCCGCATGACGTCCTTCTGCCATTGCGGTAACCACCAGATCCGCACCACGCACGGCGTCACCACCAGCg                                   <  1:3924902/70‑1 (MQ=255)
                                      ccGCCTGCCTTCCTGCGGCCATTGCAGTGACAACCAGATCCGCGCCATGAACTGCATCACCACCAGca                                   >  1:1818578/1‑67 (MQ=255)
                                         cATGACGTCCTTCTGCCATTGCGGTAACCACCAGATCCGCACCACGCACGGCGTCACCACCAGCGAAGAtt                             <  1:3212266/71‑1 (MQ=255)
                                          aTGACGTCCTTCTGCCATTGCGGTAACCACCAGATCCGCACCACGCACGGCGTCACCACCAGCGAAGAttt                            <  1:2672164/71‑1 (MQ=255)
                                             aCGTCCTTCTGCCATTGCGGTAACCACCAGATCCGCACCACGCACGGCGTCACCACCAGCGAAGATTTTc                          <  1:2910727/70‑1 (MQ=255)
                                                tCCTGCGGCCATTGCAGTGACAACCAGATCCGCGCCATGAACTGCATca                                            >  1:476990/1‑49 (MQ=255)
                                                 ccTTCTGCCATTGCGGTAACCACCAGATCCGTACCACGCACGGCGTCACCACCAGCGAAGATTTTCGGAtt                     >  1:3917505/1‑71 (MQ=255)
                                                 ccTTCTGCCATTGCGGTAACCACCAGATCCGCACCACGCACGGCGTCACCACCAGCGAAGAtt                             <  1:2607303/63‑1 (MQ=255)
                                                 ccTTCTGCCATTGCGGTAACCACCAGATCCGCACCACGCACGGCGTCACCACCAGCGAAGATTTTCGGAtt                     >  1:5446307/1‑71 (MQ=255)
                                                  cTTCTGCCATTGCGGTAACCACCAGATCCGCACCACGCACGGCGTCACCACCAGCGAAGATTTTCgg                        >  1:4583245/1‑67 (MQ=255)
                                                  cTTCTGCCATTGCGGTAACCACCAGATCCGCACCACGCACGGCGTCACCACCAGCGAAGATTTTCGGATtg                    <  1:5380573/71‑1 (MQ=255)
                                                  cTTCTGCCATTGCGGTAACCACCAGATCCGCACCACGCACGGCGTCACCACCAGCGAAGATTTTCGGATtg                    <  1:5371234/71‑1 (MQ=255)
                                                  cTTCTGCCATTGCGGTAACCACCAGATCCGCACCACGCACGGCGTCACCACCAGCGAAGATTTTCGGATtg                    >  1:2885266/1‑71 (MQ=255)
                                                    tCTGCCATTGCGGTAACCACCAGATCCGCACCACGCACGGCGTCACCACCAGCGAAGATTTTCGGATtggt                  >  1:2667186/1‑71 (MQ=255)
                                                    tCTGCCATTGCGGTAACCACCAGATCCGCACCACGCACGGCGTCACCACCAGCGAAGATTTTCGGATtggt                  >  1:2673643/1‑71 (MQ=255)
                                                     cTGCCATTGCGGTAACCACCAGATCCGCACCACGCACGGCGTCACCACCAGCGAAGATTTTCGGa                       >  1:868906/1‑65 (MQ=255)
                                                     cTGCCATTGCGGTAACCACCAGATCCGCACCACGCACGGCGTCACCACCAGCGAAGATTTTCGGa                       >  1:2007297/1‑65 (MQ=255)
                                                        ccATTGCGGTAACCACCAGATCCGCACCACGCACGGCGTCACCACCAGCGAAGATTTTCGGATtggtg                 <  1:4466455/68‑1 (MQ=255)
                                                         cATTGCGGTAACCACCAGATCCGCACCACGCACGGCGTCACCACCAGCGAAGATTTTCGGATTGGTGGTCt             >  1:3575533/1‑71 (MQ=255)
                                                          aTTGCGGTAACCACCAGATCCGCACCACGCACGGCGTCACCACCAGCGAAGATTTTc                          <  1:4542820/57‑1 (MQ=255)
                                                          aTTGCAGTGACAACCAGATCCGCGCCATGAACTGCATCACCACCAGCaaa                                 <  1:1409255/50‑1 (MQ=255)
                                                          aTTGCAGTGACAACCAGATCCGCGCCATGAACTGCATCACCACCAGCaaa                                 <  1:1465982/50‑1 (MQ=255)
                                                          aTTGCAGTGACAACCAGATCCGCGCCATGAACTGCATCACCACCAGCAAAGActtt                           >  1:5448334/1‑56 (MQ=255)
                                                           ttGCGGTAACCACCAGATCCGCACCACGCACGGCGTCACCACCAGCGAAGATTTTCg                         <  1:4321572/57‑1 (MQ=255)
                                                           ttGCAGTGACAACCAGATCCGCGCCATGAACTGCATCACCACCAGCAAAGActttt                          <  1:1510369/56‑2 (MQ=255)
                                                           ttGCAGTGACAACCAGATCCGCGCCATGAACTGCATCACCACCAGCAAAGActttt                          <  1:2644609/56‑2 (MQ=255)
                                                            tGCAGTGACAACCAGATCCGCGCCATGAACTGCATCACCACCAGCAAAGActttt                          >  1:4026992/1‑54 (MQ=255)
                                                             gCGGTAACCACCAGATCCGCACCACGCACGGCGTCACCACCAGCGAAGATTTTCGTATTGGTGGTCTCGt          >  1:787550/1‑70 (MQ=255)
                                                             gCGGTAACCACCAGATCCGCACCACGCACGGCGTCACCACCAGCGAAGATTTTCGGATTGGTGGTCTGGt          <  1:4658968/70‑1 (MQ=255)
                                                             gCGGTAACCACCAGATCCGCACCACGCACGGCGTCACCACCAGCGAAGATTTTCGGATTGGTGGTCTGGt          <  1:5056635/70‑1 (MQ=255)
                                                             gCGGTAACCACCAGATCCGCACCACGCACGGCGTCACCACCAGCGAAGATTTTCGGATTGGTGGTCTGGt          >  1:5336268/1‑70 (MQ=255)
                                                                gTAACCACCAGATCCGCACCACGCACGGCGTcaccac                                        <  1:3927757/37‑1 (MQ=255)
                                                                gTAACCACCAGATCCGCACCACGCACGGCGTcaccac                                        <  1:3093120/37‑1 (MQ=255)
                                                                gTAACCACCAGATCCGCACCACGCACGGCGTcaccac                                        <  1:1387995/37‑1 (MQ=255)
                                                                gTAACCACCAGATCCGCACCACGCACGGCGTCACCACCAGCGAAGa                               <  1:5371855/46‑1 (MQ=255)
                                                                gTAACCACCAGATCCGCACCACGCACGGCGTCACCACCAGCGAAGa                               <  1:206108/46‑1 (MQ=255)
                                                                 tAACCACCAGATCCGCACCACGCACGGCGTCACCACCAGCGAAGATTTTCGGATTGGTGGTCTGGTAACgg     <  1:5342025/71‑1 (MQ=255)
                                                                  aaCCACCAGATCCGCACCACGCACGGCGTCACCACCAGCGAAGATTTTc                          >  1:3175616/1‑49 (MQ=255)
                                                                  aaCCACCAGATCCGCACCACGCACGGCGTCACCACCAGCGAAGATTTTc                          >  1:304173/1‑49 (MQ=255)
                                                                   accaccAGATCCGCACCACGCACGGCGTCACCACCAGCGAAGATTTTCGGATTGGTGGTCTGGTAACGGTa   >  1:1746014/1‑71 (MQ=255)
                                                                   accaccAGATCCGCACCACGCACGGCGTCACCACCAGCGAAGATTTTCGGATTGGTGGTCTGGTAACGGTa   >  1:5180842/1‑71 (MQ=255)
                                                                    ccaccaGATCCGCACCACGCACGGCGTCACCACCAGCGAAGATTTTCGGATTGGTGGTCTGGTAACGGTAc  <  1:5418865/71‑1 (MQ=255)
                                                                    ccaccaGATCCGCACCACGCACGGCGTCACCACCAGCGAAGATTTTCGGATTGGTGGTCTGGTAACGGTAc  <  1:1456046/71‑1 (MQ=255)
                                                                    ccaccaGATCCGCACCACGCACGGCGTCACCACCAGCGAAGATTTTCGGATTGGTGGTCTGGTAACGGTAc  <  1:81339/71‑1 (MQ=255)
                                                                     |                                                                     
GACTGATTTTACCCCCAGCCAGTCAATAATCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAACCACCAGATCCGCACCACGCACGGCGTCACCACCAGCGAAGATTTTCGGATTGGTGGTCTGGTAACGGTAC  >  minE/1629826‑1629964

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: