Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,512,230 A→G 16.1% W371R (TGG→CGG)  ypdG ← predicted enzyme IIC component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,512,2300AG16.1% 70.5 / 8.2 31W371R (TGG→CGG) ypdGpredicted enzyme IIC component of PTS
Reads supporting (aligned to +/- strand):  ref base A (16/10);  new base G (2/3);  total (18/13)
Fisher's exact test for biased strand distribution p-value = 6.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

CAGTTATGACTGCCCCCAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGAC  >  W3110S.gb/2512161‑2512290
                                                                     |                                                            
cAGTTATGACTGCCCCCAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCa                                                              <  1:3157124/70‑1 (MQ=255)
  gTTATGACTGCCCCCAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGAt                                                           <  1:3473844/71‑1 (MQ=255)
    tATGACTGCCCCCAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAg                                                         >  1:1025745/1‑71 (MQ=255)
    tATGACTGCCCCCAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAg                                                         >  1:2875017/1‑71 (MQ=255)
          tGCCCCCAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACt                                                   >  1:1913573/1‑71 (MQ=255)
                 aGCGAGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGaa                                            >  1:2112104/1‑71 (MQ=255)
                    gcgATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATa                                                          >  1:1159279/1‑54 (MQ=255)
                    gcgATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAg                                                         >  1:1395473/1‑55 (MQ=255)
                    gcgATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAg                                                         >  1:1692363/1‑55 (MQ=255)
                                 tataGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCggg                                              <  1:3071075/53‑1 (MQ=255)
                                   taGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCg                                                       <  1:2930707/42‑1 (MQ=255)
                                   taGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCg                                                       <  1:3115338/42‑1 (MQ=255)
                                   taGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGaa                                            >  1:2701845/1‑53 (MQ=255)
                                   taGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCa                                    >  1:2969096/1‑61 (MQ=255)
                                   taGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCccag                          >  1:1312669/1‑71 (MQ=255)
                                   taGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCccag                          >  1:1308380/1‑71 (MQ=255)
                                   taGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCccag                          >  1:3623178/1‑71 (MQ=255)
                                   taGACGCCCAGGTTAGTGACCAGGGCCATGCCCAGATAGCGGACTCCgg                                               >  1:1288315/1‑49 (MQ=39)
                                       cGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCg                                                       <  1:847720/38‑1 (MQ=255)
                                         cccAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCggg                                              <  1:3074538/45‑1 (MQ=255)
                                             ggTTAGTGACCAGCGGCCATGCCCGGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAgcg                >  1:3209924/1‑71 (MQ=255)
                                             ggTTAGTGACCAGCGGCCATGCCCGGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAgcg                >  1:2452084/1‑71 (MQ=255)
                                                tAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAgcggcg             >  1:1918425/1‑71 (MQ=255)
                                                 aGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGTAACCATTGCACTGCGCCCAGCCAGACAgcggcg             <  1:3378874/70‑1 (MQ=255)
                                                 aGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAgcggcg             <  1:1073250/70‑1 (MQ=255)
                                                     aCCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGtt          >  1:306256/1‑69 (MQ=255)
                                                     aCCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGtt          >  1:727461/1‑69 (MQ=255)
                                                        aGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGcc     <  1:138161/71‑1 (MQ=255)
                                                            gCCATGCCGGGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGAc  <  1:618572/70‑1 (MQ=255)
                                                            gCCATGCCCGGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGAc  <  1:3457781/70‑1 (MQ=255)
                                                            gCCATGCCCGGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGAc  <  1:3045715/70‑1 (MQ=255)
                                                                tGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAgcggc              <  1:850379/54‑1 (MQ=255)
                                                                     |                                                            
CAGTTATGACTGCCCCCAGCGCGATCCCCGCCATATAGACGCCCAGGTTAGTGACCAGCGGCCATGCCCAGATAGCGGACTCCGGGAACCATTGCACTGCGCCCAGCCAGACAGCGGCGGTTGAGCCGAC  >  W3110S.gb/2512161‑2512290

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: