Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,035,613 G→A 9.1% P384S (CCG→TCG)  recJ ← ssDNA exonuclease, 5' ‑‑> 3'‑specific

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,035,6130GA 9.1% 122.4 / 2.7 44P384S (CCG→TCG) recJssDNA exonuclease, 5' ‑‑> 3'‑specific
Reads supporting (aligned to +/- strand):  ref base G (19/21);  new base A (2/2);  total (21/23)
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.30e-01

CCCCTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATA  >  W3110S.gb/3035548‑3035682
                                                                 |                                                                     
ccccTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGt                                                                  <  1:2381457/71‑1 (MQ=255)
   cTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGaa                                                               <  1:1136361/71‑1 (MQ=255)
   cTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGaa                                                               <  1:867043/71‑1 (MQ=255)
   cTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGaa                                                               <  1:2220617/71‑1 (MQ=255)
     gAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGaaaa                                                             >  1:3078817/1‑71 (MQ=255)
     gAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGaaaa                                                             >  1:3667987/1‑71 (MQ=255)
     gAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGaaaa                                                             >  1:175440/1‑71 (MQ=255)
             cGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCttt                                                     <  1:2115375/71‑1 (MQ=255)
             cGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCttt                                                     <  1:410297/71‑1 (MQ=255)
             cGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCGCttt                                                     <  1:1474792/71‑1 (MQ=255)
              gACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTg                                                    <  1:258184/71‑1 (MQ=255)
               aCCGGAACCTTACAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGa                                                   <  1:2302241/71‑1 (MQ=255)
                ccGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGAt                                                  <  1:3242865/71‑1 (MQ=255)
                    aaCCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATgcgc                                              >  1:2648819/1‑71 (MQ=255)
                        tttCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAg                                          >  1:3160700/1‑71 (MQ=255)
                          tCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGc                                         >  1:3799758/1‑70 (MQ=255)
                            aGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAg                                      <  1:64385/71‑1 (MQ=255)
                                  ccGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGctc                                                        <  1:2588222/47‑1 (MQ=255)
                                  ccGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATAcc                                <  1:2997074/71‑1 (MQ=255)
                                  ccGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATAc                                 >  1:780643/1‑70 (MQ=255)
                                   cGTCACCTGCTGGCGCAAACGCGATAACCGAACGGTGAAAACGCTCTTTGATGCGCGAAGCCAg                                      <  1:939814/64‑1 (MQ=255)
                                   cGTCACCTGCTGGCGCAAACGCGATAACCGAACGGTGAAAACGCTCTTTGATGCGCGAAGCCAg                                      <  1:2304177/64‑1 (MQ=255)
                                    gTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCttt                                                     <  1:3160654/48‑1 (MQ=255)
                                     tCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGAc                             <  1:935521/71‑1 (MQ=255)
                                     tCACCTGCTGGCGCAAACGCGATAACCGAACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGaa                                    >  1:3856372/1‑64 (MQ=255)
                                     tCACCTGCTGGCGCAAACGCGATAACCGAACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGaa                                    >  1:3852201/1‑64 (MQ=255)
                                       aCCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACaa                           >  1:463/1‑71 (MQ=255)
                                         ctgctgGCGCAAACGCGATAACCGGACGGTGAAAACGCTCttt                                                     <  1:3699711/43‑1 (MQ=255)
                                         ctgctgGCGCAAACGCGATAACCGGACGGTGAAAACGCTCttt                                                     <  1:486394/43‑1 (MQ=255)
                                            ctgGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGccc                      >  1:374937/1‑71 (MQ=255)
                                            ctgGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGccc                      >  1:3192384/1‑71 (MQ=255)
                                               gcgcAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCt                     <  1:3307577/69‑1 (MQ=255)
                                               gcgcAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCt                     <  1:2591772/69‑1 (MQ=255)
                                                 gcAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATg                 >  1:2304006/1‑71 (MQ=255)
                                                    aaCGCGATAACCGGACGGTGAAAACGCTCTTTGATgcgc                                              >  1:2407613/1‑39 (MQ=255)
                                                    aaCGCGATAACCGGACGGTGAAAACGCTCTTTGATgcgc                                              >  1:3757905/1‑39 (MQ=255)
                                                     aCGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGc                                         <  1:366316/43‑1 (MQ=255)
                                                     aCGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCa              >  1:2085952/1‑70 (MQ=255)
                                                      cgcgATAAGCGGACGGTGAAAACGCTCTTTGATgc                                                >  1:2509786/1‑35 (MQ=255)
                                                      cgcgATAACCGGACGGTGAAAACGCTCTTTGATgc                                                >  1:1161457/1‑35 (MQ=255)
                                                      cgcgATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCAtt            >  1:3861490/1‑71 (MQ=255)
                                                           tAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACg                         >  1:3782788/1‑53 (MQ=255)
                                                           tAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGcc                       >  1:3533129/1‑55 (MQ=255)
                                                              ccGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGa                                     <  1:472317/38‑1 (MQ=255)
                                                                ggACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGaa                                    >  1:1486090/1‑37 (MQ=255)
                                                                ggACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATa  <  1:3054691/71‑1 (MQ=255)
                                                                ggACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATa  <  1:1483746/71‑1 (MQ=255)
                                                                ggACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATa  <  1:3082141/71‑1 (MQ=255)
                                                                 gACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATa  <  1:3757565/70‑1 (MQ=255)
                                                                 |                                                                     
CCCCTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATA  >  W3110S.gb/3035548‑3035682

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: