Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 214,925 A→C 45.1% 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,9250AC45.1% 55.3 / 315.3 182D12A (GAT→GCT) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base A (62/38);  new base C (39/43);  total (101/81)
Fisher's exact test for biased strand distribution p-value = 5.36e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 188 total aligned reads displayed.

AACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGC  >  minE/214855‑214994
                                                                      |                                                                     
aaCTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGa                                                                       <  1:2080065/71‑1 (MQ=255)
     tAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCg                                                                  <  1:4087201/71‑1 (MQ=255)
        aaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGaaa                                                               >  1:3812500/1‑71 (MQ=255)
        aaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGaaa                                                               >  1:3424792/1‑71 (MQ=255)
          cGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAac                                                              >  1:3377298/1‑70 (MQ=255)
          cGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAac                                                              >  1:1723462/1‑70 (MQ=255)
          cGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAac                                                              >  1:115467/1‑70 (MQ=255)
          cGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAac                                                              >  1:922039/1‑70 (MQ=255)
             cAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAaccacc                                                          <  1:2153528/71‑1 (MQ=255)
              aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCg                                                         <  1:2492574/71‑1 (MQ=255)
                gTCTGACATAAATGACCGCTATGAGCACTGGAATTACACGCCAGATCGTTCTCGATACCGAAACCACCgg                                                        >  1:3221532/1‑70 (MQ=255)
                   tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCg                                                                  <  1:3079984/57‑1 (MQ=255)
                   tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                     >  1:3789386/1‑70 (MQ=255)
                     aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                     <  1:3631898/68‑1 (MQ=255)
                     aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATg                                                    <  1:3705916/69‑1 (MQ=255)
                       aTAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATa                                                                     >  1:3781495/1‑50 (MQ=255)
                        tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCgca                                                                <  1:3513928/54‑3 (MQ=255)
                        tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                               <  1:2128412/71‑1 (MQ=255)
                         aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGaa                                                  >  1:2360856/1‑67 (MQ=255)
                         aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGaa                                                  >  1:3510070/1‑67 (MQ=255)
                         aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                                 >  1:2840969/1‑68 (MQ=255)
                          aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAac                                                              >  1:2545592/1‑54 (MQ=255)
                          aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTa                                                      >  1:2966818/1‑62 (MQ=255)
                           aTGACCTCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAt                                            <  1:1009127/71‑1 (MQ=255)
                           aTGACCGCTATGAGCACTGCAATTTCACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                        <  1:119135/59‑1 (MQ=255)
                           aTGACCGCTATGAGCACTGCAATTACACTCCAGATCGTTCTCGCTACCGCAACCACCgg                                                        <  1:4045414/59‑1 (MQ=255)
                           aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAt                                            >  1:31411/1‑71 (MQ=255)
                           aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAt                                            >  1:1567237/1‑71 (MQ=255)
                           aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAt                                            <  1:2880265/71‑1 (MQ=255)
                           aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAt                                            <  1:2915213/71‑1 (MQ=255)
                           aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAt                                            >  1:2449916/1‑71 (MQ=255)
                           aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAt                                            <  1:2414989/71‑1 (MQ=255)
                              aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                          <  1:112392/70‑1 (MQ=255)
                                cGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTg                                       <  1:2165184/71‑1 (MQ=255)
                                cGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGa                                             >  1:1371130/1‑65 (MQ=255)
                                 gCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTgc                                      >  1:1902022/1‑71 (MQ=255)
                                  cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                     >  1:3171573/1‑71 (MQ=255)
                                  cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTgcg                                     >  1:735084/1‑71 (MQ=255)
                                   tATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                     >  1:2030707/1‑70 (MQ=255)
                                    aTGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCa                                   >  1:3860506/1‑71 (MQ=255)
                                    aTGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCa                                   >  1:669868/1‑71 (MQ=255)
                                    aTGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCa                                   >  1:3269394/1‑71 (MQ=255)
                                        gCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTAt                               <  1:3869370/71‑1 (MQ=255)
                                        gCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTAt                               <  1:4047990/71‑1 (MQ=255)
                                         cACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                              <  1:2241099/71‑1 (MQ=255)
                                         cACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGAATGAAACACATTGGTGCGCACTAt                               >  1:1685542/1‑70 (MQ=255)
                                          aCTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                             <  1:953597/71‑1 (MQ=255)
                                          aCTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                             <  1:948838/71‑1 (MQ=255)
                                          aCTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                              >  1:3244556/1‑70 (MQ=255)
                                          aCTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                             <  1:3736525/71‑1 (MQ=255)
                                           cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                     <  1:1442894/62‑1 (MQ=255)
                                           cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                              >  1:1851653/1‑69 (MQ=255)
                                           cTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAg                                              <  1:2011097/53‑1 (MQ=255)
                                           cTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                              >  1:95507/1‑69 (MQ=255)
                                           cTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCAACGGTATGAACCAGATTGGTg                                       >  1:1741805/1‑60 (MQ=255)
                                            tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg                           >  1:852748/1‑71 (MQ=255)
                                            tGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGa                                                   <  1:2770827/47‑1 (MQ=255)
                                            tGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg                           >  1:749496/1‑71 (MQ=255)
                                            tGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg                           >  1:685006/1‑71 (MQ=255)
                                              cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                 >  1:4119954/1‑47 (MQ=39)
                                              cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                         >  1:1975389/1‑71 (MQ=255)
                                              cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCCCACTATGAAgg                          <  1:2312947/70‑1 (MQ=255)
                                               aaTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                         >  1:162587/1‑70 (MQ=255)
                                                 ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                      <  1:1786871/71‑1 (MQ=255)
                                                   acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                           >  1:3160677/1‑48 (MQ=39)
                                                   acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                            >  1:2557188/1‑63 (MQ=255)
                                                   acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                      >  1:2654405/1‑69 (MQ=255)
                                                   acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCca                       >  1:1130001/1‑68 (MQ=255)
                                                   acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                                 <  1:2276653/42‑1 (MQ=255)
                                                   acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCAc                                  >  1:581693/1‑57 (MQ=255)
                                                   acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCAc                                  >  1:1428092/1‑57 (MQ=255)
                                                   acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcaca                     <  1:1478826/70‑1 (MQ=255)
                                                    cacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                      >  1:2209623/1‑68 (MQ=255)
                                                    cacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAg                   <  1:4022486/71‑1 (MQ=255)
                                                     acGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAg                   <  1:470463/70‑1 (MQ=255)
                                                     acGCCAGAGCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAg                   <  1:3807945/70‑1 (MQ=255)
                                                      cGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt                 <  1:3552063/71‑1 (MQ=255)
                                                       gCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATc                <  1:1652303/71‑1 (MQ=255)
                                                          aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg                           >  1:177462/1‑57 (MQ=255)
                                                          aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt             >  1:1083241/1‑71 (MQ=255)
                                                          aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAt              >  1:1837360/1‑70 (MQ=255)
                                                            aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                           <  1:2409863/39‑1 (MQ=37)
                                                            aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt             >  1:3480568/1‑69 (MQ=255)
                                                            aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga           <  1:3892604/71‑1 (MQ=255)
                                                            aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga           <  1:1754827/71‑1 (MQ=255)
                                                             tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga           <  1:1173453/70‑1 (MQ=255)
                                                             tCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                            >  1:1151316/1‑53 (MQ=255)
                                                              cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt             >  1:3504368/1‑67 (MQ=255)
                                                              cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGt                                        >  1:1577855/1‑40 (MQ=255)
                                                               gTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt        <  1:3886545/71‑1 (MQ=255)
                                                                ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGt                                        <  1:3471051/38‑1 (MQ=37)
                                                                ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGt                                        <  1:529121/38‑1 (MQ=37)
                                                                ttCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTAt                               <  1:393874/47‑1 (MQ=255)
                                                                 tctcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                              <  1:1917225/47‑1 (MQ=255)
                                                                 tctcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTg      >  1:2756164/1‑71 (MQ=255)
                                                                  ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTg      <  1:649934/70‑1 (MQ=255)
                                                                   tcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTg      <  1:333172/69‑1 (MQ=255)
                                                                     gATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTg   >  1:1003198/1‑70 (MQ=255)
                                                                     gATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGc  >  1:2339531/1‑71 (MQ=255)
                                                                      aTACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                      <  1:1228779/50‑1 (MQ=255)
                                                                      |                                                                     
AACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGC  >  minE/214855‑214994

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: