Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,629,889 G→A 12.3% T635T (ACC→ACT 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,8890GA12.3% 188.4 / 13.6 73T635T (ACC→ACTaegAfused predicted FeS binding subunit and predicted NAD/FAD‑binding subunit of oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (32/32);  new base A (4/5);  total (36/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

AGATTTGACTGATTTTACCCCCAGCCAGTCAATAATCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAACCACCAGATCCGCACCACGCACGGCGTCACCACCAGCGAAGATTTTCGGATTGGTGGTCTGGT  >  minE/1629820‑1629956
                                                                     |                                                                   
agaTTTGACTGATTTTACCCCCAGCCAGTCAATAATCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCgg                                                                    >  1:810388/1‑71 (MQ=255)
 gaTTTGACTGATTTTACCCCCAGCCAGTCAATAATCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGt                                                                   <  1:1386993/71‑1 (MQ=255)
    ttGACTGATTTTACCCCCAGCCAGTCAATAATCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAac                                                                >  1:660831/1‑71 (MQ=255)
      gACTGATTTTACCCCCAGCCAGTCAATAATCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAacca                                                              >  1:2190950/1‑71 (MQ=255)
        cTGATTTTACCCCCAGCCAGTCAATAATCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAaccacc                                                            <  1:889878/71‑1 (MQ=255)
        cTGATTTTACCCCCAGCCAGTCAATAATCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAaccacc                                                            <  1:5292382/71‑1 (MQ=255)
        cTGATTTTACCCCCAGCCAGTCAATAATCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAaccacc                                                            <  1:4151787/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)
              ttACCCCCAGCCAGTCAATAATCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTa                                                                  <  1:1065820/59‑1 (MQ=255)
              ttACCCCCAGCCAGTCAATAATCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTa                                                                  <  1:3513813/59‑1 (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:3935993/1‑71 (MQ=255)
                    ccagccagTCAATAATCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAACCACCAGATCCGcacca                                                >  1:2145963/1‑71 (MQ=255)
                          agTCAATAATCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAAccacca                                                           >  1:320343/1‑54 (MQ=255)
                           gTCAATAATCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTa                                                                  <  1:825756/46‑1 (MQ=255)
                              aatgatCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAACCACCAGa                                                         <  1:2454803/52‑1 (MQ=255)
                              aataatCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAACCACCAGAt                                                        >  1:4407690/1‑53 (MQ=255)
                              aataatCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAACCACCAGATCCGCACCACGCACgg                                         >  1:4137203/1‑68 (MQ=255)
                              aataatCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAACCACCAGATCCGCACCACGCACgg                                         >  1:3839548/1‑68 (MQ=255)
                              aataatCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAACCACCAGATCCGCACCACGCACGGCGt                                      <  1:1828145/71‑1 (MQ=255)
                              aataatCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAACCACCAGATCCGCACCACGCACGGCGt                                      >  1:2434561/1‑71 (MQ=255)
                              aataatCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAACCACCAGATCCGCACCACGCACGGCGt                                      <  1:5431478/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:1565966/1‑48 (MQ=255)
                                 aatCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAACCACCAg                                                          >  1:1523683/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:885890/1‑65 (MQ=255)
                                          tgCCGCATGACGTCCTTCTGCCATTGCGGTAACCACCAGATCCGCACCACGCACGGCGTcaccac                                >  1:2781339/1‑65 (MQ=255)
                                          tgCCGCATGACGTCCTTCTGCCATTGCGGTAACCACCAGATCCGCACCACGCACGGCGTCACCACCAGCg                           <  1:3924902/70‑1 (MQ=255)
                                          tgCCGCATGACGTCCTTCTGCCATTGCGGTAACCACCAGATCCGCACCACGCACGGCGTCACCACCAGCg                           <  1:2499490/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:2885266/1‑71 (MQ=255)
                                                        cTTCTGCCATTGCGGTAACCACCAGATCCGCACCACGCACGGCGTCACCACCAGCGAAGATTTTCGGATtg            <  1:5380573/71‑1 (MQ=255)
                                                        cTTCTGCCATTGCGGTAACCACCAGATCCGCACCACGCACGGCGTCACCACCAGCGAAGATTTTCGGATtg            <  1:5371234/71‑1 (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:1465982/50‑1 (MQ=255)
                                                                aTTGCAGTGACAACCAGATCCGCGCCATGAACTGCATCACCACCAGCaaa                         <  1:1409255/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:5336268/1‑70 (MQ=255)
                                                                   gCGGTAACCACCAGATCCGCACCACGCACGGCGTCACCACCAGCGAAGATTTTCGGATTGGTGGTCTGGt  <  1:5056635/70‑1 (MQ=255)
                                                                   gCGGTAACCACCAGATCCGCACCACGCACGGCGTCACCACCAGCGAAGATTTTCGGATTGGTGGTCTGGt  <  1:4658968/70‑1 (MQ=255)
                                                                     |                                                                   
AGATTTGACTGATTTTACCCCCAGCCAGTCAATAATCCCCTGTGCCGCATGACGTCCTTCTGCCATTGCGGTAACCACCAGATCCGCACCACGCACGGCGTCACCACCAGCGAAGATTTTCGGATTGGTGGTCTGGT  >  minE/1629820‑1629956

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