Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,747,003 T→C 12.8% 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,0030TC12.8% 122.2 / 8.9 47intergenic (+87/+122)srmB/yfiEATP‑dependent RNA helicase/predicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (21/20);  new base C (2/4);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.99e-01

AGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAT  >  minE/1746934‑1747068
                                                                     |                                                                 
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)
                         aCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAatatat                                                       >  1:1611487/1‑57 (MQ=255)
                            cAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGaa                                                   >  1:3935573/1‑58 (MQ=255)
                            cAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGaa                                                   >  1:752795/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)
                                                                  ggCTCATGCAATATATTGAATTTGCACGATCTTGTAGGc                                >  1:1938790/1‑39 (MQ=255)
                                                                    cTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCAc                <  1:2101473/53‑1 (MQ=255)
                                                                     |                                                                 
AGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAT  >  minE/1746934‑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: