Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,587,198 C→T 100% 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,1980CT94.3% 94.5 / ‑1.2 35E840K (GAA→AAA) polAfused DNA polymerase I 5'‑>3' exonuclease, 3'‑>5' polymerase and 3'‑>5' exonuclease
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (9/24);  total (9/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCAATCATCGCCCGTTTGATAATGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATACAGACGGCGTCCGTCC  >  W3110S.gb/3587131‑3587262
                                                                   |                                                                
gCAATCATCGCCCGTTTGATAATGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGc                                                               <  1:1419511/71‑1 (MQ=255)
      atcGCCCGTTTGATAATGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTg                                                          <  1:987543/70‑1 (MQ=255)
      atcGCCCGTTTGATAATGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTg                                                          >  1:3612331/1‑70 (MQ=255)
      atcGCCCGTTTGATAATGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTg                                                          >  1:194660/1‑70 (MQ=255)
      atcGCCAGTTTGATAATGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTg                                                          >  1:1323511/1‑70 (MQ=255)
              tttGATAATGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCAcgacga                                                  >  1:274302/1‑70 (MQ=255)
              tttGATAATGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAg                                                 >  1:3259340/1‑71 (MQ=255)
              tttGATAATGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAg                                                 >  1:3238037/1‑71 (MQ=255)
                      tGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCAtt                                          <  1:3413126/70‑1 (MQ=255)
                          ggcggcggTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTgg                                     <  1:1910400/71‑1 (MQ=255)
                              gcggTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGAttt                                 <  1:204286/71‑1 (MQ=255)
                              gcggTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGAttt                                 <  1:701851/71‑1 (MQ=255)
                               cggTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGAttt                                 <  1:2710792/70‑1 (MQ=255)
                                ggTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGATTTGa                               >  1:1314363/1‑71 (MQ=255)
                                 gTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGATTTGa                               <  1:1806869/70‑1 (MQ=255)
                                 gTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGATTTGa                               <  1:1376925/70‑1 (MQ=255)
                                 gTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGATTTGa                               <  1:2742050/70‑1 (MQ=255)
                                        gCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCa                                            >  1:3679510/1‑50 (MQ=255)
                                           ttGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGAttt                                 <  1:3107562/58‑1 (MQ=255)
                                           ttGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGAttt                                 <  1:1701710/58‑1 (MQ=255)
                                           ttGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGAttt                                 <  1:884894/58‑1 (MQ=255)
                                           ttGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGAttt                                 <  1:2360861/58‑1 (MQ=255)
                                           ttGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGAttt                                 <  1:3586903/58‑1 (MQ=255)
                                           ttGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGAttt                                 <  1:677945/58‑1 (MQ=255)
                                           ttGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGa                    <  1:3141645/71‑1 (MQ=255)
                                           ttGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGa                    <  1:2874805/71‑1 (MQ=255)
                                           ttGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGa                    <  1:854633/71‑1 (MQ=255)
                                           ttGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGa                    <  1:803635/71‑1 (MQ=255)
                                                          cTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATACAGACGGcgtccg     <  1:2048098/71‑1 (MQ=255)
                                                           tGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGATTTGa                               <  1:3495827/44‑1 (MQ=255)
                                                             cACGTTTAGCCGCTGCACGACTAGCACCATTGCTGGATTTGATATCCGGCAGATACAGACGGCgtccgtcc  <  1:3006798/71‑1 (MQ=255)
                                                             cACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATACAGACGGcgtccgtc   >  1:558494/1‑70 (MQ=255)
                                                             cACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATACAGACGGCgtccgtcc  <  1:3642916/71‑1 (MQ=255)
                                                             cACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATACAGACGGCgtccgtcc  <  1:3321155/71‑1 (MQ=255)
                                                             cACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATACAGACGGCgtccgtcc  <  1:2637076/71‑1 (MQ=255)
                                                                   |                                                                
GCAATCATCGCCCGTTTGATAATGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATACAGACGGCGTCCGTCC  >  W3110S.gb/3587131‑3587262

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: