Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,747,012 A→G 7.6% intergenic (+96/+113) srmB → / ← yfiE ATP‑dependent RNA helicase/predicted DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,747,0120AG 7.6% 187.3 / 6.4 66intergenic (+96/+113)srmB/yfiEATP‑dependent RNA helicase/predicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (26/35);  new base G (2/3);  total (28/38)
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.90e-01

TTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCG  >  minE/1746943‑1747078
                                                                     |                                                                  
ttAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCAAGCAata                                                                    >  1:2458556/1‑70 (MQ=255)
  aaaaGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAatatat                                                                 >  1:1931771/1‑71 (MQ=255)
   aaaGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATAtt                                                                >  1:3546306/1‑71 (MQ=255)
     aGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTg                                                               >  1:2733293/1‑70 (MQ=255)
       tGCTGGCAGACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAAt                                                            >  1:137012/1‑71 (MQ=255)
       tGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAAt                                                            >  1:1851349/1‑71 (MQ=255)
         cTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAAtt                                                           <  1:122205/70‑1 (MQ=255)
            gCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTg                                                         <  1:746896/69‑1 (MQ=255)
            gCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGc                                                        <  1:3160789/70‑1 (MQ=255)
            gCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGc                                                        <  1:2413217/70‑1 (MQ=255)
             cAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAAttt                                                          >  1:3294067/1‑67 (MQ=255)
              aaaCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCAc                                                      >  1:4054176/1‑70 (MQ=255)
              aaaCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACg                                                     <  1:934481/71‑1 (MQ=255)
                aCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGa                                                    >  1:1688455/1‑70 (MQ=255)
                aCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGa                                                    >  1:3920529/1‑70 (MQ=255)
                    aaaaCTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTg                                               <  1:154733/71‑1 (MQ=255)
                    aaaaCTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTg                                               <  1:1705711/71‑1 (MQ=255)
                     aaaCTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTg                                               <  1:996861/70‑1 (MQ=255)
                     aaaCTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTg                                               <  1:2087123/70‑1 (MQ=255)
                       aCTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAg                                            >  1:3672688/1‑71 (MQ=255)
                        cTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAgg                                           >  1:2490902/1‑71 (MQ=255)
                         tGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTg                                               <  1:763461/66‑1 (MQ=255)
                         tGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTa                                             >  1:1160001/1‑68 (MQ=255)
                          gCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTg                                                         <  1:3119280/55‑1 (MQ=255)
                          gCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGcc                                         >  1:1007775/1‑71 (MQ=255)
                              gATGCGCTACGCTTATCAGGCCTACATGGCCCATGCAATGTATTGAATTTg                                                         <  1:2628304/51‑1 (MQ=14)
                               aTGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGAt                                     <  1:2462306/70‑1 (MQ=255)
                               aTGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGAt                                     <  1:3605692/70‑1 (MQ=255)
                               aTGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGAt                                     <  1:671991/70‑1 (MQ=255)
                               aTGCGCTACGCTTATCAGGCCTACATGGCCCATGCAATGTATTGAATTTGCACGCTTTTGTAGGCCGGAt                                     <  1:3721914/70‑1 (MQ=14)
                                tGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACg                                                     <  1:2074871/53‑1 (MQ=255)
                                tGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGAt                                                   >  1:1513567/1‑55 (MQ=255)
                                    cTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCg                               <  1:3979899/71‑1 (MQ=255)
                                     tACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGt                              >  1:428956/1‑71 (MQ=255)
                                      aCGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTg                                               >  1:3517398/1‑53 (MQ=255)
                                      aCGCTTATCAGGCCTACATGGCCCATGCAATGTATTGAATTTGCACGCTTTTGTAGGCCGGATAAGGCGtt                             >  1:2921837/1‑71 (MQ=14)
                                       cGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTc                            <  1:3459371/71‑1 (MQ=255)
                                       cGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTc                            <  1:1198617/71‑1 (MQ=255)
                                       cGCTTATCAGGCCTACATGGCCCATGCAATGTATTGAATTTGCACGCTTTTGTAGGCCGGATAAGGCGtt                             <  1:2649708/70‑1 (MQ=14)
                                        gCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTa                                             >  1:1034879/1‑53 (MQ=255)
                                         cTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTg                                               <  1:3115032/50‑1 (MQ=255)
                                            aTCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCAcg                         <  1:3878535/69‑1 (MQ=255)
                                             tCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCAcgccg                      <  1:1596953/71‑1 (MQ=255)
                                               aGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGt                              <  1:3085821/61‑1 (MQ=255)
                                               aGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCa                    >  1:1144103/1‑71 (MQ=255)
                                               aGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCa                    <  1:2810920/71‑1 (MQ=255)
                                               aGGCCTACATGGCCCATGCAATGTATTGAATTTGCACGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCa                    >  1:3702629/1‑71 (MQ=11)
                                                    tACGTGGCTCATGCAATATATTGAATTTGCACGATCtt                                                <  1:1999201/38‑1 (MQ=255)
                                                    tACGTGGCTCATGCAATATATTGAATTTGCACGATCtt                                                <  1:2325422/38‑1 (MQ=255)
                                                    tACGTGGCTCATGCAATATATTGAATTTGCACGATCt                                                 <  1:1617068/37‑1 (MQ=255)
                                                    tACGTGGCTCATGCAATATATTGAATTTGCACGATCt                                                 <  1:2057335/37‑1 (MQ=255)
                                                    tACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGAt                                     <  1:2649184/49‑1 (MQ=255)
                                                      cGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATaa                                   >  1:2978402/1‑49 (MQ=255)
                                                      cGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATaa                                   >  1:2896352/1‑49 (MQ=255)
                                                      cGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg               <  1:3837988/69‑1 (MQ=255)
                                                      cGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc              <  1:580446/70‑1 (MQ=255)
                                                      cGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa             >  1:3792565/1‑71 (MQ=255)
                                                       gTGGCTCATGCAATATATTGAATTTGCACGATCTTGTa                                             >  1:1910464/1‑38 (MQ=255)
                                                        tGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCg                                        <  1:759210/42‑1 (MQ=255)
                                                          gCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCAcgccg                      >  1:792362/1‑58 (MQ=255)
                                                          gCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaa         >  1:1902471/1‑71 (MQ=255)
                                                              aTGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATa       <  1:2955016/69‑1 (MQ=255)
                                                              aTGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATa       >  1:3990894/1‑69 (MQ=255)
                                                               tGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAg    <  1:3406973/71‑1 (MQ=255)
                                                                 cAATATATTGAATTTGCACGATCTTGTAGGCCGGATaa                                   >  1:421252/1‑38 (MQ=255)
                                                                 cAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATc                  >  1:2937799/1‑55 (MQ=255)
                                                                 cAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa             <  1:2650320/60‑1 (MQ=255)
                                                                 cAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAgcg  <  1:261330/71‑1 (MQ=255)
                                                                   atatatTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATaa      <  1:2485958/65‑1 (MQ=255)
                                                                   atatatTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAgcg  <  1:2130416/69‑1 (MQ=255)
                                                                    tatatTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATaa      <  1:1206774/64‑1 (MQ=255)
                                                                     atatTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATc                  >  1:800510/1‑51 (MQ=38)
                                                                     atatAGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCa                           >  1:1405506/1‑42 (MQ=255)
                                                                     |                                                                  
TTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCG  >  minE/1746943‑1747078

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: