Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 214,925 A→C 91.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,9250AC91.5% 2690.1 / 137.9 1116D12A (GAT→GCT) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base A (65/31);  new base C (586/434);  total (651/465)
Fisher's exact test for biased strand distribution p-value = 5.21e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

Only 100 of 1116 total aligned reads displayed.

GACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATT  >  minE/214866‑214989
                                                           |                                                                
ctcAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAac                                                         <  1:2474477/67‑1 (MQ=255)
  cAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAacc                                                        >  1:486814/1‑68 (MQ=255)
    aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                                     >  1:435472/1‑69 (MQ=255)
     gTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccac                                                      <  1:2316234/67‑1 (MQ=255)
     gTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccac                                                      <  1:1223072/67‑1 (MQ=255)
      tCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                   <  1:2951227/69‑1 (MQ=255)
      tCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                   <  1:2111843/69‑1 (MQ=255)
      tCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                   <  1:590287/69‑1 (MQ=255)
      tCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                   <  1:2365224/69‑1 (MQ=255)
      tCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                   <  1:2457284/69‑1 (MQ=255)
       cTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGt                                                  >  1:2646931/1‑69 (MQ=255)
       cTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGt                                                  >  1:1198372/1‑69 (MQ=255)
       cTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGt                                                  >  1:279709/1‑69 (MQ=255)
       cTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGt                                                  >  1:3420733/1‑69 (MQ=255)
            aTAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                              <  1:3122682/68‑1 (MQ=255)
            aTAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                              >  1:1776273/1‑68 (MQ=255)
            aTAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATg                                               <  1:841480/67‑1 (MQ=255)
             tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                              >  1:3067592/1‑67 (MQ=255)
             tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                              >  1:1850284/1‑67 (MQ=255)
               aaTGACCGCTATGAGCACTTCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                          <  1:1112030/69‑1 (MQ=255)
               aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGt                                                  <  1:3044481/61‑1 (MQ=255)
               aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATg                                               >  1:1692021/1‑64 (MQ=255)
               aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                              >  1:1000043/1‑65 (MQ=255)
               aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                              >  1:2778132/1‑65 (MQ=255)
               aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                              >  1:2915491/1‑65 (MQ=255)
               aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                              >  1:981542/1‑65 (MQ=255)
               aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                              >  1:791176/1‑65 (MQ=255)
               aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                              >  1:3329296/1‑65 (MQ=255)
               aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                          <  1:1528699/69‑1 (MQ=255)
                 tGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGa                                        >  1:2040950/1‑69 (MQ=255)
                     cGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTgg                                    >  1:875360/1‑69 (MQ=255)
                         aTGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAg                                         >  1:397835/1‑60 (MQ=255)
                            aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATg                                               >  1:2076331/1‑51 (MQ=255)
                            aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAg                                         >  1:16299/1‑57 (MQ=255)
                            aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                             >  1:1654945/1‑69 (MQ=255)
                            aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                             <  1:1140281/69‑1 (MQ=255)
                              cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACt                            <  1:1326800/68‑1 (MQ=255)
                               aCTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                   >  1:2873217/1‑44 (MQ=255)
                                cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                           >  1:1598501/1‑67 (MQ=255)
                                  gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                   >  1:2986416/1‑41 (MQ=255)
                                  gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                   >  1:945325/1‑41 (MQ=255)
                                  gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                   >  1:1275406/1‑41 (MQ=255)
                                  gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                       >  1:1375044/1‑69 (MQ=255)
                                  gCAATTACACGCCAGATAGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                       >  1:965059/1‑69 (MQ=255)
                                   cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg                      >  1:1945086/1‑69 (MQ=255)
                                   cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg                      >  1:3383907/1‑69 (MQ=255)
                                   cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg                      >  1:216033/1‑69 (MQ=255)
                                       tACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGcc                   >  1:2739265/1‑68 (MQ=255)
                                       tACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGcc                   >  1:2254765/1‑68 (MQ=255)
                                       tACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGcc                   >  1:516133/1‑68 (MQ=255)
                                       tACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGcc                   >  1:2501436/1‑68 (MQ=255)
                                        acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                      >  1:648477/1‑48 (MQ=39)
                                        acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                     <  1:36110/49‑1 (MQ=39)
                                        acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                <  1:2094289/54‑1 (MQ=255)
                                        acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGcc                   >  1:3299890/1‑67 (MQ=255)
                                        acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGcc                   >  1:325964/1‑67 (MQ=255)
                                        acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGcc                   >  1:271036/1‑67 (MQ=255)
                                        acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGcc                   >  1:2187078/1‑67 (MQ=255)
                                        acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGcc                   >  1:1910279/1‑67 (MQ=255)
                                        acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTg                                     <  1:2416848/49‑1 (MQ=255)
                                        acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTgcgc                               <  1:106117/55‑1 (MQ=255)
                                           cGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                         <  1:1563754/58‑1 (MQ=255)
                                           cGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAg              <  1:1078952/69‑1 (MQ=255)
                                           cGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAg              <  1:677040/69‑1 (MQ=255)
                                              cAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt            <  1:1342258/68‑1 (MQ=255)
                                               aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa          <  1:2613924/69‑1 (MQ=255)
                                               aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa          <  1:309158/69‑1 (MQ=255)
                                                gATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa          <  1:742396/68‑1 (MQ=255)
                                                gATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTg                                     <  1:2847052/41‑1 (MQ=255)
                                                  tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCTCACTATGAAGGCCACAAGATCa          >  1:1253999/1‑66 (MQ=255)
                                                  tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg                      >  1:380190/1‑54 (MQ=255)
                                                  tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATc           >  1:1482720/1‑65 (MQ=255)
                                                  tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa          >  1:3363363/1‑66 (MQ=255)
                                                  tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa          >  1:2571520/1‑66 (MQ=255)
                                                  tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa          >  1:3141543/1‑66 (MQ=255)
                                                  tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa          >  1:3019072/1‑66 (MQ=255)
                                                  tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa          >  1:2824590/1‑66 (MQ=255)
                                                  tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa          >  1:2681200/1‑66 (MQ=255)
                                                  tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa          >  1:2400161/1‑66 (MQ=255)
                                                  tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa          >  1:2335834/1‑66 (MQ=255)
                                                  tCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa          >  1:553490/1‑66 (MQ=255)
                                                   cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                         <  1:2529809/50‑1 (MQ=255)
                                                   cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa          >  1:1037444/1‑65 (MQ=255)
                                                   cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa          >  1:1409651/1‑65 (MQ=255)
                                                   cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa          >  1:1807150/1‑65 (MQ=255)
                                                   cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa          >  1:631042/1‑65 (MQ=255)
                                                   cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt        <  1:2293387/67‑1 (MQ=255)
                                                   cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga      <  1:3180163/69‑1 (MQ=255)
                                                   cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga      <  1:697419/69‑1 (MQ=255)
                                                   cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga      <  1:1872629/69‑1 (MQ=255)
                                                   cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga      <  1:2591469/69‑1 (MQ=255)
                                                   cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga      <  1:3209821/69‑1 (MQ=255)
                                                    gTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgag     <  1:1828704/69‑1 (MQ=255)
                                                    gTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgag     <  1:1460034/69‑1 (MQ=255)
                                                    gTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgag     <  1:1174427/69‑1 (MQ=255)
                                                      tctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt   <  1:2224215/69‑1 (MQ=255)
                                                      tctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt   <  1:457621/69‑1 (MQ=255)
                                                      tctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt   >  1:144551/1‑69 (MQ=255)
                                                      tctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt   <  1:911149/69‑1 (MQ=255)
                                                        tcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa          <  1:1736716/60‑1 (MQ=255)
                                                           |                                                                
GACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATT  >  minE/214866‑214989

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