Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 236,101 A→C 46.4% 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,1010AC46.4% 78.6 / inf 181D12A (GAT→GCT) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base A (54/43);  new base C (45/39);  total (99/82)
Fisher's exact test for biased strand distribution p-value = 8.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 188 total aligned reads displayed.

AACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGC  >  W3110S.gb/236031‑236170
                                                                      |                                                                     
aaCTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGa                                                                       <  1:2802552/71‑1 (MQ=255)
  cTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATa                                                                     >  1:4142555/1‑71 (MQ=255)
   gcgtacaactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTAc                                                                    >  1:4886509‑M2/11‑71 (MQ=255)
      tacaactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACcgc                                                                 <  1:2673899‑M2/64‑2 (MQ=255)
      tacaactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCg                                                                  >  1:3621451‑M2/8‑70 (MQ=255)
      aaaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGa                                                                       >  1:3276291/1‑65 (MQ=255)
       acaactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACcgc                                                                 >  1:2787218‑M2/7‑69 (MQ=255)
         aCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAac                                                              >  1:100163/1‑71 (MQ=255)
           ctCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAacca                                                            <  1:4186270‑M2/69‑1 (MQ=255)
            aCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAaccac                                                           <  1:1308127/71‑1 (MQ=255)
             cAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAaccacc                                                          >  1:4643843/1‑71 (MQ=255)
              aaGTGTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCg                                                         >  1:4877928/1‑71 (MQ=255)
              aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccac                                                           <  1:1585044/69‑1 (MQ=255)
              aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCg                                                         <  1:2278783/71‑1 (MQ=255)
               aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                        <  1:4031518/71‑1 (MQ=255)
                 tCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTa                                                      <  1:3746027/71‑1 (MQ=255)
                   tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAac                                                              >  1:2297121/1‑61 (MQ=255)
                   tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTAt                                                     >  1:4350701/1‑70 (MQ=255)
                   tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTAt                                                     >  1:2653224/1‑70 (MQ=255)
                     aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACTACCGGTATg                                                    <  1:4462072/69‑1 (MQ=255)
                     aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATg                                                    <  1:1600369/69‑1 (MQ=255)
                     aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATg                                                    <  1:1759926/69‑1 (MQ=255)
                     aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATg                                                    <  1:2449038/69‑1 (MQ=255)
                     aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATg                                                    <  1:4711429/69‑1 (MQ=255)
                      cATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                 <  1:893459/71‑1 (MQ=255)
                      cATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGaaa                                                               >  1:4296/1‑57 (MQ=255)
                       aTAAATGCCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                 >  1:4223910/1‑70 (MQ=255)
                       aTAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAcc                                                <  1:4680825/71‑1 (MQ=255)
                       aTAAATGACCGCTATGAGCACTGCAATTACACGACAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                 <  1:1804312/70‑1 (MQ=255)
                        tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                               <  1:154481/71‑1 (MQ=255)
                        tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                               <  1:2404001/71‑1 (MQ=255)
                        tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCa                                               <  1:2038083/71‑1 (MQ=255)
                         aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAacca                                                            >  1:299968/1‑57 (MQ=255)
                          aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                 >  1:2266129/1‑67 (MQ=255)
                          aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                 >  1:5135315/1‑67 (MQ=255)
                           aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAt                                            <  1:1846563/71‑1 (MQ=255)
                           aTGACCGCTACGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAac                                                              >  1:551965/1‑53 (MQ=255)
                            tGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                           >  1:5097137/1‑71 (MQ=255)
                            tGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                           <  1:4317765/71‑1 (MQ=255)
                            tGACCGCTATGAGCACTGCAATTACACACCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                           <  1:3685464/71‑1 (MQ=255)
                             gACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                          >  1:1053199/1‑71 (MQ=255)
                             gACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                          >  1:2647296/1‑71 (MQ=255)
                                cGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTg                                       >  1:1074981/1‑71 (MQ=255)
                                  cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                     >  1:4814754/1‑71 (MQ=255)
                                  cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTgcg                                     <  1:1926363/71‑1 (MQ=255)
                                    aTGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGa                                             >  1:3596417/1‑61 (MQ=255)
                                    aTGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGa                                             >  1:3334341/1‑61 (MQ=255)
                                    aTGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCa                                   >  1:4084066/1‑71 (MQ=255)
                                     tGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                   >  1:3412550/1‑54 (MQ=255)
                                       aGCACTGCAATTACACGCCAGATCGTTCTCGATACCGaaa                                                               >  1:2233879/1‑40 (MQ=255)
                                       aGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGt                                        >  1:1343470/1‑63 (MQ=255)
                                          aCTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                             <  1:3316757/71‑1 (MQ=255)
                                          aCTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                             <  1:2419581/71‑1 (MQ=255)
                                           cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                              <  1:3073929/69‑1 (MQ=255)
                                           cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                              <  1:3063450/69‑1 (MQ=255)
                                           cTGCAATTACACGCCAGATCGTTCTCGATACCGAAacca                                                            <  1:4572049/39‑1 (MQ=255)
                                           cTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                              >  1:2533904/1‑69 (MQ=255)
                                            tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAg                                              <  1:4521797/52‑1 (MQ=255)
                                             gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                           <  1:4761003/54‑1 (MQ=255)
                                             gCAATTACACGCCAGATCGTTCTCGATACCGAAaccacc                                                          >  1:3994626/1‑39 (MQ=255)
                                             gCAATTACACGCCAGATCGTTCTCGATACCGAAaccacc                                                          >  1:2590305/1‑39 (MQ=255)
                                              cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGa                                             >  1:952792/1‑51 (MQ=255)
                                              cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                         >  1:246431/1‑71 (MQ=255)
                                              cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                         >  1:1500900/1‑71 (MQ=255)
                                              cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                         >  1:1191971/1‑71 (MQ=255)
                                               aaTTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                         >  1:291419/1‑70 (MQ=255)
                                                 ttACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                        <  1:763002/37‑1 (MQ=38)
                                                 ttACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                        <  1:4654366/37‑1 (MQ=38)
                                                 ttACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCa                                               <  1:4375882/46‑1 (MQ=255)
                                                  tACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTgcg                                     >  1:1015061/1‑55 (MQ=255)
                                                  tACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTgcg                                     >  1:3781297/1‑55 (MQ=255)
                                                  tACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTgc                                      <  1:1565194/54‑1 (MQ=255)
                                                   acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCca                       >  1:1812946/1‑68 (MQ=255)
                                                   acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCa                                   >  1:2168819/1‑56 (MQ=255)
                                                      cGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                                >  1:3271909/1‑56 (MQ=255)
                                                           gATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt             >  1:1529664/1‑70 (MQ=255)
                                                           gATCGTTCTCGATACCGAAACCACCGGTATGAACCAGa                                             >  1:2178305/1‑38 (MQ=255)
                                                            aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt             >  1:598083/1‑69 (MQ=255)
                                                            aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga           <  1:1282973/71‑1 (MQ=255)
                                                            aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga           <  1:4519895/71‑1 (MQ=255)
                                                            aTCGTTCTCGCTACCGCAACCACGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt             >  1:859248/1‑68 (MQ=255)
                                                             tCGTTCTCGATACCGAAACCACCGGTATGAACCAGAt                                            >  1:4627794/1‑37 (MQ=255)
                                                             tCGTTCTCGATACAGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATAAtt             >  1:3713523/1‑68 (MQ=255)
                                                              cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg                           >  1:3209365/1‑53 (MQ=255)
                                                              cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgag          <  1:810107/70‑1 (MQ=255)
                                                              cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTgcg                                     >  1:668525/1‑43 (MQ=255)
                                                              cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                            >  1:1681825/1‑52 (MQ=255)
                                                              cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt             >  1:3096068/1‑67 (MQ=255)
                                                              cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgag          >  1:3635428/1‑70 (MQ=255)
                                                               gTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgaga         <  1:3475389/70‑1 (MQ=255)
                                                               gTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt        <  1:261680/71‑1 (MQ=255)
                                                                ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                          <  1:596764/52‑1 (MQ=255)
                                                                ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga           <  1:4986628/67‑1 (MQ=255)
                                                                ttCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                              <  1:705681/48‑1 (MQ=255)
                                                                ttCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                              <  1:1210723/48‑1 (MQ=255)
                                                                 tctcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                             <  1:1606982/48‑1 (MQ=255)
                                                                 tctcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                             <  1:1449551/48‑1 (MQ=255)
                                                                   tcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTgtt    >  1:3112171/1‑69 (MQ=255)
                                                                    cGATACCGAAACCACAGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt        >  1:406638/1‑66 (MQ=255)
                                                                      aTACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga           <  1:1389752/61‑1 (MQ=255)
                                                                      |                                                                     
AACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGC  >  W3110S.gb/236031‑236170

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: