Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,951,877 G→A 28.3% A928V (GCT→GTT)  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,951,8770GA28.3% 86.9 / 37.3 53A928V (GCT→GTT) recBexonuclease V (RecBCD complex), beta subunit
Reads supporting (aligned to +/- strand):  ref base G (20/18);  new base A (8/7);  total (28/25)
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.13e-01

CGGTGACGCACCGCGCGGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCAAACCAGAGTAGC  >  W3110S.gb/2951807‑2951946
                                                                      |                                                                     
cGGTGACGCACCGCGCGGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTgcag                                                                       >  1:4699545/1‑71 (MQ=255)
cGGTGACGCACCGCGCGGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTgcag                                                                       >  1:1459244/1‑71 (MQ=255)
cGGTGACGCACCGCGCGGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTgcag                                                                       >  1:398116/1‑71 (MQ=255)
cGGTGACGCACCGCGCGGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTgcag                                                                       >  1:2376513/1‑71 (MQ=255)
 ggTGACGCACCGCGCGGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTgcag                                                                       >  1:5061671/1‑70 (MQ=255)
 ggTGACGCACCGCGCGGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTgcag                                                                       >  1:3176567/1‑70 (MQ=255)
   tGACGCACCGCGCGGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCAt                                                                    >  1:2655481/1‑71 (MQ=255)
     aCGCACCGCGCGGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATc                                                                   <  1:4602331/70‑1 (MQ=255)
     aCGCACCGCGCGGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATc                                                                   <  1:2609274/70‑1 (MQ=255)
     aCGCACCGCGCGGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCg                                                                  <  1:4818902/71‑1 (MQ=255)
     aCGCACCGCGCGGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAACATc                                                                   <  1:976131/70‑1 (MQ=255)
         aCCGCGCGGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGa                                                                 <  1:3252264/68‑1 (MQ=255)
         aCCGCGCGGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGa                                                                 <  1:2365209/68‑1 (MQ=255)
         aCCGCGCGGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACAt                                                              >  1:1311832/1‑71 (MQ=255)
          ccGCGCGGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAACATCGACATc                                                             <  1:1222869/71‑1 (MQ=255)
            gcgcgGAAACTGATGTTGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAACATCGACATCCa                                                           <  1:2538375/71‑1 (MQ=255)
            gcgcgGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATc                                                                   >  1:2826103/1‑63 (MQ=255)
            gcgcgGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATc                                                                   >  1:1690262/1‑63 (MQ=255)
            gcgcgGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCa                                                           <  1:1490761/71‑1 (MQ=255)
               cgGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGcc                                                        <  1:238952/71‑1 (MQ=255)
               cgGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAACATCGACATCCAGcc                                                        >  1:285843/1‑71 (MQ=255)
               cgGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAACATCGACATCCAGcc                                                        >  1:4787436/1‑71 (MQ=255)
                   aaCTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAgg                                                    <  1:155742/71‑1 (MQ=255)
                        aTGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCa                                               <  1:3188753/71‑1 (MQ=255)
                         tgtgGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCaa                                              <  1:1742645/71‑1 (MQ=255)
                           tggtgTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGcc                                                        <  1:1585024/59‑1 (MQ=255)
                                ttAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAACATCGACATCCAGCCGAGGCATCAAATCCt                                         <  1:1937736/69‑1 (MQ=255)
                                  aaCGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTggg                                      >  1:3369229/1‑70 (MQ=255)
                                  aaCGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTggg                                      >  1:191040/1‑70 (MQ=255)
                                  aaCGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTggg                                      >  1:2203508/1‑70 (MQ=255)
                                  aaCGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGc                                     >  1:1223435/1‑71 (MQ=255)
                                  aaCGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGc                                     >  1:1990578/1‑71 (MQ=255)
                                  aaCGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGc                                     >  1:3033210/1‑71 (MQ=255)
                                  aaCGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGc                                     >  1:2712459/1‑71 (MQ=255)
                                  aaCGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAACATCGACATCCAGCCGAGGCATCAAATCCTggg                                      <  1:4848315/70‑1 (MQ=255)
                                     gTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGAt                                  <  1:2662578/71‑1 (MQ=255)
                                      tCGGTTCTTCAACGACGCTTGCAACGCCTGCAACATCGACATCCAGCCGAGGCa                                                  >  1:4497307/1‑54 (MQ=255)
                                      tCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATa                                 >  1:2347082/1‑71 (MQ=255)
                                      tCGGTTCTTCAACGACGCTGGCAACGCCTGCAACATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATa                                 >  1:3005629/1‑71 (MQ=255)
                                          ttcttcAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGt                             <  1:1356566/71‑1 (MQ=255)
                                            cttcAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGa                           >  1:1103755/1‑71 (MQ=255)
                                            cttcAACGACGCTGGCAACGCCTGCAACATCGACATCCAGCCGAGGCATCAAATCCt                                         <  1:3586626/57‑1 (MQ=255)
                                             ttcAACGACGCTGGCAACGCCTGCAACATCGACATCCAGCCGAGGCATCAAATCCt                                         <  1:942785/56‑1 (MQ=255)
                                              tcAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATcc                                          >  1:308882/1‑54 (MQ=255)
                                              tcAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATcc                                          >  1:3236722/1‑54 (MQ=255)
                                              tcAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGAc                          >  1:4977533/1‑70 (MQ=255)
                                                 acgacgCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACg                      >  1:2917222/1‑71 (MQ=255)
                                                 acgacgCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACg                      >  1:3955861/1‑71 (MQ=255)
                                                   gacgCTGGCAACGCCTGCAACATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCt                    >  1:3150395/1‑71 (MQ=255)
                                                      gCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGtt                 >  1:3584309/1‑71 (MQ=255)
                                                      gCTGGCAACGCCTGCAACATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGtt                 >  1:63598/1‑71 (MQ=255)
                                                             aCGCCTGCAACATCGACATCCAGCCGAGGCATCAAATc                                           >  1:3302693/1‑38 (MQ=255)
                                                             aCGCCTGCAACATCGACATCCAGCCGAGGCATCAAATc                                           >  1:270840/1‑38 (MQ=255)
                                                               gCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCAAACCag        <  1:1506439/71‑1 (MQ=255)
                                                               gCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCAAACCag        <  1:1916894/71‑1 (MQ=255)
                                                               gCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCAAACCag        <  1:4771044/71‑1 (MQ=255)
                                                                 cTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCAAACCagag      <  1:4642779/71‑1 (MQ=255)
                                                                 cTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCAAACCagag      >  1:2536446/1‑71 (MQ=255)
                                                                 cTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCAAACCagag      <  1:2176746/71‑1 (MQ=255)
                                                                    cagcaTCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCAAACCAGAGTAg   >  1:4227451/1‑71 (MQ=255)
                                                                     agcaTCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCAAACCAGAGTAGc  <  1:3517067/71‑1 (MQ=255)
                                                                     agcaTCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCAAACCAGAGTAGc  <  1:1322877/71‑1 (MQ=255)
                                                                      |                                                                     
CGGTGACGCACCGCGCGGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCAAACCAGAGTAGC  >  W3110S.gb/2951807‑2951946

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: