Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 389,198 T→C 11.9% M48V (ATG→GTG)  ykiB ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb389,1980TC11.9% 150.4 / 12.3 59M48V (ATG→GTG) ykiBhypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (22/30);  new base C (3/4);  total (25/34)
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.82e-01

GGTTATTATTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATT  >  W3110S.gb/389136‑389265
                                                              |                                                                   
ggTTATTATTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCACATTCATAg                                                             >  1:4132086/1‑71 (MQ=255)
 gTTATTATTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCACATTCATAgg                                                            <  1:2296374/71‑1 (MQ=255)
  ttattaTTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCAt                                                               <  1:4478176/67‑1 (MQ=255)
   tattatTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGa                                                           >  1:1200936/1‑70 (MQ=255)
    attattTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGata                                                         >  1:2021286/1‑71 (MQ=255)
    attattTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGata                                                         >  1:4377743/1‑71 (MQ=255)
    attattTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGata                                                         >  1:3888670/1‑71 (MQ=255)
    attattTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGata                                                         >  1:3772701/1‑71 (MQ=255)
    attattTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCACATTCATAGGata                                                         >  1:2747036/1‑71 (MQ=255)
     ttattTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGatat                                                        <  1:1536358/71‑1 (MQ=255)
     ttattTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGata                                                         <  1:271199/70‑1 (MQ=255)
       attTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATAtt                                                       >  1:2569896/1‑70 (MQ=255)
        ttttttGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCACATTCATAGGATAttt                                                      <  1:1820351/70‑1 (MQ=255)
         tttttGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCata                                                                    >  1:3078012/1‑55 (MQ=255)
          ttttGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCata                                                                    >  1:2991952/1‑54 (MQ=255)
                  tttGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTc                                                                 <  1:1969254/49‑1 (MQ=255)
                  tttGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGata                                                         <  1:3794486/57‑1 (MQ=255)
                    tGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCata                                                                    >  1:3548315/1‑44 (MQ=255)
                    tGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCat                                                                     >  1:941974/1‑43 (MQ=255)
                    tGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgtttgtgt                                         >  1:2477247/1‑71 (MQ=255)
                    tGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgtttgtgt                                         >  1:1933487/1‑71 (MQ=255)
                     ggTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCAt                                                               <  1:2761235/48‑1 (MQ=255)
                     ggTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgtttgtgt                                         <  1:1652109/70‑1 (MQ=255)
                      gTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCACATTCATAGGATATTTATCTGgtttgtgtt                                        <  1:783049/70‑1 (MQ=255)
                       ttGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCatat                                                                   <  1:1804620/42‑1 (MQ=255)
                         gcgcTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgtttgtgtttgt                                     >  1:2189006/1‑70 (MQ=255)
                         gcgcTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgtttgtgtttgt                                     >  1:3840812/1‑70 (MQ=255)
                          ccctgctATAAATTGAATAATTAAAACTTTGTTGCACATTCATAGGATATTTATCTGgtt                                              <  1:1798917/58‑1 (MQ=255)
                            ctgctATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGata                                 >  1:1141211/1‑71 (MQ=255)
                            ctgctATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGata                                 >  1:3187923/1‑71 (MQ=255)
                              gctATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGatat                                >  1:340535/1‑70 (MQ=255)
                               ctATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGAtattt                              <  1:1771067/71‑1 (MQ=255)
                               ctATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGAtatt                               <  1:3310298/70‑1 (MQ=255)
                               ctATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGAtatt                               <  1:3120613/70‑1 (MQ=255)
                                 ataAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGatat                                >  1:1265624/1‑67 (MQ=255)
                                  taAATTGAATAATTAAAACTTTGTTGCATATTCATAGGatat                                                        <  1:1833850/42‑1 (MQ=255)
                                  taAATTGAATAATTAAAACTTTGTTGCATATTCATAGGata                                                         <  1:737278/41‑1 (MQ=255)
                                  taAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGAtatttat                            <  1:1847618/70‑1 (MQ=255)
                                  taAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGAtatt                               <  1:1562151/67‑1 (MQ=255)
                                   aaaTTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTAt                                                    >  1:149781/1‑45 (MQ=255)
                                    aaTTGAATAATTAAAACTTTGTTGCATATTCATAgg                                                            >  1:2049525/1‑36 (MQ=255)
                                    aaTTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgtttgtgt                                         >  1:1373620/1‑55 (MQ=255)
                                    aaTTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGAtatttat                            <  1:2568011/68‑1 (MQ=255)
                                     aTTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTtattat                         <  1:1950207/70‑1 (MQ=255)
                                     aTTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATg                        >  1:2111163/1‑71 (MQ=255)
                                      ttGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGc                       >  1:2535098/1‑71 (MQ=255)
                                        gAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgtt                                              <  1:1052062/46‑1 (MQ=255)
                                        gAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCg                      <  1:2990162/70‑1 (MQ=255)
                                        gAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCg                      <  1:694177/70‑1 (MQ=255)
                                        gAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCg                      <  1:1295996/70‑1 (MQ=255)
                                        gAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCg                      <  1:4231296/70‑1 (MQ=255)
                                        gAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCg                      <  1:2586948/70‑1 (MQ=255)
                                         aataatTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCg                      <  1:337254/69‑1 (MQ=255)
                                         aataatTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTa                    <  1:4410183/71‑1 (MQ=255)
                                              ttAAAACTTTGTTGCATATTCATAGGATATTTATCTGgtttgt                                           <  1:1545770/43‑1 (MQ=255)
                                              ttAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGAtatttat                            <  1:3515003/58‑1 (MQ=255)
                                               tAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGttat                                <  1:3645795/53‑1 (MQ=255)
                                                   aCTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGAtatt                               <  1:1881447/50‑1 (MQ=255)
                                                      ttgttgCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGAtatt                               >  1:4410151/1‑47 (MQ=255)
                                                      ttgttgCACATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATttatta                          >  1:1971326/1‑52 (MQ=255)
                                                        gttgCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACaa      >  1:519591/1‑70 (MQ=255)
                                                        gttgCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACaa      >  1:432680/1‑70 (MQ=255)
                                                           gCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACaaaa    <  1:589449/69‑1 (MQ=255)
                                                            cATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAAtt  <  1:4177646/70‑1 (MQ=255)
                                                              tatTCATAGGATATTTATCTGGTTTGTGTTTGTGAtatttat                            <  1:3755782/42‑1 (MQ=255)
                                                              tatTCATAGGATATTTATCTGGTTTGTGTTTGTGAtattt                              <  1:3998141/40‑1 (MQ=255)
                                                              tatTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCaaa           <  1:681632/59‑1 (MQ=255)
                                                              tatTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACaa      <  1:680803/64‑1 (MQ=255)
                                                              |                                                                   
GGTTATTATTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATT  >  W3110S.gb/389136‑389265

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: ‑