Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 236,107 A→C 52.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,1070AC52.3% 51.0 / 76.9 67E14A (GAA→GCA) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base A (18/14);  new base C (25/10);  total (43/24)
Fisher's exact test for biased strand distribution p-value = 2.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

CGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTG  >  W3110S.gb/236041‑236175
                                                                  |                                                                    
cGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAac                                                                   >  1:494688/1‑70 (MQ=255)
   cAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                                    >  1:1628993/1‑66 (MQ=255)
   cAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                                    >  1:884852/1‑66 (MQ=255)
   cAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                                    >  1:501194/1‑66 (MQ=255)
    aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAac                                                                   >  1:208141/1‑66 (MQ=255)
     aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                             >  1:731991/1‑71 (MQ=255)
        cTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTAGATACCGAAACCACCGGTAt                                                          >  1:888141/1‑71 (MQ=255)
         tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTAt                                                          >  1:599394/1‑70 (MQ=255)
             aTAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                      >  1:1483172/1‑70 (MQ=255)
                    aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTgg                                              <  1:1481463/71‑1 (MQ=255)
                        cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                          <  1:1235363/71‑1 (MQ=255)
                        cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTgg                                              <  1:703961/67‑1 (MQ=255)
                           tGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCa                                        <  1:159805/70‑1 (MQ=255)
                           tGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCa                                        <  1:967285/70‑1 (MQ=255)
                           tGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCa                                        <  1:1004747/70‑1 (MQ=255)
                            gAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGACACCACCGGTATGAACCAGATTGGTGCGCa                                        <  1:869124/69‑1 (MQ=255)
                             aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                                     >  1:497534/1‑71 (MQ=255)
                             aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                                     >  1:996129/1‑71 (MQ=255)
                             aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                                     >  1:1087883/1‑71 (MQ=255)
                              gCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                               <  1:815331/60‑1 (MQ=255)
                                aCTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                                  <  1:877908/71‑1 (MQ=255)
                                aCTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                                  <  1:688119/71‑1 (MQ=255)
                                 cTGCAATTACACGCCAGATCGTTCTTGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                                 <  1:518816/71‑1 (MQ=255)
                                 cTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                   <  1:640324/69‑1 (MQ=255)
                                  tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                                       <  1:1546564/64‑1 (MQ=255)
                                  tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                                       <  1:1127632/64‑1 (MQ=255)
                                   gCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                               >  1:304755/1‑71 (MQ=255)
                                    cAATTACACGCCAGGTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                              >  1:1711069/1‑71 (MQ=255)
                                    cAATTACACGCCAGATCGTTCTTGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                              >  1:1563419/1‑71 (MQ=255)
                                    cAATTACACGCCAGATCGTTCTTGCTACCGCAACCACCGGTATGAACAAGATTGGTGCGCACTATGAAGGc                              >  1:436675/1‑71 (MQ=255)
                                    cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                                       >  1:1281927/1‑62 (MQ=255)
                                    cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                                       >  1:778720/1‑62 (MQ=255)
                                    cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                                     >  1:202239/1‑64 (MQ=255)
                                    cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                                     >  1:1821145/1‑64 (MQ=255)
                                    cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                              >  1:1759989/1‑71 (MQ=255)
                                    cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                              >  1:350589/1‑71 (MQ=255)
                                    cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                              >  1:816723/1‑71 (MQ=255)
                                    cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCACATTGGTGCGCACTATGAAgg                               >  1:125011/1‑70 (MQ=255)
                                    cAATTACACGCCAGATCGTTCTCGCTACAGCAACCACAGGTATGAACCAGATTGGTGCGCACTATGAAGGc                              >  1:359115/1‑71 (MQ=255)
                                    cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGaa                                                       >  1:1276839/1‑46 (MQ=255)
                                      aTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                      <  1:598452/45‑1 (MQ=38)
                                      aTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                      <  1:200/45‑1 (MQ=38)
                                         acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                                 >  1:1130581/1‑63 (MQ=255)
                                           acGCCAGATCGTTCTTGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                           >  1:660609/1‑67 (MQ=255)
                                               cAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATAGGTGCGCACt                                      >  1:872467/1‑52 (MQ=255)
                                                  aTCGTTCTCTATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga                <  1:497454/71‑1 (MQ=255)
                                                  aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                  >  1:1298383/1‑69 (MQ=255)
                                                  aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga                <  1:443928/71‑1 (MQ=255)
                                                   tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                  >  1:474418/1‑68 (MQ=255)
                                                   tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                  >  1:640063/1‑68 (MQ=255)
                                                   tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                  >  1:1234040/1‑68 (MQ=255)
                                                   tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAt                   >  1:137290/1‑67 (MQ=255)
                                                   tCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgag               <  1:360027/71‑1 (MQ=255)
                                                    cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACaa                         >  1:749358/1‑60 (MQ=255)
                                                    cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCa                                        >  1:465799/1‑45 (MQ=255)
                                                    cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCa                                        >  1:440303/1‑45 (MQ=255)
                                                    cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCa                                        >  1:908944/1‑45 (MQ=255)
                                                    cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                  >  1:1616076/1‑67 (MQ=255)
                                                    cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                  >  1:553259/1‑67 (MQ=255)
                                                    cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                  >  1:1332441/1‑67 (MQ=255)
                                                    cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                  >  1:471565/1‑67 (MQ=255)
                                                    cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                  >  1:446787/1‑67 (MQ=255)
                                                    cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                  >  1:407711/1‑67 (MQ=255)
                                                    cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                  >  1:324347/1‑67 (MQ=255)
                                                     gTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt             <  1:147362/71‑1 (MQ=255)
                                                      ttCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt             <  1:1534543/70‑1 (MQ=255)
                                                      ttCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt             <  1:110810/70‑1 (MQ=255)
                                                        ctcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTgg          >  1:963755/1‑71 (MQ=255)
                                                             tACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                   <  1:214736/41‑1 (MQ=255)
                                                                cGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTg  <  1:1751859/71‑1 (MQ=255)
                                                                  |                                                                    
CGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTG  >  W3110S.gb/236041‑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: