Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 214,925 A→C 70.7% 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,9250AC70.7% 547.9 / 180.6 259D12A (GAT→GCT) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base A (46/30);  new base C (105/78);  total (151/108)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.46e-01

Only 100 of 261 total aligned reads displayed.

ACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTG  >  minE/214856‑214993
                                                                     |                                                                    
aCTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGAt                                                                     <  1:1490733/71‑1 (MQ=255)
    tAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATAcc                                                                  >  1:2807213/1‑70 (MQ=255)
    tAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCg                                                                 <  1:3415628/71‑1 (MQ=255)
      aaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                                >  1:3022583/1‑70 (MQ=255)
         cGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAacc                                                            <  1:1426611/71‑1 (MQ=255)
           aCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAaccac                                                          <  1:2005728/71‑1 (MQ=255)
            cAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                                         >  1:3380534/1‑71 (MQ=255)
            cAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAacca                                                           >  1:1486282/1‑69 (MQ=255)
             aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccac                                                          <  1:169756/69‑1 (MQ=255)
             aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccac                                                          <  1:432702/69‑1 (MQ=255)
                tCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGa                                                                      >  1:1109400/1‑54 (MQ=255)
                tCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                                >  1:2641339/1‑60 (MQ=255)
                 cTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTAt                                                    >  1:3481670/1‑71 (MQ=255)
                  tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTa                                                     >  1:922177/1‑69 (MQ=255)
                  tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                    >  1:2279659/1‑70 (MQ=255)
                  tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                    >  1:1189199/1‑70 (MQ=255)
                  tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                    >  1:1180179/1‑70 (MQ=255)
                  tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                    >  1:1806175/1‑70 (MQ=255)
                    aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATg                                                   <  1:1327691/69‑1 (MQ=255)
                    aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATg                                                   <  1:1583745/69‑1 (MQ=255)
                    aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATg                                                   <  1:266476/69‑1 (MQ=255)
                    aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                  <  1:2704676/70‑1 (MQ=255)
                    aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                  <  1:1011496/70‑1 (MQ=255)
                       tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAac                                                             <  1:1123343/56‑1 (MQ=255)
                       tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                              <  1:2875862/71‑1 (MQ=255)
                          aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAt                                           <  1:3327086/71‑1 (MQ=255)
                          aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAt                                           >  1:791460/1‑71 (MQ=255)
                           tGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                          <  1:3447435/71‑1 (MQ=255)
                           tGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                          <  1:2832734/71‑1 (MQ=255)
                           tGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                          <  1:2152516/71‑1 (MQ=255)
                           tGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                          <  1:2099808/71‑1 (MQ=255)
                             aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                >  1:3440058/1‑63 (MQ=255)
                             aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTgg                                        <  1:1788026/71‑1 (MQ=255)
                             aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTgg                                        <  1:1592841/71‑1 (MQ=255)
                             aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTgg                                        <  1:747421/71‑1 (MQ=255)
                              ccGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGt                                       >  1:683799/1‑71 (MQ=255)
                              ccGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGt                                       >  1:446497/1‑71 (MQ=255)
                               cGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGt                                       <  1:3351373/70‑1 (MQ=255)
                                  tATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                    <  1:1607121/70‑1 (MQ=255)
                                    tGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCa                                              <  1:1401040/58‑1 (MQ=255)
                                       gCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTAt                              <  1:2574102/71‑1 (MQ=255)
                                       gCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTAt                              <  1:2381620/71‑1 (MQ=255)
                                        cACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTAt                              >  1:1377776/1‑70 (MQ=255)
                                         aCTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                             >  1:770860/1‑70 (MQ=255)
                                         aCTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                             <  1:2326374/70‑1 (MQ=255)
                                          cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                             >  1:3277402/1‑69 (MQ=255)
                                          cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                             >  1:2918964/1‑69 (MQ=255)
                                           tGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg                          >  1:3154400/1‑71 (MQ=255)
                                            gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                         >  1:2238403/1‑71 (MQ=255)
                                             cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGt                                       <  1:1862404/56‑1 (MQ=255)
                                             cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgga                        >  1:2512260/1‑70 (MQ=255)
                                             cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                        >  1:1906146/1‑71 (MQ=255)
                                             cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                        >  1:877114/1‑71 (MQ=255)
                                             cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                        >  1:1655429/1‑71 (MQ=255)
                                             cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACAGGTATGAACCAGATTGGTGCGCACTATGAAGGc                        >  1:481238/1‑71 (MQ=255)
                                             cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGa                                            >  1:2931513/1‑51 (MQ=255)
                                             cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                        >  1:1053076/1‑71 (MQ=255)
                                             cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                        >  1:1095036/1‑71 (MQ=255)
                                             cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                        >  1:234560/1‑71 (MQ=255)
                                               aTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                          <  1:673999/51‑1 (MQ=255)
                                               aTTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                           >  1:633106/1‑66 (MQ=255)
                                                ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                              <  1:535054/46‑1 (MQ=39)
                                                ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGt                                       <  1:3230891/53‑1 (MQ=255)
                                                ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGt                                       <  1:31366/53‑1 (MQ=255)
                                                ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCca                      >  1:2447661/1‑70 (MQ=255)
                                                  acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAcc                                               >  1:1869546/1‑43 (MQ=38)
                                                  acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                              >  1:2139694/1‑44 (MQ=38)
                                                  acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                              >  1:2431377/1‑44 (MQ=38)
                                                  acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                              >  1:1745558/1‑44 (MQ=38)
                                                  acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgc                                     >  1:1358509/1‑53 (MQ=255)
                                                  acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgc                                     >  1:2205662/1‑53 (MQ=255)
                                                  acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                     >  1:598736/1‑69 (MQ=255)
                                                    acGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGa                 <  1:3072209/71‑1 (MQ=255)
                                                    acGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                     >  1:183566/1‑67 (MQ=255)
                                                    acGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                     >  1:1650865/1‑67 (MQ=255)
                                                     cGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt                <  1:252415/71‑1 (MQ=255)
                                                         aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                               >  1:147442/1‑52 (MQ=255)
                                                         aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt            <  1:1994015/71‑1 (MQ=255)
                                                         aGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt            >  1:2071752/1‑71 (MQ=255)
                                                           aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt            >  1:584507/1‑69 (MQ=255)
                                                           aTCGTTCTCGATACCGAAACCACCGGTATGAACCa                                              <  1:2734472/35‑1 (MQ=255)
                                                           aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGAtt                                          <  1:1707438/39‑1 (MQ=255)
                                                           aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTg                                         <  1:1447547/40‑1 (MQ=255)
                                                             cGTTCTCGCTACCGCAACCACCGGTATGAACCAGa                                            >  1:56205/1‑35 (MQ=37)
                                                             cGTTCTCGCTACCGCAACCACCGGTATGAACCAGa                                            >  1:2600735/1‑35 (MQ=37)
                                                              gTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgaga        >  1:2961280/1‑70 (MQ=255)
                                                               ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                             <  1:372567/48‑1 (MQ=39)
                                                               ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga          <  1:3048365/67‑1 (MQ=255)
                                                                tctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt       <  1:1962102/69‑1 (MQ=255)
                                                                 ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt            >  1:3090645/1‑63 (MQ=255)
                                                                 ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt            >  1:885909/1‑63 (MQ=255)
                                                                 ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga          <  1:3206889/65‑1 (MQ=255)
                                                                 ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTgg    >  1:858132/1‑71 (MQ=255)
                                                                 ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTgg    >  1:1397882/1‑71 (MQ=255)
                                                                 ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTgg    >  1:2628477/1‑71 (MQ=255)
                                                                 ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTgg    >  1:3036919/1‑71 (MQ=255)
                                                                 ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTg     >  1:2559035/1‑70 (MQ=255)
                                                                 ctcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTg     <  1:3429697/70‑1 (MQ=255)
                                                                  tcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGt   >  1:1733912/1‑71 (MQ=255)
                                                                   cGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTg  >  1:1714920/1‑71 (MQ=255)
                                                                     |                                                                    
ACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTG  >  minE/214856‑214993

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: