Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,713,297 A→G 12.8% T21A (ACC→GCC)  nth → DNA glycosylase and apyrimidinic (AP) lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,713,2970AG12.8% 118.0 / 10.0 47T21A (ACC→GCC) nthDNA glycosylase and apyrimidinic (AP) lyase
Reads supporting (aligned to +/- strand):  ref base A (8/33);  new base G (4/2);  total (12/35)
Fisher's exact test for biased strand distribution p-value = 3.01e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.56e-01

TCTGATGAATAAAGCAAAACGCCTGGAGATCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGAATTGCTGATTGCCGTACTGCTTTCCGCTCAGGCGACC  >  W3110S.gb/1713233‑1713365
                                                                |                                                                    
tCTGATGAATAAAGCAAAACGCCTGGAGATCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCaccacc                                                                 >  1:3676162/1‑70 (MQ=255)
       aaTAAAGCAAAACGCCTGGAGATCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTa                                                          <  1:2255979/70‑1 (MQ=255)
       aaTAAAGCAAAACGCCTGGAGATCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTa                                                          <  1:4026602/70‑1 (MQ=255)
        aTAAAGCAAAACGCCTGGAGATCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTaa                                                         <  1:2476468/70‑1 (MQ=255)
        aTAAAGCAAAACGCCTGGAGATCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTa                                                          <  1:3678711/69‑1 (MQ=255)
            aGCAAAACGCCTGGAGATCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCAGCTTAATTTCAg                                                   <  1:1891284/71‑1 (MQ=255)
              cAAAACGCCTGGAGATCCTCACTCGCCTTCGTGAGAACAATCCTCATCCCACCACCGa                                                               <  1:3501180/58‑1 (MQ=255)
              cAAAACGCCTGGAGATCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGa                                                               <  1:561929/58‑1 (MQ=255)
                      cTGGAGATCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCtt                                           <  1:3976402/70‑1 (MQ=255)
                         gagaTCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTg                                        <  1:1751006/70‑1 (MQ=255)
                            aTCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTg                                        <  1:109830/67‑1 (MQ=255)
                            aTCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTg                                        <  1:3705180/67‑1 (MQ=255)
                              ccTCACTCGCCTGCGTGAGAACAATCCTCATCCCGCCACCGAGCTTAATTTCAGTTCGCCTTTTGAATTGc                                  >  1:419510/1‑71 (MQ=255)
                              ccTCACTCGCCTGCGTGAGAACAATCCTCATCCCGCCACCGAGCTTAATTTCAGTTCGCCTTTTGAATTGc                                  >  1:2668684/1‑71 (MQ=255)
                               cTCACTCGCCTGCGTGAGAACAATCCTCATCCcaccac                                                                  >  1:2476669/1‑38 (MQ=255)
                               cTCACTCGCCTGCGTGAGAACAATCCTCATCCcaccac                                                                  >  1:2589436/1‑38 (MQ=255)
                               cTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGa                                                               <  1:2320926/41‑1 (MQ=255)
                               cTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGAATTGCt                                 <  1:1780330/71‑1 (MQ=255)
                               cTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGAATTGCt                                 <  1:2764197/71‑1 (MQ=255)
                               cTCACTCGCCTGAGTGAGAACAATCCTCATCCCACCACCGa                                                               <  1:1090553/41‑1 (MQ=255)
                                tCACTCGCCTGCGTGAGAACAATCCTCATCCCGCCACCGAGCTTAATTTCAGTTCGCCTTTTGAATTGCt                                 >  1:467531/1‑70 (MQ=255)
                                  aCTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGAATTGCTGAt                              >  1:3525188/1‑71 (MQ=255)
                                   cTCGCCTGCGTGAGAACAATCCTCATCCCGCCACCGAGCTTAAttt                                                      <  1:1291921/46‑1 (MQ=255)
                                   cTCGCCTGCGTGAGAACAATCCTCATCCCACCACCTAGCtt                                                           <  1:1511741/41‑1 (MQ=255)
                                   cTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCtt                                                           <  1:3469773/41‑1 (MQ=255)
                                   cTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCtt                                                           <  1:2875122/41‑1 (MQ=255)
                                    tCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCtt                                                           <  1:3571200/40‑1 (MQ=255)
                                     cGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAg                                                   >  1:2525327/1‑47 (MQ=255)
                                     cGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGaa                                      >  1:272657/1‑60 (MQ=255)
                                     cGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGaa                                      >  1:4080017/1‑60 (MQ=255)
                                        cTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGa                                       <  1:302632/56‑1 (MQ=255)
                                        cTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGa                                       <  1:2115548/56‑1 (MQ=255)
                                        cTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGAATTGCTGATTGCCg                         <  1:3579549/70‑1 (MQ=255)
                                        cTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGAATTGCTGATTGCCGt                        <  1:1247185/71‑1 (MQ=255)
                                        cTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGAATTGCTGATTGCCGt                        <  1:1723052/71‑1 (MQ=255)
                                         tGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGAATTGCt                                 <  1:498676/61‑1 (MQ=255)
                                         tGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGAATTGCt                                 <  1:720895/61‑1 (MQ=255)
                                         tGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGAATTGCt                                 <  1:1295955/61‑1 (MQ=255)
                                           cGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGAAt                                     >  1:190576/1‑55 (MQ=255)
                                                 aaCAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTg                                        <  1:2594646/46‑1 (MQ=255)
                                                 aaCAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTg                                        <  1:1831055/46‑1 (MQ=255)
                                                 aaCAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTg                                        <  1:1405532/46‑1 (MQ=255)
                                                 aaCAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTg                                        <  1:477791/46‑1 (MQ=255)
                                                 aaCAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTg                                        <  1:1250244/46‑1 (MQ=255)
                                                 aaCAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGAATTGCTGAt                              >  1:2423427/1‑56 (MQ=255)
                                                 aaCAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGAATTGCTGAt                              >  1:1008241/1‑56 (MQ=255)
                                                        cTCATCCCGCCACCGAGCTTAATTTCAGTTCGCCTTTTg                                        <  1:150333/39‑1 (MQ=255)
                                                          cATCCCGCCACCGAGCTTAATTTCAGTTCGCCTTTTGaa                                      >  1:2758271/1‑39 (MQ=255)
                                                          cATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGaa                                      >  1:3142999/1‑39 (MQ=255)
                                                            tCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGAATTGCTGATTGCCGt                        <  1:575875/51‑1 (MQ=255)
                                                              ccaccaCCGAGCTTAATTTCAGTTCGCCTTTTGAATTGCTGATTGCCGTACTGCTTTCCGCTCAGGCGAcc  >  1:1012810/1‑71 (MQ=255)
                                                              ccaccaCCGAGCTTAATTTCAGTTCGCCTTTTGAATTGCTGATTGCCGTACTGCTTTCCGCTCAGGCGAcc  <  1:654633/71‑1 (MQ=255)
                                                                |                                                                    
TCTGATGAATAAAGCAAAACGCCTGGAGATCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGAATTGCTGATTGCCGTACTGCTTTCCGCTCAGGCGACC  >  W3110S.gb/1713233‑1713365

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: