Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,227,651 C→G 21.1% G36G (GGC→GGG) ‡ ygjK → predicted glycosyl hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,227,6510CG21.1% 143.9 / 33.4 71G36G (GGC→GGG) ‡ygjKpredicted glycosyl hydrolase
Reads supporting (aligned to +/- strand):  ref base C (16/40);  new base G (9/6);  total (25/46)
Fisher's exact test for biased strand distribution p-value = 3.38e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTTTTTCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGG  >  W3110S.gb/3227592‑3227719
                                                           |                                                                    
aGTTTTTCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCGCAGTa                                                           >  1:3616519/1‑71 (MQ=255)
     ttCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGCGGCGCGCCAGTAcg                                                         <  1:3285005/68‑2 (MQ=37)
     ttCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGCGGCGCGCCAGTAcg                                                         <  1:5180493/68‑2 (MQ=37)
     ttCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGCGGCGCGCCAGTAcg                                                         <  1:1263489/68‑2 (MQ=37)
     ttCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGCGGCGCGCCAGTAcg                                                         <  1:4010414/68‑2 (MQ=37)
     ttCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGCGGCGCGCCAGTAcg                                                         <  1:1846295/68‑2 (MQ=37)
     ttCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGCGGCGCGCCAGTAcg                                                         <  1:1803750/68‑2 (MQ=37)
       ccGCCCATGCCGCTAACGCTGACAATTATAAAAACGTGATTAACCGTACTGCGGCGCGCCAGTAcggt                                                       >  1:1262586/1‑65 (MQ=25)
       ccGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCGCAGTACATGaa                                                     <  1:2621785/70‑1 (MQ=255)
       ccGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGCGGCGCGCCAGTAcggt                                                       >  1:1857517/1‑65 (MQ=37)
       ccGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGCGGCGCGCCAGTAcggt                                                       >  1:3191528/1‑65 (MQ=37)
       ccGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGCGGCGCGCCAGTAcggt                                                       >  1:3174358/1‑65 (MQ=37)
       ccGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGCGGCGCGCCAGTAcggt                                                       >  1:3332819/1‑65 (MQ=37)
       ccGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGCGGCGCGCCAGTAcggt                                                       >  1:197431/1‑65 (MQ=37)
       ccGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGCGGCGCGCCAGTAcggt                                                       >  1:4523540/1‑65 (MQ=37)
       ccGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGCGGCGCGCCAGTAcggt                                                       >  1:1060512/1‑65 (MQ=37)
       ccGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGCGGCGCGCCAGTAcggt                                                       >  1:4926167/1‑65 (MQ=37)
                ccGCTAACGCCGACAATTATAAAAACGTTATTAACCGTACTGGCGCGCCGCAGTACATGAAGgattacgat                                           >  1:4629915/1‑71 (MQ=255)
                ccGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCGCAGTACATGAAGgattacgat                                           >  1:5176388/1‑71 (MQ=255)
                   cTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCGCAGTACATGAAGgattacg                                             <  1:1585350/66‑1 (MQ=255)
                   cTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCGCAGTACATGAAGgattacg                                             <  1:3586652/66‑1 (MQ=255)
                   cTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCGCAGTACATGAAGgattacg                                             <  1:1293413/66‑1 (MQ=255)
                     aaCGCCGACAATTATAAAAACGTGATTAACCGTACTGgcg                                                                     >  1:1370540/1‑40 (MQ=255)
                     aaCGCCGACAATTATAAAAACGTGATTAACCGTACTGgcg                                                                     >  1:1241596/1‑40 (MQ=255)
                     aaCGCCGACAATTATAAAAACGTGATTAACCGTACTGgcg                                                                     >  1:2075389/1‑40 (MQ=255)
                                ttATAAAAACGTGATTAACCGTACTGGCGcgccgc                                                               <  1:2464103/35‑1 (MQ=255)
                                ttATAAAAACGTGATTAACCGTACTGGCGcgccgc                                                               <  1:3161597/35‑1 (MQ=255)
                                ttATAAAAACGTGATTAACCGTACTGGCGcgccgc                                                               <  1:330412/35‑1 (MQ=255)
                                ttATAAAAACGTGATTAACCGTACTGGCGcgccgc                                                               <  1:3590445/35‑1 (MQ=255)
                                ttATAAAAACGTGATTAACCGTACTGGCGcgccgc                                                               <  1:4253942/35‑1 (MQ=255)
                                ttATAAAAACGTGATTAACCGTACTGGCGcgccgc                                                               <  1:1707313/35‑1 (MQ=255)
                                ttATAAAAACGTGATTAACCGTACTGGCGcgccgc                                                               <  1:4595558/35‑1 (MQ=255)
                                ttATAAAAACGTGATTAACCGTACTGGCGcgccgc                                                               <  1:1356880/35‑1 (MQ=255)
                                ttATAAAAACGTGATTAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGc                           <  1:656806/71‑1 (MQ=255)
                                ttATAAAAACGTGATTAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGc                           <  1:2467448/71‑1 (MQ=255)
                                ttATAAAAACGTGATTAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGc                           <  1:2515227/71‑1 (MQ=255)
                                ttATAAAAACGTGATTAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGc                           <  1:4725515/71‑1 (MQ=255)
                                ttATAAAAACGTGATTAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGc                           <  1:1997490/71‑1 (MQ=255)
                                ttATAAAAACGTGATTAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGc                           <  1:3287632/71‑1 (MQ=255)
                                ttATAAAAACGTGATTAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGc                           <  1:1016704/71‑1 (MQ=255)
                                ttATAAAAACGTGATTAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGc                           <  1:1064458/71‑1 (MQ=255)
                                ttATAAAAACCTGATTAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGc                           <  1:3244944/71‑1 (MQ=255)
                                       aaCGTGATTAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTaa                    <  1:980321/71‑1 (MQ=255)
                                              ttAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTaa                    <  1:2622447/64‑1 (MQ=255)
                                              ttAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTaa                    <  1:2650900/64‑1 (MQ=255)
                                              ttAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTaa                    <  1:1324293/64‑1 (MQ=255)
                                              ttAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTaa                    <  1:2601004/64‑1 (MQ=255)
                                              ttAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTaa                    <  1:4235892/64‑1 (MQ=255)
                                              ttAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTa                     <  1:4697185/63‑1 (MQ=255)
                                              ttAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTa                     <  1:2660527/63‑1 (MQ=255)
                                              ttAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTa                     <  1:3116871/63‑1 (MQ=255)
                                              ttAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTa                     <  1:2566745/63‑1 (MQ=255)
                                                aaCCGTACTGGCGCGCCGCAGTACATGAAGTATTACGATTACGACGATCACCAGCGTTTTAATc                  >  1:3310775/1‑64 (MQ=255)
                                                aaCCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTaa                    >  1:3909829/1‑62 (MQ=255)
                                                aaCCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAAt                   >  1:1737061/1‑63 (MQ=255)
                                                aaCCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAAt                   >  1:1364112/1‑63 (MQ=255)
                                                aaCCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAAt                   >  1:1778352/1‑63 (MQ=255)
                                                aaCCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAAt                   >  1:3970477/1‑63 (MQ=255)
                                                aaCCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATc                  >  1:1979777/1‑64 (MQ=255)
                                                aaCCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATc                  >  1:948254/1‑64 (MQ=255)
                                                aaCCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATc                  >  1:1823851/1‑64 (MQ=255)
                                                aaCCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATc                  >  1:1299525/1‑64 (MQ=255)
                                                 aCCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTc          <  1:1837735/71‑1 (MQ=255)
                                                  ccGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCg         <  1:5178824/71‑1 (MQ=255)
                                                       cTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGAtctc    <  1:2530456/71‑1 (MQ=255)
                                                       cTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGAtctc    <  1:3122429/71‑1 (MQ=255)
                                                       cTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGAtctc    <  1:1309012/71‑1 (MQ=255)
                                                       cTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGAtctc    <  1:834273/71‑1 (MQ=255)
                                                       cTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGAtctc    <  1:1171832/71‑1 (MQ=255)
                                                       cTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATACGTTTTTCGAtctc    <  1:589014/71‑1 (MQ=255)
                                                         ggCGCGCCGCAGTACATGAAGGATTACGATTACGATGATCACCAGCGTTTTAATCCGTTTTTCGATCTCg   >  1:4363482/1‑70 (MQ=255)
                                                         ggCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGAtctc    <  1:5575533/69‑1 (MQ=255)
                                                         ggCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCg   >  1:2767853/1‑70 (MQ=255)
                                                         ggCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCg   >  1:800835/1‑70 (MQ=255)
                                                          gcgcgcCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCgg  <  1:693678/70‑1 (MQ=255)
                                                          gcgcgcCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCgg  <  1:1252220/70‑1 (MQ=255)
                                                          gcgcgcCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCgg  <  1:1047218/70‑1 (MQ=255)
                                                          gcgcgcCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCgg  <  1:664052/70‑1 (MQ=255)
                                                          gcgcgcCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCgg  <  1:1054746/70‑1 (MQ=255)
                                                          gcgcgcCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCgg  <  1:2333159/70‑1 (MQ=255)
                                                          gcgcgcCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCgg  <  1:4462256/70‑1 (MQ=255)
                                                          gcgcgcCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCgg  <  1:1996065/70‑1 (MQ=255)
                                                           |                                                                    
AGTTTTTCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGG  >  W3110S.gb/3227592‑3227719

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: