Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,285,704 A→G 21.4% S124S (TCA→TCG 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,7040AG21.4% 53.0 / 10.2 28S124S (TCA→TCGagaDN‑acetylgalactosamine‑specific enzyme IID component of PTS
Reads supporting (aligned to +/- strand):  ref base A (12/10);  new base G (3/3);  total (15/13)
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

CAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGC  >  W3110S.gb/3285634‑3285774
                                                                      |                                                                      
cAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCtcatca                                                                        <  1:3723649/71‑1 (MQ=255)
   tCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTcatcat                                                                       <  1:669878/69‑1 (MQ=255)
   tCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTcatcat                                                                       <  1:2999605/69‑1 (MQ=255)
   tCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCAttt                                                                     <  1:3671605/71‑1 (MQ=255)
   tCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCAtt                                                                      >  1:3390849/1‑70 (MQ=255)
    cGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTg                                                                    >  1:3828651/1‑71 (MQ=255)
     gCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTgc                                                                   >  1:2591348/1‑71 (MQ=255)
         ggATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTcatcat                                                                       <  1:3553489/63‑1 (MQ=255)
          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)
                                                                      |                                                                      
CAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGC  >  W3110S.gb/3285634‑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: