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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,227,6560CG21.1% 142.8 / 34.2 71P38R (CCG→CGG) ‡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 = 9.34e-01

AGTTTTTCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTG  >  W3110S.gb/3227592‑3227725
                                                                |                                                                     
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:2650900/64‑1 (MQ=255)
                                              ttAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTaa                          <  1:2622447/64‑1 (MQ=255)
                                              ttAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTaa                          <  1:2601004/64‑1 (MQ=255)
                                              ttAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTaa                          <  1:1324293/64‑1 (MQ=255)
                                              ttAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTaa                          <  1:4235892/64‑1 (MQ=255)
                                              ttAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTa                           <  1:3116871/63‑1 (MQ=255)
                                              ttAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTa                           <  1:4697185/63‑1 (MQ=255)
                                              ttAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTa                           <  1:2660527/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:3970477/1‑63 (MQ=255)
                                                aaCCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAAt                         >  1:1778352/1‑63 (MQ=255)
                                                aaCCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAAt                         >  1:1737061/1‑63 (MQ=255)
                                                aaCCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAAt                         >  1:1364112/1‑63 (MQ=255)
                                                aaCCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATc                        >  1:948254/1‑64 (MQ=255)
                                                aaCCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATc                        >  1:1979777/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:834273/71‑1 (MQ=255)
                                                       cTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGAtctc          <  1:1171832/71‑1 (MQ=255)
                                                       cTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGAtctc          <  1:1309012/71‑1 (MQ=255)
                                                       cTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGAtctc          <  1:3122429/71‑1 (MQ=255)
                                                       cTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGAtctc          <  1:2530456/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:664052/70‑1 (MQ=255)
                                                          gcgcgcCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCgg        <  1:693678/70‑1 (MQ=255)
                                                          gcgcgcCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCgg        <  1:1054746/70‑1 (MQ=255)
                                                          gcgcgcCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCgg        <  1:1047218/70‑1 (MQ=255)
                                                          gcgcgcCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCgg        <  1:1252220/70‑1 (MQ=255)
                                                          gcgcgcCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCgg        <  1:4462256/70‑1 (MQ=255)
                                                          gcgcgcCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCgg        <  1:1996065/70‑1 (MQ=255)
                                                          gcgcgcCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCgg        <  1:2333159/70‑1 (MQ=255)
                                                               ccgcAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTg  <  1:5465399/71‑1 (MQ=255)
                                                               ccgcAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTg  <  1:5158607/71‑1 (MQ=255)
                                                               ccgcAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTg  <  1:514521/71‑1 (MQ=255)
                                                               ccgcAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTg  <  1:160947/71‑1 (MQ=255)
                                                               ccgcAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTg  <  1:3469850/71‑1 (MQ=255)
                                                               ccgcAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTg  <  1:882307/71‑1 (MQ=255)
                                                                cgcAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTg  <  1:3757235/70‑1 (MQ=255)
                                                                |                                                                     
AGTTTTTCCGCCCATGCCGCTAACGCCGACAATTATAAAAACGTGATTAACCGTACTGGCGCGCCGCAGTACATGAAGGATTACGATTACGACGATCACCAGCGTTTTAATCCGTTTTTCGATCTCGGAGCCTG  >  W3110S.gb/3227592‑3227725

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: