Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 236,107 A→C 54.3% 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,1070AC54.3% 163.6 / 187.5 166E14A (GAA→GCA) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base A (36/40);  new base C (47/43);  total (83/83)
Fisher's exact test for biased strand distribution p-value = 6.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

Only 100 of 180 total aligned reads displayed.

AAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTGAA  >  W3110S.gb/236038‑236177
                                                                     |                                                                      
aaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGaa                                                                       <  1:3636488/71‑1 (MQ=255)
aaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                                        <  1:3649088/70‑1 (MQ=255)
 caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                                      >  1:1113820‑M2/6‑71 (MQ=255)
 caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                                      >  1:2992054‑M2/6‑71 (MQ=255)
 aaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGaaa                                                                      >  1:2747492/1‑71 (MQ=255)
   cGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAacc                                                                    <  1:4626504/71‑1 (MQ=255)
      cAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                                                 >  1:4228261/1‑71 (MQ=255)
      cAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                                                 >  1:3418743/1‑71 (MQ=255)
        aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                               <  1:1202272/71‑1 (MQ=255)
        aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCgg                                                               >  1:3122144/1‑71 (MQ=255)
         gTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCg                                                                >  1:4254490/1‑69 (MQ=255)
           cTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTAt                                                            >  1:1805847/1‑71 (MQ=255)
           cTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTAt                                                            >  1:4672548/1‑71 (MQ=255)
              aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                                         >  1:1150974/1‑71 (MQ=255)
              aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                          <  1:1574862/70‑1 (MQ=255)
              aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                          <  1:2682757/70‑1 (MQ=255)
              aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAac                                                                     <  1:146346/59‑1 (MQ=255)
              aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGaa                                                         >  1:2719949/1‑71 (MQ=255)
              aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGa                                                          <  1:2501731/70‑1 (MQ=255)
                aTAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                        >  1:505475/1‑70 (MQ=255)
                 tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCa                                                      <  1:705295/71‑1 (MQ=255)
                  aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCg                                                                >  1:4079212/1‑60 (MQ=255)
                  aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCg                                                                >  1:1563414/1‑60 (MQ=255)
                   aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATg                                                           >  1:741472/1‑64 (MQ=255)
                   aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGa                                                    <  1:2911570/71‑1 (MQ=255)
                     tGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                                  <  1:2872179/71‑1 (MQ=255)
                        ccGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTgg                                                <  1:207230/70‑1 (MQ=255)
                        ccGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGt                                               >  1:2397781/1‑71 (MQ=255)
                           cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCg                                                                <  1:3210479/51‑1 (MQ=255)
                            tATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATg                                                           <  1:380443/55‑1 (MQ=255)
                              tGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCAAACCACCGGTATGAACCAGAt                                                   <  1:2281671/60‑1 (MQ=255)
                                aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                                       >  1:3173510/1‑71 (MQ=255)
                                aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                                       >  1:1004069/1‑71 (MQ=255)
                                  cACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                     <  1:3586474/71‑1 (MQ=255)
                                   aCTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                     <  1:1477909/70‑1 (MQ=255)
                                   aCTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                                    <  1:176907/71‑1 (MQ=255)
                                    cTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                     >  1:588660/1‑69 (MQ=255)
                                     tGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                                                 <  1:1437343/40‑1 (MQ=255)
                                     tGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccac                                                                  <  1:3768351/39‑1 (MQ=38)
                                     tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                                    <  1:4591856/69‑1 (MQ=255)
                                     tGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCgg                                                               <  1:2127903/42‑1 (MQ=255)
                                      gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGt                                               >  1:947090/1‑57 (MQ=255)
                                       cAATTACACGCCAGATCGTTCTCGCTACCGCAaccac                                                                  >  1:2776130/1‑37 (MQ=38)
                                       cAATTACACGCCAGATCGTTCTCGCTACCGCAaccac                                                                  >  1:4415381/1‑37 (MQ=38)
                                       cAATTACACGCCAGATCGTTCTCGCTACCGCAaccac                                                                  >  1:464248/1‑37 (MQ=38)
                                       cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                               >  1:2298044/1‑40 (MQ=255)
                                       cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAcc                                                       <  1:1947885/48‑1 (MQ=39)
                                       cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAcc                                                       <  1:844499/48‑1 (MQ=39)
                                       cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                                >  1:1412939/1‑71 (MQ=255)
                                       cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                                >  1:4007542/1‑71 (MQ=255)
                                          ttACACGCCAGATCGTTCTCGATACCGAAACCACCGGTAt                                                            <  1:4476593/40‑1 (MQ=255)
                                          ttACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                                    <  1:616363/64‑1 (MQ=255)
                                          ttACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                             <  1:1922553/71‑1 (MQ=255)
                                          ttACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCca                              >  1:2542920/1‑70 (MQ=255)
                                           tACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGt                                               >  1:1625861/1‑52 (MQ=255)
                                            acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                                  <  1:3669763/48‑1 (MQ=39)
                                            acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcaca                            <  1:2933097/70‑1 (MQ=255)
                                            acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                             >  1:124083/1‑69 (MQ=255)
                                            acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATg                                                           >  1:3176747/1‑39 (MQ=255)
                                            acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGa                                                    >  1:3267349/1‑46 (MQ=255)
                                            acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGa                                                    >  1:436897/1‑46 (MQ=255)
                                            acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                             >  1:4372900/1‑69 (MQ=255)
                                            acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                             >  1:312823/1‑69 (MQ=255)
                                               cGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt                        <  1:308889/71‑1 (MQ=255)
                                                gCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                             >  1:2846036/1‑65 (MQ=255)
                                                 ccAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa                      >  1:2967536/1‑71 (MQ=255)
                                                  cAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                                         <  1:203686/51‑1 (MQ=255)
                                                  cAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAt                     >  1:3197548/1‑71 (MQ=255)
                                                   aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                    <  1:4603498/71‑1 (MQ=255)
                                                   aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                    <  1:676075/71‑1 (MQ=255)
                                                    gATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                                       >  1:1991060/1‑51 (MQ=255)
                                                    gATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                    >  1:800992/1‑70 (MQ=255)
                                                     aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                    >  1:2662580/1‑69 (MQ=255)
                                                      tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAt                     >  1:3443034/1‑67 (MQ=255)
                                                       cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTgg                                                >  1:1940933/1‑39 (MQ=255)
                                                       cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                    >  1:4281171/1‑67 (MQ=255)
                                                        gTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                     >  1:1877815/1‑49 (MQ=255)
                                                         ttCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                     <  1:1748616/48‑1 (MQ=255)
                                                         ttCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                     <  1:1072824/48‑1 (MQ=255)
                                                         ttCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAtt              <  1:769867/71‑1 (MQ=255)
                                                         ttCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAtt              <  1:1159819/71‑1 (MQ=255)
                                                          tctcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTg             >  1:2614883/1‑71 (MQ=255)
                                                           ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGTCCACAAGATCATTGAGAtt              >  1:2141195/1‑69 (MQ=255)
                                                            tcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGt           >  1:1906321/1‑71 (MQ=255)
                                                              gctACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                    <  1:3692954/58‑1 (MQ=255)
                                                               aTACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                     <  1:3019337/42‑1 (MQ=255)
                                                               aTACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAtt              <  1:1182512/65‑1 (MQ=255)
                                                                tACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCg       >  1:4101883/1‑71 (MQ=255)
                                                                  ccGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt               <  1:2886011/61‑1 (MQ=255)
                                                                  ccGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt               <  1:388854/61‑1 (MQ=255)
                                                                  ccGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt               <  1:348959/61‑1 (MQ=255)
                                                                  ccGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCg       <  1:210310/69‑1 (MQ=255)
                                                                  ccGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCg       <  1:2340138/69‑1 (MQ=255)
                                                                  ccGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGtt     <  1:2720167/71‑1 (MQ=255)
                                                                  ccGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                     <  1:4549789/39‑1 (MQ=255)
                                                                  ccGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGtt     <  1:3960808/71‑1 (MQ=255)
                                                                   cGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGtt     >  1:2185848/1‑70 (MQ=255)
                                                                   cGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTg    <  1:2084750/71‑1 (MQ=255)
                                                                    gAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTGa   <  1:3825581/71‑1 (MQ=255)
                                                                     aaaCCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTGaa  <  1:414649/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: