Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 236,101 A→C 45.9% D12A (GAT→GCT)  dnaQ → DNA polymerase III epsilon subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb236,1010AC45.9% 20.9 / 145.7 83D12A (GAT→GCT) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base A (11/34);  new base C (22/16);  total (33/50)
Fisher's exact test for biased strand distribution p-value = 3.19e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01

AACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGG  >  W3110S.gb/236031‑236167
                                                                      |                                                                  
aaCTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGa                                                                    >  1:2248333/1‑71 (MQ=255)
aaCTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGa                                                                    <  1:133328/71‑1 (MQ=255)
 aCTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGa                                                                    <  1:2131729/70‑1 (MQ=255)
  cgcgtacaactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTa                                                                  >  1:1546287‑M2/12‑71 (MQ=255)
   gcgtacaactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTAc                                                                 >  1:2885711‑M2/11‑71 (MQ=255)
     tAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCg                                                               <  1:2252811/71‑1 (MQ=255)
       acaactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACcgc                                                              >  1:516653‑M2/7‑69 (MQ=255)
       aaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATAc                                                                 <  1:223028/67‑1 (MQ=255)
       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)
           gACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGAt                                                                   <  1:2134711/61‑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:1572169/70‑1 (MQ=255)
                     aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGa                                                <  1:65551/70‑1 (MQ=255)
                       aTAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGa                                                                    <  1:952154/48‑1 (MQ=255)
                       aTAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTAcc                                                                <  1:3145324/52‑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:1719885/1‑71 (MQ=255)
                                              cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                      >  1:968896/1‑71 (MQ=255)
                                              cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                      >  1:3013051/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:933621/39‑1 (MQ=255)
                                                            aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGAtt                                        <  1:1712743/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)
                                                                    cGATACCGAAACCACCGGTATGAACCAGATTGGTgcg                                  <  1:762956/37‑1 (MQ=255)
                                                                     gATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTgg  <  1:3069685/68‑1 (MQ=255)
                                                                      |                                                                  
AACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGG  >  W3110S.gb/236031‑236167

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: