Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,963,700 C→T 8.0% A108T (GCC→ACC)  yggD ← predicted DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,963,7000CT 8.0% 221.5 / 7.2 75A108T (GCC→ACC) yggDpredicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (19/50);  new base T (3/3);  total (22/53)
Fisher's exact test for biased strand distribution p-value = 3.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

CGTCATAGAGCATCGCATTGCGGGTTATGCTATTGACGTTGCTGATAAAATCCCAGGTCATGTCATCGGCAAATCCGGGCGTTTCATCCTGTTCGTTCAGGTAATGCCAGAACTGGGAAAAATGCGTGATATCGGCGTA  >  minE/1963631‑1963769
                                                                     |                                                                     
cgTCATAGAGCATCGCATTGCGGGTTATGCTATTGACGTTGCTGATAAAATCCCAGGTCATGTCATCGGc                                                                       <  1:3852976/70‑1 (MQ=255)
 gTCATAGAGCATCGCATTGCGGGTTATGCTATTGACGTTGCTGATAAAATCCCAGGTCATGTCATCGGCaa                                                                     <  1:34411/71‑1 (MQ=255)
  tCATAGAGCATCGCATTGCGGGTTATGCTATTGACGTTGCTGATAAAATCCCAGGTCATGTCATCGGCaaa                                                                    <  1:3375317/71‑1 (MQ=255)
    aTAGAGCATCGCATTGCGGGTTATGCTATTGACGTTGCTGATAAAATCCCAGGTCATGTCATCGGCAAATc                                                                  <  1:258728/71‑1 (MQ=255)
     tAGAGCATCGCATTGCGGGTTATGCTATTGACGTTGCTGATAAAATCCCAGGTCATGTCATCGGCAAATcc                                                                 <  1:673926/71‑1 (MQ=255)
     tAGAGCATCGCATTGCGGGTTATGCTATTGACGTTGCTGATAAAATCCCAGGTCATGTCATCGGCAAATcc                                                                 <  1:1474884/71‑1 (MQ=255)
      agagCATCGCATTGCGGGTTATGCTATTGACGTTGCTGATAAAATCCCAGGTCATGTCATCGGCAAATcc                                                                 >  1:1253237/1‑70 (MQ=255)
        agCATCGCATTGCGGGTTATGCTATTGACGTTGCTGATAAAATCCCAGGTCATGTCATCGGCAAATCCggg                                                              >  1:3128717/1‑71 (MQ=255)
          cATCGCATTGCGGGTTATGCTATTGACGTTGCTGATAAAATCCCAGGTCATGTCATCGGCAAATc                                                                  >  1:959717/1‑65 (MQ=255)
            tCGCATTGCGGGTTATGCTATTGACGTTGCTGATAAAATCCCAGGTCATGTCATCGGCAAATCCGGGAGtt                                                          >  1:666216/1‑71 (MQ=255)
            tCGCATTGCGGGTTATGCTATTGACGTTGCTGATAAAATCCCAGGTCATGTCATCGGCAAATCCGGGAGt                                                           <  1:2317433/70‑1 (MQ=255)
                 ttGCGGGTTATGCTATTGACGTTGCTGATAAAATCCCAGGTCATGTCATCGGCAAATCCGGGCGTTTCATc                                                     >  1:335311/1‑71 (MQ=255)
                      ggTTATGCTATTGACGTTGCTGATAAAATCCCAGGTCATGTCATCGGCAAATCCGGGCGTTTCATCCTgtt                                                >  1:1604698/1‑71 (MQ=255)
                      ggTTATGCTATTGACGTTGCTGATAAAATCCCAGGTCATGTCATCGGCAAATCCGGGCGTTTCATCCTgt                                                 >  1:767515/1‑70 (MQ=255)
                        ttATGCTATTGACGTTGCTGATAAAATCCCAGGTCATGTCATCGGCAAATc                                                                  <  1:2746269/51‑1 (MQ=255)
                        ttATGCTATTGACGTTGCTGATAAAATCCCAGGTCATGTCATCGGCAAATCCGGGCGTTTCATCCTgttcg                                              <  1:3817167/71‑1 (MQ=255)
                           tGCTATTGACGTTGCTGATAAAATCCCAGGTCATGTCATCGGCAAATCCGGGCGTTTCATCCTgttcgttc                                           <  1:942210/71‑1 (MQ=255)
                             cTATTGACGTTGCTGATAAAATCCCAGGTCATGTCATCGGCAAATCCGGGCGTTTCATCCTGTTCGTTCa                                          <  1:3570027/70‑1 (MQ=255)
                             cTATTGACGTTGCTGATAAAATCCCAGGTCATGTCATCGGCAAATCCGGGCGTTTCATCCTGTTCGTTCa                                          <  1:419218/70‑1 (MQ=255)
                             cTATTGACGTTGCTGATAAAATCCCAGGTCATGTCATCGGCAAATCCGGGCGTTTCATCCTGTTCGTTCAg                                         <  1:3186307/71‑1 (MQ=255)
                              tATTGACGTTGCTGATAAAATCCCAGGTCATGTCATCGGCAAATCCGGGCGTTTCATCCTGTTCGTTCa                                          <  1:2827323/69‑1 (MQ=255)
                              tATTGACGTTGCTGATAAAATCCCAGGTCATGTCATCGGCAAATCCGGGCGTTTCATCCTGTTCGTTCa                                          <  1:2822803/69‑1 (MQ=255)
                               aTTGACGTTGCTGATAAAATCCCAGGTCATGTCATCGGCAAATCCGGGCGTTTCATCCTGTTCGTTCAgg                                        >  1:762154/1‑70 (MQ=255)
                               aTTGACCTTGCTGATAAAATCCCAGGTCATGTCATCGGCAAATCCGGGCGTTTCATCCTGATCGTTCAgg                                        >  1:766809/1‑70 (MQ=255)
                                ttGACGTTGCTGATAAAATCCCAGGTCATGTCATCGGCAAATCCGGGCGt                                                           <  1:2807129/50‑1 (MQ=255)
                                 tGACGTTGCTGATAAAATCCCAGGTCATGTCATCGGCAAATCCGGGCGttt                                                         <  1:513643/51‑1 (MQ=255)
                                  gACGTTGCTGATAAAATCCCAGGTCATGTCATCGGCAAATCCGGGCGTTTCATCCTGTTCGTTCAg                                         >  1:3539354/1‑66 (MQ=255)
                                   aCGTTGCTGATAAAATCCCAGGTCATGTCATCGGCAAATCCGGGCGTTTCATcc                                                    <  1:2498111/54‑1 (MQ=255)
                                   aCGTTGCTGATAAAATCCCAGGTCATGTCATCGGCAAATCCGGGCGTTTCATCCTGTTCGTTCAGGTAAt                                    <  1:3751539/70‑1 (MQ=255)
                                     gTTGCTGATAAAATCCCAGGTCATGTCATCGGCAAATCCGGGCGTTTCATCCTgttcgttc                                           <  1:3665195/61‑1 (MQ=255)
                                     gTTGCTGATAAAATCCCAGGTCATGTCATCGGCAAATCCGGGCGTTTCATCCTGTTCGTTCAGGTAATGcc                                 <  1:540453/71‑1 (MQ=255)
                                       tGCTGATAAAATCCCAGGTCATGTCATCGGCAAATCCGGGCGTTTCATc                                                     >  1:1244797/1‑49 (MQ=255)
                                        gCTGATAAAATCCCAGGTCATGTCATCGGCAAATCCGGGCGTTTCa                                                       >  1:3721264/1‑46 (MQ=255)
                                         cTGATAAAATCCCAGGTCATGTCATCGGTAAATCCGGGCGTTTCATCCTGTTCGTTCAGGTAATGCCAGaa                             <  1:378287/71‑1 (MQ=255)
                                         cTGATAAAATCCCAGGTCATGTCATCGGCAAATCCGGGCGTTTCATCCTGTTCGTTCAGGTAATGCCAGaa                             <  1:2638948/71‑1 (MQ=255)
                                           gATAAAATCCCAGGTCATGTCATCGGCAAATCCGGGCGTTTCATCCTGTTCGTTCAGGTa                                      <  1:609926/60‑1 (MQ=255)
                                            aTAAAATCCCAGGTCATGTCATCGGCAAATCCGGGCGtt                                                          >  1:3215399/1‑39 (MQ=255)
                                            aTAAAATCCCAGGTCATGTCATCGGCAAATCCGGGCGTTTCATCCTgttc                                               <  1:2356740/50‑1 (MQ=255)
                                            aTAAAATCCCAGGTCATGTCATCGGCAAATCCGGGCGTTTCATCCTgttc                                               <  1:2274384/50‑1 (MQ=255)
                                            aTAAAATCCCAGGTCATGTCATCGGCAAATCCGGGCGTTTCATCCTGTTCGTTCAGGTAATGCCAGAAc                            <  1:1656927/69‑1 (MQ=255)
                                            aTAAAATCCCAGGTCATGTCATCGGCAAATCCGGGCGTTTCATCCTGTTCGTTCAGGTAATGCCAGAAc                            <  1:3587071/69‑1 (MQ=255)
                                             tAAAATCCCAGGTCATGTCATCGGTAAATCCGGGCGt                                                           <  1:1858338/37‑1 (MQ=255)
                                               aaaTCCCAGGTCATGTCATCGGTAAATCCGGGCGttt                                                         >  1:1970418/1‑37 (MQ=255)
                                               aaaTCCCAGGTCATGTCATCGGCAAATCCGGGCGTTTCATCCTGTTCGTTCAg                                         >  1:1706389/1‑53 (MQ=255)
                                               aaaTCCCAGGTCATGTCATCGGCAAATCCGGGCGTTTCATCCTGTTCGTTCAGGTAATGCCAGAACTggg                        <  1:194239/70‑1 (MQ=255)
                                               aaaTCCCAGGTCATGTCATCGGCAAATCCGGGCGTTTCATCCTGTTCGTTCAGGTAATGCCAGAACTggg                        <  1:3583405/70‑1 (MQ=255)
                                                aaTCCCAGGTCATGTCATCGGCAAATCCGGGCGttt                                                         >  1:913086/1‑36 (MQ=255)
                                                 aTCCCAGGTCATGTCATCGGCAAATCCGGGCGTTTCATCCTGTTCGTTCAGGTAATGCCAGAACTGGGaaa                     <  1:3684159/71‑1 (MQ=255)
                                                 aTCCCAGGTCATGTCATCGGCAAATCCGGGCGTTTCATCCTGTTCGTTCAGGTAATGCCAGAACTGGGaa                      >  1:2675730/1‑70 (MQ=255)
                                                  tCCCAGGTCATGTCATCGGCAAATCCGGGCGTTTCATCCTGTTCGTTCAGGTAATGCCAGAACTGGGaaaa                    <  1:2105668/71‑1 (MQ=255)
                                                  tCCCAGGTCATGTCATCGGCAAATCCGGGCGTTTCATCCTGTTCGTTCAGGTAAAGCCAGAACTGGGaaaa                    <  1:1255528/71‑1 (MQ=255)
                                                    ccAGGTCATGTCATCGGTAAATCCGGGCGTTTCATCCTGTTCGTTCAGGTAATGCCAGAACTGGGAAAAAt                  <  1:3532042/71‑1 (MQ=255)
                                                    ccAGGTCATGTCATCGGCAAATCCGGGCGTTTCATCCTGTTCGTTCAGGTAATGCCAGAACTGGGAAAAAt                  <  1:2350700/71‑1 (MQ=255)
                                                     cAGGTCATGTCATCGGCAAATCCGGGCGTTTCATCCTGTTCGTTCAGGTAATGCCAGAACTGGGAAAAATg                 >  1:2885701/1‑71 (MQ=255)
                                                       ggTCATGTCATCGGCAAATCCGGGCGTTTCATCCTGTTCGTTCAGGTAATGCCAGAACTGGGAAAAATGCg               <  1:3842758/71‑1 (MQ=255)
                                                       ggTCATGTCATCGGCAAATCCGGGCGTTTCATCCTGTTCGTTCAGGTAATGCCAGAACTGGGAAAAATGCg               <  1:1447015/71‑1 (MQ=255)
                                                       ggTCATGTCATCGGCAAATCCGGGCGTTTCATCCTGTTCGTTCAGGTAATGCCAGAACTGGGAAAAATGCg               >  1:532817/1‑71 (MQ=255)
                                                        gtcatgtcatCGGCAAATCCGGGCGTTTCATCCTGTTCGTTCAGGTAATGCCAGAACTGGGAAAAATGCGt              <  1:349033/71‑1 (MQ=255)
                                                         tcatgtcatCGGCAAATCCGGGCGTTTCATCCTGTTCGTTCAGGTAATGCCAGAACTGGGAAAAATGCGt              <  1:2706963/70‑1 (MQ=255)
                                                           atgtcatCGGCAAATCCGGGCGTTTCATCCTgttcgt                                             <  1:485416/37‑1 (MQ=255)
                                                           atgtcatCGGCAAATCCGGGCGTTTCATCCTgttcgt                                             <  1:875547/37‑1 (MQ=255)
                                                           atgtcatCGGCAAATCCGGGCGTTTCATCCTGTTCGTTCAGGTAATGCCAGAACTGGGAAAAATGCGTGa            <  1:3035662/70‑1 (MQ=255)
                                                            tgtcatCGGCAAATCCGGGCGTTTCATCCTGTTCGTTCAGGTAATGCCAGAACTGGGAAAAATg                 <  1:511820/64‑1 (MQ=255)
                                                            tgtcatCGGCAAATCCGGGCGTTTCATCCTGTTCGTTCAGGTAATGCCAGAACTGGGAAAAATGCGTGata          <  1:3054990/71‑1 (MQ=255)
                                                             gtcatCGGTAAATCCGGGCGTTTCATCCTGTTCGTTCAGGTAAt                                    >  1:2564093/1‑44 (MQ=255)
                                                             gtcatCGGTAAATCCGGGCGTTTCATCCTGTTCGTTCAGGTAAt                                    >  1:1839993/1‑44 (MQ=255)
                                                              tcatCGGCAAATCCGGGCGTTTCATCCTGTTCGTTCAGGTAATGCCAGAACTGGGAAAAATGCGTGatat         <  1:2658681/70‑1 (MQ=255)
                                                              tcatCGGCAAATCCGGGCGTTTCATCCTGTTCGTTCAGGTAATGCCAGAACTGGGAAAAATGCGTGATATc        <  1:2183907/71‑1 (MQ=255)
                                                              tcatCGGCAAATCCGGGCGTTTCATCCTGTTCGTTCAGGTAATGCCAGAACTGGGAAAAATGCGTGATATc        <  1:2103992/71‑1 (MQ=255)
                                                                atCGGCAAATCCGGGCGTTTCATCCTGTTCGTTCAGGTAATGCCAGAACTGGGAAAAATGCGTGATATCg       <  1:1182503/70‑1 (MQ=255)
                                                                 tCGGCAAATCCGGGCGTTTCATCCTGTTCGTTCAGGt                                       <  1:608384/37‑1 (MQ=255)
                                                                 tCGGCAAATCCGGGCGTTTCATCCTGTTCGTTCAGGTAATGCCAGAACTGGGaaa                     <  1:621017/55‑1 (MQ=255)
                                                                 tCGGCAAATCCGGGCGTTTCATCCTGTTCGTTCAGGTAATGCCAGAACTGGGaaa                     <  1:3477840/55‑1 (MQ=255)
                                                                   ggCAAATCCGGGCGTTTCATCCTGTTCGTTCAGGTAATGCCAGAACTGGGAAAAATGCGTGATATCGGCg    >  1:365244/1‑70 (MQ=255)
                                                                   ggCAAATCCGGGCGTTTCATCCTGTTCGTTCAGGTAATGCCAGAACTGGGAAAAATGCGTGATATCGGCGt   <  1:627211/71‑1 (MQ=255)
                                                                    gCAAATCCGGGCGTTTCATCCTGTTCGTTCAGGTAATGCCAGAACTGGGAAAAATGCGTGATATCGGCGTa  <  1:817889/71‑1 (MQ=255)
                                                                     |                                                                     
CGTCATAGAGCATCGCATTGCGGGTTATGCTATTGACGTTGCTGATAAAATCCCAGGTCATGTCATCGGCAAATCCGGGCGTTTCATCCTGTTCGTTCAGGTAATGCCAGAACTGGGAAAAATGCGTGATATCGGCGTA  >  minE/1963631‑1963769

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: