Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,237,438 G→A 24.5% P461S (CCA→TCA)  ycgB ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,237,4380GA24.5% 90.8 / 28.5 49P461S (CCA→TCA) ycgBconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (24/13);  new base A (9/3);  total (33/16)
Fisher's exact test for biased strand distribution p-value = 7.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.09e-01

CCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTGAATAT  >  W3110S.gb/1237371‑1237497
                                                                   |                                                           
ccccACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGaaa                                                           >  1:3111914/1‑70 (MQ=255)
ccccACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGaaa                                                           >  1:3182712/1‑70 (MQ=255)
 cccACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAata                                                         <  1:3242203/71‑1 (MQ=255)
  ccACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAatat                                                        <  1:3043460/71‑1 (MQ=255)
  ccACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAatat                                                        <  1:4350523/71‑1 (MQ=255)
  ccACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAatat                                                        <  1:637103/71‑1 (MQ=255)
  ccACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAatat                                                        <  1:2178892/71‑1 (MQ=255)
     caGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAAc                                                     >  1:1146805/1‑71 (MQ=255)
     caGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAAc                                                     >  1:3225830/1‑71 (MQ=255)
     caGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGAAATATAAc                                                     >  1:730161/1‑71 (MQ=255)
          gATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGc                                                >  1:3047434/1‑71 (MQ=255)
             gCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGAAATATAACGCAGCGt                                              >  1:4751258/1‑70 (MQ=255)
              cacaTGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTc                                             <  1:911214/70‑1 (MQ=255)
              cacaTGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTc                                             <  1:1439113/70‑1 (MQ=255)
              cacaTGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGAAat                                                          >  1:2904124/1‑57 (MQ=255)
              cacaTGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGAAATATAACGCAGCGTc                                             <  1:4148464/70‑1 (MQ=255)
                 aTGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCg                                         <  1:2499067/71‑1 (MQ=255)
                  tGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGa                                                            >  1:707524/1‑51 (MQ=255)
                  tGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGAAATATAACGCAGCGTCAGCg                                         >  1:1760346/1‑70 (MQ=255)
                   gCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGAAATATAACGCAGCGTCAGCGaa                                       >  1:2653027/1‑71 (MQ=255)
                     ttCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATaa                                                      >  1:3238689/1‑54 (MQ=255)
                       cAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACgg                                    >  1:4808907/1‑70 (MQ=255)
                         ggACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGaa                                       >  1:4724935/1‑65 (MQ=255)
                         ggACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGaa                                       >  1:3432141/1‑65 (MQ=255)
                         ggACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAAAGCAGCGACAGCGaa                                       >  1:204012/1‑65 (MQ=255)
                             ttcttTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGaa                                                           >  1:369578/1‑41 (MQ=255)
                             ttcttTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGaa                                                           >  1:3766446/1‑41 (MQ=255)
                               cttTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCcgcg                            >  1:4150383/1‑70 (MQ=255)
                               cttTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCcgcg                            >  1:2134846/1‑70 (MQ=255)
                               cttTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCcgcg                            >  1:1726444/1‑70 (MQ=255)
                                  tGCGCCCCCGATCCAGCGGTGCGCGATTATGTGAAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCaa                         <  1:1069648/70‑1 (MQ=255)
                                          cGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCaaa                        >  1:4819115/1‑63 (MQ=255)
                                          cGATCCAGCGGTGCGCGATTATGTGAAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGtt                >  1:4126497/1‑71 (MQ=255)
                                          cGATCCAGCGGTGCGCGATTATGTGAAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGtt                >  1:1058804/1‑71 (MQ=255)
                                          cGATCCAGCGGTGCGCGATTATGTGAAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGtt                >  1:4262759/1‑71 (MQ=255)
                                               cAGCGGTGCGCGATTATGTGAAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGa           <  1:4050919/71‑1 (MQ=255)
                                                aGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTc               >  1:3195740/1‑66 (MQ=255)
                                                  cGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCt        <  1:379859/71‑1 (MQ=255)
                                                   ggTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCt        >  1:4401166/1‑70 (MQ=255)
                                                   ggTGCGCGATTATGTGAAATATAACGCAGCGTCAGCGAACGGTcgc                                >  1:911713/1‑46 (MQ=255)
                                                    gTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCa             >  1:1251788/1‑64 (MQ=255)
                                                    gTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCa             >  1:3081629/1‑64 (MQ=255)
                                                      gcgcGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTc               <  1:2229251/60‑1 (MQ=255)
                                                      gcgcGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTc               <  1:1934720/60‑1 (MQ=255)
                                                       cgcgATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTGAata   <  1:2225988/71‑1 (MQ=255)
                                                        gcgATTATGTGGAATATAACGCAGCGTCAGCGAACg                                     <  1:4155991/36‑1 (MQ=255)
                                                        gcgATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTGAatat  >  1:436093/1‑71 (MQ=255)
                                                        gcgATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTGAatat  >  1:4893909/1‑71 (MQ=255)
                                                        gcgATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTGAatat  >  1:152409/1‑71 (MQ=255)
                                                         cgATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTGAatat  >  1:2877686/1‑70 (MQ=255)
                                                           aTTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTGAata   >  1:2704355/1‑67 (MQ=255)
                                                                   |                                                           
CCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTGAATAT  >  W3110S.gb/1237371‑1237497

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: