Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,427,890 T→G 100% I246S (ATT→AGT)  yfcH → conserved hypothetical protein with NAD(P)‑binding Rossmann‑fold domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,427,8900TG82.7% 65.9 / 2.7 23I246S (ATT→AGT) yfcHconserved hypothetical protein with NAD(P)‑binding Rossmann‑fold domain
Reads supporting (aligned to +/- strand):  ref base T (3/1);  new base G (4/15);  total (7/16)
Fisher's exact test for biased strand distribution p-value = 6.72e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CCCGGTACGCAATGAACAATTTGCCCATGCGCTCGGTCATGCGCTGCATCGCCCGGCCATTTTGCGCGTCCCCGCGACCGCCATTCGACTGTTAATGGGCGAATCTTCAGTACT  >  W3110S.gb/2427831‑2427944
                                                           |                                                      
ccgggTACGCAATGAACAATTTGCCCATGCGCTCGGTCATGCGCTGCATCGCCCGGCCATTTTGCGCGt                                               >  1:1593599/4‑69 (MQ=255)
cccGGTACGCAATGAACAATTTGCCCATGCACTCGGTCATGCGCTGCATCGCCCGGCCAGTTTGcgcg                                                >  1:1968716/1‑68 (MQ=255)
            tGAACAATTTGCCCATGCGCTCGGTCATGCGCTGCATCGCCCGGCCAGTTTGCGCGTCCCCGCGACCGc                                   <  1:2094257/69‑1 (MQ=255)
                 aaTTTGCCCATGCGCTCGGTCATGCGCTGCATCGCCCGGCCATTTTGCGCGTCCCCGCGACCGCCATTc                              >  1:2580196/1‑69 (MQ=255)
                 aaTTTGCCCATGCGCTCGGTCATGCGCTGCATCGCCCAGCCATTTTGCGCGTCCCCGCGACCGCCATTc                              >  1:668393/1‑69 (MQ=255)
                         cATGCGCTCGGTCATGCGCTGCATCGCCCGGCCAGTTTGCGCGTCCGCGCGACCGCCATTCGACTGtt                       <  1:2925475/68‑1 (MQ=255)
                         cATGCGCTCGGTCATGCGCTGCATCGCCCGGCCAGTTTGCGCGTCCCCGCGACCGCCATTCGACTGtt                       <  1:2674304/68‑1 (MQ=255)
                         cATGCGCTCGGTCATGCGCTGCATCGCCCGGCCAGTTTGCGCGTCCCCGCGACCGCCATTCGACTGtt                       <  1:1199300/68‑1 (MQ=255)
                         cATGCGCTCGGTCATGCGCTGCATCGCCCGGCCAGTTTGCGCGTCCCCGCGACCGCCATTCGACTGTTa                      <  1:195441/69‑1 (MQ=255)
                          aTGCGCTCGGTCATGCGCTGCATCGCCCGGCCAGTTTGCGCGTCCCCGCGACCGCCATTCGACTGTTaa                     <  1:629788/69‑1 (MQ=255)
                          aTGCGCTCGGTCATGCGCTGCATCGCCCGGCCAGTTTGCGCGTCCCCGCGACCGCCATTCGACTGTTaa                     <  1:253016/69‑1 (MQ=255)
                             cgcTCGGTCATGCGCTGCATCGCCCGGCCAGTTTGCGCGTCCCCGCGACCGCCATTCGACTGtt                       <  1:1412144/64‑1 (MQ=255)
                             cgcTCGGTCATGCGCTGCATCGCCCGGCCAGTTTGCGCGTCCCCGCGACCGCCATTCGACTGtt                       <  1:1377392/64‑1 (MQ=255)
                             cgcTCGGTCATGCGCTGCATCGCCCGGCCAGTTTGCGCGTCCCCGCGACCGCCATTCGACTGTTAATgg                  >  1:2047107/1‑69 (MQ=255)
                                     cATGCGCTGCATCGCCCGGCCAGTTTGCGCGTCCCCGCGACCGCCATTCGACTGTTAATGGGCGAATc           >  1:1001002/1‑68 (MQ=255)
                                     cATGCGCTGCATCGCCCGGCCAGTTTGCGCGTCCCCGCGACCGCCATTCGACTGTTAATGGGCGAATc           >  1:2560499/1‑68 (MQ=255)
                                     cATGCGCTGCATCGCCCGGCCAGTTTGCGCGTCCCCGCGACCGCCATTCGACTTTAATGGGCGAATCtt         <  1:1586715/69‑1 (MQ=255)
                                      aTGCGCTGCATCGCCCGGCCAGTTTGCGCGTCCCCGCGACCGCCTTTCGACTGTTAATGGGCGAATCtt         <  1:650643/69‑1 (MQ=255)
                                      aTGCGCTGCATCGCCCGGCCAGTTTGCGCGTCCCCGCGACCGCCATTCGACTGTTAATGGGCGAATCtt         <  1:1792446/69‑1 (MQ=255)
                                      aTGCGCTGCATCGCCCGGCCAGTTTGCGCGTCCCCGCGACCGCCATTCGACTGTTAATGGGCGAATCtt         <  1:3318522/69‑1 (MQ=255)
                                      aTGCGCTGCATCGCCCGGCCAGTTTGCGCGTCCCCGCGACCGCCATTCGACTGTTAATGGGCGAATCtt         <  1:1609678/69‑1 (MQ=255)
                                         cgcTGCATCGCCCGGCCATTTTGCGCGTCCCCGCGACCGCCATTCGACTGTTAATGGGCGAATCTTCAg      <  1:567517/69‑1 (MQ=255)
                                             gCATCGCCCGGCCAGTTTGCGCGTCCCCGCGACCGCCATTCGACTGTTAATGGGCGAATCTTCAGTACt  <  1:325732/69‑1 (MQ=255)
                                                           |                                                      
CCCGGTACGCAATGAACAATTTGCCCATGCGCTCGGTCATGCGCTGCATCGCCCGGCCATTTTGCGCGTCCCCGCGACCGCCATTCGACTGTTAATGGGCGAATCTTCAGTACT  >  W3110S.gb/2427831‑2427944

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: