Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,688,001 A→G 66.7% S440G (AGT→GGT)  yifB → predicted bifunctional protein, enzyme and transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,688,0010AG66.7% 24.5 / 18.4 24S440G (AGT→GGT) yifBpredicted bifunctional protein, enzyme and transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (4/4);  new base G (14/2);  total (18/6)
Fisher's exact test for biased strand distribution p-value = 1.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAACAACGCGTAATGGCCGCCAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGTTGGAAGGAACACTGA  >  W3110S.gb/3687932‑3688067
                                                                     |                                                                  
aaaCAACGCGTAATGGCCGCCAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAg                                                                   >  1:1424368/1‑71 (MQ=255)
aaaCAACGCGTAATGGCCGCCAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAg                                                                   >  1:1612439/1‑71 (MQ=255)
aaaCAACGCGTAATGGCCGCCAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAg                                                                   >  1:211907/1‑71 (MQ=255)
               gccgccAGAGAGCGCCAATTTAAGAGGCAGAATAAACTGAATGCCTGGCTGGATGGTCCGGAAATACGcc                                                     >  1:1402103/1‑70 (MQ=255)
                     agagagCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATGGTCCGGAAATACGCCAATTCt                                               >  1:596344/1‑70 (MQ=255)
                              cAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATGGTCCGGaaa                                                           <  1:802048/49‑1 (MQ=255)
                                aTTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATGGTCCGGAAATa                                                         >  1:1319466/1‑49 (MQ=255)
                                   tAAGCGGCAGAATAAACTGAATGCCTGGCTGGATGGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGc                                 >  1:1815780/1‑70 (MQ=255)
                                   tAAGCGGCAGAATAAACTGAATGCCTGGCTGGATGGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGc                                 >  1:467368/1‑70 (MQ=255)
                                        ggCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGa                               <  1:1830340/67‑1 (MQ=255)
                                          cAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGa                               <  1:1160658/65‑1 (MQ=255)
                                               tAAACTGAATGCCTGGCTGGATGGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGAtgt                    >  1:859580/1‑71 (MQ=255)
                                               tAAACTGAATGCCTGGCTGGATGGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGAtgt                    >  1:451577/1‑71 (MQ=255)
                                                 aaCTGAATGCCTGGCTGGATGGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGAtgtg                   <  1:1788864/70‑1 (MQ=255)
                                                    tGAATGCCTGGCTGGATGGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGc                         >  1:972233/1‑61 (MQ=255)
                                                    tGAATGCCTGGCTGGATGGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGc                         >  1:927267/1‑61 (MQ=255)
                                                    tGAATGCCTGGCTGGATGGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGc                         >  1:584207/1‑61 (MQ=255)
                                                      aaTGCCTGGCTGGATGGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGaa                              >  1:1748204/1‑54 (MQ=255)
                                                      aaTGCCTGGCTGGATGGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGa                            >  1:679766/1‑56 (MQ=255)
                                                      aaTGCCTGGCTGGATGGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGa                            >  1:1326338/1‑56 (MQ=255)
                                                       aTGCCTGGCTGGATGGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGa                            >  1:1545771/1‑55 (MQ=255)
                                                          ccTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGa                               <  1:1402727/49‑1 (MQ=255)
                                                              gctggATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGTTGGAAGGAacac     >  1:1233388/1‑71 (MQ=255)
                                                                  gATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGTTGGAAGGAACACTGa  <  1:695967/70‑1 (MQ=255)
                                                                     |                                                                  
AAACAACGCGTAATGGCCGCCAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGTTGGAAGGAACACTGA  >  W3110S.gb/3687932‑3688067

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: