Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,747,003 T→C 8.0% intergenic (+87/+122) 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,0030TC 8.0% 172.8 / 5.1 62intergenic (+87/+122)srmB/yfiEATP‑dependent RNA helicase/predicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (25/32);  new base C (2/3);  total (27/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAA  >  minE/1746933‑1747071
                                                                      |                                                                    
aaGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCt                                                                      >  1:1985982/1‑71 (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:2413217/70‑1 (MQ=255)
                      gCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGc                                                 <  1:3160789/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: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)
                                    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: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)
                                                                      |                                                                    
AAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAA  >  minE/1746933‑1747071

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: