Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,285,713 C→T 24.2% S127S (AGC→AGT agaD → N‑acetylgalactosamine‑specific enzyme IID component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,285,7130CT24.2% 42.2 / 8.9 25S127S (AGC→AGTagaDN‑acetylgalactosamine‑specific enzyme IID component of PTS
Reads supporting (aligned to +/- strand):  ref base C (10/8);  new base T (3/3);  total (13/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGC  >  W3110S.gb/3285644‑3285774
                                                                     |                                                             
gATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTgccagcc                                                              <  1:3009389/71‑1 (MQ=255)
    ggCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAgg                                                           >  1:257320/1‑70 (MQ=255)
         tgcgatTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCATCCAGGGAAAcc                                                     <  1:2556018/71‑1 (MQ=255)
               tttCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAAc                                                      >  1:2268544/1‑64 (MQ=255)
                ttCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTgg                                               <  1:1930675/70‑1 (MQ=255)
                ttCTGGTTTACCTTGTTGCCAATTATGGCGGGGATTTGCTCATCGTTTGCCAGTCAGGGGAACctgttg                                                >  1:2277954/1‑69 (MQ=255)
                  cTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGaaa                                                       >  1:2818685/1‑60 (MQ=255)
                  cTGGTTTACTTTGTTGACGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGaaa                                                       >  1:980420/1‑60 (MQ=255)
                        tACCTTTTTGCCAATTATGGCGGGGATTTGCTCATCGTTTGCCAGTCAgggg                                                         <  1:3050803/52‑2 (MQ=19)
                        tACCTTGTTGCCTATTATGGCGGGGATTTGCTCATCGTTTGCCAGTCAgggg                                                         <  1:435324/52‑2 (MQ=25)
                        tACCTTGTTGCCAATTATGGCGGGGATTTGCTCATCGTTTGCCAGTCAgggg                                                         <  1:1839672/52‑2 (MQ=25)
                           tttGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGccagccag                                                            >  1:1397923/1‑46 (MQ=255)
                           tttGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTgg                                               >  1:490224/1‑59 (MQ=255)
                           tttGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTgg                                               >  1:213939/1‑59 (MQ=255)
                                       tATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAAc                                                      >  1:2197917/1‑40 (MQ=255)
                                       tATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACct                                                    <  1:1144370/42‑1 (MQ=255)
                                        aTGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGttt                       <  1:102421/70‑1 (MQ=255)
                                          ggCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTa                      <  1:623632/69‑1 (MQ=255)
                                           gCGGGGATTTGCTCATCGTTTGCCAGTCAGGGGAACctgttg                                                >  1:727635/1‑42 (MQ=19)
                                           gCGGGGATTTGCTCATCGTTTGCCAGTCAGGGGAACCTGTTgg                                               >  1:3743939/1‑43 (MQ=19)
                                                 aTTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTAt                                   <  1:3097571/49‑1 (MQ=255)
                                                    tGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCtt          >  1:3189565/1‑71 (MQ=255)
                                                    tGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCt           <  1:1559815/70‑1 (MQ=255)
                                                      cTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCtttt        <  1:2351797/71‑1 (MQ=255)
                                                        catcatTTGCCAGGCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGcc                           <  1:3132766/50‑1 (MQ=255)
                                                           catTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTg   >  1:2227638/1‑71 (MQ=255)
                                                            atTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTgc  >  1:2581569/1‑71 (MQ=255)
                                                                     |                                                             
GATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGC  >  W3110S.gb/3285644‑3285774

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: