Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 415,058 T→C 9.8% E373E (GAA→GAG sbcD ← exonuclease, dsDNA, ATP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb415,0580TC 9.8% 160.8 / 8.9 61E373E (GAA→GAGsbcDexonuclease, dsDNA, ATP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (30/25);  new base C (3/3);  total (33/28)
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.68e-01

CGGCGAGGGTATGCAACGTCGTGGTGAAAAGATGCTGCAGACGTTGCTGCTGCGATTCATCCAGTTCTTCCAGTGCCAGACGGCGATTGAACACCTCTTCGACGCTGAGTTCGCTGAGGGTTTCACGCTGTTGGCTG  >  W3110S.gb/414991‑415127
                                                                   |                                                                     
cGGCGAGGGTATGCAACGTCGTGGTGAAAAGATGCTGCAGACGTTGCTGCTGCGATTCATCCAGttct                                                                       <  1:249266/68‑1 (MQ=255)
    gAGGGTATGCAACGTCGTGGTGAAAAGATGCTGCAGACGTTGCTGCTGCGATTCATCCAGTTCTTcc                                                                    <  1:2092339/67‑1 (MQ=255)
        gTATGCAACGTCGTGGTGAAAAGATGCTGCAGACGTTGCTGCTGCGATTCATCCAGTTCTTCCAGTGc                                                               >  1:776716/1‑68 (MQ=255)
           tGCAACGTCGTGGTGAAAAGATGCTGCAGACGTTGCTGCTGCGATTCATCCAGttcttc                                                                     <  1:4280059/59‑1 (MQ=255)
           tGCAACGTCGTGGTGAAAAGATGCTGCAGACGTTGCTGCTGCGATTCATCCAGTTCTTCCAGTGCCAGAc                                                          <  1:1777543/70‑1 (MQ=255)
           tGCAACGTCGTGGTGAAAAGATGCTGCAGACGTTGCTGCTGCGATTCATCCAGTTCTTCCAGTGCCAGACg                                                         >  1:255916/1‑71 (MQ=255)
               aCGTCGTGGTGAAAAGATGCTGCAGACGTTGCTGCTGCGATTCATCCAGttcttc                                                                     <  1:1717493/55‑1 (MQ=255)
               aCGTCGTGGTGAAAAGATGCTGCAGACGTTGCTGCTGCGATTCATCCAGTTCTTCCAGTGCCa                                                             <  1:311108/63‑1 (MQ=255)
               aCGTCGTGGTGAAAAGATGCTGCAGACGTTGCTGCTGCGATTCATCCAGTTCTTCCAGTGCCAGACGGCGa                                                     >  1:1031741/1‑71 (MQ=255)
               aCGTCGTGGTGAAAAGATGCTGCAGACGTTGCTGCTGCGATTCATCCAGTTCTTCCAGTGCCAGACGGCGa                                                     >  1:2295547/1‑71 (MQ=255)
               aCGTCGTGGTGAAAAGATGCTGCAGACGTTGCTGCTGCGATTCATCCAGTTCTTCCAGTGCCAGACGGCGa                                                     >  1:3737196/1‑71 (MQ=255)
               aCGTCGTGGTGAAAAGATGCTGCAGACGTTGCTGCTGCGATTCATCCAGTTCCTCCAGTGc                                                               <  1:2117386/61‑1 (MQ=255)
                cgtcgtGGTGAAAAGATGCTGCAGACGTTGCTGCTGCGATTCATCCAGTTCTTCCAGTGc                                                               <  1:569225/60‑1 (MQ=255)
                 gtcgtGGTGAAAAGATGCTGCAGACGTTGCTGCTGCGATTCATCCAGTTCTTCCAGTGCCAGa                                                           >  1:900490/1‑63 (MQ=255)
                  tcgtGGTGAAAAGATGCTGCAGACGTTGCTGCTGCGATTCATCCAGTTCTTCCAGTGCCAGACGGCGATTg                                                  <  1:1186217/71‑1 (MQ=255)
                  tcgtGGTGAAAAGATGCTGCAGACGTTGCTGCTGCGATTCATCCAGTTCTTCCAGTGCCAGACGGCGATTg                                                  <  1:387370/71‑1 (MQ=255)
                    gtggtgAAAAGATGCTGCAGACGTTGCTGCTGCGATTCATCCAGTTCTTCCAGTGCCAGACGGCGATTGaa                                                >  1:2688962/1‑71 (MQ=255)
                     tggtgAAAAGATGCTGCAGACGTTGCTGCTGCGATTCATCCAGTTCCTCCAGTGCCAGACg                                                         <  1:2874322/61‑1 (MQ=255)
                       gtgAAAAGATGCTGCAGACGTTGCTGCTGCGATTCATCCAGTTCTTCCAGTGCCAGACGGCGAtt                                                   >  1:3531855/1‑65 (MQ=255)
                       gtgAAAAGATGCTGCAGACGTTGCTGCTGCGATTCATCCAGTTCTTCCAGTGCCAGACGGCGATTGAaca                                              <  1:2402828/70‑1 (MQ=255)
                       gtgAAAAGATGCTGCAGACGTTGCTGCTGCGATTCATCCAGTTCCTCCAGTGCCAGACGGc                                                       >  1:2735893/1‑61 (MQ=255)
                        tgAAAAGATGCTGCAGACGTTGCTGCTGCGATTCATCCAGttcttc                                                                     <  1:1054423/46‑1 (MQ=255)
                         gAAAAGATGCTGCAGACGTTGCTGCTGCGATTCATCCAGTTCTTCCAg                                                                  >  1:3630111/1‑48 (MQ=255)
                         gAAAAGATGCTGCAGACGTTGCTGCTGCGATTCATCCAGTTCTTCCAGt                                                                 >  1:1062339/1‑49 (MQ=255)
                         gAAAAGATGCTGCAGACGTTGCTGCTGCGATTCATCCAGTTCTTCCAGt                                                                 >  1:2002691/1‑49 (MQ=255)
                          aaaaGATGCTGCAGACGTTGCTGCTGCGATTCATCCAGTTCTTCCAGTGCCAGACGGCGATTGAACACctc                                          <  1:1467354/71‑1 (MQ=255)
                            aaGATGCTGCAGACGTTGCTGCTGCGATTCATCCAGttct                                                                       <  1:2244569/40‑1 (MQ=255)
                            aaGATGCTGCAGACGTTGCTGCTGCGATTCATCCAGttct                                                                       >  1:3108698/1‑40 (MQ=255)
                            aaGATGCTGCAGACGTTGCTGCTGCGATTCATCCAGTTCTTCCAGTGCCAGACGGCGATTGaa                                                >  1:309334/1‑63 (MQ=255)
                            aaGATGCTGCAGACGTTGCTGCTGCGATTCATCCAGTTCTTCCAGTGCCAGACGGCGATTGaa                                                >  1:303655/1‑63 (MQ=255)
                             aGATGCTGCAGACGTTGCTGCTGCGATTCATCCAGTTCTTCCAGTGCCAGACGGCGATTGAACACCTCTTc                                       <  1:3062252/71‑1 (MQ=255)
                               aTGCTGCAGACGTTGCTGCTGCGATTCATCCAGttcttc                                                                     <  1:233608/39‑1 (MQ=255)
                               aTGCTGCAGACGTTGCTGCTGCGATTCATCCAGTTCTTCCAGTGCCAGACGGCGATTGaa                                                >  1:1938098/1‑60 (MQ=255)
                               aTGCTGCAGACGTTGCTGCTGCGATTCATCCAGTTCCTCCAGTGCCAGACGGCGATTGAACACCTCTTCGa                                     >  1:4245890/1‑71 (MQ=255)
                               aTGCTGCAGACGTTGCTGCTGCGATTCATCCAGTTCCTCCAGTGCCAGACGGCGATTGAACACCTCTTCGa                                     >  1:4251066/1‑71 (MQ=255)
                                 gctgcAGACGTTGCTGCTGCGATTCATCCAGTTCTTCCAGTGCCAGACGGCGa                                                     >  1:2482205/1‑53 (MQ=255)
                                 gctgcAGACGTTGCTGCTGCGATTCATCCAGTTCTTCCAGTGCCAGACGGCGATTg                                                  <  1:4591081/56‑1 (MQ=255)
                                    gcAGACGTTGCTGCTGCGATTCATCCAGTTCTTCCAGTGCCAGACGGCGATTGAACACCTCTTCGACGCt                                 >  1:909082/1‑70 (MQ=255)
                                      aGACGTTGCTGCTGCGATTCATCCAGTTCTTCCAGTGCCa                                                             <  1:1954865/40‑1 (MQ=255)
                                      aGACGTTGCTGCTGCGATTCATCCAGTTCTTCCAGTGCCAGACGGCGATTGAACACCTCTTCGACGCTGAg                              <  1:2123259/71‑1 (MQ=255)
                                      aGACGTTGCTGCTGCGATTCATCCAGTTCTTCCAGTGCCAGACGGCGATTGAACACCTCTTCGACGCTGAg                              <  1:4146320/71‑1 (MQ=255)
                                        aCGTTGCTGCTGCGATTCATCCAGTTCTTCCAGTGCCAg                                                            <  1:2446301/39‑1 (MQ=255)
                                         cGTTGCTGCTGCGATTCATCCAGTTCTTCCAGTGCCAGACGGCGATTGAaca                                              >  1:4035792/1‑52 (MQ=255)
                                         cGTTGCTGCTGCGATTCATCCAGTTCTTCCAGTGCCAGACGGCGATTGAaca                                              >  1:3131069/1‑52 (MQ=255)
                                          gTTGCTGCTGCGATTCATCCAGTTCTTCCAGTGCCAGACGGCGATTGAACACCTCTTCGACGCTGAGTTCg                          <  1:262301/71‑1 (MQ=255)
                                          gTTGCTGCTGCGATTCATCCAGTTCTTCCAGTGCCAGACGGCGATTGAACACCTCTTCGACGCTGAGTTCg                          >  1:3365155/1‑71 (MQ=255)
                                             gctgctgcGATTCATCCAGTTCTTCCAGTGCCAGACg                                                         >  1:2384021/1‑37 (MQ=255)
                                             gctgctgcGATTCATCCAGTTCTTCCAGTGCCAGACGGc                                                       <  1:1036028/39‑1 (MQ=255)
                                               tgctgcGATTCATCCAGTTCTTCCAGTGCCAGACGGCGa                                                     >  1:648007/1‑39 (MQ=255)
                                               tgctgcGATTCATCCAGTTCTTCCAGTGCCAGACGGCGAttaa                                                 >  1:4094667/1‑41 (MQ=255)
                                               tgctgcGATTCATCCAGTTCTTCCAGTGCCAGACGGCGATTGAACACCTCTTc                                       <  1:931007/53‑1 (MQ=255)
                                               tgctgcGATTCATCCAGTTCTTCCAGTGCCAGACGGCGATTCAACAcc                                            <  1:591837/48‑1 (MQ=255)
                                                 ctgcGATTCATCCAGTTCTTCCAGTGCCAGACGGCGATTGAACACCTCTTCGACGCTGAGTTCGCTGAggg                   >  1:968297/1‑71 (MQ=255)
                                                   gcGATTCATCCAGTTCTTCCAGTGCCAGACGGCGATTGAACACCTCTTCGACGCTGAGTTCGCTg                       <  1:909560/65‑1 (MQ=255)
                                                    cGATTCATCCAGTTCTTCCAGTGCCAGACGGCGATTGAACACCTCTTCGACGCTGAGTTCGCTGAGGGttt                >  1:3819700/1‑71 (MQ=255)
                                                    cGATTAATCCAGTTCTTCCAGTGCCAGACGGCGATTGAACACCTCTTCGACGCTGAGTTCGCTGAGGGttt                >  1:1265156/1‑71 (MQ=255)
                                                         cATCCAGTTCTTCCAGTGCCAGACGGCGATTGAACACCTCTTCGACGCTGAGTTCGCTGAGGGTTTCACGc           >  1:4360699/1‑71 (MQ=255)
                                                          aTCCAGTTCTTCCGGTGCCAGACGGCGATTGAACACCTCTTCGACGCTGAGTTCGCTGAGGGTTTCACGCt          <  1:424811/71‑1 (MQ=255)
                                                           tCCAGTTCTTCCAGTGCCAGACGGCGATTGAACACCTCTTc                                       <  1:3380044/41‑1 (MQ=255)
                                                            ccAGTTCTTCCAGTGCCAGACGGCGATTGAACACCTCtt                                        >  1:1858587/1‑39 (MQ=255)
                                                            ccAGTTCTTCCAGTGCCAGACGGCGATTGAACACCTCTTc                                       >  1:3208731/1‑40 (MQ=255)
                                                            ccAGTTCTTCCAGTGCCAGACGGCGATTGAACACCTCTTCGa                                     >  1:2413127/1‑42 (MQ=255)
                                                             cAGTTCTTCCAGTGCCAGACGGCGATTGAACACCTCTTCGACGCTGAGTTCGCTGAGGGttt                >  1:699762/1‑62 (MQ=255)
                                                             cAGTTCTTCCAGTGCCAGACGGCGATTGAACACCTCTTCGACGCTGAGTTCGCTGAGGGTTTCACGCTGtt       <  1:4592364/71‑1 (MQ=255)
                                                             cAGTTCCTCCAGTGCCAGACGGCGATTGAACACCTCTTCGACGCTGAGTTCGCTGAGGGTTTCACGCTGtt       <  1:3109484/71‑1 (MQ=255)
                                                              aGTTCTTCCAGTGCCAGACGGCGATTGAACACCTCTTCGACGCTGAGTTCGCTGAGTGTTTCACGCTGTtg      <  1:2359014/71‑1 (MQ=255)
                                                               gTTCTTCCAGTGCCAGACGGCGATTGAACACCTCTTc                                       <  1:1269452/37‑1 (MQ=255)
                                                               gTTCTTCCAGTGCCAGACGGCGATTGAACACCTCTTCGACGCTGAGTTCGCTGAGGGTTTCACGCTGTtgg     <  1:781603/71‑1 (MQ=255)
                                                                ttcttcCAGTGCCAGACGGCGATTGAACACCTCTTCGACGCTGAGTTCGCTGAGGGTTTCACg            <  1:628010/63‑1 (MQ=255)
                                                                ttcctcCAGTGCCAGACGGCGATTGAACACCTCTTc                                       <  1:1791653/36‑1 (MQ=255)
                                                                 tcttcCAGTGCCAGACGGCGATTGAACACCTCTTCGACGCTGAg                              >  1:2187677/1‑44 (MQ=255)
                                                                  cttcCAGTGCCAGACGGCGATTGAACACCTCTTCGa                                     >  1:2410877/1‑36 (MQ=255)
                                                                  cttcCAGTGCCAGACGGCGATTGAACACCTCTTCGACGc                                  >  1:3404925/1‑39 (MQ=255)
                                                                  cttcCAGTGCCAGACGGCGATTGAACACCTCTTCGACGCTGAGTTCGCTGa                      >  1:2149596/1‑51 (MQ=255)
                                                                  cttcCAGTGCCAGACGGCGATTGAACACCTCTTCGACGCTGAGTTCGCTGAGGGTTTCa              >  1:390727/1‑59 (MQ=255)
                                                                  cttcCAGTGCCAGACGGCGATTGAACACCTCTTCGACGCTGAGTTCGCTGAGGGTTTCAc             >  1:3594341/1‑60 (MQ=255)
                                                                  cttcCAGTGCCAGACGGCGATTGAACACCTCTTCGACGCTGAGTTCGCTGAGGGTTTCACGCTGTtggctg  >  1:4583055/1‑71 (MQ=255)
                                                                  cttcCAGTGCCAGACGGCGATTGAACACCTCTTCGACGCTGAGTTCGCTGAGGGTTTCACGCTGTtggct   <  1:1330875/70‑1 (MQ=255)
                                                                  cttcCAGTGCCAGACGGCGATTGAACACCTCTTCGACGCTGAGTTCGCTGAGGGTTTCACGCTGTtgg     >  1:192486/1‑68 (MQ=255)
                                                                  cttcCAGTGCCAGACGGCGATTGAACACCACTTCGACGCTGAGTTCGCTGAGGGTTTCACGCTGTtggctg  >  1:547791/1‑71 (MQ=255)
                                                                   |                                                                     
CGGCGAGGGTATGCAACGTCGTGGTGAAAAGATGCTGCAGACGTTGCTGCTGCGATTCATCCAGTTCTTCCAGTGCCAGACGGCGATTGAACACCTCTTCGACGCTGAGTTCGCTGAGGGTTTCACGCTGTTGGCTG  >  W3110S.gb/414991‑415127

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: