Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,587,198 C→T 29.6% 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,1980CT29.6% 76.3 / 39.9 54E840K (GAA→AAA) polAfused DNA polymerase I 5'‑>3' exonuclease, 3'‑>5' polymerase and 3'‑>5' exonuclease
Reads supporting (aligned to +/- strand):  ref base C (16/22);  new base T (13/3);  total (29/25)
Fisher's exact test for biased strand distribution p-value = 1.54e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

CGTTTGATAATGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATACAGACGG  >  W3110S.gb/3587143‑3587253
                                                       |                                                       
cGTTTGATAATGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTCAGCCGCTGCacgacg                                          <  1:2655107/71‑1 (MQ=255)
cGTTTGATAATGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTCAGCCGCTGCacgacg                                          <  1:3165804/71‑1 (MQ=255)
cGTTTGATAATGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTCAGCCGCTGCacgacg                                          <  1:148984/71‑1 (MQ=255)
        aaTGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCAt                                  >  1:96731/1‑71 (MQ=255)
        aaTGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCAt                                  >  1:2420021/1‑71 (MQ=255)
        aaTGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCAt                                  >  1:2388737/1‑71 (MQ=255)
        aaTGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCAt                                  >  1:2812382/1‑71 (MQ=255)
        aaTGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACCTTTAGCCGCTGCACGACGAGCACCAt                                  >  1:1878073/1‑71 (MQ=255)
            tCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCt                              >  1:1713566/1‑71 (MQ=255)
                  gcggTTCCCTGCATTGGCGCGTTAATGGCTGCATGTTTAGCCGCTGCACGACGAGCACCATTGCTGGAttt                        <  1:3331224/71‑1 (MQ=255)
                   cggTTCCCTGCATTGGCGCGTTAATGGCTGCATGTTTAGCCGCTGCACGACGAGCACCATTGCTGGAttt                        <  1:3254718/70‑1 (MQ=255)
                   cggTTCCCTGCATTGGCGCGTTAATGGCTGCATGTTTAGCCGCTGCACGACGAGCACCATTGCTGGAttt                        <  1:3106709/70‑1 (MQ=255)
                     gTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGat                     <  1:3078920/71‑1 (MQ=255)
                     gTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGat                     <  1:1374536/71‑1 (MQ=255)
                     gTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGat                     <  1:371677/71‑1 (MQ=255)
                     gTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGat                     <  1:835900/71‑1 (MQ=255)
                            gCATTGGCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCgg               >  1:986708/1‑70 (MQ=255)
                            gCATTGGCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCgg               >  1:3354573/1‑70 (MQ=255)
                            gCATTGGCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCgg               >  1:2980757/1‑70 (MQ=255)
                            gCATTGGCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCgg               >  1:2974520/1‑70 (MQ=255)
                            gCATTGGCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCgg               >  1:2695523/1‑70 (MQ=255)
                            gCATTGGCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCgg               >  1:2689880/1‑70 (MQ=255)
                            gCATTGGCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCgg               >  1:261498/1‑70 (MQ=255)
                            gCATTGGCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCgg               >  1:389216/1‑70 (MQ=255)
                            gCATTGGCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCgg               >  1:2493317/1‑70 (MQ=255)
                            gCATTGGCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCgg               >  1:2449954/1‑70 (MQ=255)
                            gCATTGGCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCgg               >  1:2283367/1‑70 (MQ=255)
                            gCATTGGCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCgg               >  1:2159046/1‑70 (MQ=255)
                            gCATTGGCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCgg               >  1:1353648/1‑70 (MQ=255)
                            gCATTGGCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCgg               >  1:2139533/1‑70 (MQ=255)
                            gCATTGGCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCgg               >  1:1475278/1‑70 (MQ=255)
                                 ggCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGATTTGata                    >  1:48861/1‑60 (MQ=255)
                                 ggCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGATTTGata                    >  1:512304/1‑60 (MQ=255)
                                 ggCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGATTTGata                    >  1:556275/1‑60 (MQ=255)
                                 ggCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGATTTGata                    >  1:827062/1‑60 (MQ=255)
                                 ggCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGATTTGata                    >  1:1691730/1‑60 (MQ=255)
                                 ggCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGATTTGata                    >  1:894478/1‑60 (MQ=255)
                                 ggCGCGTTAATGGCTGCACGTTTAGCCGCTGCACGACGAGCACCATTGCTGGATTTGata                    >  1:2114523/1‑60 (MQ=255)
                                 ggCGCGTTAATGGCTGCACGTTTAGCCGCAGCacgacg                                          >  1:1986342/1‑38 (MQ=37)
                                   cgcgTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTTGATTTGATATCCGGCAGATACa       <  1:1194502/71‑1 (MQ=255)
                                   cgcgTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATACa       <  1:939067/71‑1 (MQ=255)
                                   cgcgTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATACa       <  1:1111391/71‑1 (MQ=255)
                                   cgcgTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATACa       <  1:2317816/71‑1 (MQ=255)
                                   cgcgTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATACa       <  1:1424557/71‑1 (MQ=255)
                                   cgcgTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATACa       <  1:1471950/71‑1 (MQ=255)
                                   cgcgTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATACa       <  1:3063704/71‑1 (MQ=255)
                                   cgcgTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATACa       <  1:214721/71‑1 (MQ=255)
                                   cgcgTTAATGGCGGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATACa       <  1:1986372/71‑1 (MQ=255)
                                    gcgTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATACa       <  1:67593/70‑1 (MQ=255)
                                         aaTGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATACAGACgg  <  1:1350500/70‑1 (MQ=255)
                                         aaTGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATACAGACgg  <  1:539913/70‑1 (MQ=255)
                                         aaTGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATACAGACgg  <  1:1666412/70‑1 (MQ=255)
                                         aaTGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATACAGACgg  <  1:1010906/70‑1 (MQ=255)
                                          aTGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATACAGACgg  <  1:2948281/69‑1 (MQ=255)
                                                       |                                                       
CGTTTGATAATGTCGGCGGCGGTTCCCTGCATTGGCGCGTTAATGGCTGCACGTTCAGCCGCTGCACGACGAGCACCATTGCTGGATTTGATATCCGGCAGATACAGACGG  >  W3110S.gb/3587143‑3587253

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: