Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,713,297 A→G 15.4% 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,2970AG15.4% 125.8 / 15.9 52T21A (ACC→GCC) nthDNA glycosylase and apyrimidinic (AP) lyase
Reads supporting (aligned to +/- strand):  ref base A (24/20);  new base G (4/4);  total (28/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.29e-01

TCTGATGAATAAAGCAAAACGCCTGGAGATCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGAATTGCTGATTGCCGTACTGCTTTCCGCTCAGGCGACC  >  W3110S.gb/1713233‑1713365
                                                                |                                                                    
tCTGATGAATAAAGCAAAACGCCTGGAGATCCTCACTCGCCTGCGTGAGAACAATCCTCATCccac                                                                     >  1:4187436/1‑66 (MQ=255)
 cTGATGAATAAAGCAAAACGCCTGGAGATCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCaccacc                                                                 >  1:3901584/1‑69 (MQ=255)
  tgatgaATAAAGCAAAACGCCTGGAGATCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAg                                                              >  1:3625414/1‑71 (MQ=255)
    atgaATAAAGCAAAACGCCTGGAGATCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCt                                                            <  1:2480039/71‑1 (MQ=255)
     tgaATAAAGCAAAACGCCTGGAGATCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCt                                                            >  1:240823/1‑70 (MQ=255)
         tAAAGCAAAACGCCTGGAGATCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAAtt                                                       >  1:4533889/1‑71 (MQ=255)
           aaGCAAAACGCCTGGAGATCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAAttt                                                      >  1:3293050/1‑70 (MQ=255)
              cAAAACGCCTGGAGATCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAAttt                                                      <  1:3305215/67‑1 (MQ=255)
               aaaaCGCCTGGAGATCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGtt                                                 <  1:4630221/71‑1 (MQ=255)
                 aaCGCCTGGAGATCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCg                                               <  1:1712514/71‑1 (MQ=255)
                       tGGAGATCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCtttt                                         >  1:135047/1‑71 (MQ=255)
                        ggAGATCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAAttt                                                      >  1:2721366/1‑57 (MQ=255)
                         gagaTCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCGCCACCGAGCTTAATTTCAGTTCGCCTTTTg                                        >  1:2159544/1‑70 (MQ=255)
                         gagaTCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTc                                                     <  1:2184082/57‑1 (MQ=255)
                         gagaTCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTc                                                     <  1:4537692/57‑1 (MQ=255)
                         gagaTCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTgt                                       <  1:4298980/71‑2 (MQ=255)
                         gagaTCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTg                                        <  1:4384249/70‑1 (MQ=255)
                         gagaTCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTg                                        <  1:920463/70‑1 (MQ=255)
                            aTCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCGCCACCGAGCTTAATTTCAGTTCGCCTTTTGAAtt                                    <  1:4129614/71‑1 (MQ=255)
                               cTCACTCGCCTGCGTGAGAACAATCCTCATCCCGCCACCGAGCTTAATTTCAGTTCGCCTTTTGAATTGCt                                 <  1:694475/71‑1 (MQ=255)
                               cTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTg                                        <  1:4389746/64‑1 (MQ=255)
                                  aCTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGAATTGCTg                                <  1:397188/69‑1 (MQ=255)
                                   cTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGa                                                               <  1:2313742/37‑1 (MQ=255)
                                   cTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCg                                               <  1:1718003/53‑1 (MQ=255)
                                    tCGCCTGCGTGAGAACAATCCTCATCCCACCACCg                                                                >  1:3410599/1‑35 (MQ=255)
                                    tCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGAATTGCt                                 <  1:3652485/66‑1 (MQ=255)
                                      gCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTaa                                                         >  1:1355263/1‑40 (MQ=255)
                                      gCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGAATTGCTGa                               >  1:2311167/1‑66 (MQ=255)
                                      gCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGAATTGCTGATTGcc                          >  1:905341/1‑71 (MQ=255)
                                       ccTGCGTGAGAACAATCCTCATCCCGCCACCGAGCTTAATTTCAGTTCGCCTTTTGAATTGCTGATTGCCg                         <  1:4514258/71‑1 (MQ=255)
                                         tGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGAATTGCTGATTGCCGTa                       >  1:4371109/1‑71 (MQ=255)
                                           cGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGAATTGCTGATTGCCGTACt                     >  1:2162677/1‑71 (MQ=255)
                                           cGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGAATTGCTGATTGCCGTACt                     >  1:684914/1‑71 (MQ=255)
                                            gTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGAATTGCTGATTGCCGTACTg                    >  1:2305726/1‑71 (MQ=255)
                                               agaACAATCCTCATCCCGCCACCGAGCTTAATTTCAGTTCGCCTTTTGAATTGCTGATTGCCGTACTGCtt                 <  1:1920954/71‑1 (MQ=255)
                                               agaACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGc                                              >  1:3365057/1‑42 (MQ=255)
                                               agaACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTg                                        <  1:2953750/48‑1 (MQ=255)
                                               agaACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGAAt                                     >  1:1959932/1‑51 (MQ=255)
                                                gaACAATCCTCATCCCACCACCGAGCTTAATTTCAGtt                                                 >  1:1322668/1‑38 (MQ=255)
                                                gaACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGc                                              >  1:3778660/1‑41 (MQ=255)
                                                gaACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGaa                                      >  1:3452230/1‑49 (MQ=255)
                                                gaACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGAATTGCTGATTGCCGTACTGCttt                <  1:2294114/71‑1 (MQ=255)
                                                 aaCAATCCTCATCCCACCACCGAGCTTAATTTCAGtt                                                 >  1:1743094/1‑37 (MQ=255)
                                                     aTCCTCATCCCGCCACCGAGCTTAATTTCAGTTCGCCTTTTGAATTGCTGATTGCCGTACTGCTTTCCGCt           >  1:1097878/1‑71 (MQ=255)
                                                     aTCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGAAt                                     <  1:211329/45‑1 (MQ=255)
                                                     aTCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGAATTGCTGATTGCCGTACTGCTTTCCGCt           >  1:3455966/1‑71 (MQ=255)
                                                      tCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTg                                        <  1:3225344/41‑1 (MQ=255)
                                                      tCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTg                                        <  1:1353817/41‑1 (MQ=255)
                                                       ccTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGAATTg                                   >  1:3918474/1‑45 (MQ=255)
                                                         tCATCCCACCACCGAGCTTAATTTCAGTTCCCCTTTTGAATTGCt                                 <  1:2839547/45‑1 (MQ=255)
                                                          cATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGAATTGCTGa                               >  1:896474/1‑46 (MQ=255)
                                                           aTCCCGCCACCGAGCTTAATTTCAGTTCGCCTTTTg                                        >  1:392438/1‑36 (MQ=255)
                                                           aTCCCGCCACCGAGCTTAATTTCAGTTCGCCTTTTg                                        >  1:4596509/1‑36 (MQ=255)
                                                           aTCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGAATTGCTGATTGCCGTACTGCTTTCCGCTCAGGc      >  1:3640503/1‑70 (MQ=255)
                                                              ccaccaCCGAGCTTAATTTCAGTTCGCCTTTTGAATTGc                                  >  1:1108542/1‑39 (MQ=255)
                                                              ccaccaCCGAGCTTAATTTCAGTTCGCCTTTTGAATTGCTGATTGCCGTACTGCTTTCCGCTCAGGCGAcc  >  1:847265/1‑71 (MQ=255)
                                                                accaccGAGCTTAATTTCAGTTCGCCTTTTGAATTGCt                                 <  1:899471/38‑1 (MQ=255)
                                                                |                                                                    
TCTGATGAATAAAGCAAAACGCCTGGAGATCCTCACTCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGAATTGCTGATTGCCGTACTGCTTTCCGCTCAGGCGACC  >  W3110S.gb/1713233‑1713365

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: