Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,953,568 G→A 7.7% S364S (AGC→AGT recB ← exonuclease V (RecBCD complex), beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,953,5680GA 7.7% 228.8 / 6.4 78S364S (AGC→AGTrecBexonuclease V (RecBCD complex), beta subunit
Reads supporting (aligned to +/- strand):  ref base G (43/29);  new base A (3/3);  total (46/32)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

CTGAAATTCATCGATCATTGCCACCGGGAATCGCGTACGGATCGCCGCTGCCAACACCTCACCGCTTTCGCTACGCAGCGCGGAATCGAGCCGACTTAACATGTCATCAAAACCCAATTCGCCACGGCGGCGTTTTT  >  W3110S.gb/2953499‑2953635
                                                                     |                                                                   
cTGAAATTCATCGATCATTGCCACCGGGAATCGCGTACGGATCGCCGCTGCCAACACCTCACCGCTTTCGc                                                                    >  1:1300208/1‑71 (MQ=255)
   aaaTTCATCGATCATTGCCACCGGGAATCGCGTACGGATCGCCGCTGCCAACACCTCACCGCTTTCGCTAc                                                                 >  1:3115952/1‑71 (MQ=255)
     aTTCATCGATCATTGCCACCGGGAATCGCGTACGGATCGCCGCTGCCAACACCTCACCGCTTTCGCTACg                                                                >  1:2629917/1‑70 (MQ=255)
     aTTCATCGATCATTGCCACCGGGAATCGCGTACGGATCGCCGCTGCCAACACCTCACCGCTTTCGCTACg                                                                >  1:1560338/1‑70 (MQ=255)
        cATCGATCATTGCCACCGGGAATCGCGTACGGATCGCCGCTGCCAACACCTCACCGCTTTCGCTACGCAgc                                                            >  1:2753595/1‑71 (MQ=255)
        cATCGATCATTGCCACCGGGAATCGCGTACGGATCGCCGCTGCCAACACCTCACCGCTTTCGCTACGCAgc                                                            >  1:1212159/1‑71 (MQ=255)
         aTCGATCATTGCCACCGGGAATCGCGTACGGATCGCCGCTGCCAATACCTCGCCGCTTTCACTACGCAgc                                                            <  1:2626844/70‑1 (MQ=255)
            gATCATTGCCACCGGGAATCGCGTACGGATCGCCGCTGCCAACACCTCCCCGCTTTCGCTACGCAGcgcg                                                         >  1:2083044/1‑70 (MQ=255)
            gATCATTGCCACCGGGAATCGCGTACGGATCGCCGCTGCCAACACCTCACCGCTTTCGCTACGCAGcgcg                                                         >  1:2423125/1‑70 (MQ=255)
            gATCATTGCCACCGGGAATCGCGTACGGATCGCCGCTGCCAACACCTCACCGCTTTCGCTACGCAGcgcg                                                         >  1:2159143/1‑70 (MQ=255)
            gATCATTGCCACCGGGAATCGCGTACGGATCGCCGCTGCCAACACCTCACCGCTTTCGCTACGCAGCGCgg                                                        <  1:2736937/71‑1 (MQ=255)
               cATTGCCACCGGGAATCGCGTACGGATCGCCGCTGCCAACACCTCACCGCTTTCGCTACGCAGCGCGGAAt                                                     >  1:1556706/1‑71 (MQ=255)
               cATTGCCACCGGGAATCGCGTACGGATCGCCGCTGCCAACACCTCACCGCTTTCGCTACGCAGCGCGGAAt                                                     >  1:3701044/1‑71 (MQ=255)
                aTTGCCACCGGGAATCGCGTACGGATCGCCGCTGCCAACACCTCACCGCTTTCGCTACGCAGCGCGGa                                                       <  1:2993884/68‑1 (MQ=255)
                aTTGCCACCGGGAATCGCGTACGGATCGCCGCTGCCAACACCTCACCGCTTTCGCTACGCAGCGCGGAATc                                                    >  1:2177681/1‑71 (MQ=255)
                aTTGCCACCGGGAATCGCGTACGGATCGCCGCTGCCAACACCTCACCGCTTTCGCTACGCAGCGCGGAATc                                                    <  1:503982/71‑1 (MQ=255)
                  tGCCACCGGGAATCGCGTACGGATCGCCGCTGCCAACACCTCACCGCTTTCGCTACGCAGCGCGGAAt                                                     >  1:822470/1‑68 (MQ=255)
                  tGCCACCGGGAATCGCGTACGGATCGCCGCTGCCAACACCTCACCGCTTTCGCTACGCAGCGCGGAATc                                                    >  1:407915/1‑69 (MQ=255)
                  tGCCACCGGGAATCGCGTACGGATCGCCGCTGCCAACACCTCACCGCTTTCGCTACGCAGCGCGGAATc                                                    >  1:2564844/1‑69 (MQ=255)
                  tGCCACCGGGAATCGCGTACGGATCGCCGCTGCCAACACCTCACCGCTTTCGCTACGCAGCGCGGAATc                                                    >  1:3419518/1‑69 (MQ=255)
                       ccGGGAATCGCGTACGGATCGCCGCTGCCAACACCTCACCGCTTTCGCTACGCAGCGCGGAATCGAGCCGa                                             >  1:3324484/1‑71 (MQ=255)
                         gggAATCGCGTACGGATCGCCGCTGCCAACACCTCACCGCTTTCGCTACGCAGCGCGGAATCGAGCCGACt                                           >  1:3239112/1‑71 (MQ=255)
                         gggAATCGCGTACGGATCGCCGCTGCCAACACCTCACCGCTTTCGCTACGCAGCGCGGAATCGAGCCGACt                                           >  1:3243677/1‑71 (MQ=255)
                         gggAATCGCGTACGGATCGCCGCTGCCAACACCTCACCGCTTTCGCTACGCAGCGCGGAATCGAGCCGACt                                           >  1:2133479/1‑71 (MQ=255)
                          ggAATCGCGTACGGATCGCCGCTGCCAACACCTCACCGCTTTCGCTACGCAGCGCGGAATCGAGCCGACtt                                          <  1:2856445/71‑1 (MQ=255)
                           gAATCGCGTACGGATCGCCGCTGCCAACACCTCACCGCTTTCGCTACGCAGCGCGGAATCGAGCCGACTTa                                         <  1:2597860/71‑1 (MQ=255)
                           gAATCGCGTACGGATCGCCGCTGCCAACACCTCACCGCTTTCGCTACGCAGCGCGGAATCGAGCCGACTTa                                         <  1:789434/71‑1 (MQ=255)
                             aTCGCGTACGGATCGCCGCTGCCAATACCTCGCCGCTTTCACTACGCAGCGCGGAATCGAg                                                 >  1:537272/1‑61 (MQ=255)
                             aTCGCGTACGGATCGCCGCTGCCAATACCTCGCCGCTTTCACTACGCAGCGCGGAATCGAGCCGACTTaa                                        >  1:1768972/1‑70 (MQ=255)
                             aTCGCGTACGGATCGCCGCTGCCAATACCTCGCCGCTTTCACTACGCAGCGCGGAATCGAGCCGACTTa                                         >  1:3149450/1‑69 (MQ=255)
                             aTCGCGTACGGATCGCAGCTGCCAACACCTCACCGCTTTCGCTACGCAGCGCGGAATCGAGCCGACTTaa                                        >  1:3550209/1‑70 (MQ=255)
                              tCGCGTACGGATCGCCGCTGCCAACACCTCACCGCTTTCGCTACGCAGCGCgg                                                        <  1:1972288/53‑1 (MQ=255)
                              tCGCGTACGGATCGCCGCTGCCAACACCTCACCGCTTTCGCTACGCAGCGCgg                                                        <  1:2720511/53‑1 (MQ=255)
                               cgcgTACGGATCGCCGCTGCCAACACCTCACCGCTTTCGCTACGCAGCGCGGAATCGAGCCGACTTAACAt                                     >  1:2976430/1‑71 (MQ=255)
                                gcgTACGGATCGCCGCTGCCAACACCTCACCGCTTTCGCTa                                                                  >  1:2278075/1‑41 (MQ=255)
                                gcgTACGGATCGCCGCTGCCAACACCTCACCGCTTTCGCTa                                                                  >  1:2325876/1‑41 (MQ=255)
                                gcgTACGGATCGCCGCTGCCAACACCTCACCGCTTTCGCTa                                                                  >  1:442646/1‑41 (MQ=255)
                                   tACGGATCGCCGCTGCCAATACCTCGCCGCTTTCACTACGCAGCGCGGAATCGAGCCGACTTAACATGtca                                 <  1:3024148/71‑1 (MQ=255)
                                   tACGGATCGCCGCTGCCAATACCTCGCCGCTTTCACTACGCAGCGCGGAATCGAGCCGACTTAACATGtca                                 <  1:903590/71‑1 (MQ=255)
                                     cGGATCGCCGCTGCCAACACCTCACCGCTTTCGCTACGCAGCGCGGAATCGAGCCGACTTAACATGtcat                                >  1:3506040/1‑70 (MQ=255)
                                     cGGATCGCCGCTGCCAACACCTCACCGCTTTCGCTACGCAGCGCGGAATCGAGCCGACTTAACATGtcat                                >  1:3499900/1‑70 (MQ=255)
                                       gATCGCCGCTGCCAACACCTCACCGCTTTCGCTACGCAGCGCGGAATCGAGCCGACTTAACATGTCATCaa                             <  1:3017828/71‑1 (MQ=255)
                                       gATCGCCGCTGCCAACACCTCACCGCTTTCGCTACGCAGCGCGGAATCGAGCCGACTTAACATGTCATCaa                             <  1:3315415/71‑1 (MQ=255)
                                         tCGCCGCTGCCAACACCTCACCGCTTTCGCTACGCAGCGCGGAATCGAGCCGACTTAACATGTCATCaaaa                           >  1:1645351/1‑71 (MQ=255)
                                            ccgcTGCCAACACCTCACCGCTTTCGCTACGCAGCGCGGAATCGAGCCGACTTAACATGTCATCAAAAccc                        <  1:3739362/71‑1 (MQ=255)
                                            ccgcTGCCAACACCTCACCGCTTTCGCTACGCAGCGCGGAATCGAGCCGACTTAACATGTCATCAAAAccc                        <  1:187824/71‑1 (MQ=255)
                                            ccgcTGCCAACACCTCACCGCTTTCGCTACGCAGCGCGGAATCGAGCCGACTTAACATGTCATCAAAAccc                        <  1:76332/71‑1 (MQ=255)
                                              gcTGCCAACACCTCACCGCTTTCGCTACGCAGCGCGGAATCGAGCCGACTTAACATGTCATCAAAACCCaa                      <  1:584845/71‑1 (MQ=255)
                                               cTGCCAACACCTCACCGCTTTCGCTACGCAGCGCGGAATCGAGCCGACTTAACATGTCATCAAAACCCAAt                     <  1:2712990/71‑1 (MQ=255)
                                               cTGCCAACACCTCACCGCTTTCGCTACGCAGCGCGGAATCGAGCCGACTTAACATGTCATCAAAACCCAAt                     <  1:797647/71‑1 (MQ=255)
                                                 gCCAACACCTCACCGCTTTCGCTACGCAGCGCGGAATCGAGCCGACTTAACATGTCATCAAAACCCAAt                     >  1:3266603/1‑69 (MQ=255)
                                                  ccAACAGCTCACCGCTTTCGCTACGCAGCGCGGAATCGAGCCGACTTAACATGTCATCAAAACCCaa                      >  1:497688/1‑67 (MQ=255)
                                                  ccAACACCTCACCGCTTTCGCTACGCAGCGCGGAATCGAGCCGACTTAACATGTCATCAAAACCCAATTc                   <  1:3260674/70‑1 (MQ=255)
                                                  ccAACACCTCACCGCTTTCGCTACGCAGCGCGGAATCGAGCCGACTTAACATGTCATCAAAACCCAATTCg                  <  1:940513/71‑1 (MQ=255)
                                                   cAACACCTCACCGCTTTCGCTACGCAGCGCGGAATCGAGCCGACTTAACATGTCATCAAAAc                          >  1:1545422/1‑62 (MQ=255)
                                                   cAACACCTCACCGCTTTCGCTACGCAGCGCGGAATCGAGCCGACTTAACATGTCATCAAAACCCAATTCGc                 <  1:1179994/71‑1 (MQ=255)
                                                   cAACACCTCACCGCTTTCGCTACGCAGCGCGGAATCGAGCCGACTTAACATGTCATCAAAACCCAATTCGc                 <  1:2643265/71‑1 (MQ=255)
                                                   cAACACCTCACCGCTTTCGCTACGCAGCGCGGAATCGAGCCGACTTAACATGTCATCAAAACCCAATTCGc                 <  1:707511/71‑1 (MQ=255)
                                                   cAACACCTCACCGCTTTCGCTACGCAGCGCGGAATCGAGCCGACTTAACATGTCATCAAAACCCAATTCGc                 <  1:412326/71‑1 (MQ=255)
                                                       acCTCACCGCTTTCGCTACGCAGCGCGGAATCGAGCCGAc                                            >  1:3288361/1‑40 (MQ=255)
                                                        ccTCACCGCTTTCGCTACGCAGCGCGGAATCGAGCCGACTTAACATGTCATCAAAACCCAATTCGCCAcgg            >  1:1439902/1‑71 (MQ=255)
                                                         cTCACCGCTTTCGCTACGCAGCGCGGAATCGAGCCGACTTAACATGTCATCAAAACCCAATTCGc                 <  1:3185423/65‑1 (MQ=255)
                                                         cTCACCGCTTTCGCTACGCAGCGCGGAATCGAGCCGACTTAACATGTCATCAAAACCCAATTCGc                 <  1:1792203/65‑1 (MQ=255)
                                                          tCACCGCTTTCGCTACGCAGCGCGGAATCGAGCCGACTTAACATGtc                                  <  1:3831407/47‑1 (MQ=255)
                                                          tCACCGCTTTCGCTACGCAGCGCGGAATCGAGCCGACTTAACATGtc                                  <  1:556394/47‑1 (MQ=255)
                                                          tCACCGCTTTCGCTACGCAGCGCGGAATCGAGCCGACTTAACATGTCATCAAAAccc                        <  1:1907152/57‑1 (MQ=255)
                                                            aCCGCTTTCGCTACGCAGCGCGGAATCGAGCCGACTTAACa                                      >  1:204873/1‑41 (MQ=255)
                                                            aCCGCTTTCGCTACGCAGCGCGGAATCGAGCCGACTTAACATGtcat                                >  1:103516/1‑47 (MQ=255)
                                                            aCCGCTTTCGCTACGCAGCGCGGAATCGAGCCGACTTAACATGtcat                                >  1:2017791/1‑47 (MQ=255)
                                                            aCCGCTTTCGCTACGCAGCGCGGAATCGAGCCGACTTAACATGtcat                                >  1:3268266/1‑47 (MQ=255)
                                                              cGCTTTCGCTACGCAGCGCGGAATCGAGCCGACTTaa                                        >  1:1615184/1‑37 (MQ=255)
                                                              cGCTTTCGCTACGCAGCGCGGAATCGAGCCGACTTAACATGTCATCAAAACCCAATTCGCCTcggcggg        >  1:2280115/1‑68 (MQ=255)
                                                                cTTTCGCTACGCAGCGCGGAATCGAGCCGACTTAACATGTCATCaaaa                           >  1:1655224/1‑48 (MQ=255)
                                                                cTTTCGCTACGCAGCGCGGAATCGAGCCGACTTAACATGTCATCaaaa                           >  1:3540307/1‑48 (MQ=255)
                                                                cTTTCGCTACGCAGCGCGGAATCGAGCCGACTTAACATGTCATCAAAACCCAATTCGCCACGGCGGCGttt    <  1:2153598/71‑1 (MQ=255)
                                                                cTTTCGCTACGCAGCGCGGAATCGAGCCGACTTAACATGTCATCAAAACCCAATTCGCCACGGCGGCGt      >  1:2749503/1‑69 (MQ=255)
                                                                   tCGCTACGCAGCGCGGAATCGAGCCGACTTAACATGTCATCAAAACCCAATTCGCCACGGCGGCGttttt  <  1:1821598/70‑1 (MQ=255)
                                                                   tCGCTACGCAGCGCGGAATCGAGCCGACTTAACATGTCATCAAAACCCAATTCGCCACGGCGGCGttt    >  1:1570550/1‑68 (MQ=255)
                                                                     |                                                                   
CTGAAATTCATCGATCATTGCCACCGGGAATCGCGTACGGATCGCCGCTGCCAACACCTCACCGCTTTCGCTACGCAGCGCGGAATCGAGCCGACTTAACATGTCATCAAAACCCAATTCGCCACGGCGGCGTTTTT  >  W3110S.gb/2953499‑2953635

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: