Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,823,089 A→G 9.6% L320L (TTA→TTG norW → NADH:flavorubredoxin oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,823,0890AG 9.6% 136.5 / 6.6 52L320L (TTA→TTGnorWNADH:flavorubredoxin oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (28/19);  new base G (2/3);  total (30/22)
Fisher's exact test for biased strand distribution p-value = 6.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCTTCTCGGCAATAACACGCCGCTGAAACTCCCGGCGATGCTGGTGAAAATCAAAACGCCGGAATTACCGCTGCATCTGGCAGGCGAAACCCAGCGTCAGGATTTACGCTGGCAAATTAATACCGAACGC  >  minE/1823022‑1823152
                                                                   |                                                               
aTCTTCTCGGCAATAACACGCCGCTGAAACTCCCGGCGATGCTGGTGAAAATCAAAACGCCGGAATTACCg                                                              >  1:1269668/1‑71 (MQ=255)
 tcttctCGGCAATAACACGCCGCTGAAACTCCCGGCGATGCTGGTGAAAATCAAAACGCCGGAATTGCCGc                                                             >  1:202922/1‑71 (MQ=255)
   ttctCGGCAATAACACGCCGCTGAAACTCCCGGCGATGCTGGTGAAAATCAAAACGCCGGAATTACCGCTg                                                           >  1:93463/1‑71 (MQ=255)
   ttctCGGCAATAACACGCCGCTGAAACTCCCGGCGATGCTGGTGAAAATCAAAACGCCGGAATTACCGCTg                                                           >  1:1808349/1‑71 (MQ=255)
   ttctCGGCAATAACACGCCGCTGAAACTCCCGGCGATGCTGGTGAAAATCAAAACGCCGGAATTACCGCTg                                                           >  1:1589819/1‑71 (MQ=255)
    tctcGGCAATAACACGCCGCTGAAACTCCCGGCGATGCTGGTGAAAATCAAAACGCCGGAATTAc                                                                <  1:1360882/65‑1 (MQ=255)
             tAACACGCCGCTGAAACTCCCGGCGATGCTGGTGAAAATCAAAACGCCGGAATTACCGCTGCATc                                                       <  1:1060730/65‑1 (MQ=255)
              aaCACGCCGCTGAAACTCCCGGCGATGCTGGTGAAAATCAAAACGCCGGAATTACCGCTGCATCTGGCa                                                  >  1:2234597/1‑69 (MQ=255)
               acacGCCGCTGAAACTCCCGGCGATGCTGGTGAAAATCAAAACGCCGGAATTACCGCTGCATCTGGCAGGc                                               >  1:2070355/1‑71 (MQ=255)
                 acGCCGCTGAAACTCCCGGCGATGCTGGTGAAAATCAAAACGCCGGAATTACCGCTGCATCTGGCAGGcta                                             <  1:2585401/71‑3 (MQ=255)
                 acGCCGCTGAAACTCCCGGCGATGCTGGTGAAAATCAAAACGCCGGAATTACCGCTGCATCTGGCAGGCGa                                             <  1:1507083/71‑1 (MQ=255)
                     cgcTGAAACTCCCGGCGATGCTGGTGAAAATCAAAACGCCGGAATTACCGCTGCATCTGGCAGGCGAAAcc                                         >  1:781522/1‑71 (MQ=255)
                     cgcTGAAACTCCCGGCGATGCTGGTGAAAATCAAAACGCCGGAATTACCGCTGCATCTGGCAGGCGAAAcc                                         <  1:2075290/71‑1 (MQ=255)
                     cgcTGAAACTCCCGGCGATGCTGGTGAAAATCAAAACGCCGGAATTACCGCTGCATCTGGCAGGCGAAAcc                                         >  1:2099300/1‑71 (MQ=255)
                     cgcTGAAACTCCCGGCGATGCTGGTGAAAATCAAAACGCCGGAATTACCGCTGCATCTGGCAGGCGAAAc                                          <  1:323172/70‑1 (MQ=255)
                         gAAACTCCCGGCGATGCTGGTGAAAATCAAAACGCCGGAATTACCGCTGCATCTGGCAGGCGAAACCCAGc                                     >  1:1591159/1‑71 (MQ=255)
                         gAAACTCCCGGCGATGCTGGTGAAAATCAAAACGCCGGAATTACCGCTGCATCTGGCAGGCGAAACCCAGc                                     >  1:2108872/1‑71 (MQ=255)
                          aaaCTCCCGGCGATGCTGGTGAAAATCAAAACGCCGGAATTACCGCTGCATCTGGCAGGCGAAACCCAGCg                                    >  1:3647007/1‑71 (MQ=255)
                          aaaCTCCCGGCGATGCTGGTGAAAATCAAAACGCCGGAATTACCCCTGCATCTGGCAGGCGAAACCCAGCg                                    >  1:4119880/1‑71 (MQ=255)
                           aaCTCCCGGCGATGCTGGTGAAAATCAAAACGCCGGAATTGCCGCTGCATCTGGCAGGCGAAACCCAGCGt                                   >  1:175388/1‑71 (MQ=255)
                               cccGGCGATGCTGGTGAAAATCAAAACGCCGGAATTACCGCTGCATCTGGCATGCGAAACCCAGCGTCAg                                >  1:235925/1‑70 (MQ=255)
                                ccGGCGATGCTGGTGAAAATCAAAACGCCGGAATTGCCGCTGCATCTGGCAGGCGAAACCCAGCGTCAg                                <  1:967395/69‑1 (MQ=255)
                                ccGGCGATGCTGGTGAAAATCAAAACGCCGGAATTACCGCTGCATCTGGCAGGCGAAACCCAGCGTCAgg                               <  1:2705618/70‑1 (MQ=255)
                                 cGGCGATGCTGGTGAAAATCAAAACGCCGGAATTACCGCTGCATCTGGCAGGCGAAACCCAGCGTCAGGAt                             >  1:3805540/1‑71 (MQ=255)
                                 cGGCGATGCTGGTGAAAATCAAAACGCCGGAATTACCGCTGCATCTGGCAGGCGAAACCCAGCGTCAGGAt                             >  1:3417501/1‑71 (MQ=255)
                                  ggCGATGCTGGTGAAAATCAAAACGCCGGAATTACCGCTGCATCTGGCAGGCGAAACCCAGCGTCAGGAtt                            <  1:4024323/71‑1 (MQ=255)
                                  ggCGATGCTGGTGAAAATCAAAACGCCGGAATTACCGCTGCATCTGGCAGGCGAAACCCAGCGTCAGGAtt                            >  1:78357/1‑71 (MQ=255)
                                  ggCGATGCTGGTGAAAATCAAAACGCCGGAATTACCGCTGCATCTGGCAGGCGAAACCCAGCGTCAGGAtt                            <  1:1167939/71‑1 (MQ=255)
                                    cGATGCTGGTGAAAATCAAAACGCCGGAATTACCGCTGCATCTGGCAGGCGAAACCCAGCGTCAGGATTTa                          >  1:3915324/1‑71 (MQ=255)
                                       tGCTGGTGAAAATCAAAACGCCGGAATTACCGCTGCATCTGGCAGGCGAAACCCAGCGTCAGGATTTACGc                       <  1:629891/71‑1 (MQ=255)
                                        gCTGGTGAAAATCAAAACGCCGGAATTACCGCTGCATCTGGCAGGAGAAACCCAGCGTCAGGATTTACGCt                      >  1:3146888/1‑71 (MQ=255)
                                         cTGGTGAAAATCAAAACGCCGGAATTACCGCTGCATCTGGCAGGCGAAACCCAGCGTCAGGATTTACGCTg                     >  1:1437981/1‑71 (MQ=255)
                                         cTGGTGAAAATCAAAACGCCGGAATTACCGCTGCATCTGGCAGGCGAAACCCAGCGTCAGGATTTACGCTg                     >  1:465125/1‑71 (MQ=255)
                                           ggTGAAAATCAAAACGCCGGAATTACCGCTGCATCTGGCAGGCGAAACCCAGCGTCAGGATTTACGCTgg                    >  1:3112475/1‑70 (MQ=255)
                                            gTGAAAATCAAAACGCCGGAATTACCGCTGCATCTGGCAGGCGAAACCCAGCGTCAGGATTTACGCTgg                    <  1:1410300/69‑1 (MQ=255)
                                            gTGAAAATCAAAACGCCGGAATTACCGCTGCATCTGGCAGGCGAAACCCAGCGTCAGGATTTACGCTgg                    <  1:3155043/69‑1 (MQ=255)
                                            gTGAAAATCAAAACGCCGGAATTACCGCTGCATCTGGCAGGCGAAACCCAGCGTCAGGATTTACGCTGGc                   >  1:579475/1‑70 (MQ=255)
                                             tGAAAATCAAAACGCCGGAATTACCGCTGCATCTGGCAGGCGAAACCCAGCGTCAGGATTTACGCt                      >  1:2983975/1‑66 (MQ=255)
                                                 aaTCAAAACGCCGGAATTGCCGCTGCATCTGGCAGGCGAAACCCAGc                                     <  1:2584354/47‑1 (MQ=255)
                                                   tCAAAACGCCGGAATTACCGCTGCATCTGGCAGGCGAAACCCAGCGTCAGGATTTACGCTGGCAAATTaa            >  1:1821486/1‑70 (MQ=255)
                                                   tCAAAACGCCGGAATTACCGCTGCATCTGGCAGGCGAAACCCAGCGTCAGGATTTACGCTGGCAAATTAAt           <  1:3202004/71‑1 (MQ=255)
                                                    cAAAACGCCGGAATTACCGCTGCATCTGGCAGGCGAAAccc                                        >  1:1052171/1‑41 (MQ=255)
                                                     aaaaCGCCGGAATTACCGCTGCATCTGGCAGGCGAAACCCAGCGTCAGGATTTACGCTGGCAAATTAATAc         <  1:3132585/71‑1 (MQ=255)
                                                       aaCGCCGGAATTACCGCTGCATCTGGCATGCGAAACCCAGCGTCAGGATTTACGCTGGCAAATTAATACCg       <  1:540771/71‑1 (MQ=255)
                                                       aaCGCCGGAATTACCGCTGCATCTGGCAGGCGAAACCCAGCGTCAGGATTTACGCTGGCAAATTAATACCg       <  1:2105474/71‑1 (MQ=255)
                                                            cGGAATTGCCGCTGCATCTGGCAGGCGAAACCCAGCGTCAGGATTTACGCTGGCAAATTAATACCGAACGc  <  1:1692123/71‑1 (MQ=255)
                                                              gAATTACCGCTGCATCTGGCAGGCGAAACCCAGCGTCAGGAtt                            <  1:3229901/43‑1 (MQ=255)
                                                              gAATTACCGCTGCATCTGGCAGGCGAAACCCAGCGTCAGGAtt                            <  1:316312/43‑1 (MQ=255)
                                                              gAATTACCGCTGCATCTGGCAGGCGAAACCCAGCGTCAGGAtt                            <  1:2392666/43‑1 (MQ=255)
                                                              gAATTACCGCTGCATCTGGCAGGCGAAACCCAGCGTCAGGATTTACGCTGGCAAATTAATAc         >  1:2209607/1‑62 (MQ=255)
                                                              gAATTACCGCTGCATCTGGCAGGCGAAACCCAGCGTCAGGATTTACGCTGGCAAATTAATAc         >  1:2418536/1‑62 (MQ=255)
                                                              gAATTACCGCTGCATCTGGCAGGCGAAACCCAGCGTCAGGATTTACGCTGGCAAATTAATAc         >  1:2352524/1‑62 (MQ=255)
                                                                   aCCGCTGCATCTGGCAGGCGAAACCCAGCGTCAGGATTTACGCTGGCAAATTAAt           >  1:1945300/1‑55 (MQ=255)
                                                                   |                                                               
ATCTTCTCGGCAATAACACGCCGCTGAAACTCCCGGCGATGCTGGTGAAAATCAAAACGCCGGAATTACCGCTGCATCTGGCAGGCGAAACCCAGCGTCAGGATTTACGCTGGCAAATTAATACCGAACGC  >  minE/1823022‑1823152

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: