Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 236,107 A→C 100% E14A (GAA→GCA)  dnaQ → DNA polymerase III epsilon subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb236,1070AC95.8% 3203.8 / 51.7 887E14A (GAA→GCA) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base A (12/26);  new base C (468/380);  total (481/406)
Fisher's exact test for biased strand distribution p-value = 4.60e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 905 total aligned reads displayed.

AACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTGA  >  W3110S.gb/236039‑236176
                                                                    |                                                                     
caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                                     >  1:2909474‑M2/6‑71 (MQ=255)
     cAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                                                >  1:870197/1‑71 (MQ=255)
       aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                              <  1:3036703/71‑1 (MQ=255)
       aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                              <  1:1794703/71‑1 (MQ=255)
       aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                              <  1:2667813/71‑1 (MQ=255)
             aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATg                                                          <  1:3125811/69‑1 (MQ=255)
                  aaTGACCGCTATGAGCACTGCAATTCCACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                                        >  1:1562808/1‑66 (MQ=255)
                  aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                       >  1:1001019/1‑67 (MQ=255)
                  aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                       >  1:2601731/1‑67 (MQ=255)
                  aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                       >  1:3432/1‑67 (MQ=255)
                  aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                       >  1:2105790/1‑67 (MQ=255)
                  aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                       >  1:1350737/1‑67 (MQ=255)
                  aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                       >  1:239320/1‑67 (MQ=255)
                  aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                       >  1:1279790/1‑67 (MQ=255)
                  aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                       >  1:1261356/1‑67 (MQ=255)
                  aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                       >  1:247472/1‑67 (MQ=255)
                  aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                       >  1:2511384/1‑67 (MQ=255)
                  aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTA‑‑GCAACCACCGGTATGAACCAGAtt                                                 >  1:1439841/1‑71 (MQ=255)
                  aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                                       >  1:226343/1‑67 (MQ=255)
                  aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                                       >  1:77593/1‑67 (MQ=255)
                    tGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                                 <  1:1772608/71‑1 (MQ=255)
                          cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                              >  1:55294/1‑52 (MQ=255)
                          cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                           <  1:1402081/71‑1 (MQ=255)
                           tATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                                                >  1:1665060/1‑49 (MQ=255)
                             tGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcgc                                          <  1:2201381/69‑1 (MQ=255)
                             tGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcgc                                          <  1:2170576/69‑1 (MQ=255)
                              gAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                                        >  1:2766842/1‑70 (MQ=255)
                               aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                                      >  1:3348480/1‑71 (MQ=255)
                               aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                                      >  1:2792930/1‑71 (MQ=255)
                               aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                                      <  1:233293/71‑1 (MQ=255)
                               aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                                      >  1:1535270/1‑71 (MQ=255)
                                 cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                    <  1:2840484/71‑1 (MQ=255)
                                 cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                    <  1:930622/71‑1 (MQ=255)
                                  aCTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                    <  1:2941118/70‑1 (MQ=255)
                                   cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                                        >  1:1092683/1‑49 (MQ=39)
                                   cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                         <  1:3218619/48‑1 (MQ=39)
                                      cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                               >  1:901551/1‑71 (MQ=255)
                                      cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                               >  1:3009312/1‑71 (MQ=255)
                                      cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                               >  1:1177727/1‑71 (MQ=255)
                                      cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                               >  1:1371179/1‑71 (MQ=255)
                                      cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                               >  1:2225524/1‑71 (MQ=255)
                                      cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                               >  1:40522/1‑71 (MQ=255)
                                      cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                               >  1:1733646/1‑71 (MQ=255)
                                        aTTACACGCCATATGGTTGTCGCTACCGCAACCACCGGTATGa                                                         <  1:1593211/43‑1 (MQ=16)
                                        aTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                         <  1:444679/43‑1 (MQ=38)
                                        aTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                         <  1:700127/43‑1 (MQ=38)
                                         ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                    <  1:2882674/63‑1 (MQ=255)
                                          tACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                           >  1:3167856/1‑55 (MQ=255)
                                          tACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                           >  1:2719566/1‑55 (MQ=255)
                                           acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                                  >  1:3313911/1‑63 (MQ=255)
                                           acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                            >  1:1880381/1‑69 (MQ=255)
                                           acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                            >  1:2429790/1‑69 (MQ=255)
                                           acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                            >  1:2136294/1‑69 (MQ=255)
                                           acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                            >  1:20708/1‑69 (MQ=255)
                                            cacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                            >  1:3052999/1‑68 (MQ=255)
                                            cacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                            >  1:3248633/1‑68 (MQ=255)
                                            cacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACaa                          <  1:73679/70‑1 (MQ=255)
                                            cacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACaa                          <  1:1044631/70‑1 (MQ=255)
                                            cacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACaa                          <  1:2621815/70‑1 (MQ=255)
                                              cGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                                 >  1:1123484/1‑45 (MQ=38)
                                                  aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                   >  1:2710826/1‑71 (MQ=255)
                                                    aTCGTTCTCGCTACCGCAACCACCGGTTTGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga                 <  1:255034/71‑1 (MQ=255)
                                                    aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga                 <  1:652207/71‑1 (MQ=255)
                                                    aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga                 <  1:2564548/71‑1 (MQ=255)
                                                    aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga                 <  1:2374034/71‑1 (MQ=255)
                                                    aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga                 <  1:612141/71‑1 (MQ=255)
                                                    aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga                 <  1:1151150/71‑1 (MQ=255)
                                                    aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga                 <  1:830442/71‑1 (MQ=255)
                                                    aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga                 <  1:799499/71‑1 (MQ=255)
                                                    aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga                 <  1:1991707/71‑1 (MQ=255)
                                                    aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga                 <  1:1694183/71‑1 (MQ=255)
                                                    aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga                 <  1:1847734/71‑1 (MQ=255)
                                                    aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga                 <  1:1906359/71‑1 (MQ=255)
                                                    aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga                 <  1:2690209/71‑1 (MQ=255)
                                                     tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                   >  1:597685/1‑68 (MQ=255)
                                                      cGTTCTCGCTACCGCAACCACCGGTATGAACCAGa                                                   >  1:1963455/1‑35 (MQ=37)
                                                      cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATc                      >  1:476954/1‑64 (MQ=255)
                                                      cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATc                      >  1:678367/1‑64 (MQ=255)
                                                      cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                   >  1:496501/1‑67 (MQ=255)
                                                      cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                   >  1:1200645/1‑67 (MQ=255)
                                                      cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                   >  1:1928265/1‑67 (MQ=255)
                                                      cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                   >  1:2291660/1‑67 (MQ=255)
                                                       gTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgaga               <  1:1317364/70‑1 (MQ=255)
                                                       gTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgaga               <  1:571777/70‑1 (MQ=255)
                                                        ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                   <  1:162013/65‑1 (MQ=255)
                                                        ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                   <  1:522675/65‑1 (MQ=255)
                                                        ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                   <  1:1467583/65‑1 (MQ=255)
                                                        ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga                 <  1:631081/67‑1 (MQ=255)
                                                        ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAtt             <  1:973302/71‑1 (MQ=255)
                                                        ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAtt             <  1:3289030/71‑1 (MQ=255)
                                                        ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAtt             <  1:2646634/71‑1 (MQ=255)
                                                          ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGGACTATGAAGGCCACAAGATCATTGAGATTgg           >  1:372496/1‑71 (MQ=255)
                                                          ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAt                    >  1:1234593/1‑62 (MQ=255)
                                                          ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAt                    >  1:2037132/1‑62 (MQ=255)
                                                          ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga                 >  1:1507945/1‑65 (MQ=255)
                                                              ctACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGcc       <  1:42783/70‑1 (MQ=255)
                                                                   gAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTGa  <  1:3096070/71‑1 (MQ=255)
                                                                   gAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTGa  <  1:2972004/71‑1 (MQ=255)
                                                                   gAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTGa  <  1:1064168/71‑1 (MQ=255)
                                                                    aaaCCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTGa  <  1:3146308/70‑1 (MQ=255)
                                                                    |                                                                     
AACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTGA  >  W3110S.gb/236039‑236176

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: