Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,746,998 G→A 35.4% 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,9980GA35.4% 121.6 / 83.7 96intergenic (+82/+127)srmB/yfiEATP‑dependent RNA helicase/predicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (22/40);  new base A (15/19);  total (37/59)
Fisher's exact test for biased strand distribution p-value = 5.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCAT  >  minE/1746929‑1747061
                                                                     |                                                               
gTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGt                                                                <  1:3926903/71‑1 (MQ=255)
  ttAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTgg                                                              >  1:975318/1‑71 (MQ=255)
  ttAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTgg                                                              >  1:216265/1‑71 (MQ=255)
    aaGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGc                                                             <  1:2335264/70‑1 (MQ=255)
     aGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCt                                                            >  1:3889713/1‑70 (MQ=255)
     aGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTc                                                           <  1:3397245/71‑1 (MQ=255)
         gcggcttgagactgatgacaaaggtaaaatTGCCTGATGCGCTACGCTTATCAGGCCTACATGGCCCATGc                                                       <  1:4351900‑M2/41‑1 (MQ=255)
         gcggcttgagactgatgacaaaggtaaaatTGCCTGATGCGCTACGCTTATCAGGCCTACATGGCCCATGc                                                       <  1:3883591‑M2/41‑1 (MQ=255)
         gcggcttgagactgatgacaaaggtaaaatTGCCTGATGCGCTACGCTTATCAGGCCTACATGGCCCATGc                                                       <  1:2736229‑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:2818381‑M2/51‑1 (MQ=255)
                   actgatgacaaaggtaaaatTGCCTGATGCGCTACGCTTATCAGGCCTACATGGCCCATGCAATGTATTGa                                             <  1:2292355‑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:4283624/1‑71 (MQ=255)
                              aCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGAt                                  >  1:1474693/1‑71 (MQ=255)
                              aCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGAt                                  >  1:1613768/1‑71 (MQ=255)
                               cGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATc                                 <  1:947765/71‑1 (MQ=255)
                               cGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATc                                 <  1:5576898/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:4336819/70‑1 (MQ=255)
                                             aTGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGAt                    <  1:4543787/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:1601060/57‑1 (MQ=14)
                                                cgcTACGCTTATCAGGCCTACATGGCCCATGCAATGTATTGAATTTGCACGCTTTTg                              <  1:2101676/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:1918981/1‑35 (MQ=255)
                                                    aCGCTTATCAGGCCTACGTGGCTCATGCAatatat                                                >  1:2099620/1‑35 (MQ=255)
                                                    aCGCTTATCAGGCCTACGTGGCTCATGCAatatat                                                >  1:3085331/1‑35 (MQ=255)
                                                    aCGCTTATCAGGCCTACGTGGCTCATGCAatatat                                                >  1:3487589/1‑35 (MQ=255)
                                                    aCGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAAtt                                          <  1:1440307/41‑1 (MQ=255)
                                                    aCGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAAtt                                          <  1:5037317/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:1157651/1‑41 (MQ=255)
                                                      gCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTg                                        >  1:325515/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:5315591/55‑1 (MQ=255)
                                                                ccTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAgg                <  1:877077/55‑1 (MQ=255)
                                                                 cTACGTGGCTCATGCAATATATTGAATTTGCACGATc                                 >  1:2937576/1‑37 (MQ=255)
                                                                    cGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGc               >  1:3182663/1‑52 (MQ=255)
                                                                     gTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATaa                  >  1:4465721/1‑48 (MQ=255)
                                                                     gTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATaa                  >  1:4243018/1‑48 (MQ=255)
                                                                     gTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCAt  <  1:3191174/64‑1 (MQ=255)
                                                                     gTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCAt  <  1:4993722/64‑1 (MQ=255)
                                                                     |                                                               
GTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCAT  >  minE/1746929‑1747061

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: