Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,489,255 G→A 12.1% A103T (GCG→ACG)  ydcF → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,489,2550GA12.1% 90.4 / 3.6 33A103T (GCG→ACG) ydcFconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (18/11);  new base A (2/2);  total (20/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

TCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGATTGAAGACCAGTCAACAAA  >  W3110S.gb/1489193‑1489322
                                                              |                                                                   
tCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGa                                                                 <  1:1860960/67‑1 (MQ=255)
  gCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGc                                                           >  1:2329129/1‑71 (MQ=255)
       gCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCtcatc                                                      <  1:2456140/71‑1 (MQ=255)
       gCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCtcat                                                       >  1:2473706/1‑70 (MQ=255)
            cGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTc                                                 >  1:1652289/1‑71 (MQ=255)
               aCTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTgg                                              <  1:2567930/71‑1 (MQ=255)
                   caacaCTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGcaca                                          >  1:1690895/1‑71 (MQ=255)
                   caacaCTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGcaca                                          <  1:2671654/71‑1 (MQ=255)
                   caacaCTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGACGGATATCGCACATCAGTTCTGGcaca                                          >  1:1095799/1‑71 (MQ=255)
                    aacaCTATCCGCACCACTGGCAGAGCAGAAGCGAACATCCTGGCGGATATCGCTCATCAGTTCTGGCACAt                                         >  1:2704403/1‑71 (MQ=255)
                         tATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCg                                     >  1:1558501/1‑70 (MQ=255)
                           tCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCa                                   >  1:2346358/1‑70 (MQ=255)
                           tCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCa                                   >  1:2470282/1‑70 (MQ=255)
                           tCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCa                                   >  1:2473502/1‑70 (MQ=255)
                           tCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCAt                                  >  1:2565878/1‑71 (MQ=255)
                             cGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATg                                 <  1:2771454/70‑1 (MQ=255)
                                       gcagagcagaAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTg                       <  1:2121350/70‑1 (MQ=255)
                                       gcagagcagaAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTg                       <  1:2447550/70‑1 (MQ=255)
                                           agcataAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGATTg                  <  1:601688/71‑1 (MQ=255)
                                            gcagaAGCGACCATCCTGACGGATATCGCTCAGCAGTTCTGGCACATTCCGCATGAAAAAATCTGGAGTg                  >  1:673160/1‑70 (MQ=255)
                                             cagaAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGATTg                  >  1:1833634/1‑69 (MQ=255)
                                                aaGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGATTGAAGAc             <  1:82935/71‑1 (MQ=255)
                                                  gCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAAt                          >  1:1123342/1‑56 (MQ=255)
                                                  gCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAAt                          >  1:1802891/1‑56 (MQ=255)
                                                  gCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGATTGAAGAc             >  1:1783121/1‑69 (MQ=255)
                                                  gCGACCATCCTGACGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGATTGAAGACCa           <  1:751481/71‑1 (MQ=255)
                                                  gCGACCATCCTGACGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGATTGAAGACCa           <  1:1466059/71‑1 (MQ=255)
                                                      ccATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGATTg                  <  1:2197598/60‑1 (MQ=255)
                                                         tCCTGGCGGATATCGCTCATCAGTTCTTGCACATTCCGCATGAAAAAATCTGGATTGAAGAc             <  1:1515297/62‑1 (MQ=255)
                                                           cTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGATTGAAGACCAGTcaacaa   >  1:655746/1‑70 (MQ=255)
                                                           cTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGATTGAAGACCAGTCAACaaa  >  1:2253260/1‑71 (MQ=255)
                                                           cTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGATTGAAGACCAGTCAACaaa  >  1:59222/1‑71 (MQ=255)
                                                           cTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGATTGAAGACCAGTCAACaaa  >  1:962958/1‑71 (MQ=255)
                                                              |                                                                   
TCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGATTGAAGACCAGTCAACAAA  >  W3110S.gb/1489193‑1489322

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: