Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,285,692 A→G 17.6% G120G (GGA→GGG 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,6920AG17.6% 73.6 / 9.4 34G120G (GGA→GGGagaDN‑acetylgalactosamine‑specific enzyme IID component of PTS
Reads supporting (aligned to +/- strand):  ref base A (13/15);  new base G (3/3);  total (16/18)
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.51e-01

ACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTA  >  W3110S.gb/3285623‑3285754
                                                                     |                                                              
aCTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGaa                                                               <  1:1647203/71‑1 (MQ=255)
   gTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAAtt                                                             >  1:3364990/1‑70 (MQ=255)
           cAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCtcatca                                                    <  1:3723649/71‑1 (MQ=255)
              tCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAAt                                                              <  1:1513757/58‑1 (MQ=255)
              tCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAAt                                                              <  1:2274355/58‑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)
                              tgcgatTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCtcatc                                                     <  1:1988291/51‑1 (MQ=255)
                              tgcgatTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCtcatc                                                     <  1:135075/51‑1 (MQ=255)
                              tgcgatTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCATCCAGGGAAAcc                                 <  1:2556018/71‑1 (MQ=255)
                                cgatTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGatt                                                              >  1:2698072/1‑38 (MQ=255)
                                    tttCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCtc                                                        >  1:1059585/1‑42 (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:213939/1‑59 (MQ=255)
                                                tttGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTgg                           >  1:490224/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)
                                                                     |                                                              
ACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTA  >  W3110S.gb/3285623‑3285754

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: