Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 214,925 A→C 55.5% D12A (GAT→GCT)  dnaQ → DNA polymerase III epsilon subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE214,9250AC55.5% 186.4 / 202.8 182D12A (GAT→GCT) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base A (44/37);  new base C (57/44);  total (101/81)
Fisher's exact test for biased strand distribution p-value = 8.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 185 total aligned reads displayed.

TAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGG  >  minE/214860‑214991
                                                                 |                                                                  
tAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGa                                                                    >  1:344069/1‑66 (MQ=255)
tAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATAc                                                                 >  1:3165669/1‑69 (MQ=255)
tAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCg                                                               <  1:1169068/71‑1 (MQ=255)
  aaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGaa                                                             <  1:651895/71‑1 (MQ=255)
  aaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                              >  1:1111766/1‑70 (MQ=255)
  aaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                              <  1:3487274/70‑1 (MQ=255)
  aaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                              >  1:2284414/1‑70 (MQ=255)
     actcAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAac                                                           >  1:323615/4‑70 (MQ=255)
     actcAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAac                                                           >  1:3754532/4‑70 (MQ=255)
        cAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                                       >  1:3828932/1‑71 (MQ=255)
          aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                     <  1:3568046/71‑1 (MQ=255)
          aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCgg                                                     <  1:1604579/71‑1 (MQ=255)
            tCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                              <  1:2535304/60‑1 (MQ=255)
             cTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                  >  1:2104254/1‑71 (MQ=255)
             cTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                  >  1:2569026/1‑71 (MQ=255)
             cTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                  >  1:2394724/1‑71 (MQ=255)
              tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                  >  1:951067/1‑70 (MQ=255)
              tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                  >  1:926086/1‑70 (MQ=255)
              tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTAt                                                  >  1:212133/1‑70 (MQ=255)
               gACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCgca                                                             <  1:13803/58‑3 (MQ=255)
               gACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCgca                                                             <  1:75863/58‑3 (MQ=255)
               gACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATg                                                 >  1:1463658/1‑70 (MQ=255)
                aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATg                                                 <  1:961907/69‑1 (MQ=255)
                aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                               >  1:2611657/1‑71 (MQ=255)
                aCATAAATGACCGCTATCAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGaa                                               >  1:3060425/1‑71 (MQ=255)
                 cATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                              >  1:2267009/1‑71 (MQ=255)
                   tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCa                                            <  1:629233/71‑1 (MQ=255)
                    aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                              >  1:2366385/1‑68 (MQ=255)
                     aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                              >  1:3764344/1‑67 (MQ=255)
                     aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAg                                           <  1:861723/70‑1 (MQ=255)
                      aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAt                                         >  1:1915870/1‑71 (MQ=255)
                         aCCGCTATGAGCACTCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGt                                     >  1:3410912/1‑71 (MQ=255)
                          ccGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATg                                                 <  1:1570897/59‑1 (MQ=255)
                          ccGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATg                                                 <  1:1091951/59‑1 (MQ=255)
                             cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTgcg                                  >  1:2430405/1‑71 (MQ=255)
                             cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTgcg                                  >  1:476613/1‑71 (MQ=255)
                               aTGAGCACTGCAATTACACGCCAGATCGTTCTCGa                                                                    >  1:3162889/1‑35 (MQ=255)
                               aTGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCa                                >  1:22319/1‑71 (MQ=255)
                               aTGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCa                                >  1:3513056/1‑71 (MQ=255)
                                 gAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                <  1:704247/53‑1 (MQ=255)
                                   gCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTAt                            >  1:3780218/1‑71 (MQ=255)
                                     aCTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                          <  1:3176102/71‑1 (MQ=255)
                                     aCTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                          <  1:2525916/71‑1 (MQ=255)
                                       tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                     <  1:1289852/42‑1 (MQ=255)
                                        gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                       >  1:26844/1‑71 (MQ=255)
                                         cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                       >  1:2062236/1‑70 (MQ=255)
                                         cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                      >  1:2793690/1‑71 (MQ=255)
                                         cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                      >  1:2805829/1‑71 (MQ=255)
                                           aTTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGa                                          >  1:3846715/1‑49 (MQ=255)
                                           aTTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCa                                >  1:176094/1‑59 (MQ=255)
                                            ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATg                                                 <  1:1216136/41‑1 (MQ=38)
                                            ttACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCa                                            <  1:2888039/46‑1 (MQ=255)
                                            ttACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCa                                            <  1:1154206/46‑1 (MQ=255)
                                            ttACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCa                                            <  1:2177354/46‑1 (MQ=255)
                                            ctACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATg                                                 <  1:2644892/40‑1 (MQ=38)
                                             tACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                   >  1:814204/1‑70 (MQ=255)
                                             tACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                   >  1:1680840/1‑70 (MQ=255)
                                             tACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTCGTGCGCACTATGAAGGCcac                   >  1:3072485/1‑70 (MQ=255)
                                              acacGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                     >  1:1664280/1‑35 (MQ=38)
                                              acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGa                                          >  1:1106066/1‑46 (MQ=39)
                                              acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                       >  1:2387973/1‑49 (MQ=39)
                                              acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGa                                          >  1:387369/1‑46 (MQ=255)
                                              acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                   >  1:281241/1‑69 (MQ=255)
                                               cacGCCAGATCGTTCTCGCTACCGCAACCACCGGTa                                                   <  1:2747909/36‑1 (MQ=37)
                                               cacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                   >  1:3392773/1‑68 (MQ=255)
                                               cacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                   >  1:981899/1‑68 (MQ=255)
                                               cacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                   >  1:311043/1‑68 (MQ=255)
                                               cacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                   >  1:2115310/1‑68 (MQ=255)
                                               cacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACaa                 <  1:2711557/70‑1 (MQ=255)
                                                acGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGa               <  1:2948754/71‑1 (MQ=255)
                                                 cGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt              <  1:3084819/71‑1 (MQ=255)
                                                 cGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                   >  1:2708166/1‑66 (MQ=255)
                                                 cGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                   >  1:3016273/1‑66 (MQ=255)
                                                  gCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                           <  1:1984919/57‑1 (MQ=255)
                                                  gCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATc             <  1:3291871/71‑1 (MQ=255)
                                                    cAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                          >  1:2124365/1‑56 (MQ=255)
                                                     aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAt           >  1:2495641/1‑70 (MQ=255)
                                                     aGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          <  1:3777104/71‑1 (MQ=255)
                                                      gATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                               <  1:2979531/49‑1 (MQ=39)
                                                      gATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                               <  1:303488/49‑1 (MQ=39)
                                                      gATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          >  1:1741491/1‑70 (MQ=255)
                                                      gATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTgg                                      <  1:973439/42‑1 (MQ=255)
                                                       aTCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                            <  1:3404282/35‑1 (MQ=37)
                                                       aTCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                            <  1:1327606/35‑1 (MQ=37)
                                                       aTCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                            <  1:1897701/35‑1 (MQ=37)
                                                       aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga        <  1:1720890/71‑1 (MQ=255)
                                                        tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                             >  1:1048496/1‑49 (MQ=39)
                                                        tCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga        <  1:3629231/70‑1 (MQ=255)
                                                         cGTTCTCGCTACCGCAACCACCGGTATGAACCAGa                                          >  1:3338960/1‑35 (MQ=37)
                                                         cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          >  1:100565/1‑67 (MQ=255)
                                                         cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          >  1:2963051/1‑67 (MQ=255)
                                                         cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          >  1:1657426/1‑67 (MQ=255)
                                                          gTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt     >  1:459079/1‑71 (MQ=255)
                                                          gTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt     <  1:1139455/71‑1 (MQ=255)
                                                           ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTAt                            <  1:2047881/47‑1 (MQ=39)
                                                           ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTAt                            <  1:1585866/47‑1 (MQ=39)
                                                           ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTAt                            <  1:3171696/47‑1 (MQ=39)
                                                             ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                          >  1:3040722/1‑47 (MQ=39)
                                                             ctcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          >  1:2056558/1‑63 (MQ=255)
                                                              tcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTgg  >  1:364652/1‑70 (MQ=255)
                                                                 |                                                                  
TAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGG  >  minE/214860‑214991

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: