Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 3,587,198 C→T 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,1980CT95.8% 67.3 / ‑3.6 24E840K (GAA→AAA) polAfused DNA polymerase I 5'‑>3' exonuclease, 3'‑>5' polymerase and 3'‑>5' exonuclease
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (11/12);  total (12/12)
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.86e-01

CGCCCGTTTGATAATGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATACAGACGGCGTCCG  >  W3110S.gb/3587139‑3587259
                                                           |                                                             
cGCCCGTTTGATAATGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCac                                                    >  1:1087466/1‑71 (MQ=255)
cGCCCGTTTGATAATGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCac                                                    >  1:1743911/1‑71 (MQ=255)
   ccGTTTGATAATGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCacgac                                                 >  1:824245/1‑71 (MQ=255)
   ccGTTTGATAATGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCacgac                                                 >  1:3775124/1‑71 (MQ=255)
   ccGTTTGATAATGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCacgac                                                 >  1:1886529/1‑71 (MQ=255)
   ccGTTTGATAATGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCacga                                                  >  1:865131/1‑70 (MQ=255)
   ccGTTTGATAATGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCacga                                                  >  1:4025966/1‑70 (MQ=255)
   ccGTTTGATAATGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCacga                                                  >  1:3039663/1‑70 (MQ=255)
            aaTGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCacgac                                                 >  1:174103/1‑62 (MQ=255)
            aaTGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCacgac                                                 >  1:3510950/1‑62 (MQ=255)
                tCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCt                                    >  1:644331/1‑71 (MQ=255)
                         gTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGAt                                <  1:22494/66‑1 (MQ=255)
                         gTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGAt                                <  1:3841912/66‑1 (MQ=255)
                         gTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGAt                                <  1:1343189/66‑1 (MQ=255)
                         gTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGAt                                <  1:350048/66‑1 (MQ=255)
                         gTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGAt                                <  1:3405371/66‑1 (MQ=255)
                         gTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGAt                                <  1:3054044/66‑1 (MQ=255)
                         gTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGAt                                <  1:2868489/66‑1 (MQ=255)
                         gTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGAt                                <  1:2518940/66‑1 (MQ=255)
                         gTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGAt                                <  1:2344411/66‑1 (MQ=255)
                         gTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGATTTGat                           <  1:617390/71‑1 (MQ=255)
                          ttCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGATTTGata                          <  1:89034/71‑1 (MQ=255)
                                                   tGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATACAGACGGcgtccg  >  1:1116835/1‑70 (MQ=255)
                                                       cGTTTAGCCGCTGCACGACGAGCACCATTCTGGAt                                <  1:2467269/35‑1 (MQ=25)
                                                           |                                                             
CGCCCGTTTGATAATGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATACAGACGGCGTCCG  >  W3110S.gb/3587139‑3587259

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: