Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,587,198 C→T 83.9% E840K (GAA→AAA)  polA ← fused DNA polymerase I 5'‑>3' exonuclease, 3'‑>5' polymerase and 3'‑>5' exonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,587,1980CT83.9% 71.1 / 8.3 31E840K (GAA→AAA) polAfused DNA polymerase I 5'‑>3' exonuclease, 3'‑>5' polymerase and 3'‑>5' exonuclease
Reads supporting (aligned to +/- strand):  ref base C (3/2);  new base T (12/14);  total (15/16)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

TCATCGCCCGTTTGATAATGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATACAGACGGCGTCCGTCCA  >  W3110S.gb/3587135‑3587263
                                                               |                                                                 
tcatcGCCCGTTTGATAATGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCg                                                              >  1:933466/1‑69 (MQ=255)
tcatcGCCCGTTTGATAATGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCg                                                              >  1:898350/1‑69 (MQ=255)
  atcGCCCGTTTGATAATGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGc                                                             >  1:220724/1‑68 (MQ=255)
     gCCCGTTTGATAATGTCGGCGGCGCTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGc                                                          >  1:2107069/1‑68 (MQ=255)
       ccGTTTGATAATGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCac                                                        <  1:1773101/68‑1 (MQ=255)
              aTAATGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGc                                                 <  1:3148596/68‑1 (MQ=255)
                aaTGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCAc                                               <  1:1910194/68‑1 (MQ=255)
                  tGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCAcc                                              >  1:3530367/1‑67 (MQ=255)
                   gTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCAtt                                           >  1:695305/1‑69 (MQ=255)
                   gTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCAtt                                           >  1:1745268/1‑69 (MQ=255)
                   gTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCAt                                            >  1:1263844/1‑68 (MQ=255)
                   gTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCAt                                            >  1:197000/1‑68 (MQ=255)
                    tCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTg                                          >  1:3097491/1‑69 (MQ=255)
                    tCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTg                                          >  1:2352748/1‑69 (MQ=255)
                          gcggTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCAc                                               >  1:443676/1‑58 (MQ=255)
                             gTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGATTTg                                 <  1:1943381/69‑1 (MQ=255)
                             gTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGATTTg                                 <  1:1492261/69‑1 (MQ=255)
                               tCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGATTTGat                               <  1:1064538/69‑1 (MQ=255)
                                        tGGCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCa                       >  1:413267/1‑68 (MQ=255)
                                         ggCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTg                                       <  1:2644392/51‑1 (MQ=255)
                                         ggCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTg                                       <  1:3565451/51‑1 (MQ=255)
                                         ggCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTg                                       <  1:102587/51‑1 (MQ=255)
                                             cgTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATACa                 <  1:827060/69‑1 (MQ=255)
                                                 aaTGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATACAGAc              <  1:3049240/68‑1 (MQ=255)
                                                          aCGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATACAGACGGcgtccg      >  1:1292375/1‑67 (MQ=255)
                                                           cGTTTAGCCGCTGCACGACGAGCACCATTGCTGGAttt                                  <  1:2931972/38‑1 (MQ=255)
                                                           cGTTTAGCCGCTGCACGACGAGCACCATTGCTGGAttt                                  <  1:2248125/38‑1 (MQ=255)
                                                           cGTTTAGCCGCTGCACGACGAGCACCATTGCTGGAttt                                  <  1:1422759/38‑1 (MQ=255)
                                                           cGTTTAGCCGCTGCACGACGAGCACCATTGCTGGAttt                                  <  1:861376/38‑1 (MQ=255)
                                                           cGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATACAGACGGcgtccgt     <  1:3550688/67‑1 (MQ=255)
                                                            gTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATACAGACGGCGTCCGTCCa  >  1:3435104/1‑69 (MQ=255)
                                                               |                                                                 
TCATCGCCCGTTTGATAATGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATACAGACGGCGTCCGTCCA  >  W3110S.gb/3587135‑3587263

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: