Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 230,040 C→A 28.5% 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,0400CA28.5% 53.7 / 13.9 28R281L (CGC→CTC) yafCpredicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (14/6);  new base A (4/4);  total (18/10)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

AGCTCCTCCGGGAGCTGTTTTTACATGCTCGCTAAGGAAATCGATAAAAGCCCGGATGCGCGTACTTACCGCACGGTCGCTGTAATAGACGGCACTGAATGGCATTTC  >  W3110S.gb/229982‑230089
                                                          |                                                 
aGCTCCTCCGGGAGCTGTTTTTACATGCTCGCTAAGGAAATCGATAAAAGCCCGGATGCGCGTACTTACCg                                       >  1:198016/1‑71 (MQ=255)
              cTGTTTTTACATGCTCGCTAAGGAAATCGATAAAAGCCCGGATGCGCGTACTTACCGCACGGTCGCTGt                           >  1:2424133/1‑69 (MQ=255)
              cTGTTTTTACATGCTCGCTAAGGAAATCGATAAAAGCCCGGATGCGCGTACTTACCGCACGGTCGCTGTaa                         >  1:889612/1‑71 (MQ=255)
              cTGTTTTTACATGCTCGCTAAGGAAATCGATAAAAGCCCGGATGCGCGTACTTACCGCACGGTCGCTGTaa                         >  1:3041321/1‑71 (MQ=255)
              cTGTTTTTACATGCTCGCTAAGGAAATCGATAAAAGCCCGGATGCGCGTACTTACCGCACGGTCGCTGTaa                         >  1:2593620/1‑71 (MQ=255)
              cTGTTTTTACATGCTCGCTAAGGAAATCGATAAAAGCCCGGATGCGCGTACTTACCGCACGGTCGCTGTaa                         >  1:1480273/1‑71 (MQ=255)
              cTGTTTTTACATGCTCGCTAAGGAAATCGATAAAAGCCCGGATGCGCGTACTTACCGCACGGTCGCTGTaa                         >  1:2414052/1‑71 (MQ=255)
              cTGTTTTTACATGCTCGCTAAGGAAATCGATAAAAGCCCGGATGCGCGTACTTACCGCACGGTCGCTGTaa                         >  1:1676815/1‑71 (MQ=255)
              cTGTTTTTACATGCTCGCTAAGGAAATCGATAAAAGCCCGGATGCGCGTACTTACCGCACGGTCGCTGTaa                         >  1:1902373/1‑71 (MQ=255)
              cTGTTTTTACATGCTCGCTAAGGAAATCGATAAAAGCCCGGATGCGCGTACTTACCGCACGGTCGCTGTaa                         >  1:1790737/1‑71 (MQ=255)
              cTGTTTTTACATGCTCGCTAAGGAAATCGATAAAAGCCCGGATGCGCGTACTTACCGCACGGTCGCTGTa                          >  1:943582/1‑70 (MQ=255)
                 tttttACATGCTCGCTAAGGAAATCGATAAAAGCCCGGATGCGCGTACTTACCGCACGGTCGCTGTaa                         >  1:3139006/1‑68 (MQ=255)
                 tttttACATGCTCGCTAAGGAAATCGATAAAAGCCCGGATGCGCGTACTTACCGCACGGTCGCTGTaa                         >  1:2810734/1‑68 (MQ=255)
                      aCATGCTCGCTAAGGAAATCGATAAAAGCCCGGATGAGCGTACTTACCGCACGGTCGCTGTAATAGACgg                  <  1:104491/70‑1 (MQ=255)
                           cTCGCTAAGGAAATCGATAAAAGCCCGGATGAGCGTACTTACCGCACGGTCGCTGTAATAGACGGCACTg             >  1:1377758/1‑70 (MQ=255)
                           cTCGCTAAGGAAATCGATAAAAGCCCGGATGAGCGTACTTACCGCACGGTCGCTGTAATAGACGGCACTg             >  1:1252751/1‑70 (MQ=255)
                           cTCGCTAAGGAAATCGATAAAAGCCCGGATGAGCGTACTTACCGCACGGTCGCTGTAATAGACGGCACTg             >  1:934709/1‑70 (MQ=255)
                           cTCGCTAAGGAAATCGATAAAAGCCCGGATGAGCGTACTAACCGCACGGTCGCTGTAAt                        >  1:2475762/1‑59 (MQ=255)
                               cTAAGGAAATCGATAAAAGCCCGGATGCGCGTACTTACCGCACGGTCGCTGTAATAGACGGCACTGAATgg        <  1:1503330/71‑1 (MQ=255)
                               cTAAGGAAATCGATAAAAGCCCGGATGCGCGTACTTACCGCACGGTCGCTGTAATAGACGGCACTGAATgg        <  1:2892814/71‑1 (MQ=255)
                               cTAAGGAAATCGATAAAAGCCCGGATGCGCGTACTTACCGCACGGTCGCTGTAATAGACGGCACTGAATgg        <  1:3263400/71‑1 (MQ=255)
                               cTAAGGAAATCGATAAAAGCCCGGATGCGCGTACTTACCGCACGGTCGCTGTAATAGACGGCACTGAATgg        <  1:734657/71‑1 (MQ=255)
                               cTAAGGAAATCGATAAAAGCCCGGATGCGCGTACTTACCGCACGGTCGCTGTAATAGACGGCACTGAATgg        <  1:1736187/71‑1 (MQ=255)
                                tAAGGAAATCGATAAAAGCCCGGATGCGCGTACTTACCGCACGGTCGCTGTAATAGACGGCACTGAATgg        <  1:1625301/70‑1 (MQ=255)
                                     aaaTCGATAAAAGCCCGGATGAGCGTACTTACCGCACGGTCGCTGTAATAGACGGCACTGAATGGCATTTc  <  1:2550166/71‑1 (MQ=255)
                                     aaaTCGATAAAAGCCCGGATGAGCGTACTTACCGCACGGTCGCTGTAATAGACGGCACTGAATGGCATTTc  <  1:234360/71‑1 (MQ=255)
                                     aaaTCGATAAAAGCCCGGATGAGCGTACTTACCGCACGGTCGCTGTAATAGACGGCACTGAATGGCATTTc  <  1:1124052/71‑1 (MQ=255)
                                             aaaaGCCCGGATGCGCGTACTTACCGCACGGTCGCTGTAATAGa                     >  1:1084508/1‑44 (MQ=255)
                                                          |                                                 
AGCTCCTCCGGGAGCTGTTTTTACATGCTCGCTAAGGAAATCGATAAAAGCCCGGATGCGCGTACTTACCGCACGGTCGCTGTAATAGACGGCACTGAATGGCATTTC  >  W3110S.gb/229982‑230089

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: