Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 236,107 A→C 63.1% 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,1070AC63.1% 314.5 / 150.0 181E14A (GAA→GCA) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base A (30/37);  new base C (65/49);  total (95/86)
Fisher's exact test for biased strand distribution p-value = 1.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 196 total aligned reads displayed.

AAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTGAA  >  W3110S.gb/236038‑236177
                                                                     |                                                                      
aaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGaa                                                                       <  1:1539063/71‑1 (MQ=255)
aaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                                        <  1:2288900/70‑1 (MQ=255)
 caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                                      >  1:3922132‑M2/6‑71 (MQ=255)
 caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                                      >  1:2084728‑M2/6‑71 (MQ=255)
 caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                                      >  1:1358630‑M2/6‑71 (MQ=255)
    gACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAacca                                                                   <  1:1239674/71‑1 (MQ=255)
      cAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                                                 >  1:4788626/1‑71 (MQ=255)
       aaGTCTGACATAAATGACCGCTATGAGCACTGAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                               >  1:1043925/1‑71 (MQ=255)
        aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                               <  1:3075521/71‑1 (MQ=255)
        aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCgg                                                               <  1:4145326/71‑1 (MQ=255)
           cTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                            >  1:3025529/1‑71 (MQ=255)
            tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                            >  1:1124790/1‑70 (MQ=255)
            tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTAt                                                            >  1:4081491/1‑70 (MQ=255)
            tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTAt                                                            >  1:4355575/1‑70 (MQ=255)
              aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                            <  1:3686647/68‑1 (MQ=255)
              aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATg                                                           >  1:380147/1‑69 (MQ=255)
              aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                                         <  1:24970/71‑1 (MQ=255)
              aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                          <  1:2446062/70‑1 (MQ=255)
              aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                                        <  1:1246516/56‑1 (MQ=255)
              aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATg                                                           <  1:2321075/69‑1 (MQ=255)
              aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGa                                                          <  1:2736707/70‑1 (MQ=255)
                  aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                        >  1:3542819/1‑68 (MQ=255)
                  aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                        >  1:2549623/1‑68 (MQ=255)
                   aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                        >  1:1308454/1‑67 (MQ=255)
                   aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGaa                                                                       >  1:4617353/1‑52 (MQ=255)
                   aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                                        >  1:1393894/1‑67 (MQ=255)
                    aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAt                                                   <  1:2424717/71‑1 (MQ=255)
                     tGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                                  <  1:2391214/71‑1 (MQ=255)
                     tGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAtt                                                  >  1:4217966/1‑71 (MQ=255)
                      gACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                                 >  1:1899324/1‑71 (MQ=255)
                      gACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACAGGTATGAACCAGATTg                                                 >  1:4754433/1‑71 (MQ=255)
                         cGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAac                                                                     >  1:737641/1‑48 (MQ=255)
                           cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                                      <  1:1851824/61‑1 (MQ=255)
                           cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                            >  1:1033280/1‑71 (MQ=255)
                              tGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                                         <  1:576620/71‑1 (MQ=255)
                              tGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCAc                                         <  1:4174845/71‑1 (MQ=255)
                               gAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGaa                                                         >  1:3843687/1‑54 (MQ=255)
                                aGCACTGCAATTACACGCCAGCTCGTTCTCGCTACCGcaa                                                                      >  1:4293870/1‑40 (MQ=255)
                                aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                                      >  1:1554863/1‑40 (MQ=255)
                                aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                                       >  1:349807/1‑71 (MQ=255)
                                aGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGa                                                          >  1:4033359/1‑52 (MQ=255)
                                     tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                                      <  1:4476240/51‑1 (MQ=255)
                                     tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                                 <  1:2844911/56‑1 (MQ=255)
                                      gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                                   >  1:3580077/1‑69 (MQ=255)
                                      gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg                                  >  1:3758437/1‑70 (MQ=255)
                                      gCAATTACACGCCAGATCGTTCTCGATACCGAAac                                                                     >  1:610244/1‑35 (MQ=255)
                                      gCAATTACACGCCAGATCGTTCTCGATACCGAAac                                                                     >  1:4572064/1‑35 (MQ=255)
                                      gCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg                                  >  1:3757188/1‑70 (MQ=255)
                                       cAATTACACGCCAGATCGTTCTCGCTACCGCAacca                                                                   >  1:513326/1‑36 (MQ=38)
                                       cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                                >  1:4245761/1‑71 (MQ=255)
                                       cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                                >  1:1636560/1‑71 (MQ=255)
                                       cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                                >  1:4687226/1‑71 (MQ=255)
                                       cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                                >  1:1185659/1‑71 (MQ=255)
                                         aTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGa                                                    >  1:3323035/1‑49 (MQ=39)
                                          ttACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCa                                                      <  1:4092709/46‑1 (MQ=255)
                                          ttACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAt                                                   <  1:1754146/49‑1 (MQ=255)
                                            acacGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                               >  1:3227525/1‑35 (MQ=38)
                                            acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                             >  1:1159769/1‑69 (MQ=255)
                                            acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                             >  1:1130054/1‑69 (MQ=255)
                                            acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                             >  1:3174151/1‑69 (MQ=255)
                                            acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                             >  1:2138281/1‑69 (MQ=255)
                                            acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                             >  1:519927/1‑69 (MQ=255)
                                            acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCca                              >  1:2221134/1‑68 (MQ=255)
                                             cacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAg                          <  1:279774/71‑1 (MQ=255)
                                               cGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt                        <  1:2366760/71‑1 (MQ=255)
                                                gCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt                        >  1:1980734/1‑70 (MQ=255)
                                                  cAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                                   >  1:2251411/1‑57 (MQ=255)
                                                   aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                                  <  1:867012/41‑1 (MQ=38)
                                                   aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                    >  1:4534571/1‑71 (MQ=255)
                                                   aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                    <  1:656279/71‑1 (MQ=255)
                                                     aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                                  <  1:1485050/39‑1 (MQ=37)
                                                     aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                     <  1:3379935/52‑1 (MQ=255)
                                                     aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                    <  1:3212186/69‑1 (MQ=255)
                                                     aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                    <  1:3912922/69‑1 (MQ=255)
                                                     aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga                  <  1:4199038/71‑1 (MQ=255)
                                                     aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGAtt                                                  <  1:265106/39‑1 (MQ=255)
                                                     aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTg                                                 <  1:188245/40‑1 (MQ=255)
                                                     aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTg                                                 <  1:4449259/40‑1 (MQ=255)
                                                     aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga                  <  1:2690386/71‑1 (MQ=255)
                                                      tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGa                                                    >  1:2732780/1‑36 (MQ=37)
                                                      tCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                    >  1:3096372/1‑68 (MQ=255)
                                                       cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                    >  1:693036/1‑67 (MQ=255)
                                                        gTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt               >  1:3320199/1‑71 (MQ=255)
                                                         ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                             <  1:3849997/56‑1 (MQ=255)
                                                         ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                             <  1:2207428/56‑1 (MQ=255)
                                                         ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGATGCGCACTATGAAGGCcac                             <  1:2943102/56‑1 (MQ=255)
                                                         ttCTCGATACCGAAACCACCGGTATGAACCAGATTg                                                 >  1:2225386/1‑36 (MQ=255)
                                                         ttCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                    <  1:2341564/65‑1 (MQ=255)
                                                         ttCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga                  <  1:266376/67‑1 (MQ=255)
                                                          tctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTg             >  1:1456965/1‑71 (MQ=255)
                                                            tcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgag                 <  1:3734092/65‑1 (MQ=255)
                                                            tcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACt                                        <  1:9340/42‑1 (MQ=255)
                                                              gATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGc         >  1:3620005/1‑71 (MQ=255)
                                                              gATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGc         >  1:695602/1‑71 (MQ=255)
                                                                 aCCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGt      >  1:3506328/1‑71 (MQ=255)
                                                                 aCCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAg                          <  1:962794/51‑1 (MQ=255)
                                                                  ccGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGtt     <  1:1578319/71‑1 (MQ=255)
                                                                   cGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTg    <  1:4575409/71‑1 (MQ=255)
                                                                    gAAACCACCGGTATGAACCAGATTGGTGCGCACTAt                                      <  1:2999639/36‑1 (MQ=255)
                                                                     aaaCCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTGaa  <  1:419100/71‑1 (MQ=255)
                                                                     |                                                                      
AAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTGAA  >  W3110S.gb/236038‑236177

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: