Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 671,132 C→T 9.7% R299H (CGT→CAT)  holA ← DNA polymerase III, delta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb671,1320CT 9.7% 289.3 / 16.8 93R299H (CGT→CAT) holADNA polymerase III, delta subunit
Reads supporting (aligned to +/- strand):  ref base C (34/50);  new base T (4/5);  total (38/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.14e-01

CCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATACCC  >  W3110S.gb/671064‑671198
                                                                    |                                                                  
cccACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGt                                                                  <  1:2292839/71‑1 (MQ=255)
cccACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGt                                                                  <  1:5119690/71‑1 (MQ=255)
cccACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGt                                                                  <  1:3331405/71‑1 (MQ=255)
cccACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGt                                                                  <  1:3366257/71‑1 (MQ=255)
cccACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGt                                                                  <  1:3578545/71‑1 (MQ=255)
cccACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGt                                                                  <  1:3595167/71‑1 (MQ=255)
 ccACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGt                                                                  <  1:2437067/70‑1 (MQ=255)
   acacTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAAc                                                               >  1:724623/1‑71 (MQ=255)
   acacTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAAc                                                               >  1:2529986/1‑71 (MQ=255)
   acacTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAAc                                                               >  1:2646525/1‑71 (MQ=255)
   acacTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAAc                                                               >  1:2687347/1‑71 (MQ=255)
   acacTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAAc                                                               >  1:1995348/1‑71 (MQ=255)
   acacTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAAc                                                               >  1:1745781/1‑71 (MQ=255)
   acacTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAAc                                                               >  1:4076636/1‑71 (MQ=255)
     actgactgACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACt                                                              >  1:5235341/1‑70 (MQ=255)
         actgacCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCg                                                           >  1:2032286/1‑69 (MQ=255)
         actgacCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGt                                                          >  1:3134803/1‑70 (MQ=255)
         actgacCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGt                                                          >  1:1076179/1‑70 (MQ=255)
         actgacCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGt                                                          >  1:2208361/1‑70 (MQ=255)
         actgacCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTc                                                         >  1:4731983/1‑71 (MQ=255)
         actgacCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTc                                                         >  1:5077401/1‑71 (MQ=255)
         actgacCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTc                                                         >  1:4616054/1‑71 (MQ=255)
         actgacCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTc                                                         >  1:1831829/1‑71 (MQ=255)
         actgacCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTc                                                         >  1:3495839/1‑71 (MQ=255)
         actgacCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGACTGACGTAACTGCGTc                                                         >  1:1048012/1‑71 (MQ=255)
           tgacCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCt                                                        >  1:4811048/1‑70 (MQ=255)
             acCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGAc                                                     >  1:3793363/1‑71 (MQ=255)
             acCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGAc                                                     >  1:2181258/1‑71 (MQ=255)
             acCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAAATGCGTCTGAc                                                     >  1:2164055/1‑71 (MQ=255)
                      ttGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGAt                                            >  1:1007750/1‑71 (MQ=255)
                         tttGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGAttt                                          <  1:4329183/70‑1 (MQ=255)
                         tttGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGAttt                                          <  1:1644682/70‑1 (MQ=255)
                         tttGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGAttt                                          <  1:235965/70‑1 (MQ=255)
                         tttGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGAttt                                          <  1:2270666/70‑1 (MQ=255)
                         tttGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGAttt                                          <  1:2661982/70‑1 (MQ=255)
                         tttGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGAttt                                          <  1:1651494/70‑1 (MQ=255)
                         tttGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGAttt                                          <  1:2833064/70‑1 (MQ=255)
                         tttGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGAttt                                          <  1:3387281/70‑1 (MQ=255)
                          ttGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGAttt                                          <  1:2364276/69‑1 (MQ=255)
                             aGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACg                                      <  1:3387509/70‑1 (MQ=255)
                             aGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACg                                      <  1:5379679/70‑1 (MQ=255)
                               ggTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGAttt                                          <  1:5309051/64‑1 (MQ=255)
                                       ccGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGa                                                      >  1:4790842/1‑44 (MQ=255)
                                       ccGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGa                                                      >  1:5576730/1‑44 (MQ=255)
                                       ccGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCcat                            >  1:3221896/1‑70 (MQ=255)
                                                tCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCca                             <  1:3377818/60‑1 (MQ=255)
                                                tCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCca                             <  1:3224855/60‑1 (MQ=255)
                                                tCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCca                             <  1:828333/60‑1 (MQ=255)
                                                tCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCca                             <  1:985287/60‑1 (MQ=255)
                                                  aGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGAttt                                          <  1:65046/45‑1 (MQ=255)
                                                  aGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGAttt                                          <  1:734317/45‑1 (MQ=255)
                                                    gAGTTGCACGGCCTGATGTAACTGCGTCTGACTTAAGCGAttt                                          <  1:1003541/43‑1 (MQ=255)
                                                    gAGTTGCACGGCCTGATGTAACTGCGTCTGACTTAAGCGAttt                                          <  1:5346173/43‑1 (MQ=255)
                                                    gAGTTGCACGGCCTGATGTAACTGCGTCTGACTTAAGCGAttt                                          <  1:3523723/43‑1 (MQ=255)
                                                    gAGTTGCACGGCCTGATGTAACTGCGTCTGACTTAAGCGAttt                                          <  1:5061771/43‑1 (MQ=255)
                                                      gTTGCACGGCCTGATGTAACTGCGTCTGACTTAAGCGtttt                                          <  1:5017883/41‑1 (MQ=38)
                                                      gTTGCACGGCCTGATGTAACTGCGTCTGACTTAAGCGATTTaa                                        >  1:2157582/1‑43 (MQ=255)
                                                      gTTGCACGGCCTGATGTAACTGCGTCTGACTTAAGCGATTTaa                                        >  1:5059191/1‑43 (MQ=255)
                                                      gTTGCACGGCCTGATGTAACTGCGTCTGACTTAAGCGATTTaa                                        >  1:405024/1‑43 (MQ=255)
                                                      gTTGCACGGCCTGATGTAACTGCGTCTGACTTAAGCGATTTaa                                        >  1:3005599/1‑43 (MQ=255)
                                                       ttGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGc                                <  1:5280669/50‑1 (MQ=255)
                                                       ttGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGc                                <  1:2017861/50‑1 (MQ=255)
                                                       ttGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGc                                <  1:4532361/50‑1 (MQ=255)
                                                       ttGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGc                                <  1:1428357/50‑1 (MQ=255)
                                                       ttGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGc                                <  1:702585/50‑1 (MQ=255)
                                                       ttGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGc                                <  1:685039/50‑1 (MQ=255)
                                                        tGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGc                                <  1:150592/49‑1 (MQ=255)
                                                         gCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGc         >  1:2474628/1‑71 (MQ=255)
                                                                ctaaCGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATAccc  <  1:1753276/68‑1 (MQ=255)
                                                                cTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCt           <  1:1437721/62‑1 (MQ=255)
                                                                cTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCt           <  1:3612110/62‑1 (MQ=255)
                                                                cTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCt           <  1:383907/62‑1 (MQ=255)
                                                                cTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCt           <  1:138717/62‑1 (MQ=255)
                                                                cTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCt           <  1:5432388/62‑1 (MQ=255)
                                                                cTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCt           <  1:4375500/62‑1 (MQ=255)
                                                                cTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATAccc  <  1:2097588/71‑1 (MQ=255)
                                                                cTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATAccc  <  1:3144268/71‑1 (MQ=255)
                                                                cTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATAccc  <  1:772100/71‑1 (MQ=255)
                                                                cTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATAccc  <  1:95804/71‑1 (MQ=255)
                                                                cTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATAccc  <  1:3084606/71‑1 (MQ=255)
                                                                cTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATAccc  <  1:219437/71‑1 (MQ=255)
                                                                cTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATAccc  <  1:522890/71‑1 (MQ=255)
                                                                cTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATAccc  <  1:1484870/71‑1 (MQ=255)
                                                                cTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATAccc  <  1:1492946/71‑1 (MQ=255)
                                                                cTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATAccc  <  1:4493354/71‑1 (MQ=255)
                                                                cTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATAccc  <  1:4124246/71‑1 (MQ=255)
                                                                   aCGTAACTGCGTCTGACTTAATCGATTTAACGCCTCGCCcatca                          >  1:2396561/1‑44 (MQ=255)
                                                                   aCGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCcatca                          >  1:1366144/1‑44 (MQ=255)
                                                                   aCGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCcatca                          >  1:1604260/1‑44 (MQ=255)
                                                                   aCGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCcatca                          >  1:4783417/1‑44 (MQ=255)
                                                                   aCGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCcatca                          >  1:322547/1‑44 (MQ=255)
                                                                   aCGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCcatca                          >  1:3138433/1‑44 (MQ=255)
                                                                   aCGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCcatca                          >  1:87479/1‑44 (MQ=255)
                                                                    |                                                                  
CCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATACCC  >  W3110S.gb/671064‑671198

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: