Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 236,107 A→C 43.4% 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,1070AC43.4% 7.3 / 150.7 90E14A (GAA→GCA) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base A (19/32);  new base C (21/18);  total (40/50)
Fisher's exact test for biased strand distribution p-value = 1.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.11e-01

AAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTGAA  >  W3110S.gb/236038‑236177
                                                                     |                                                                      
aaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                                        <  1:1810031/70‑1 (MQ=255)
aaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                                        <  1:1533224/70‑1 (MQ=255)
 caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                                      >  1:2070409‑M2/6‑71 (MQ=255)
 aaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGaaa                                                                      <  1:2543879/71‑1 (MQ=255)
  aCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                                        <  1:2025366/68‑1 (MQ=255)
      cAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAaccacc                                                                 >  1:2310338/1‑71 (MQ=255)
       aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                                                 <  1:558937/70‑1 (MQ=255)
       aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCg                                                                >  1:67679/1‑71 (MQ=255)
        aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCg                                                                <  1:2046015/70‑1 (MQ=255)
              aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATg                                                           <  1:2299628/69‑1 (MQ=255)
              aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                          <  1:258238/70‑1 (MQ=255)
              aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTAt                                                            <  1:2269842/68‑1 (MQ=255)
              aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGaa                                                         >  1:1021938/1‑71 (MQ=255)
              aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGa                                                          <  1:65551/70‑1 (MQ=255)
              aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGa                                                          <  1:1572169/70‑1 (MQ=255)
                  aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                        >  1:2494316/1‑68 (MQ=255)
                  aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                                        >  1:1329915/1‑68 (MQ=255)
                   aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGa                                                    <  1:972228/71‑1 (MQ=255)
                    aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAt                                                   <  1:2050382/71‑1 (MQ=255)
                    aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAt                                                   <  1:1407353/71‑1 (MQ=255)
                        ccGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCCGAtt                                                  <  1:481808/68‑1 (MQ=255)
                           cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                                 <  1:632408/66‑1 (MQ=255)
                            tATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                            >  1:3071852/1‑70 (MQ=255)
                             aTGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCa                                          >  1:1202460/1‑71 (MQ=255)
                              tGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCAc                                         <  1:2075131/71‑1 (MQ=255)
                                aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                                         >  1:372468/1‑69 (MQ=255)
                                aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                                       <  1:1532090/71‑1 (MQ=255)
                                aGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACt                                        <  1:136456/70‑1 (MQ=255)
                                aGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACt                                        <  1:2185141/70‑1 (MQ=255)
                                  cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACt                                        <  1:1224876/68‑1 (MQ=255)
                                  cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                     <  1:2685888/71‑1 (MQ=255)
                                   aCTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGt                                                              >  1:1311583/1‑45 (MQ=255)
                                   aCTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGt                                               <  1:615569/60‑1 (MQ=255)
                                   aCTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                                    <  1:304135/71‑1 (MQ=255)
                                    cTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCa                                                      <  1:3055312/52‑1 (MQ=255)
                                    cTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                     >  1:1333040/1‑69 (MQ=255)
                                    cTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                                   <  1:1841496/71‑1 (MQ=255)
                                     tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                               <  1:2331010/42‑1 (MQ=255)
                                      gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                                 >  1:2648597/1‑71 (MQ=255)
                                       cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg                                  <  1:174383/69‑1 (MQ=255)
                                       cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                                >  1:968896/1‑71 (MQ=255)
                                       cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                                >  1:3013051/1‑71 (MQ=255)
                                       cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                                >  1:1719885/1‑71 (MQ=255)
                                       cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                                >  1:2018684/1‑71 (MQ=255)
                                         aTTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAg                                                     <  1:2906011/48‑1 (MQ=255)
                                          ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                             <  1:1746334/71‑1 (MQ=255)
                                          ttACACGCCAGATCGTTCTCGATACCGAAACTACCGGt                                                              <  1:3650/38‑1 (MQ=255)
                                          ttACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                             <  1:1540281/71‑1 (MQ=255)
                                           tACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGa                                                    >  1:2120146/1‑47 (MQ=39)
                                           tACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                                >  1:1711725/1‑67 (MQ=255)
                                           tACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                             >  1:741591/1‑70 (MQ=255)
                                           tACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAtt                                                  <  1:3067891/49‑1 (MQ=255)
                                            acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                                 >  1:2765583/1‑49 (MQ=39)
                                            acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGaa                                                         >  1:1712307/1‑41 (MQ=255)
                                            acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                             >  1:2342085/1‑69 (MQ=255)
                                             cacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACaa                           <  1:2765454/70‑1 (MQ=255)
                                                gCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                          <  1:1577191/36‑1 (MQ=37)
                                                 ccAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa                      <  1:220652/71‑1 (MQ=255)
                                                   aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCa                                          >  1:1473374/1‑49 (MQ=39)
                                                    gATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                    >  1:134097/1‑70 (MQ=255)
                                                     aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATTAAGGCCACAAGATCAtt                    >  1:1084783/1‑69 (MQ=255)
                                                     aTCGTTCTCGATACCGAAACCACCGGTATGAACCAg                                                     <  1:2901390/36‑1 (MQ=255)
                                                     aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGAtt                                                  <  1:1712743/39‑1 (MQ=255)
                                                     aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGAtt                                                  <  1:933621/39‑1 (MQ=255)
                                                     aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTg                                                 <  1:45467/40‑1 (MQ=255)
                                                     aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                    >  1:843164/1‑69 (MQ=255)
                                                      tCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgag                 <  1:1309152/71‑1 (MQ=255)
                                                       cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTgg                                                >  1:2888790/1‑39 (MQ=255)
                                                         ttCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                     <  1:1288669/48‑1 (MQ=255)
                                                         ttCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                     <  1:1720805/48‑1 (MQ=255)
                                                         ttCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga                  <  1:833213/67‑1 (MQ=255)
                                                         ttCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAtt              >  1:1198447/1‑71 (MQ=255)
                                                           ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTgg            >  1:310778/1‑71 (MQ=255)
                                                           ctcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                                 >  1:2754342/1‑50 (MQ=255)
                                                             cgctACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTg          >  1:1002793/4‑71 (MQ=255)
                                                             cGATACCGAAACCACCGGTATGAACCAGATTGGTgcg                                            <  1:762956/37‑1 (MQ=255)
                                                              gctACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGt           <  1:1717316/67‑1 (MQ=255)
                                                              gATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTgg            <  1:3069685/68‑1 (MQ=255)
                                                                tACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgaga                <  1:931410/62‑1 (MQ=255)
                                                                tACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCg       >  1:2140872/1‑71 (MQ=255)
                                                                tACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                                    >  1:1388671/1‑42 (MQ=255)
                                                                tACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCg       <  1:1825236/71‑1 (MQ=255)
                                                                tACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCg       <  1:2440185/71‑1 (MQ=255)
                                                                 aCCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                                 <  1:1419097/44‑1 (MQ=255)
                                                                 aCCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGt      >  1:2083594/1‑71 (MQ=255)
                                                                 aCCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGt      >  1:675265/1‑71 (MQ=255)
                                                                 aCCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGt      >  1:902696/1‑71 (MQ=255)
                                                                 aCCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGt      >  1:258796/1‑71 (MQ=255)
                                                                  ccGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTg                   <  1:1619773/57‑1 (MQ=255)
                                                                   cGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGttc    <  1:2345517/71‑2 (MQ=255)
                                                                    gAAACCACCGGTATGAACCAGATTGGTGCGCACTa                                       >  1:612045/1‑35 (MQ=255)
                                                                    gAAACCACCGGTATGAACCAGATTGGTGCGCACTa                                       >  1:453327/1‑35 (MQ=255)
                                                                    gAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                     <  1:3083630/37‑1 (MQ=255)
                                                                    gAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACaa                           >  1:2671927/1‑47 (MQ=255)
                                                                    gAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAt                     >  1:452913/1‑53 (MQ=255)
                                                                     aaaCCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTg          >  1:1228977/1‑63 (MQ=255)
                                                                     aaaCCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTGaa  <  1:2098420/71‑1 (MQ=255)
                                                                     |                                                                      
AAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTGAA  >  W3110S.gb/236038‑236177

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