Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 230,040 C→A R281L (CGC→CTC)  yafC ← predicted DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb230,0400CA93.9% 85.0 / 1.0 33R281L (CGC→CTC) yafCpredicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base A (13/18);  total (15/18)
Fisher's exact test for biased strand distribution p-value = 1.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

TGACAGCTCCTCCGGGAGCTGTTTTTACATGCTCGCTAAGGAAATCGATAAAAGCCCGGATGCGCGTACTTACCGCACGGTCGCTGTAATAGACGGCACTGAATGGCATTTCCACTGGCAACACTTTATCT  >  W3110S.gb/229978‑230108
                                                              |                                                                    
tGACAGCTCCTCCGGGAGCTGTTTTTACATGCTCGCTAAGGAAATCGATAAAAGCCCGGATGAGCGTa                                                                 <  1:1872123/68‑1 (MQ=255)
  aCAGCTCCTCCGGGAGCTGTTTTTACATGCTCGCTAAGGAAATCGATAAAAGCCCGGATGAGCGTACTTa                                                             <  1:1325228/70‑1 (MQ=255)
  aCAGCTCCTCCGGGAGCTGTTTTTACATGCTCGCTAAGGAAATCGATAAAAGCCCGGATGAGCGTACTTa                                                             <  1:1349200/70‑1 (MQ=255)
  aCAGCTCCTCCGGGAGCTGTTTTTACATGCTCGCTAAGGAAATCGATAAAAGCCCGGATGAGCGTACTTa                                                             <  1:210598/70‑1 (MQ=255)
                aGCTGTTTTTACATGCTCGCTAAGGAAATCGATAAAAGCCCGGATGAGCGTACTTACCGCACGGTCGCTg                                               >  1:234968/1‑70 (MQ=255)
                aGCTGTTTTTACATGCTCGCTAAGGAAATCGATAAAAGCCCGGATGAGCGTACTTACCGCACGGTCGCTg                                               >  1:1075365/1‑70 (MQ=255)
                aGCTGTTTTTACATGCTCGCTAAGGAAATCGATAAAAGCCCGGATGAGCGTACTTACCGCACGGTCGCTg                                               >  1:3659519/1‑70 (MQ=255)
                aGCTGTTTTTACATGCTCGCTAAGGAAATCGATAAAAGCCCGGATGAGCGTACTTACCGCACGGTCGCTg                                               >  1:3788155/1‑70 (MQ=255)
                aGCTGTTTTTACATGCTCGCTAAGGAAATCGATAAAAGCCCGGATGAGCGTACTTACCGCACGGTCGCTg                                               >  1:2039201/1‑70 (MQ=255)
                aGCTGTTTTTACATGCTCGCTAAGGAAATCGATAAAAGCCCGGATGAGCGTACTTACCGCACGGTCGCTGt                                              <  1:268356/71‑1 (MQ=255)
                aGCTGTTTTTACATGCTCGCTAAGGAAATCGATAAAAGCCCGGATGAGCGTACTTACCGCACGGTCGCTGt                                              <  1:2518865/71‑1 (MQ=255)
                aGCTGTTTTTACATGCTCGCTAAGGAAATCGATAAAAGCCCGGATGAGCGTACTTACCGCACGGTCGCTGt                                              <  1:1944577/71‑1 (MQ=255)
                aGCTGTTTTTACATGCTCGCTAAGGAAATCGATAAAAGCCCGGATGAGCGTACTTACCGCACGGTCGCTGt                                              <  1:1035937/71‑1 (MQ=255)
                aGCTGTTTTTACATGCTCGCTAAGGAAATCGATAAAAGCCCGGATGAGCGTACTTACCGCACGGTCGCTGt                                              <  1:1581589/71‑1 (MQ=255)
                aGCTGTTTTTACATGCTCGCTAAGGAAATCGATAAAAGCCCGGATGAGCGTACTTACCGCACGGTCGCTGt                                              <  1:2582043/71‑1 (MQ=255)
                aGCTGATTTTACATGCTCGCTAAGGAAATCGATAAAAGCCCGGATGAGCGTACTTACCGCACGGACGCtc                                               >  1:1407371/1‑69 (MQ=255)
                               cTCGCTAAGGAAATCGATAAAAGCCCGGATGAGCGTACTTAc                                                            >  1:1068307/1‑42 (MQ=255)
                               cTCGCTAAGGAAATCGATAAAAGCCCGGATGAGCGTACTTAc                                                            >  1:2619039/1‑42 (MQ=255)
                               cTCGCTAAGGAAATCGATAAAAGCCCGGATGAGCGTACTTACCGGACGGTCGCTGTAATAGACGGCACTg                                <  1:3831630/70‑1 (MQ=255)
                               cTCGCTAAGGAAATCGATAAAAGCCCGGATGAGCGTACTTACAGCACGGTCGCTGTAATAGACGGCACTg                                <  1:868738/70‑1 (MQ=255)
                                    tAAGGAAATCGATAAAAGCCCGGATGAGCGTACTTACCGCACGGTCGCTGTAATAGACGGCACTg                                <  1:2624974/65‑1 (MQ=255)
                                    tAAGGAAATCGATAAAAGCCCGGATGAGCGTACTTACCGCACGGTCGCTGTAATAGACGGCACTg                                <  1:3483759/65‑1 (MQ=255)
                                    tAAGGAAATCGATAAAAGCCCGGATGAGCGTACTTACCGCACGGTCGCTGTAATAGACGGCACTg                                <  1:3676545/65‑1 (MQ=255)
                                    tAAGGAAATCGATAAAAGCCCGGATGAGCGTACTTACCGCACGGTCGCTGTAATAGACGGCACTg                                <  1:4059543/65‑1 (MQ=255)
                                    tAAGGAAATCGATAAAAGCCCGGATGAGCGTACTTACCGCACGGTCGCTGTAATAGACGGCACTg                                <  1:1733993/65‑1 (MQ=255)
                                      aGGAAATCGATAAAAGCCCGGATGAGCGTACTTACCGCACGGTCGCTGTAATAGACGGCACTGaa                              >  1:251500/1‑65 (MQ=255)
                                      aGGAAATCGATAAAAGCCCGGATGAGCGTACTTACCGCACGGTCGCTGTAATAGACGGCACTGaa                              >  1:3675826/1‑65 (MQ=255)
                                      aGGAAATCGATAAAAGCCCGGATGAGCGTACTTACCGCACGGTCGCTGTAATAGACGGCACTGaa                              >  1:1965808/1‑65 (MQ=255)
                                      aGGAAATCGATAAAAGCCCGGATGAGCGTACTTACCGCACGGTCGCTGTAATAGACGGCACTGaa                              >  1:497870/1‑65 (MQ=255)
                                       ggAAATCGATAAAAGCCCGGATGAGCGTACTTACCGCACGGTCGCTg                                               >  1:2632077/1‑47 (MQ=255)
                                       ggAAATCGATAAAAGCCCGGATGAGCGTACTTACCGCACGGTCGCTGTAATAGACGGCACTg                                <  1:233/62‑1 (MQ=255)
                                                            tGCGCGTACTTACCGCACGGTCGCTGTAATAGACGGCACTGAATGGCATTTCCACTGGCAACACTTTATCt  >  1:2498046/1‑71 (MQ=255)
                                                            tGCGCGTACTTACCGCACGGTCGCTGTAATAGACGGCACTGAATGGCATTTCCACTGGCAACACTTTATCt  >  1:2296268/1‑71 (MQ=255)
                                                              |                                                                    
TGACAGCTCCTCCGGGAGCTGTTTTTACATGCTCGCTAAGGAAATCGATAAAAGCCCGGATGCGCGTACTTACCGCACGGTCGCTGTAATAGACGGCACTGAATGGCATTTCCACTGGCAACACTTTATCT  >  W3110S.gb/229978‑230108

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: