Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 236,101 A→C 54.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
*W3110S.gb236,1010AC54.1% 163.4 / 157.2 150D12A (GAT→GCT) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base A (49/20);  new base C (45/36);  total (94/56)
Fisher's exact test for biased strand distribution p-value = 6.29e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 154 total aligned reads displayed.

CTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGT  >  W3110S.gb/236033‑236168
                                                                    |                                                                   
cTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGa                                                                     >  1:2559962/1‑69 (MQ=255)
cTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGa                                                                     >  1:736272/1‑69 (MQ=255)
cTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGa                                                                     >  1:3663535/1‑69 (MQ=255)
   tAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATa                                                                   >  1:870383/1‑68 (MQ=255)
   tAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATa                                                                   >  1:1888548/1‑68 (MQ=255)
   tAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATa                                                                   >  1:1810571/1‑68 (MQ=255)
   gtacaactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTAc                                                                  <  1:2591558‑M2/61‑1 (MQ=255)
      caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACcgc                                                               >  1:2545064‑M2/6‑68 (MQ=255)
      caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACcgc                                                               >  1:2704774‑M2/6‑68 (MQ=255)
      caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACcgc                                                               <  1:3344441‑M2/64‑2 (MQ=255)
       aCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGaa                                                              >  1:107593/1‑69 (MQ=255)
       aCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGaa                                                              >  1:1413561/1‑69 (MQ=255)
           cAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAacca                                                          >  1:1981684/1‑69 (MQ=255)
           cAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAacca                                                          >  1:3695156/1‑69 (MQ=255)
            aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAaccac                                                         >  1:3128832/1‑69 (MQ=255)
            aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAaccac                                                         >  1:3298321/1‑69 (MQ=255)
                        aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                 >  1:144536/1‑65 (MQ=255)
                        aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                             <  1:1747720/69‑1 (MQ=255)
                        aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                             <  1:2467139/69‑1 (MQ=255)
                        aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGa                                                 >  1:584992/1‑65 (MQ=255)
                        aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGa                                                 >  1:1242257/1‑65 (MQ=255)
                        aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCa                                             <  1:3344565/69‑1 (MQ=255)
                           gACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAt                                          >  1:2280036/1‑69 (MQ=255)
                             ccGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTg                                        >  1:2962867/1‑69 (MQ=255)
                                cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTg                                     <  1:1051605/69‑1 (MQ=255)
                                 tATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTg                                     >  1:3558630/1‑68 (MQ=255)
                                 tATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTg                                     >  1:1122370/1‑68 (MQ=255)
                                     aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                      <  1:2154716/47‑1 (MQ=255)
                                     aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                      <  1:2302443/47‑1 (MQ=255)
                                     aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCa                                 <  1:409497/68‑1 (MQ=255)
                                     aGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCAc                                >  1:3067210/1‑69 (MQ=255)
                                     aGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCAc                                >  1:2110450/1‑69 (MQ=255)
                                     aGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCAc                                >  1:3071372/1‑69 (MQ=255)
                                      gCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACt                               >  1:3498605/1‑69 (MQ=255)
                                      gCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACt                               >  1:1648426/1‑69 (MQ=255)
                                         cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                            <  1:2160331/69‑1 (MQ=255)
                                         cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                            <  1:3000520/69‑1 (MQ=255)
                                         cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                            <  1:2422836/69‑1 (MQ=255)
                                         cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                            <  1:737278/69‑1 (MQ=255)
                                         cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                            <  1:87742/69‑1 (MQ=255)
                                         cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                            <  1:92129/69‑1 (MQ=255)
                                          tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCTCCGGTATGAACCAGATTGGTGCGCACTATGa                           <  1:2957916/69‑1 (MQ=255)
                                          tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                            <  1:1596566/68‑1 (MQ=255)
                                          tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                           <  1:1586981/69‑1 (MQ=255)
                                           gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATg                                                  >  1:3085739/1‑45 (MQ=38)
                                           gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                 >  1:359687/1‑46 (MQ=39)
                                           gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                 >  1:1932187/1‑46 (MQ=39)
                                           gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcgc                                  <  1:3017307/61‑1 (MQ=255)
                                           gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                          <  1:3078235/69‑1 (MQ=255)
                                           gCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                          >  1:942136/1‑69 (MQ=255)
                                           gCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                          >  1:2701207/1‑69 (MQ=255)
                                            cAATTACACGCCATATCGTTCTCGATACCGAAACCACCGGTATg                                                  <  1:138922/44‑1 (MQ=255)
                                            cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg                         >  1:1428566/1‑69 (MQ=255)
                                            cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg                         >  1:2013725/1‑69 (MQ=255)
                                            cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg                         >  1:318919/1‑69 (MQ=255)
                                            cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg                         >  1:3090507/1‑69 (MQ=255)
                                             aaTTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATg                                                  <  1:2527769/43‑1 (MQ=255)
                                                tACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAtt                                         >  1:2889193/1‑49 (MQ=255)
                                                tACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAt                                          >  1:2365547/1‑48 (MQ=255)
                                                 acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAcc                                              >  1:2745267/1‑43 (MQ=38)
                                                 acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAcc                                              >  1:1833942/1‑43 (MQ=38)
                                                 acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                            >  1:1349508/1‑61 (MQ=255)
                                                 acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                            >  1:3153695/1‑61 (MQ=255)
                                                 acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGcc                      >  1:2387550/1‑67 (MQ=255)
                                                 acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGcc                      >  1:2391838/1‑67 (MQ=255)
                                                 acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGcc                      >  1:2352181/1‑67 (MQ=255)
                                                 acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGcc                      >  1:2496134/1‑67 (MQ=255)
                                                 acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                            >  1:1933431/1‑61 (MQ=255)
                                                 acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                            >  1:3526036/1‑61 (MQ=255)
                                                 acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                            >  1:1120060/1‑61 (MQ=255)
                                                 acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                            >  1:1852498/1‑61 (MQ=255)
                                                 acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                            >  1:2869047/1‑61 (MQ=255)
                                                 acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                            >  1:2773201/1‑61 (MQ=255)
                                                 acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGcc                      >  1:601957/1‑67 (MQ=255)
                                                 acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGcc                      >  1:2169406/1‑67 (MQ=255)
                                                 acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                    <  1:529995/69‑1 (MQ=255)
                                                 acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                    <  1:2013471/69‑1 (MQ=255)
                                                 acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                    <  1:768694/69‑1 (MQ=255)
                                                 acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                    <  1:1911443/69‑1 (MQ=255)
                                                 acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                    <  1:3455320/69‑1 (MQ=255)
                                                 acacGCAAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                            >  1:2086726/1‑61 (MQ=255)
                                                          aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa             >  1:688755/1‑67 (MQ=255)
                                                           tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                            <  1:616464/51‑1 (MQ=255)
                                                           tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa             >  1:435282/1‑66 (MQ=255)
                                                           tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa             >  1:3294571/1‑66 (MQ=255)
                                                            cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                            <  1:3432525/50‑1 (MQ=255)
                                                            cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                            <  1:1240556/50‑1 (MQ=255)
                                                            cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACt                               >  1:1821767/1‑47 (MQ=255)
                                                            cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACt                               >  1:3209875/1‑47 (MQ=255)
                                                            cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                            <  1:829698/50‑1 (MQ=255)
                                                            cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                    <  1:2117898/58‑1 (MQ=255)
                                                            cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                    <  1:760367/58‑1 (MQ=255)
                                                             gTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgag        <  1:2374312/69‑1 (MQ=255)
                                                               tctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt      >  1:105902/1‑69 (MQ=255)
                                                               tctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt      <  1:2167592/69‑1 (MQ=255)
                                                                ctcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt               <  1:925680/59‑1 (MQ=255)
                                                                ctcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAtt     >  1:842401/1‑69 (MQ=255)
                                                                ctcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAtt     >  1:1108024/1‑69 (MQ=255)
                                                                   gATACCGAAACCACCGGTATGAACCAGATTGGTGCGCTCTATGAAgg                        >  1:2057118/1‑47 (MQ=255)
                                                                   gATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGt  >  1:3600526/1‑69 (MQ=255)
                                                                    |                                                                   
CTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGT  >  W3110S.gb/236033‑236168

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: