Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,587,198 C→T 53.5% 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,1980CT53.5% 2.2 / 61.5 43E840K (GAA→AAA) polAfused DNA polymerase I 5'‑>3' exonuclease, 3'‑>5' polymerase and 3'‑>5' exonuclease
Reads supporting (aligned to +/- strand):  ref base C (18/2);  new base T (13/10);  total (31/12)
Fisher's exact test for biased strand distribution p-value = 1.93e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCGCCCGTTTGATAATGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATACA  >  W3110S.gb/3587137‑3587248
                                                             |                                                  
atcGCCCGTTTGATAATGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGc                                              >  1:1640231/1‑68 (MQ=255)
atcGCCCGTTTGATAATGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTg                                            >  1:850103/1‑70 (MQ=255)
atcGCCCGTTTGATAATGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTg                                            >  1:811548/1‑70 (MQ=255)
atcGCCCGTTTGATAATGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTg                                            >  1:310835/1‑70 (MQ=255)
atcGCCCGTTTGATAATGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTg                                            >  1:221094/1‑70 (MQ=255)
atcGCCCGTTTGATAATGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTg                                            >  1:1436987/1‑70 (MQ=255)
atcGCCCGTTTGATAATGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTg                                            >  1:1962874/1‑70 (MQ=255)
atcGCCCGTTTGATAATGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTg                                            >  1:1162709/1‑70 (MQ=255)
  cGCCCGTTTGATAATGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCa                                          >  1:536035/1‑70 (MQ=255)
     ccGTTTGATAATGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCacgac                                      <  1:148154/71‑1 (MQ=255)
        tttGATAATGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAg                                   <  1:1827805/71‑1 (MQ=255)
        tttGATAATGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAg                                   <  1:693069/71‑1 (MQ=255)
        tttGATAATGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAg                                   <  1:688824/71‑1 (MQ=255)
          tGATAATGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCa                                 >  1:326186/1‑71 (MQ=255)
                                   cATTGGCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGata               >  1:1219720/1‑64 (MQ=255)
                                   cATTGGCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGata               >  1:1877540/1‑64 (MQ=255)
                                   cATTGGCGCGTTAATGGATGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGata               >  1:1367405/1‑64 (MQ=255)
                                    aTTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAg       <  1:206700/71‑1 (MQ=255)
                                    aTTGGCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCgg          >  1:1745157/1‑68 (MQ=255)
                                    aTTGGCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCgg          >  1:1703753/1‑68 (MQ=255)
                                    aTTGGCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCgg          >  1:1507422/1‑68 (MQ=255)
                                    aTTGGCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCgg          >  1:468564/1‑68 (MQ=255)
                                    aTTGGCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCgg          >  1:453689/1‑68 (MQ=255)
                                    aTTGGCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCgg          >  1:407455/1‑68 (MQ=255)
                                     ttGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAg       <  1:1901738/70‑1 (MQ=255)
                                       ggCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATa    >  1:1611827/1‑71 (MQ=255)
                                       ggCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATa    >  1:599620/1‑71 (MQ=255)
                                       ggCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATa    >  1:555250/1‑71 (MQ=255)
                                       ggCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATa    >  1:36197/1‑71 (MQ=255)
                                       ggCGCGTTAATGGCTGCACGTTCAGCCGCTGCacgac                                      >  1:1377432/1‑37 (MQ=255)
                                       ggCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGAt     >  1:1289430/1‑70 (MQ=255)
                                       ggCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGAt     >  1:1696073/1‑70 (MQ=255)
                                       ggCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGAt     >  1:835138/1‑70 (MQ=255)
                                       ggCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATa    >  1:1494457/1‑71 (MQ=255)
                                       ggCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATa    >  1:528289/1‑71 (MQ=255)
                                       ggCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATa    >  1:1247248/1‑71 (MQ=255)
                                       ggCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATa    >  1:1751417/1‑71 (MQ=255)
                                         cgcgTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATACa  <  1:1544441/71‑1 (MQ=255)
                                         cgcgTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATACa  <  1:1920576/71‑1 (MQ=255)
                                         cgcgTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATACa  <  1:1275763/71‑1 (MQ=255)
                                         cgcgTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATACa  <  1:1273485/71‑1 (MQ=255)
                                         cgcgTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATACa  <  1:158396/71‑1 (MQ=255)
                                         cgcgTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATACa  <  1:1756897/71‑1 (MQ=255)
                                                             |                                                  
ATCGCCCGTTTGATAATGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATACA  >  W3110S.gb/3587137‑3587248

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: