Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 236,101 A→C 60.0% 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,1010AC60.0% 103.8 / 69.2 75D12A (GAT→GCT) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base A (18/12);  new base C (32/13);  total (50/25)
Fisher's exact test for biased strand distribution p-value = 3.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

TGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGG  >  W3110S.gb/236034‑236167
                                                                   |                                                                  
tGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGa                                                                    <  1:1523982/68‑1 (MQ=255)
tGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGa                                                                    <  1:1333635/68‑1 (MQ=255)
tGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGa                                                                    <  1:1331508/68‑1 (MQ=255)
gcgtacaactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTTGCTAc                                                                 >  1:442343‑M2/11‑71 (MQ=255)
   tacaactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCg                                                               >  1:1753381‑M2/8‑70 (MQ=255)
   tacaactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCg                                                               >  1:373268‑M2/8‑70 (MQ=255)
   tacaactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCg                                                               >  1:657119‑M2/8‑70 (MQ=255)
     caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTAccgca                                                             <  1:1699646‑M2/65‑3 (MQ=255)
     caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTAccgca                                                             <  1:1228524‑M2/65‑3 (MQ=255)
     caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTAccgca                                                             <  1:463176‑M2/65‑3 (MQ=255)
       cGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAac                                                           >  1:494688/1‑70 (MQ=255)
          cAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                            >  1:1628993/1‑66 (MQ=255)
          cAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                            >  1:501194/1‑66 (MQ=255)
          cAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                            >  1:884852/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:967285/70‑1 (MQ=255)
                                  tGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCa                                <  1:1004747/70‑1 (MQ=255)
                                  tGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCa                                <  1:159805/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:688119/71‑1 (MQ=255)
                                       aCTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                          <  1:877908/71‑1 (MQ=255)
                                        cTGCAATTACACGCCAGATCGTTCTTGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                         <  1:518816/71‑1 (MQ=255)
                                        cTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                           <  1:640324/69‑1 (MQ=255)
                                         tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                               <  1:1127632/64‑1 (MQ=255)
                                         tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                               <  1:1546564/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:778720/1‑62 (MQ=255)
                                           cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                               >  1:1281927/1‑62 (MQ=255)
                                           cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                             >  1:202239/1‑64 (MQ=255)
                                           cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                             >  1:1821145/1‑64 (MQ=255)
                                           cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                      >  1:816723/1‑71 (MQ=255)
                                           cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                      >  1:1759989/1‑71 (MQ=255)
                                           cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                      >  1:350589/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:200/45‑1 (MQ=38)
                                             aTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                              <  1:598452/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:1234040/1‑68 (MQ=255)
                                                          tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          >  1:640063/1‑68 (MQ=255)
                                                          tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          >  1:474418/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:908944/1‑45 (MQ=255)
                                                           cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCa                                >  1:465799/1‑45 (MQ=255)
                                                           cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCa                                >  1:440303/1‑45 (MQ=255)
                                                           cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          >  1:407711/1‑67 (MQ=255)
                                                           cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          >  1:553259/1‑67 (MQ=255)
                                                           cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          >  1:446787/1‑67 (MQ=255)
                                                           cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          >  1:1332441/1‑67 (MQ=255)
                                                           cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          >  1:324347/1‑67 (MQ=255)
                                                           cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          >  1:1616076/1‑67 (MQ=255)
                                                           cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          >  1:471565/1‑67 (MQ=255)
                                                            gTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt     <  1:147362/71‑1 (MQ=255)
                                                             ttCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt     <  1:110810/70‑1 (MQ=255)
                                                             ttCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt     <  1:1534543/70‑1 (MQ=255)
                                                               ctcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTgg  >  1:963755/1‑71 (MQ=255)
                                                                   |                                                                  
TGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGG  >  W3110S.gb/236034‑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: