Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,567,710 T→C 35.0% R102R (CGT→CGC yihQ → alpha‑glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,567,7100TC35.0% 14.9 / 16.5 20R102R (CGT→CGCyihQalpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base T (8/5);  new base C (4/3);  total (12/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCAGCCGTGGTTCTGACATTAGCGCCACGCTGAATATCTCTGCCGACGATCAGGGGCGTTTATTGCTGGAACTACAAAACGACAACCTTAACCACAACCGTATCTGGCTGCGCCTTGCCGCTCA  >  W3110S.gb/3567651‑3567775
                                                           |                                                                 
ttCAGCCGTGGTTCTGACATTAGCGCCACGCTGAATATCTCTGCCGACGATCAGGGGCGtt                                                                  <  1:1827303/61‑1 (MQ=255)
 tCAGCCGTGGTTCTGACATTAGCGCCACGCTGAATATCTCTGCCGACGATCAGGGGCGTTTATTGCTgg                                                         >  1:807165/1‑69 (MQ=255)
 tCAGCCGTGGTTCTGACATTAGCGCCACGCTGAATATCTCTGCCGACGATCAGGGGCGTTTATTGCTgg                                                         >  1:3565610/1‑69 (MQ=255)
 tCAGCCGTGGTTCTGACATTAGCGCCACGCTGAATATCTCTGCCGACGATCAGGGGCGTTTATTGCTgg                                                         >  1:3411860/1‑69 (MQ=255)
 tCAGCCGTGGTTCTGACATTAGCGCCACGCTGAATATCTCTGCCGACGATCAGGGGCGTTTATTGCTg                                                          >  1:3191191/1‑68 (MQ=255)
            tCTGACATTAGCGCCACGCTGAATATCTCTGCCGACGATCAGGGGCGTTTATTGCTGGAACTACAAAAc                                              <  1:19042/69‑1 (MQ=255)
            tCTGACATTAGCGCCACGCTGAATATCTCTGCCGACGATCAGGGGCGTTTATTGCTGGAACTACAAAAc                                              <  1:23465/69‑1 (MQ=255)
             cTGACATTAGCGCCACGCTGAATATCTCTGCCGACGATCAGGGGCGtt                                                                  <  1:1589722/48‑1 (MQ=255)
             cTGACATTAGCGCCACGCTGAATATCTCTGCCGACGATCAGGGGCGtt                                                                  >  1:1600657/1‑48 (MQ=255)
             cTGACATTAGCGCCACGCTGAATATCTCTGCCGACGATCAGGGGCGtt                                                                  >  1:3453866/1‑48 (MQ=255)
             cTGACATTAGCGCCACGCTGAATATCTCTGCCGACGATCAGGGGCGCTTATTGCTGGAACTACAAAAc                                              >  1:1311267/1‑68 (MQ=255)
             cTGACATTAGCGCCACGCTGAATATCTCTGCCGACGATCAGGGGCGCTTATTGCTGGAACTACAAAAc                                              >  1:2789333/1‑68 (MQ=255)
                  aTTAGCGCCACGCTGAATATCTCTGCCGACGATCAGGGGCGCTTATTGCTGGAACTACAAAACGACaa                                         <  1:3429579/68‑1 (MQ=255)
                     aGCGCCACGCTGAATATCTCTGCCGACGATCAGGGGCGTTTATTGCTGGAACTACAAAACGACAACCtt                                     <  1:2209319/69‑1 (MQ=255)
                                   tatCTCTGCCGACGATCAGGGGCGTTTATTGCTGGAACTACAAAACGACAACCtt                                     <  1:3312738/55‑1 (MQ=255)
                                   tatCTCTGCCGACGATCAGGGGCGCTTATTGCTGGAACTACAAAACGACAACCTTAACCACAAGCGt                         <  1:2227041/67‑1 (MQ=255)
                                     tctcTGCCGACGATCAGGGGCGTTTATTGCTGGAACTACAAAACGACAAc                                        >  1:28821/1‑50 (MQ=255)
                                     tctcTGCCGACGATCAGGGGCGTTTATTGCTGGAACTACAAAACGACAAc                                        >  1:1848970/1‑50 (MQ=255)
                                       tctGCCGACGATCAGGGGCGTTTATTGCTGGAACTACAAAACGACAACCTTAACCACAACCGTATCTg                    <  1:3232321/68‑1 (MQ=255)
                                                   cAGGGGCGTTTATTGCTGGAACTACAAAACGACAACCTTAACCACAACCGTATCTGGCTGCGCCTTGcc       >  1:988929/1‑69 (MQ=255)
                                                   cAGGGGCGTTTATTGCTGGAACTACAAAACGACAACCTTAACCACAACCGTATCTGGCTGCGCCTTGc        >  1:2689719/1‑68 (MQ=255)
                                                       ggCGCTTATTGCTGGAACTACAAAACGACAACCTTAACCACAACCGTATCTGGCTGCGCCTTGCCGCt    >  1:1234312/1‑68 (MQ=255)
                                                       ggCGCTTATTGCTGGAACTACAAAACGACAACCTTAACCACAACCGTATCTGGCTGCGCCTTGCCGCTc   <  1:2673722/69‑1 (MQ=255)
                                                       ggCGCTTATTGCTGGAACTACAAAACGACAACCTTAACCACAACCGTATCTGGCTGCGCCTTGCCGCTc   >  1:67390/1‑69 (MQ=255)
                                                           tttGTTGCTGGAACTACAAAACGACAACCTTAACCACAACCGTATCTGGCTGCGCCTTGCCGCTCa  >  1:989506/1‑66 (MQ=255)
                                                           |                                                                 
TTCAGCCGTGGTTCTGACATTAGCGCCACGCTGAATATCTCTGCCGACGATCAGGGGCGTTTATTGCTGGAACTACAAAACGACAACCTTAACCACAACCGTATCTGGCTGCGCCTTGCCGCTCA  >  W3110S.gb/3567651‑3567775

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: