Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 286,496 T→C 27.8% 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,4960TC27.8% 52.0 / 23.8 36F162L (TTT→CTT) yagHpredicted xylosidase/arabinosidase
Reads supporting (aligned to +/- strand):  ref base T (14/12);  new base C (8/2);  total (22/14)
Fisher's exact test for biased strand distribution p-value = 2.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTACCGCCCGTGGGGGCCGCGCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTA  >  W3110S.gb/286434‑286563
                                                              |                                                                   
tCTACCGCCCGTGGGGGCCGCGCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGCTTGa                                                                 >  1:3138544/1‑67 (MQ=255)
           tGGGGGCCGCGCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCAc                                                  <  1:1888594/71‑1 (MQ=255)
              gggCCGCGCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGct                                               <  1:1758623/71‑1 (MQ=255)
              gggCCGCGCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGct                                               <  1:2378023/71‑1 (MQ=255)
              gggCCGCGCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGCTTGACCCGCAGACCGGCACGct                                               <  1:2476811/71‑1 (MQ=255)
                  cgcgCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCa                                                           >  1:3352818/1‑55 (MQ=255)
                   gcgcCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGCTTGACCCGCAGACCGGCACGCTCTCGc                                          >  1:1439920/1‑71 (MQ=255)
                   gcgcCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGCTTGACCCGCAGACCGGCACGCTCTCGc                                          >  1:428298/1‑71 (MQ=255)
                           acaGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACTCTCTCGccc                                        >  1:3130909/1‑65 (MQ=255)
                           acaGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAgcg                                   >  1:1456238/1‑70 (MQ=255)
                           acaGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAgcg                                   >  1:571090/1‑70 (MQ=255)
                            caGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGAc                                                                >  1:1670283/1‑40 (MQ=255)
                            caGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGAc                                                                >  1:2241637/1‑40 (MQ=255)
                                 aCCCGCACAACACCATCGTGTTACAGGCGCTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAAc                             >  1:744000/1‑70 (MQ=255)
                                   ccGCACAACACCATCGTGTTACAGGCGCTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCt                          <  1:3085518/71‑1 (MQ=255)
                                    cGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTg                         >  1:1866210/1‑71 (MQ=255)
                                    cGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTg                         >  1:1503835/1‑71 (MQ=255)
                                    cGCACAACACCATCGTGTTACAGGCGCTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTg                         >  1:566660/1‑71 (MQ=255)
                                       acaacaCCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGttt                      <  1:450720/71‑1 (MQ=255)
                                       acaacaCCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGttt                      <  1:2462703/71‑1 (MQ=255)
                                       acaacaCCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGttt                      <  1:1318696/71‑1 (MQ=255)
                                        caacaCCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTa                     <  1:3778290/71‑1 (MQ=255)
                                            acCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCgg                 <  1:727212/71‑1 (MQ=255)
                                              cATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCa               <  1:3095201/71‑1 (MQ=255)
                                                  gtgtTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTa                     >  1:984581/1‑61 (MQ=255)
                                                  gtgtTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTa                     >  1:1539907/1‑61 (MQ=255)
                                                  gtgtTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCACGAGCGCAAAACGCTGTTTa                     >  1:189907/1‑61 (MQ=255)
                                                   tgtTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCAcgccg          <  1:151630/71‑1 (MQ=255)
                                                   tgtTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCAcgccg          <  1:1180781/71‑1 (MQ=255)
                                                   tgtTACAGGCGTTTGACCCGCAGACCGGCACCCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCAcgccg          <  1:1113598/71‑1 (MQ=255)
                                                    gtTACAGGCGCTTGACCCGCAGACCGGCACGCTCTCg                                           >  1:3446387/1‑37 (MQ=255)
                                                    gtTACAGGCGCTTGACCCGCAGACCGGCACGCTCTCGcc                                         >  1:3177159/1‑39 (MQ=255)
                                                    gtTACAGGCGCTTGACCCGCAGACCGGCACGCTCTCGcc                                         >  1:626370/1‑39 (MQ=255)
                                                     ttACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGct        >  1:3546187/1‑71 (MQ=255)
                                                     ttACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGct        >  1:271847/1‑71 (MQ=255)
                                                           gCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTa  >  1:524902/1‑71 (MQ=255)
                                                              |                                                                   
TCTACCGCCCGTGGGGGCCGCGCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTA  >  W3110S.gb/286434‑286563

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: