Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,747,030 C→T 5.6% intergenic (+114/+95) 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,0300CT 5.6% 235.4 / 2.9 71intergenic (+114/+95)srmB/yfiEATP‑dependent RNA helicase/predicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (25/42);  new base T (2/2);  total (27/44)
Fisher's exact test for biased strand distribution p-value = 6.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.53e-01

AAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAAT  >  minE/1746963‑1747099
                                                                   |                                                                     
aaaaCTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTg                                                                    <  1:1705711/71‑1 (MQ=255)
aaaaCTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTg                                                                    <  1:154733/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)
      gCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGcc                                                              >  1:1007775/1‑71 (MQ=255)
           aTGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGAt                                                          <  1:671991/70‑1 (MQ=255)
           aTGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGAt                                                          <  1:3605692/70‑1 (MQ=255)
           aTGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGAt                                                          <  1:2462306/70‑1 (MQ=255)
           aTGCGCTACGCTTATCAGGCCTACATGGCCCATGCAATGTATTGAATTTGCACGCTTTTGTAGGCCGGAt                                                          <  1:3721914/70‑1 (MQ=14)
                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: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)
                                    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:3713264/71‑1 (MQ=255)
                                                               cGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTa     <  1:3575827/71‑1 (MQ=255)
                                                               cGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTa     <  1:2088843/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)
                                                                  tCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAAt  >  1:449857/1‑71 (MQ=255)
                                                                   |                                                                     
AAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAAT  >  minE/1746963‑1747099

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: