Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,746,998 G→A 8.0% intergenic (+82/+127) 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,746,9980GA 8.0% 201.9 / 5.8 62intergenic (+82/+127)srmB/yfiEATP‑dependent RNA helicase/predicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (26/31);  new base A (2/3);  total (28/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.65e-01

GTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCA  >  minE/1746929‑1747067
                                                                     |                                                                     
gTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGt                                                                      >  1:3863917/1‑71 (MQ=255)
  ttAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTgg                                                                    >  1:51145/1‑71 (MQ=255)
    aaGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCt                                                                  >  1:1985982/1‑71 (MQ=255)
      gCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACg                                                                       >  1:3586232/1‑64 (MQ=255)
        cTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACg                                                                       >  1:4084912/1‑62 (MQ=255)
        cTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACg                                                                       >  1:2015254/1‑62 (MQ=255)
              ttAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCAAGCAata                                                         >  1:2458556/1‑70 (MQ=255)
                aaaaGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAatatat                                                      >  1:1931771/1‑71 (MQ=255)
                 aaaGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATAtt                                                     >  1:3546306/1‑71 (MQ=255)
                   aGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCt                                                                  >  1:3694874/1‑56 (MQ=40)
                   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)
                            aaaCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCa                                                            <  1:3003091/53‑1 (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:2087123/70‑1 (MQ=255)
                                   aaaCTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTg                                    <  1:996861/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:671991/70‑1 (MQ=255)
                                             aTGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGAt                          <  1:3605692/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:2325422/38‑1 (MQ=255)
                                                                  tACGTGGCTCATGCAATATATTGAATTTGCACGATCtt                                     <  1:1999201/38‑1 (MQ=255)
                                                                  tACGTGGCTCATGCAATATATTGAATTTGCACGATCt                                      <  1:2057335/37‑1 (MQ=255)
                                                                  tACGTGGCTCATGCAATATATTGAATTTGCACGATCt                                      <  1:1617068/37‑1 (MQ=255)
                                                                  tACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGAt                          <  1:2649184/49‑1 (MQ=255)
                                                                    cGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATaa                        >  1:2896352/1‑49 (MQ=255)
                                                                    cGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATaa                        >  1:2978402/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)
                                                                     |                                                                     
GTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCA  >  minE/1746929‑1747067

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: