Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 236,107 A→C 56.6% 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,1070AC56.6% 192.2 / 161.7 159E14A (GAA→GCA) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base A (36/33);  new base C (46/44);  total (82/77)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

Only 100 of 161 total aligned reads displayed.

ACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTG  >  W3110S.gb/236043‑236175
                                                                |                                                                    
tCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccac                                                                >  1:131006‑M2/2‑71 (MQ=255)
tCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccac                                                                >  1:778568‑M2/2‑71 (MQ=255)
   aGTCTGACATAAATGCCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                             >  1:2687040/1‑71 (MQ=255)
   aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAacc                                                                  <  1:1569766/66‑1 (MQ=255)
   aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                             <  1:1310181/71‑1 (MQ=255)
       tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTAt                                                          >  1:1179559/1‑70 (MQ=255)
         aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATg                                                         <  1:3497549/69‑1 (MQ=255)
          cATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                                      >  1:1445659/1‑71 (MQ=255)
            tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccac                                                                <  1:2039662/59‑1 (MQ=255)
            tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccac                                                                <  1:3304728/59‑1 (MQ=255)
            tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                      >  1:710494/1‑69 (MQ=255)
            tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAcc                                                     >  1:399619/1‑70 (MQ=255)
             aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccac                                                                >  1:1287763/1‑58 (MQ=255)
             aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccac                                                                >  1:3381926/1‑58 (MQ=255)
             aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccac                                                                >  1:1508856/1‑58 (MQ=255)
             aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACAGCAaccac                                                                >  1:1781532/1‑58 (MQ=255)
             aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAac                                                                   >  1:1345628/1‑55 (MQ=255)
              aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAg                                                   <  1:2975761/70‑1 (MQ=255)
              aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                                      >  1:2357676/1‑67 (MQ=255)
              aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCa                                                    >  1:570570/1‑69 (MQ=255)
                tGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTa                                                           >  1:2519216/1‑60 (MQ=255)
                tGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAtt                                                <  1:355809/71‑1 (MQ=255)
                tGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAtt                                                <  1:308638/71‑1 (MQ=255)
                  aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                               <  1:2810635/70‑1 (MQ=255)
                  aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                               <  1:2175756/70‑1 (MQ=255)
                      cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                          >  1:1837101/1‑71 (MQ=255)
                       tATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                        <  1:1986626/56‑1 (MQ=255)
                         tGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAt                                                 <  1:298307/61‑1 (MQ=255)
                         tGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                                       >  1:1507451/1‑71 (MQ=255)
                           aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACt                                      <  1:3652424/70‑1 (MQ=255)
                           aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                                     >  1:1430097/1‑71 (MQ=255)
                            gCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGa                                                  <  1:1918796/57‑1 (MQ=255)
                            gCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTAt                                    <  1:1929646/71‑1 (MQ=255)
                             cACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                   <  1:619702/71‑1 (MQ=255)
                             cACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                   <  1:3459523/71‑1 (MQ=255)
                               cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                   >  1:2252280/1‑69 (MQ=255)
                                tGCAATTACACGCCAGATCGTTCTCGCTACCGCAac                                                                   <  1:3221797/36‑1 (MQ=38)
                                tGCAATTACACGCCAGATCGTTCTCGATACCGAAacca                                                                 >  1:367168/1‑38 (MQ=255)
                                 gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                               >  1:2928611/1‑71 (MQ=255)
                                 gCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACAAGATTGGTGCGCACTATGAAgg                               >  1:886102/1‑71 (MQ=255)
                                  cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCa                                        >  1:3372638/1‑61 (MQ=255)
                                   aaTTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                             <  1:887832/39‑1 (MQ=38)
                                   aaTTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                             <  1:1719089/39‑1 (MQ=38)
                                    aTTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTg                                               <  1:1931211/52‑1 (MQ=255)
                                    aTTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCca                            >  1:1299917/1‑71 (MQ=255)
                                     ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATg                                                         <  1:3209036/41‑1 (MQ=38)
                                     ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                                    <  1:1014464/46‑1 (MQ=39)
                                     ttACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                              <  1:2558535/68‑1 (MQ=255)
                                      tACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                                       >  1:2933083/1‑42 (MQ=38)
                                      tACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                                     >  1:3013918/1‑60 (MQ=255)
                                      tACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGa                                                  >  1:1102695/1‑47 (MQ=255)
                                       acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                                 >  1:1747825/1‑63 (MQ=255)
                                       acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                           >  1:1396514/1‑69 (MQ=255)
                                       acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                           >  1:363724/1‑69 (MQ=255)
                                       acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                           >  1:1755841/1‑69 (MQ=255)
                                       acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                           >  1:797105/1‑69 (MQ=255)
                                       acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                           >  1:2307619/1‑69 (MQ=255)
                                       acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGa                                                  >  1:2644565/1‑46 (MQ=255)
                                       acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcaca                          <  1:3519001/70‑1 (MQ=255)
                                       acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                           >  1:308840/1‑69 (MQ=255)
                                        cacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAg                        <  1:3177293/71‑1 (MQ=255)
                                        cacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAg                        <  1:1649251/71‑1 (MQ=255)
                                        cacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                   <  1:395867/60‑1 (MQ=255)
                                          cGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt                      <  1:205411/71‑1 (MQ=255)
                                           gCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATc                     <  1:1004515/71‑1 (MQ=255)
                                            ccAGATCGTTCTCGCTACCGCAACCCCCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATc                     <  1:83578/70‑1 (MQ=255)
                                            ccAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa                    <  1:1674635/71‑1 (MQ=255)
                                             cAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGa                                                  >  1:1024491/1‑40 (MQ=38)
                                             cAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAt                   >  1:637253/1‑71 (MQ=255)
                                              aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAt                                                 >  1:877412/1‑40 (MQ=38)
                                              aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                                  <  1:2743702/55‑1 (MQ=255)
                                              aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                  <  1:2266333/71‑1 (MQ=255)
                                              aGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                  <  1:3612955/71‑1 (MQ=255)
                                               gATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTg                                            <  1:3292664/44‑1 (MQ=255)
                                                aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTAt                                    <  1:2901253/51‑1 (MQ=255)
                                                aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg                                >  1:3134863/1‑55 (MQ=255)
                                                aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCca                            <  1:2393204/59‑1 (MQ=255)
                                                aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga                <  1:2396305/71‑1 (MQ=255)
                                                aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa                    <  1:137339/67‑1 (MQ=255)
                                                aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa                    <  1:517795/67‑1 (MQ=255)
                                                aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa                    <  1:1168951/67‑1 (MQ=255)
                                                 tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                   <  1:2092654/51‑1 (MQ=255)
                                                 tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                   <  1:2363089/51‑1 (MQ=255)
                                                 tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                  >  1:2581698/1‑68 (MQ=255)
                                                 tCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgag               <  1:664927/71‑1 (MQ=255)
                                                  cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                  >  1:1912814/1‑67 (MQ=255)
                                                  cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                  >  1:2150163/1‑67 (MQ=255)
                                                  cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                  >  1:2390760/1‑67 (MQ=255)
                                                   gTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgaga              <  1:3176692/70‑1 (MQ=255)
                                                    ttCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAtt            <  1:784854/71‑1 (MQ=255)
                                                      ctctATACCGAAACCACCGGTATGAACCAGATTGATGCGCAc                                       >  1:1283432/1‑42 (MQ=38)
                                                      ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                                 >  1:850571/1‑48 (MQ=39)
                                                      ctcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTg           >  1:2598438/1‑70 (MQ=255)
                                                       tcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                   <  1:84004/45‑1 (MQ=255)
                                                            aCCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga                >  1:617946/1‑59 (MQ=255)
                                                            aCCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga                >  1:2715974/1‑59 (MQ=255)
                                                              cGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTg  <  1:3417037/71‑1 (MQ=255)
                                                               gAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg                                >  1:1863063/1‑40 (MQ=255)
                                                               gAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTg  >  1:1078828/1‑70 (MQ=255)
                                                                aaaCCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGa                       <  1:3240550/48‑1 (MQ=255)
                                                                |                                                                    
ACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTG  >  W3110S.gb/236043‑236175

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: