Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,953,831 T→C 22.5% T277A (ACA→GCA)  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,8310TC22.5% 73.7 / 20.4 40T277A (ACA→GCA) recBexonuclease V (RecBCD complex), beta subunit
Reads supporting (aligned to +/- strand):  ref base T (18/13);  new base C (5/4);  total (23/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.93e-01

CCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCTGATTGCTACGGTTAAACTTGC  >  W3110S.gb/2953762‑2953896
                                                                     |                                                                 
ccTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGtc                                                                  <  1:3321100/71‑1 (MQ=255)
    cGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGtctc                                                                <  1:2477573/69‑1 (MQ=255)
    cGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGtct                                                                 >  1:1546300/1‑68 (MQ=255)
    cGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCtctt                                                              >  1:993783/1‑71 (MQ=255)
    cGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCtctt                                                              >  1:234446/1‑71 (MQ=255)
    cGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGCctct                                                               >  1:957496/1‑70 (MQ=255)
    cGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCAGCAACTGATAACTGTTTGTCtctt                                                              >  1:1038503/1‑71 (MQ=255)
     gTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCtcttc                                                             <  1:3096761/71‑1 (MQ=255)
               gtAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTCTGTCTCTTCTTCTg                                                        <  1:3802461/65‑1 (MQ=255)
               ctAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGt                                                                   >  1:3450773/1‑55 (MQ=255)
               ctAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGt                                                                   >  1:2474806/1‑55 (MQ=255)
               ctAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCtcttcttct                                                         <  1:1905535/65‑1 (MQ=255)
               ctAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAg                                                   >  1:422810/1‑71 (MQ=255)
               ctAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAg                                                   >  1:1994984/1‑71 (MQ=255)
               ctAAGAAACCCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAg                                                   >  1:378249/1‑71 (MQ=255)
                    aaaCGCTGGGAGAATTTTTCCAGCGACTGCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCt                                              <  1:1901565/71‑1 (MQ=255)
                    aaaCGCTGGGAGAATTTTTCCAGCGACCCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCt                                              <  1:3700670/71‑1 (MQ=255)
                     aaCGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCt                                              >  1:3675937/1‑70 (MQ=255)
                     aaCGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCt                                              >  1:4052942/1‑70 (MQ=255)
                     aaCGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCt                                              >  1:279004/1‑70 (MQ=255)
                       cGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCtcttcttc                                                          >  1:610191/1‑56 (MQ=255)
                       cGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCtcttcttc                                                          >  1:1397215/1‑56 (MQ=255)
                        gCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATc                                          >  1:3933464/1‑71 (MQ=255)
                                    tttCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCAt                              >  1:3574339/1‑71 (MQ=255)
                                    tttCCAGCGACTCCGGCAACTGATAACTGTTTGCCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCa                               <  1:567505/70‑1 (MQ=255)
                                    tttCCAGCGACTCCGGCAACTGATAACTGTTTGCCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCa                               <  1:251991/70‑1 (MQ=255)
                                    tttCCAGCGACTCCGGCAACTGATAACTGTTTGCCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCa                               <  1:2475048/70‑1 (MQ=255)
                                           cGACTCCGGCAACTGATAACTGTTTGCCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCt                       <  1:3214677/71‑1 (MQ=255)
                                                       cTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATcc                                <  1:4061670/50‑1 (MQ=255)
                                                       cTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATcc                                <  1:2424966/50‑1 (MQ=255)
                                                        tGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCTGATTGCTACGGt           >  1:1426152/1‑70 (MQ=255)
                                                           tAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCTGATTGCTACGGTTaaa       <  1:1964055/71‑1 (MQ=255)
                                                           tAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCTGATTGCTACGGTTaaa       <  1:2359370/71‑1 (MQ=255)
                                                             aCTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCAt                              >  1:3169102/1‑46 (MQ=255)
                                                             aCTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCAt                              >  1:2526828/1‑46 (MQ=255)
                                                             aCTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCTGATTGCTACGGTTAAACt     >  1:2031864/1‑71 (MQ=255)
                                                             aCTGTTTGCCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGa                                   >  1:1978970/1‑41 (MQ=255)
                                                             aCTGTTTGCCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCTGATTGCTa               >  1:362770/1‑61 (MQ=255)
                                                             aCTGTTTGCCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCTGATTGCTa               >  1:3390782/1‑61 (MQ=255)
                                                                gTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGt                                      >  1:1050674/1‑35 (MQ=255)
                                                                gTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCt                       <  1:555357/50‑1 (MQ=255)
                                                                gTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCt                       <  1:2988180/50‑1 (MQ=255)
                                                                gTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCt                       <  1:2476131/50‑1 (MQ=255)
                                                                gTTTGCCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCTGATTGCTACGGTTAAACTTGc  >  1:3917361/1‑71 (MQ=255)
                                                                 tttGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATcc                                <  1:2302039/40‑1 (MQ=255)
                                                                     |                                                                 
CCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCTGATTGCTACGGTTAAACTTGC  >  W3110S.gb/2953762‑2953896

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: