Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,747,003 T→C 35.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,0030TC35.8% 112.8 / 87.4 95intergenic (+87/+122)srmB/yfiEATP‑dependent RNA helicase/predicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (21/40);  new base C (15/19);  total (36/59)
Fisher's exact test for biased strand distribution p-value = 3.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAAT  >  minE/1746934‑1747072
                                                                     |                                                                     
aGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCt                                                                       >  1:3889713/1‑70 (MQ=255)
aGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTc                                                                      <  1:3397245/71‑1 (MQ=255)
    gcggcttgagactgatgacaaaggtaaaatTGCCTGATGCGCTACGCTTATCAGGCCTACATGGCCCATGc                                                                  <  1:2736229‑M2/41‑1 (MQ=255)
    gcggcttgagactgatgacaaaggtaaaatTGCCTGATGCGCTACGCTTATCAGGCCTACATGGCCCATGc                                                                  <  1:3883591‑M2/41‑1 (MQ=255)
    gcggcttgagactgatgacaaaggtaaaatTGCCTGATGCGCTACGCTTATCAGGCCTACATGGCCCATGc                                                                  <  1:4351900‑M2/41‑1 (MQ=255)
     cggcttgagactgatgacaaaggtaaaatTGCCTGATGCGCTACGCTTATCAGGCCTACATGGCCCATGCa                                                                 <  1:4894872‑M2/42‑1 (MQ=255)
      gTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCa                                                                 <  1:2254159/70‑1 (MQ=255)
      gTATTAAAAGTGCAGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCa                                                                 <  1:4626023/70‑1 (MQ=255)
        aTTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAata                                                              >  1:2821384/1‑71 (MQ=255)
         ttAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATg                                                                   <  1:5125618/65‑1 (MQ=255)
            aaaGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAatatat                                                           <  1:3474954/70‑1 (MQ=255)
             gactgatgacaaaggtaaaatTGCCTGATGCGCTACGCTTATCAGGCCTACATGGCCCATGCaa                                                                >  1:3005596‑M2/22‑64 (MQ=255)
             aaGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTg                                                         <  1:796947/71‑1 (MQ=255)
              actgatgacaaaggtaaaatTGCCTGATGCGCTACGCTTATCAGGCCTACATGGCCCATGCAATGTATTGa                                                        <  1:2292355‑M2/51‑1 (MQ=255)
              actgatgacaaaggtaaaatTGCCTGATGCGCTACGCTTATCAGGCCTACATGGCCCATGCAATGTATTGa                                                        <  1:2818381‑M2/51‑1 (MQ=255)
               gTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAatata                                                            <  1:2024131/66‑1 (MQ=255)
               gTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGaa                                                       <  1:2751358/71‑1 (MQ=255)
               gTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGa                                                        <  1:3473074/70‑1 (MQ=255)
                tGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAAt                                                      <  1:5074793/71‑1 (MQ=255)
                  cTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAAttt                                                    >  1:1283976/1‑71 (MQ=255)
                   tgacaaaggtaaaatTGCCTGATGCGCTACGCTTATCAGGCCTACATGGCCCATGCAATGTATTGAATTTg                                                   <  1:1024386‑M2/56‑1 (MQ=255)
                     gCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGc                                                  <  1:4691263/70‑1 (MQ=255)
                     gCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGc                                                  <  1:1178496/70‑1 (MQ=255)
                       aaaggtaaaatTGCCTGATGCGCTACGCTTATCAGGCCTACATGGCCCATGCaa                                                                >  1:3650503‑M2/12‑54 (MQ=255)
                         aCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGa                                              >  1:3973708/1‑70 (MQ=255)
                         aCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGAt                                             >  1:1613768/1‑71 (MQ=255)
                         aCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGAt                                             >  1:4283624/1‑71 (MQ=255)
                         aCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGAt                                             >  1:1474693/1‑71 (MQ=255)
                          cGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATc                                            <  1:5576898/71‑1 (MQ=255)
                          cGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATc                                            <  1:947765/71‑1 (MQ=255)
                           gtaaaatTGCCTGATGCGCTACGCTTATCAGGCCTACATGGCCCATGCAATGTATTGaa                                                       <  1:1234328‑M2/52‑1 (MQ=255)
                           gtaaaatTGCCTGATGCGCTACGCTTATCAGGCCTACATGGCCCATGCAATGTATTGAAttt                                                    <  1:289721‑M2/55‑1 (MQ=255)
                           gCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCt                                           <  1:3997055/71‑1 (MQ=255)
                             aaaatTGCCTGATGCGCTACGCTTATCAGGCCTACATGGCCCATGCAATGTATTGAAtt                                                     >  1:392972‑M2/6‑59 (MQ=255)
                             aaaatTGCCTGATGCGCTACGCTTATCAGGCCTACATGGCCCATGCAATGTATTGAATTTGc                                                  >  1:4729318‑M2/6‑62 (MQ=255)
                             aaaatTGCCTGATGCGCTACGCTTATCAGGCCTACATGGCCCATGCAATGTATTGAATTTGCAc                                                >  1:4841015‑M2/6‑64 (MQ=255)
                             aaaaCTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTg                                         <  1:5417556/71‑1 (MQ=255)
                              aaatTGCCTGATGCGCTACGCTTATCAGGCCTACATGGCCCATGCAATGTATTGAATTTGCACGCTTTTGt                                        <  1:3822201‑M2/67‑1 (MQ=255)
                              aaaCTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCa                                                                 <  1:3118839/46‑1 (MQ=255)
                              aaaCTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGt                                        <  1:1747399/71‑1 (MQ=255)
                              aaaCTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGt                                        <  1:5319875/71‑1 (MQ=255)
                              aaaCTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATAAATTGTATTTGCACGATCTTGt                                        <  1:2638081/71‑1 (MQ=255)
                               aaCTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTg                                         >  1:1255450/1‑69 (MQ=255)
                                  tGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGa                                              <  1:4109525/61‑1 (MQ=255)
                                  tGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCAAGa                                              <  1:1051562/61‑1 (MQ=255)
                                  tGCCTGATGCGCTACGCTTATCAGGCCTACATGGCCCATGCAATGTATTGAATTTg                                                   >  1:1216731/1‑56 (MQ=14)
                                  tGCCTGATGCGCTACGCTTATCAGGCCTACATGGCCCATGCAATGTATTGAATTTGc                                                  >  1:3313732/1‑57 (MQ=14)
                                    ccTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCg                                  <  1:3962906/71‑1 (MQ=255)
                                    ccTGATGCGCTACGCTTATCAGGCCTACATGGCCCATGCAATGTATTGAATTTGCACGCTTTTGTAGGCCg                                  <  1:337807/71‑1 (MQ=14)
                                     cTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCg                                  >  1:3841755/1‑70 (MQ=255)
                                      tGATGCGCTACGCTTATCAGGCCTACATGGCCCATGCAATGTATTGAATTTGCACGCTTTTg                                         <  1:5583573/62‑1 (MQ=14)
                                       gATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGc                                                  <  1:3236889/52‑1 (MQ=255)
                                       gATGCGCTACGCTTATCAGGCCTACATGGCCCATGCAATGTATTGAATTTGCACGCTTTTGTAGGCCGGa                                >  1:2652614/1‑70 (MQ=14)
                                        aTGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGAt                               <  1:4543787/70‑1 (MQ=255)
                                        aTGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGAt                               <  1:4336819/70‑1 (MQ=255)
                                        aTGCGCTACGCTTATCAGGCCTACATGGCCCATGCAATGTATTGAATTTGCACGCTTTTGTAGGCCGGAt                               <  1:2012413/70‑1 (MQ=14)
                                         tGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATa                              <  1:5656641/70‑1 (MQ=255)
                                         tGCGCTACGCTTATCAGGCCTACATGGCCCATGCAATGTATTGAATTTGCACGCTTTTGTAGGCCGGATaa                             <  1:3724856/71‑1 (MQ=14)
                                           cgcTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAg                            <  1:537265/70‑1 (MQ=255)
                                           cgcTACGCTTATCAGGCCTACATGGCCCATGCAATGTATTGAATTTGCACGCTTTTg                                         <  1:2101676/57‑1 (MQ=14)
                                           cgcTACGCTTATCAGGCCTACATGGCCCATGCAATGTATTGAATTTGCACGCTTTTg                                         <  1:1601060/57‑1 (MQ=14)
                                            gcTACGCTTATCAGGCCTACATGGCCCATGCAATGTATTGAATTTGCACGCTTTTGTAGGCCGGATAAGGc                          <  1:1839796/71‑1 (MQ=14)
                                             cTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCg                         <  1:3888979/71‑1 (MQ=255)
                                             cTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCg                         <  1:211605/71‑1 (MQ=255)
                                             cTACGCTTATCAGGCCTACATGGCCCATGCAATGTATTGa                                                        <  1:3492326/40‑1 (MQ=14)
                                             cTACGCTTATCAGGCCTACATGGCCCATGCAATGTATTGa                                                        <  1:5555709/40‑1 (MQ=14)
                                             cTACGCTTATCAGGCCTACATGGCCCATGCAATGTATTGAAt                                                      >  1:4402383/1‑42 (MQ=14)
                                               aCGCTTATCAGGCCTACGTGGCTCATGCAatatat                                                           >  1:627051/1‑35 (MQ=255)
                                               aCGCTTATCAGGCCTACGTGGCTCATGCAatatat                                                           >  1:1050744/1‑35 (MQ=255)
                                               aCGCTTATCAGGCCTACGTGGCTCATGCAatatat                                                           >  1:3085331/1‑35 (MQ=255)
                                               aCGCTTATCAGGCCTACGTGGCTCATGCAatatat                                                           >  1:2099620/1‑35 (MQ=255)
                                               aCGCTTATCAGGCCTACGTGGCTCATGCAatatat                                                           >  1:1918981/1‑35 (MQ=255)
                                               aCGCTTATCAGGCCTACGTGGCTCATGCAatatat                                                           >  1:3487589/1‑35 (MQ=255)
                                               aCGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAAtt                                                     <  1:5037317/41‑1 (MQ=255)
                                               aCGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAAtt                                                     <  1:1440307/41‑1 (MQ=255)
                                               aCGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAAtt                                                     <  1:3675104/41‑1 (MQ=255)
                                               aCGCTTATCAGGCCTACATGGCCCATGCAATGTATTGAATTTGCACGCTTTTGTAGGCCGGATaa                             >  1:4356994/1‑65 (MQ=14)
                                                cGCTTATCAGGCCTACATGGCCCATGCAATGTATTGAATTTGCACGCTTTTGTAGGc                                    >  1:2093406/1‑57 (MQ=14)
                                                cGCTTATCAGGCCTACATGGCCCATGCAATGTATTGAATTTGCACGCTTTTGTAGGc                                    >  1:3622613/1‑57 (MQ=14)
                                                cGCTTATCAGGCCTACATGGCCCATGCAATGTATTGAATTTGCACGCTTTTGTAGGCCGGATAAGGCGtt                       >  1:3952435/1‑70 (MQ=14)
                                                cGCTTATCAGGCCTACATGGCCCATGCAATGTATTGAATTTGCACGCTTTTGTAGGCCGGATAAGGCGtt                       >  1:3799248/1‑70 (MQ=14)
                                                 gCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTg                                                   >  1:325515/1‑41 (MQ=255)
                                                 gCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTg                                                   >  1:1157651/1‑41 (MQ=255)
                                                  cTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCAc                    <  1:116695/71‑1 (MQ=255)
                                                   ttATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGtt                       <  1:1500793/67‑1 (MQ=255)
                                                    tATCAGGCCTACATGGCCCATGCAATGTATTGAATTTGCACGCTTTTGTAGGCCGGATAAGGCGTTTAcgc                  >  1:2397133/1‑71 (MQ=11)
                                                       cAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCAcgccg                >  1:2219218/1‑70 (MQ=255)
                                                        aGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCa              <  1:3075374/71‑1 (MQ=255)
                                                           ccTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAgg                           <  1:877077/55‑1 (MQ=255)
                                                           ccTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAgg                           <  1:5315591/55‑1 (MQ=255)
                                                            cTACGTGGCTCATGCAATATATTGAATTTGCACGATc                                            >  1:2937576/1‑37 (MQ=255)
                                                               cGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGc                          >  1:3182663/1‑52 (MQ=255)
                                                                gTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATaa                             >  1:4243018/1‑48 (MQ=255)
                                                                gTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATaa                             >  1:4465721/1‑48 (MQ=255)
                                                                gTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCAt             <  1:3191174/64‑1 (MQ=255)
                                                                gTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCAt             <  1:4993722/64‑1 (MQ=255)
                                                                     tCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAAt  >  1:4811377/1‑70 (MQ=255)
                                                                     tCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAAt  >  1:2646128/1‑70 (MQ=255)
                                                                     |                                                                     
AGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAAT  >  minE/1746934‑1747072

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: