Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 412,184 G→A 14.3% L932L (CTG→TTG)  sbcC ← exonuclease, dsDNA, ATP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb412,1840GA14.3% 128.3 / 11.8 49L932L (CTG→TTG) sbcCexonuclease, dsDNA, ATP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (14/28);  new base A (3/4);  total (17/32)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

TCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTCGCTGGCTTTGCGCTGTAACAGATAGCGCCCGTGCAGCCGGGTAAGTTGCTGATTAGCGAG  >  W3110S.gb/412114‑412250
                                                                      |                                                                  
tCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAg                                                                    >  1:530423/1‑71 (MQ=255)
     gccgGAAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCt                                                               <  1:1945120/71‑1 (MQ=255)
      ccgGAAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTc                                                              <  1:3581738/71‑1 (MQ=255)
                 aCGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAACGCCTCGCTGGCTTTgc                                                   >  1:3635742/1‑71 (MQ=255)
                      tATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTCGCTGGCTTTGCGCTGt                                              <  1:2837959/71‑1 (MQ=255)
                        tCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTCGCTGGCTTTGCGCTGTaa                                            >  1:3782251/1‑71 (MQ=255)
                          gcgTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTCGCTGGCTTTGCGCTGTAACa                                          <  1:2903194/71‑1 (MQ=255)
                          gcgTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTCGCTGGCTTTGCGCTGTAACa                                          <  1:2672767/71‑1 (MQ=255)
                          gcgTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTCGCTGGCTTTGCGCTGTAACa                                          <  1:4030558/71‑1 (MQ=255)
                          gcgTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTCGCTGGCTTTGCGCTGTAACa                                          <  1:3877044/71‑1 (MQ=255)
                          gcgTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTCGCTGGCTTTGCGCTGTAACa                                          <  1:2257696/71‑1 (MQ=255)
                           cgTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAACGCCTCGCTGGCTTTGCTCTGTAACAg                                         <  1:381773/71‑1 (MQ=255)
                           cgTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAACGCCTCGCTGGCTTTGCGCTGTAACAg                                         <  1:455535/71‑1 (MQ=255)
                             tACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTCGCTGGCTTTGCGCTGTAACa                                          <  1:1933599/68‑1 (MQ=255)
                             tACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTCGCTGGCTTTGCGCTGTAACa                                          <  1:2493423/68‑1 (MQ=255)
                              aCCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTCGCTGGCTTTGCGCTGTAACAg                                         >  1:2532529/1‑68 (MQ=255)
                               ccGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTCGCTGGCTTTGCGCTGTAACAGa                                        >  1:2063462/1‑68 (MQ=255)
                               ccGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTCGCTGGCTTTGCGCTGTAACAGATAg                                     >  1:836090/1‑71 (MQ=255)
                                    tCTGCCTGCCATGTATCAACAACCTCGACTTCCAGCGCCTCGCTGGCTTTGCGCTGTAACa                                          <  1:3043864/61‑1 (MQ=255)
                                    tCTGCCTGCCAGGTATCAACAACCTCGACTTCCAACGCCTCGCTGGCTTTGCGCTGTAACAGATAGCGccc                                >  1:3828077/1‑71 (MQ=255)
                                      tgcctgccAGGTATCAACAACCTCGACTTCCAGCGCCTCGCTg                                                          >  1:3986167/1‑43 (MQ=255)
                                      tgcctgccAGGTATCAACAACCTCGACTTCCAGCGCCTCGCTg                                                          >  1:2867310/1‑43 (MQ=255)
                                      tgcctgccAGGTATCAACAACCTCGACTTCCAGCGCCTCGCTGGCTTTGCGCTGTAACAGAt                                       <  1:2275506/62‑1 (MQ=255)
                                      tgcctgccAGGTATCAACAACCTCGACTTCCAGCGCCTCGCTGGCTTTGCGCTGTAACAGAt                                       <  1:1150329/62‑1 (MQ=255)
                                      tgcctgccAGGTATCAACAACCTCGACTTCCAGCGCCTCGCTGGCTTTGCGCTGTAACAGATAGCGCCCGt                              >  1:1226934/1‑71 (MQ=255)
                                      tgcctgccAGGTATCAACAACCTCGACTTCCAGCGACTCGCTg                                                          >  1:1266612/1‑43 (MQ=255)
                                             cAGGTATCAACAACCTCGACTTCCAACGCCTCGCTGGCTTTgcgc                                                 <  1:3220437/45‑1 (MQ=255)
                                                 tATCAACAACCTCGACTTCCAGCGCCTCGCTGGCTTTGCGCTGTAACa                                          <  1:913498/48‑1 (MQ=255)
                                                 tATCAACAACCTCGACTTCCAGCGCCTCGCTGGCTTTGCGCTGTAACAGATAGCGCCCGTGCAGCCGGGTa                   <  1:955210/71‑1 (MQ=255)
                                                 tATCAACAACCTCGACTTCCAGCGCCTCGCTGGCTTTGCGCTGTAACAGATAGCGCCCGTGCAGCCGGGTa                   <  1:1091736/71‑1 (MQ=255)
                                                 tATCAACAACCTCGACTTCCAACGCCTCGCTGGCTTTgcgc                                                 <  1:1991284/41‑1 (MQ=255)
                                                  aTCAACAACCTCGACTTCCAGCGCCTCGCTGGCTTTGCGCTGTAACAg                                         >  1:1629005/1‑48 (MQ=255)
                                                  aTCAACAACCTCGACTTCCAGCGCCTCGCTGGCTTTGCGCTGTAACAGa                                        >  1:3714744/1‑49 (MQ=255)
                                                   tCAACAACCTCGACTTCCAGCGCCTCGCTGGCTTTGCGCTGTAACa                                          <  1:349086/46‑1 (MQ=255)
                                                   tCAACAACCTCGACTTCCAGCGCCTCGCTGGCTTTGCGCTGTAACa                                          <  1:453252/46‑1 (MQ=255)
                                                   tCAACAACCTCGACTTCCAGCGCCTCGCTGGCTTTGCGCTGTAACAGAt                                       >  1:2910583/1‑49 (MQ=255)
                                                   tCAACAACCTCGACTTCCAGCGCCTCGCTGGCTTTGCGCTGTAACAGAt                                       >  1:2536471/1‑49 (MQ=255)
                                                   tCAACAACCTCGACTTCCAGCGCCTCGCTGGCTTTGCGCTGTAACAGATAGCGCCCGTGCa                           <  1:12459/61‑1 (MQ=255)
                                                   tCAACAACCTCGACTTCCAGCGCCTCGCTGGCTTTGCGCTGTAACAGATAGCGCCCGTGCa                           <  1:3802972/61‑1 (MQ=255)
                                                   tCAACAACCTCGACTTCCAGCGCCTCGCTGGCTTTGCGCTGTAACAGATAGCGCCCGTGCAGCCGGGTaa                  <  1:2007508/70‑1 (MQ=255)
                                                     aacaacCTCGACTTCCAGCGCCTCGCTGGCTTTGCGCTGTaa                                            >  1:1017947/1‑42 (MQ=255)
                                                        aacCTCGACTTCCAGCGCCTCGCTGGCTTTGCGCTGTAACAGATAGCGCCCGTGCAGCCGGGTAAGTTGCt            <  1:3255533/71‑1 (MQ=255)
                                                         acCTCGACTTCCAGCGCCTCGCTGGCTTTGCGCTGTAACa                                          <  1:2386029/40‑1 (MQ=255)
                                                         acCTCGACTTCCAGCGCCTCGCTGGCTTTGCGCTGTAACAGATAGCGCCCGTGCa                           >  1:3077610/1‑55 (MQ=255)
                                                           cTCGACTTCCAACGCCTCGCTGGCTTTGCGCTGTAACAGATAGCGCCCGTGCAGCCGGGTAAGTTGCTGAt         >  1:1697417/1‑71 (MQ=255)
                                                               aCTTCCAGCGCCTCGCTGGCTTTGCGCTGTAACAGATAGCGCCCGTGCAGCCGGGTAAGTTGCt            <  1:2140802/64‑1 (MQ=255)
                                                               aCTTCCAGCGCCTCGCTGGCTTTGCGCTGTAACAGATAGCGCCCGTGCAGCAGGGTAAGTTGCt            <  1:1219001/64‑1 (MQ=255)
                                                                   ccAGCGCCTCGCTGGCTTTGCGCTGTAACAGATAGCGCCCGTGCAGCCGGGTAAGTTGCTGATTAGCgag  <  1:3911952/70‑1 (MQ=255)
                                                                   ccAGCGCCTCGCTGGCTTTGCGCTGTAACAGATAGCGCCCGTGCAGCCGGGTAAGTTGCTGATTAGCgag  <  1:2926550/70‑1 (MQ=255)
                                                                     aGCGCCTCGCTGGCTTTGCGCTGTAACAGATAGCGCCCGTGCAGCCGGGTAAGTTGCt            <  1:4030821/58‑1 (MQ=255)
                                                                      |                                                                  
TCGCCGCCGGAAAGGGTACGGGTATCGCGTACCGCATCTGCCTGCCAGGTATCAACAACCTCGACTTCCAGCGCCTCGCTGGCTTTGCGCTGTAACAGATAGCGCCCGTGCAGCCGGGTAAGTTGCTGATTAGCGAG  >  W3110S.gb/412114‑412250

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: