Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 286,496 T→C 15.1% F162L (TTT→CTT)  yagH → predicted xylosidase/arabinosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb286,4960TC15.1% 76.6 / 7.7 33F162L (TTT→CTT) yagHpredicted xylosidase/arabinosidase
Reads supporting (aligned to +/- strand):  ref base T (15/13);  new base C (2/3);  total (17/16)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01

TACTATATCTACCGCCCGTGGGGGCCGCGCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACA  >  W3110S.gb/286427‑286565
                                                                     |                                                                     
tACTATATCTACCGCCCGTGGGGGCCGCGCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGtt                                                                      >  1:3348920/1‑71 (MQ=255)
  cTATATCTACCGCCCGTGGGGGCCGCGCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGttt                                                                     <  1:1936027/70‑1 (MQ=255)
    atatCTACCGCCCGTGGGGGCCGCGCCACCACAGCAACCCGCACAACATCATCGTGTTACAGGCGTTTGAc                                                                  >  1:2958473/1‑71 (MQ=255)
       tCTACCGCCCGTGGGGGCCGCGCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCg                                                               >  1:284745/1‑71 (MQ=255)
          aCCGCCCGTGGGGGCCGCGCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCg                                                               >  1:150122/1‑68 (MQ=255)
                    ggggCCGCGCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGCTTGACCCGCAGACCGGCACg                                                   <  1:4247043/70‑1 (MQ=255)
                     gggCCGCGCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGAACGct                                                 <  1:2648254/71‑1 (MQ=255)
                        ccGCGCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTg                                                                    <  1:4436849/49‑1 (MQ=255)
                                ccacAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCa                                                             <  1:313095/48‑1 (MQ=255)
                                ccacAGCAACCCGCACAACACCATCGTGTTACAGGCGCTTg                                                                    <  1:1220441/41‑1 (MQ=255)
                                ccacAGCAACCCGCACAACACCATCGTGTTACAGGCGCTTg                                                                    <  1:2712778/41‑1 (MQ=255)
                                 cacaGCAACCCGCACAACACCATCGTGTCACAGGCGTTTGACCCGCAGACCGGCACGctct                                               <  1:657280/61‑1 (MQ=255)
                                  acaGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGAc                                                                  >  1:2278680/1‑41 (MQ=255)
                                  acaGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGAc                                                                  >  1:4005176/1‑41 (MQ=255)
                                   caGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAg                                                            <  1:3378054/46‑1 (MQ=255)
                                        aCCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCACCCGAGCGCAAAACg                              <  1:1921013/71‑1 (MQ=255)
                                          ccGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGct                                                 <  1:2316503/50‑1 (MQ=255)
                                          ccGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCt                            <  1:3278409/71‑1 (MQ=255)
                                           cGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTg                           <  1:3169384/71‑1 (MQ=255)
                                            gCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGt                          >  1:2791884/1‑71 (MQ=255)
                                             cacaACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGtt                         <  1:3461146/71‑1 (MQ=255)
                                               caacaCCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGctc                                                >  1:2943850/1‑46 (MQ=255)
                                               caacaCCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTa                       <  1:1180299/71‑1 (MQ=255)
                                                   acCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCgg                   >  1:1436740/1‑71 (MQ=255)
                                                       tCGTGTTACAGGCGCTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTa                       >  1:3193112/1‑63 (MQ=255)
                                                       tCGTGTTACAGGCGCTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTAc                      >  1:1484357/1‑64 (MQ=255)
                                                         gtgtTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCAcgcc             >  1:1965404/1‑71 (MQ=255)
                                                         gtgtTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCAcgc              >  1:864921/1‑70 (MQ=255)
                                                              aCAGGCGTTTGACCCGCAGCCCGGCACGCTCTCGCCCGAGCGCaaa                                 >  1:1582898/1‑46 (MQ=255)
                                                              aCAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCaaa                                 >  1:316716/1‑46 (MQ=255)
                                                              aCAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGctc         <  1:2465646/70‑1 (MQ=255)
                                                              aCAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGct          >  1:2201951/1‑69 (MQ=255)
                                                                   cGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTa                       >  1:4553216/1‑51 (MQ=255)
                                                                   cGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTac   <  1:2838833/71‑1 (MQ=255)
                                                                    gTTTGACCGGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTaca  >  1:4226990/1‑71 (MQ=255)
                                                                     |                                                                     
TACTATATCTACCGCCCGTGGGGGCCGCGCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACA  >  W3110S.gb/286427‑286565

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: