Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,227,657 G→C 21.1% P38P (CCG→CCC) ‡ ygjK → predicted glycosyl hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,227,6570GC21.1% 145.8 / 33.8 71P38P (CCG→CCC) ‡ygjKpredicted glycosyl hydrolase
Reads supporting (aligned to +/- strand):  ref base G (16/40);  new base C (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 = 3.84e-01

AGTTTTTCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTGGC  >  W3110S.gb/3227592‑3227727
                                                                 |                                                                      
aGTTTTTCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCGCAGTa                                                                   >  1:3616519/1‑71 (MQ=255)
     ttCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGCGGCGCGCCAGTAcg                                                                 <  1:1846295/68‑2 (MQ=37)
     ttCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGCGGCGCGCCAGTAcg                                                                 <  1:5180493/68‑2 (MQ=37)
     ttCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGCGGCGCGCCAGTAcg                                                                 <  1:1263489/68‑2 (MQ=37)
     ttCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGCGGCGCGCCAGTAcg                                                                 <  1:1803750/68‑2 (MQ=37)
     ttCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGCGGCGCGCCAGTAcg                                                                 <  1:4010414/68‑2 (MQ=37)
     ttCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGCGGCGCGCCAGTAcg                                                                 <  1:3285005/68‑2 (MQ=37)
       ccGCCCATGCCGCTAACGCTGACAATTATAAAAACGTGATTAACCGTACTGCGGCGCGCCAGTAcggt                                                               >  1:1262586/1‑65 (MQ=25)
       ccGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCGCAGTACATGaa                                                             <  1:2621785/70‑1 (MQ=255)
       ccGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGCGGCGCGCCAGTAcggt                                                               >  1:3174358/1‑65 (MQ=37)
       ccGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGCGGCGCGCCAGTAcggt                                                               >  1:1857517/1‑65 (MQ=37)
       ccGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGCGGCGCGCCAGTAcggt                                                               >  1:197431/1‑65 (MQ=37)
       ccGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGCGGCGCGCCAGTAcggt                                                               >  1:3332819/1‑65 (MQ=37)
       ccGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGCGGCGCGCCAGTAcggt                                                               >  1:4523540/1‑65 (MQ=37)
       ccGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGCGGCGCGCCAGTAcggt                                                               >  1:4926167/1‑65 (MQ=37)
       ccGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGCGGCGCGCCAGTAcggt                                                               >  1:1060512/1‑65 (MQ=37)
       ccGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGCGGCGCGCCAGTAcggt                                                               >  1:3191528/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)
                                ttATAAAAACGTGATTAACCGTACTGGCGcgccgc                                                                       <  1:330412/35‑1 (MQ=255)
                                ttATAAAAACGTGATTAACCGTACTGGCGcgccgc                                                                       <  1:2464103/35‑1 (MQ=255)
                                ttATAAAAACGTGATTAACCGTACTGGCGcgccgc                                                                       <  1:3161597/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:1997490/71‑1 (MQ=255)
                                ttATAAAAACGTGATTAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGc                                   <  1:2467448/71‑1 (MQ=255)
                                ttATAAAAACGTGATTAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGc                                   <  1:2515227/71‑1 (MQ=255)
                                ttATAAAAACGTGATTAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGc                                   <  1:3287632/71‑1 (MQ=255)
                                ttATAAAAACGTGATTAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGc                                   <  1:4725515/71‑1 (MQ=255)
                                ttATAAAAACGTGATTAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGc                                   <  1:1064458/71‑1 (MQ=255)
                                ttATAAAAACGTGATTAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGc                                   <  1:656806/71‑1 (MQ=255)
                                ttATAAAAACGTGATTAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGc                                   <  1:1016704/71‑1 (MQ=255)
                                ttATAAAAACCTGATTAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGc                                   <  1:3244944/71‑1 (MQ=255)
                                       aaCGTGATTAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTaa                            <  1:980321/71‑1 (MQ=255)
                                              ttAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTaa                            <  1:1324293/64‑1 (MQ=255)
                                              ttAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTaa                            <  1:2650900/64‑1 (MQ=255)
                                              ttAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTaa                            <  1:2622447/64‑1 (MQ=255)
                                              ttAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTaa                            <  1:2601004/64‑1 (MQ=255)
                                              ttAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTaa                            <  1:4235892/64‑1 (MQ=255)
                                              ttAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTa                             <  1:3116871/63‑1 (MQ=255)
                                              ttAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTa                             <  1:2660527/63‑1 (MQ=255)
                                              ttAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTa                             <  1:4697185/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:1364112/1‑63 (MQ=255)
                                                aaCCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAAt                           >  1:1737061/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:1823851/1‑64 (MQ=255)
                                                aaCCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATc                          >  1:1299525/1‑64 (MQ=255)
                                                aaCCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATc                          >  1:948254/1‑64 (MQ=255)
                                                 aCCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTc                  <  1:1837735/71‑1 (MQ=255)
                                                  ccGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCg                 <  1:5178824/71‑1 (MQ=255)
                                                       cTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGAtctc            <  1:1171832/71‑1 (MQ=255)
                                                       cTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGAtctc            <  1:1309012/71‑1 (MQ=255)
                                                       cTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGAtctc            <  1:834273/71‑1 (MQ=255)
                                                       cTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGAtctc            <  1:2530456/71‑1 (MQ=255)
                                                       cTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGAtctc            <  1:3122429/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:4462256/70‑1 (MQ=255)
                                                          gcgcgcCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCgg          <  1:1047218/70‑1 (MQ=255)
                                                          gcgcgcCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCgg          <  1:2333159/70‑1 (MQ=255)
                                                          gcgcgcCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCgg          <  1:693678/70‑1 (MQ=255)
                                                          gcgcgcCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCgg          <  1:664052/70‑1 (MQ=255)
                                                          gcgcgcCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCgg          <  1:1054746/70‑1 (MQ=255)
                                                          gcgcgcCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCgg          <  1:1996065/70‑1 (MQ=255)
                                                          gcgcgcCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCgg          <  1:1252220/70‑1 (MQ=255)
                                                               ccgcAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTg    <  1:5465399/71‑1 (MQ=255)
                                                               ccgcAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTg    <  1:3469850/71‑1 (MQ=255)
                                                               ccgcAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTg    <  1:5158607/71‑1 (MQ=255)
                                                               ccgcAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTg    <  1:514521/71‑1 (MQ=255)
                                                               ccgcAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTg    <  1:882307/71‑1 (MQ=255)
                                                               ccgcAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTg    <  1:160947/71‑1 (MQ=255)
                                                                cgcAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTg    <  1:3757235/70‑1 (MQ=255)
                                                                 gcAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCt     >  1:3530384/1‑68 (MQ=255)
                                                                 gcAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCt     >  1:4026043/1‑68 (MQ=255)
                                                                 gcAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTGGc  <  1:5513765/71‑1 (MQ=255)
                                                                 gcAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTGGc  <  1:3173027/71‑1 (MQ=255)
                                                                 |                                                                      
AGTTTTTCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTGGC  >  W3110S.gb/3227592‑3227727

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: