Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,077,980 C→T 69.2% Y355Y (TAC→TAT yeeR → predicted membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,077,9800CT69.2% 29.4 / 19.3 26Y355Y (TAC→TATyeeRpredicted membrane protein
Reads supporting (aligned to +/- strand):  ref base C (5/3);  new base T (9/9);  total (14/12)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTCCTGCGGGCTGACAGAAGTACATCATGCCCCGGAAAATCATAAGCCGGATCCGGAATGGTACTGCTCCTCTCTTTGTCGCGAAACAGAAACACTGTGTCAGGAAATTTATGAACGCCCTTACAACAGC  >  W3110S.gb/2077916‑2078046
                                                                |                                                                  
cctcctGCGGGCTGACAGAAGTACATCATGCCCCGGAAAATCATAAGCCGGATCCGGAATGGTATTGct                                                                >  1:291895/1‑69 (MQ=255)
cctcctGCGGGCTGACAGAAGTACATCATGCCCCGGAAAATCATAAGCCGGATCCGGAATGGTATTGct                                                                >  1:928023/1‑69 (MQ=255)
  tcctGCGGGCTGACAGAAGTACATCATGCCCCGGAAAATCATAAGCCGGATCCGGAATGGTATTGctc                                                               <  1:346830/68‑1 (MQ=255)
              aCAGAAGTACATCATGCCCCGGAAAATCATAAGCCGGATCCGGAATGGTATTGCTCCTCTCTTTGTcg                                                   >  1:3315660/1‑68 (MQ=255)
                   aGTACATCATGCCCCGGAAAATCATAAGCCGGATCCGGAATGGTATTGCTCCTCTCTTTGTCGCga                                                >  1:2343458/1‑66 (MQ=255)
                      aCATCATGCCCCGGAAAATCATAAGCCGGATCCGGAATGGTACTGCTCCTCTCTTTGTCGCgaaacaga                                          <  1:1193003/69‑1 (MQ=255)
                      aCATCATGCCCCGGAAAATCATAAGCCGGATCCGGAATGGTACTGCTCCTCTCTTTGTCGCgaaacaga                                          <  1:2213759/69‑1 (MQ=255)
                      aCATCATGCCCCGGAAAATCATAAGCCGGATCCGGAATGGTACTGCTCCTCTCTTTGTCGCgaaacaga                                          <  1:542930/69‑1 (MQ=255)
                           atGCCCCGGAAAATCATAAGCCGGATCCGGAATGGTACTGCTCCTCTCTTTGTCGCGAAACAGAAacac                                     >  1:2776941/1‑69 (MQ=255)
                            tGCCCCGGAAAATCATAAGCCGGATCCGGAATGGTATTGCTCCTCTCTTTGTCGCgaaacagaaa                                        <  1:2473229/65‑1 (MQ=255)
                              ccccGGAAAATCATAAGCCGGATCCGGAATGGTATTGCTCCTCTCTTTGTcgcg                                                 <  1:197019/54‑1 (MQ=255)
                              ccccGGAAAATCATAAGCCGGATCCGGAATGGTATTGCTCCTCTCTTTGTcgcg                                                 <  1:2486312/54‑1 (MQ=255)
                              ccccGGAAAATCATAAGCCGGATCCGGAATGGTATTGCTCCTCTCTTTGTcgcg                                                 <  1:112504/54‑1 (MQ=255)
                              ccccGGAAAATCATAAGCCGGATCCGGAATGGTATTGCTCCTCTCTTTGTcgcg                                                 <  1:301179/54‑1 (MQ=255)
                                     aaaTCATAAGCCGGATCCGGAATGGTATTGCTCCTCTCTTTGTCGCGAAACAGAAacac                                     >  1:2704656/1‑59 (MQ=255)
                                        tCATAAGCCGGATCCGGAATGGTATTGCTCCTCTCTTTGTCGCgaaacaga                                          <  1:697145/51‑1 (MQ=255)
                                        tCATAAGCCGGATCCGGAATGGTATTGCTCCTCTCTTTGTCGCgaaacaga                                          <  1:1841521/51‑1 (MQ=255)
                                        tCATAAGCCGGATCCGGAATGGTATTGCTCCTCTCTTTGTCGCga                                                >  1:2985731/1‑45 (MQ=255)
                                             aGCCGGATCCGGAATGGTATTGCTCCTCTCTTTGTCGCGAAACAGAAACACTGTGTCAGGAAATTTATg                   <  1:450297/69‑1 (MQ=255)
                                             aGCCGGATCCGGAATGGTATTGCTCCTCTCTTTGTCGCGAAACAGAAACACTGTGTCAGGAAATTTATg                   >  1:2163914/1‑69 (MQ=255)
                                                     ccGGAATGGTATTGCTCCTCTCTTTGTCGCGAAACAGAAACACTGTGTCAGGAAATTTATGAACGCCCt           >  1:3608680/1‑69 (MQ=255)
                                                     ccGGAATGGTATTCCTCCTCTCTGTGTCGCgaaaca                                            >  1:40484/1‑36 (MQ=25)
                                                     ccGGAATGGTACTGCTCCTCTCTTTGTCGCGAAACAGAAACACTGTGTCAGGAAATTTATGAACGCCCt           >  1:2245211/1‑69 (MQ=255)
                                                              tACTGCTCCTCTCTTTGTCGCGAAACAGAAACACTGTGTCAGGAAATTTATGAACGCCCTTACAACAGc  >  1:3199168/1‑69 (MQ=255)
                                                              tACTGCTCCTCTCTTTGTCGCGAAACAGAAACACTGTGTCAGGAAATTTATGAACGCCCTTACAACAGc  >  1:670300/1‑69 (MQ=255)
                                                              tACTGCTCCTCTCTTTGTCGCGAAACAGAAACACTGTGTCAGGAAATTTATGAACGCCCTTACAACAGc  >  1:1858449/1‑69 (MQ=255)
                                                                |                                                                  
CCTCCTGCGGGCTGACAGAAGTACATCATGCCCCGGAAAATCATAAGCCGGATCCGGAATGGTACTGCTCCTCTCTTTGTCGCGAAACAGAAACACTGTGTCAGGAAATTTATGAACGCCCTTACAACAGC  >  W3110S.gb/2077916‑2078046

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: