Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 641,944 G→A 11.9% E212K (GAA→AAA)  ybdR → predicted oxidoreductase, Zn‑dependent and NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb641,9440GA11.9% 157.5 / 10.2 59E212K (GAA→AAA) ybdRpredicted oxidoreductase, Zn‑dependent and NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base G (17/34);  new base A (3/4);  total (21/38)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

AAGCGTTGCAGTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACAGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGATCCC  >  W3110S.gb/641874‑642014
                                                                      |                                                                       
aaGCGTTGCAGTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCgg                                                                         <  1:2544849/71‑1 (MQ=255)
       gCAGTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACAGAt                                                                 >  1:3352780/1‑71 (MQ=255)
           tCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACAGAttttt                                                             <  1:1983004/71‑1 (MQ=255)
            cTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACAGATTTttg                                                            >  1:1374609/1‑71 (MQ=255)
             tATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACAGATTTttgt                                                           >  1:3120290/1‑71 (MQ=255)
              aTGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACa                                                                     <  1:1511969/61‑1 (MQ=255)
               tGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACAGATTTttgttg                                                         <  1:870607/71‑1 (MQ=255)
               tGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACAGATTTttgttg                                                         >  1:3620783/1‑71 (MQ=255)
               tGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACAGATTTttgttg                                                         <  1:1398869/71‑1 (MQ=255)
               tGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACAGATTTttgttg                                                         <  1:1310584/71‑1 (MQ=255)
               tGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACAGATTTttgtt                                                          <  1:2180130/70‑1 (MQ=255)
               tGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACAGATTTttgt                                                           <  1:363215/69‑1 (MQ=255)
               tGGTGCTGGTCCTGTGGGATTGTCGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACAGATTTttgttg                                                         <  1:668107/71‑1 (MQ=255)
                 gTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACAGATTTttgttgt                                                        <  1:3364558/70‑1 (MQ=255)
                 gTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACAGATTTttgttgt                                                        <  1:2686365/70‑1 (MQ=255)
                  tGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACAGAttt                                                               <  1:2429010/62‑1 (MQ=255)
                  tGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACAGATTTttgttgttg                                                      <  1:20677/71‑1 (MQ=255)
                    cTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACAGATTTttgttg                                                         <  1:2666852/66‑1 (MQ=255)
                    cTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACAGATTTTTGTTGTTGat                                                    <  1:1804336/71‑1 (MQ=255)
                    cTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACAGATTTTTGTTGTTGat                                                    >  1:3937377/1‑71 (MQ=255)
                       gTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACAGATTTttgt                                                           <  1:4075836/61‑1 (MQ=255)
                       gTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACAGATTTTTGTTGTTGatca                                                  <  1:1678242/70‑1 (MQ=255)
                            gtgGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACAGATTTttgttg                                                         >  1:4151925/1‑58 (MQ=255)
                             tgGGATTGTTGACAATCGCCTGTGGACGGTTGCTCGGTGCGGAACAGATTTTTGTTGTTGATCATCATccc                                           >  1:4285390/1‑71 (MQ=255)
                             tgGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACAGATTTTTGTTGTTGatca                                                  <  1:1709117/64‑1 (MQ=255)
                                    gttgACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACAg                                                                    <  1:2677689/40‑1 (MQ=255)
                                    gttgACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACAg                                                                    <  1:2911740/40‑1 (MQ=255)
                                       gACAATCGCCTGTGCACGGTTGCTCGGTGCGTAACAGATTTTTGTTGTTGATCATCATCACTACCGCTTg                                  >  1:2946617/1‑70 (MQ=255)
                                       gACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACAGATTTTTGTTGTTGATCATCATCCCTACCGCTTg                                  >  1:2097590/1‑70 (MQ=255)
                                        aCAATCGCCTGTGCACGGTTGCTCGGTGCGGAACAGATTTTTGTTGTTGATCATCATCCCTACCGCTTGc                                 <  1:4128438/70‑1 (MQ=255)
                                           aTCGCCTGTGCACGGTTGCTCGGTGCGGAACAGATTTttgttg                                                         <  1:395701/43‑1 (MQ=255)
                                           aTCGCCTGTGCACGGTTGCTCGGTGCGGAACAGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCAttt                             >  1:2963020/1‑71 (MQ=255)
                                            tCGCCTGTGCACGGTTGCTCGGTGCGGAACAGATTTttgttgttg                                                      <  1:1198890/45‑1 (MQ=255)
                                             cGCCTGTGCACGGTTGCTCGGTGGGGAACAGAttt                                                               <  1:259675/35‑1 (MQ=255)
                                               ccTGTGCACGGTTGCTCGGTGCGGAACAGATTTttgttg                                                         >  1:185518/1‑39 (MQ=255)
                                                cTGTGCACGGTTGCTCGGTGCGGAACAGATTTTTGTTGTTGATCATCATccc                                           <  1:1503491/52‑1 (MQ=255)
                                                 tgtgCACGGTTGCTCGGTGCGGAACAGATTTttgttgttg                                                      <  1:713989/40‑1 (MQ=255)
                                                   tgCACGGTTGCTCGGTGCGGAACAGATTTTTGTTGTTGATCATCATcc                                            >  1:2504747/1‑48 (MQ=255)
                                                   tgCACGGTTGCTCGGTGCGGAACAGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCgccgcc                      <  1:2865245/70‑1 (MQ=255)
                                                   tgCACGGTTGCTCGGTGCGGAACAGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGccgccg                     <  1:1262185/71‑1 (MQ=255)
                                                     cACGGTTGCTCGGTGCGGAACAGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGAc                   >  1:3774264/1‑71 (MQ=255)
                                                        ggTTGCTCGGTGCGGAACAGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGc                <  1:1839686/71‑1 (MQ=255)
                                                         gTTGCTCGGTGCGGAACAGATTTTTGTTGTTGATTATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCt               <  1:185090/71‑1 (MQ=255)
                                                          ttGCTCGGTGCGGAACAGATTTTTGTTGTTGATCATCATCCCTACCGCTTTCATTTCGCCGCCGACCGCTa              <  1:2549581/71‑1 (MQ=255)
                                                           tGCTCGGTGCGGAACAGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTa              <  1:1915810/70‑1 (MQ=255)
                                                            gCTCGGTGCGGAACAGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTAc             <  1:4113138/70‑1 (MQ=255)
                                                            gCTCGGTGCGGAACAGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACg            >  1:130881/1‑71 (MQ=255)
                                                             aacggTGAGAAACAGATTTTTTGTTGTTGAagaactgttttttgaag                                    >  1:2349577‑M1/3‑30 (MQ=255)
                                                               cGGTGAGAAACAGATTTTTTGTTGTTGAagaactgttttttga                                      <  1:1354047‑M1/43‑16 (MQ=255)
                                                               cGGTGAGAAACAGATTTTTTGTTGTTGAagaactgttttttg                                       <  1:2069484‑M1/42‑15 (MQ=255)
                                                               cGGTGAGAAACAGATTTTTTGTTGTTGAagaactgttttttg                                       <  1:1861547‑M1/42‑15 (MQ=255)
                                                                ggTGCGGAACAGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACGgcgc        <  1:2332878/71‑1 (MQ=255)
                                                                ggTGCGGAACAGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACGgcgc        >  1:545918/1‑71 (MQ=255)
                                                                ggTGAGAAACAGATTTTTTGTTGTTGAagaactgttttttgaa                                     >  1:2135925‑M1/1‑27 (MQ=255)
                                                                ggTGAGAAACAGATTTTTTGTTGTTGAagaactgttttttgaa                                     >  1:3912935‑M1/1‑27 (MQ=255)
                                                                ggTGAGAAACAGATTTTTTGTTGTTGAagaactgttttttg                                       <  1:2153421‑M1/41‑15 (MQ=255)
                                                                 gTGCGGAACAGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGcgcg       <  1:2490595/71‑1 (MQ=255)
                                                                  tGCGGAACAGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGa      >  1:4301808/1‑71 (MQ=255)
                                                                  tGCGGAACAGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGa      >  1:2441549/1‑71 (MQ=255)
                                                                  tGCGGAACAGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGa      >  1:2774172/1‑71 (MQ=255)
                                                                      gAACAGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGATccc  <  1:2201141/71‑1 (MQ=255)
                                                                      |                                                                       
AAGCGTTGCAGTCTATGGTGCTGGTCCTGTGGGATTGTTGACAATCGCCTGTGCACGGTTGCTCGGTGCGGAACAGATTTTTGTTGTTGATCATCATCCCTACCGCTTGCATTTCGCCGCCGACCGCTACGGCGCGATCCC  >  W3110S.gb/641874‑642014

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: