Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 236,101 A→C 100% 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,1010AC97.9% 3246.5 / 16.1 881D12A (GAT→GCT) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base A (13/6);  new base C (482/379);  total (496/385)
Fisher's exact test for biased strand distribution p-value = 3.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 881 total aligned reads displayed.

AAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGG  >  W3110S.gb/236037‑236167
                                                                |                                                                  
tacaactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACcgc                                                              <  1:2771042‑M2/64‑2 (MQ=255)
tacaactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCg                                                               >  1:1177023‑M2/8‑70 (MQ=255)
  caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                            >  1:563957‑M2/6‑71 (MQ=255)
  caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                            >  1:1948443‑M2/6‑71 (MQ=255)
  caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                            >  1:1346872‑M2/6‑71 (MQ=255)
       cAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                                       >  1:46157/1‑71 (MQ=255)
       cAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                                       >  1:747342/1‑71 (MQ=255)
         aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                     <  1:3036703/71‑1 (MQ=255)
         aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                     <  1:847266/71‑1 (MQ=255)
               aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                  <  1:1665103/68‑1 (MQ=255)
                cATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                               >  1:2527826/1‑70 (MQ=255)
                 aTAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                              >  1:804459/1‑70 (MQ=255)
                 aTAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                              >  1:2906026/1‑70 (MQ=255)
                  tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                            <  1:127643/71‑1 (MQ=255)
                  tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                            <  1:3223696/71‑1 (MQ=255)
                  tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                            <  1:2723859/71‑1 (MQ=255)
                  tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                            <  1:2174397/71‑1 (MQ=255)
                    aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                              >  1:22997/1‑67 (MQ=255)
                    aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                              >  1:1857645/1‑67 (MQ=255)
                    aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                              >  1:243692/1‑67 (MQ=255)
                    aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                              >  1:2138942/1‑67 (MQ=255)
                    aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                              >  1:2117042/1‑67 (MQ=255)
                    aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                              >  1:656846/1‑67 (MQ=255)
                    aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                              >  1:2601731/1‑67 (MQ=255)
                    aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                              >  1:2971733/1‑67 (MQ=255)
                    aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                              >  1:1001019/1‑67 (MQ=255)
                    aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAg                                           <  1:1778676/70‑1 (MQ=255)
                    aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAg                                           <  1:6829/70‑1 (MQ=255)
                    aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                              >  1:1883637/1‑67 (MQ=255)
                           gCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccac                                                        <  1:643519/50‑1 (MQ=255)
                            cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                  <  1:1919656/71‑1 (MQ=255)
                            cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                  <  1:1402081/71‑1 (MQ=255)
                            cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                  >  1:2234198/1‑71 (MQ=255)
                            cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                  >  1:1065953/1‑71 (MQ=255)
                                gAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                               >  1:3099404/1‑70 (MQ=255)
                                gAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                               >  1:358688/1‑70 (MQ=255)
                                 aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                             >  1:1972082/1‑71 (MQ=255)
                                 aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                             >  1:716026/1‑71 (MQ=255)
                                    aCTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                           <  1:2476041/70‑1 (MQ=255)
                                       gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                       >  1:877722/1‑71 (MQ=255)
                                       gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                       >  1:330041/1‑71 (MQ=255)
                                       gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                       >  1:1806478/1‑71 (MQ=255)
                                       gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                       >  1:621271/1‑71 (MQ=255)
                                        cAATTACCCGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                      >  1:579307/1‑71 (MQ=255)
                                        cAATTACACGCCAGCTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                      >  1:1602847/1‑71 (MQ=255)
                                        cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                      >  1:3009312/1‑71 (MQ=255)
                                        cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                      >  1:113551/1‑71 (MQ=255)
                                        cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                      >  1:938036/1‑71 (MQ=255)
                                        cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                      >  1:3260241/1‑71 (MQ=255)
                                        cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                      >  1:2650382/1‑71 (MQ=255)
                                        cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                      >  1:2394799/1‑71 (MQ=255)
                                        cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                      >  1:2352027/1‑71 (MQ=255)
                                        cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                      >  1:2088007/1‑71 (MQ=255)
                                        cAATTACACGCCAGATAGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                      >  1:1154277/1‑71 (MQ=255)
                                          aTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                <  1:3320123/43‑1 (MQ=38)
                                           ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                           <  1:1741978/63‑1 (MQ=255)
                                           ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                           <  1:1998906/63‑1 (MQ=255)
                                            tACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGATCCAGATTGGTgcg                                  >  1:1457909/1‑55 (MQ=255)
                                            tACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                  >  1:271551/1‑55 (MQ=255)
                                            tACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                           <  1:603287/62‑1 (MQ=255)
                                            tACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                         >  1:2383163/1‑64 (MQ=255)
                                             acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                   >  1:414826/1‑69 (MQ=255)
                                             acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                   >  1:1509516/1‑69 (MQ=255)
                                             acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                   >  1:2263481/1‑69 (MQ=255)
                                             acacGCCAGATAGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                   >  1:2039413/1‑69 (MQ=255)
                                              cacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACaa                 <  1:3150172/70‑1 (MQ=255)
                                              cacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACaa                 <  1:1044631/70‑1 (MQ=255)
                                              cacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACaa                 <  1:1495395/70‑1 (MQ=255)
                                               acGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACaa                 <  1:2624722/69‑1 (MQ=255)
                                                cGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                        >  1:382134/1‑45 (MQ=38)
                                                cGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt              <  1:2571230/71‑1 (MQ=255)
                                                cGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt              <  1:1312497/71‑1 (MQ=255)
                                                 gCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt              <  1:2834729/70‑1 (MQ=255)
                                                 gCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATc             <  1:3359909/71‑1 (MQ=255)
                                                    aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                           <  1:537229/54‑1 (MQ=255)
                                                     gATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                           <  1:2791119/53‑1 (MQ=255)
                                                      aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                           <  1:3171771/52‑1 (MQ=255)
                                                      aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                           <  1:2203336/52‑1 (MQ=255)
                                                      aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga        <  1:2680350/71‑1 (MQ=255)
                                                      aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga        <  1:2874595/71‑1 (MQ=255)
                                                      aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga        <  1:436028/71‑1 (MQ=255)
                                                      aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga        <  1:1198914/71‑1 (MQ=255)
                                                       tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgag       <  1:2694306/71‑1 (MQ=255)
                                                        cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTgg                                      >  1:483630/1‑39 (MQ=37)
                                                        cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTgg                                      >  1:1700503/1‑39 (MQ=37)
                                                        cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          >  1:3054419/1‑67 (MQ=255)
                                                        cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          >  1:1370882/1‑67 (MQ=255)
                                                        cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          >  1:1565090/1‑67 (MQ=255)
                                                        cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          >  1:1631291/1‑67 (MQ=255)
                                                        cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          >  1:1253052/1‑67 (MQ=255)
                                                        cGTTCTCGCTACCGCAACCACCGGTATGAACAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt     >  1:785385/1‑71 (MQ=255)
                                                         gTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgaga      <  1:1543039/70‑1 (MQ=255)
                                                          ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAtt    <  1:973302/71‑1 (MQ=255)
                                                          ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAtt    <  1:2938441/71‑1 (MQ=255)
                                                          ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAtt    <  1:3129800/71‑1 (MQ=255)
                                                            ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAGTATGAAGGCCACAAGATCAt           >  1:1216161/1‑62 (MQ=255)
                                                            ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAt           >  1:1094819/1‑62 (MQ=255)
                                                            ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga        >  1:2553035/1‑65 (MQ=255)
                                                            ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga        >  1:514361/1‑65 (MQ=255)
                                                            ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga        >  1:903268/1‑65 (MQ=255)
                                                                |                                                                  
AAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGG  >  W3110S.gb/236037‑236167

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: