Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,747,028 A→C 5.7% intergenic (+112/+97) 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,0280AC 5.7% 199.4 / 6.1 70intergenic (+112/+97)srmB/yfiEATP‑dependent RNA helicase/predicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (25/41);  new base C (2/2);  total (27/43)
Fisher's exact test for biased strand distribution p-value = 6.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTA  >  minE/1746959‑1747096
                                                                     |                                                                    
aCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGa                                                                      >  1:1688455/1‑70 (MQ=255)
aCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGa                                                                      >  1:3920529/1‑70 (MQ=255)
    aaaaCTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTg                                                                 <  1:1705711/71‑1 (MQ=255)
    aaaaCTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTg                                                                 <  1:154733/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)
          gCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGcc                                                           >  1:1007775/1‑71 (MQ=255)
               aTGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGAt                                                       <  1:3605692/70‑1 (MQ=255)
               aTGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGAt                                                       <  1:671991/70‑1 (MQ=255)
               aTGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGAt                                                       <  1:2462306/70‑1 (MQ=255)
               aTGCGCTACGCTTATCAGGCCTACATGGCCCATGCAATGTATTGAATTTGCACGCTTTTGTAGGCCGGAt                                                       <  1:3721914/70‑1 (MQ=14)
                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:1198617/71‑1 (MQ=255)
                       cGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTc                                              <  1:3459371/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: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)
                                        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)
                                                      tatTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACtt               <  1:971010/71‑1 (MQ=255)
                                                      tatTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACt                <  1:470337/70‑1 (MQ=255)
                                                       atTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCa                                             >  1:161987/1‑40 (MQ=40)
                                                        ttGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACttt              <  1:1492894/70‑1 (MQ=255)
                                                            aTTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGa         >  1:2924759/1‑71 (MQ=255)
                                                             tttGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGAc        <  1:1293943/71‑1 (MQ=255)
                                                              ttGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGAc        >  1:1957860/1‑70 (MQ=255)
                                                               tGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGa         >  1:3144925/1‑68 (MQ=255)
                                                                gCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACaa      <  1:1149120/70‑1 (MQ=255)
                                                                   cGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTa  <  1:2088843/71‑1 (MQ=255)
                                                                   cGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTa  <  1:3713264/71‑1 (MQ=255)
                                                                   cGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTa  <  1:3575827/71‑1 (MQ=255)
                                                                   cGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTa  <  1:3751677/71‑1 (MQ=255)
                                                                    gATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGt    >  1:1308222/1‑68 (MQ=255)
                                                                    gATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTa  <  1:4004541/70‑1 (MQ=255)
                                                                     aTCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTa  <  1:4055045/69‑1 (MQ=255)
                                                                     |                                                                    
ACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTA  >  minE/1746959‑1747096

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