Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,746,998 G→A 12.5% 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,9980GA12.5% 135.5 / 8.4 48intergenic (+82/+127)srmB/yfiEATP‑dependent RNA helicase/predicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (19/23);  new base A (2/4);  total (21/27)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

GGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAT  >  minE/1746928‑1747068
                                                                      |                                                                      
ggTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACg                                                                        >  1:3622478/1‑71 (MQ=255)
ggTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACg                                                                        >  1:1633053/1‑71 (MQ=255)
    tAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTgg                                                                     <  1:2071561/70‑1 (MQ=255)
    tAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTgg                                                                     <  1:1888175/70‑1 (MQ=255)
    tAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTgg                                                                     <  1:2834317/70‑1 (MQ=255)
      aGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTc                                                                  <  1:3535206/71‑1 (MQ=255)
                 aaaaGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAatatat                                                       >  1:867269/1‑71 (MQ=255)
                 aaaaGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAatata                                                        <  1:3585763/70‑1 (MQ=255)
                   aaGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATAtt                                                      <  1:1563487/70‑1 (MQ=255)
                      tGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAAt                                                  <  1:1155244/71‑1 (MQ=255)
                      tGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAAt                                                  <  1:939325/71‑1 (MQ=255)
                      tGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAAt                                                  <  1:1644582/71‑1 (MQ=255)
                         tGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAAttt                                                <  1:2879748/70‑1 (MQ=255)
                           gCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCa                                                             <  1:3253647/55‑1 (MQ=255)
                           gCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAAt                                                  <  1:959292/66‑1 (MQ=255)
                           gCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCa                                             >  1:104539/1‑71 (MQ=255)
                           gCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCa                                             >  1:2771160/1‑71 (MQ=255)
                             aaaCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAAttt                                                >  1:2139686/1‑66 (MQ=255)
                               aCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGc                                                                    <  1:284063/44‑1 (MQ=32)
                               aCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAatatat                                                       >  1:1611487/1‑57 (MQ=255)
                                  cAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGaa                                                   >  1:752795/1‑58 (MQ=255)
                                  cAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGaa                                                   >  1:3935573/1‑58 (MQ=255)
                                  cAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCtt                                      <  1:1695829/71‑1 (MQ=255)
                                     aaCTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCAc                                            >  1:2865907/1‑62 (MQ=255)
                                     aaCTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTg                                     >  1:2462775/1‑69 (MQ=255)
                                       cTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTg                                     <  1:3750802/67‑1 (MQ=255)
                                          ccTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCg                              >  1:917321/1‑71 (MQ=255)
                                          ccTGATGCGCTACGCTTATCAGGCCTACATGGCCCATGCAATGTATTGAATTTGCACGCTTTTGTAGGCCg                              <  1:3936422/71‑1 (MQ=14)
                                           cTGATGCGCTACGCTTATCAGGCCTACATGGCCCATGCAATGTATTGAATTTGCACGCTTTTGTAGGCCgg                             <  1:3887687/71‑1 (MQ=14)
                                            tGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGa                            >  1:2348084/1‑71 (MQ=255)
                                             gATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAatata                                                        <  1:3145273/42‑1 (MQ=255)
                                             gATGCGCTACGCTTATCAGGCCTACATGGCCCATGCaa                                                            >  1:2562083/1‑38 (MQ=14)
                                              aTGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAat                                                           >  1:2899883/1‑38 (MQ=255)
                                              aTGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGAt                                         >  1:1146702/1‑56 (MQ=255)
                                              aTGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGAt                                         >  1:3522759/1‑56 (MQ=255)
                                               tGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATa                          <  1:814088/70‑1 (MQ=255)
                                                    tACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGt                    >  1:2223844/1‑71 (MQ=255)
                                                      cGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTc                  <  1:1468605/71‑1 (MQ=255)
                                                      cGCTTATCAGGCCTACATGGCCCATGCAATGTATTGAATTTGCACGCTTTTGTAGGCCGGATAAGGCGttt                  <  1:483224/71‑2 (MQ=11)
                                                      cGCTTATCAGGCCTACATGGCCCATGCAATGTATTGAATTTGCACGCTTTTGTAGGCCGGATAAGGCGtt                   <  1:1844694/70‑1 (MQ=14)
                                                        cTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCAc                <  1:1381199/71‑1 (MQ=255)
                                                         ttATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCAcg               >  1:41877/1‑71 (MQ=255)
                                                           aTCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCAcg               <  1:806968/69‑1 (MQ=255)
                                                              aGGCCTACGTGGCTCATGCAATATATTGAATTTGCa                                             >  1:3247176/1‑36 (MQ=255)
                                                              aGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCAcgccgc           <  1:2306152/70‑1 (MQ=255)
                                                                 ccTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGc                      <  1:1517023/56‑1 (MQ=255)
                                                                  cTACATGGCCCATGCAATGTATTGAATTTGCACGCTTTTGTAGGCCGGATAAGGCGtt                   >  1:746952/1‑58 (MQ=11)
                                                                    aCGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTa                                   >  1:1853025/1‑40 (MQ=255)
                                                                    aCGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGa                            >  1:3157007/1‑47 (MQ=255)
                                                                     cGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa   <  1:1191098/71‑1 (MQ=255)
                                                                      gTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAt  >  1:1880728/1‑71 (MQ=255)
                                                                      |                                                                      
GGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAT  >  minE/1746928‑1747068

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