Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,184,035 C→A 84.6% G3* (GGA→TGA) 
R269M (AGG→ATG) 
gspD ←
gspC ←
general secretory pathway component, cryptic
general secretory pathway component, cryptic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,184,0350CA84.6% 66.3 / 14.5 39G3* (GGA→TGA) 
R269M (AGG→ATG) 
gspD
gspC
general secretory pathway component, cryptic
general secretory pathway component, cryptic
Reads supporting (aligned to +/- strand):  ref base C (3/3);  new base A (16/17);  total (19/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

AGCGTGTCCTGCACAAGGCATGAGTAGTGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGTAAGTTCCCCGACGGAAACATTTATCAAGCGGGGTACGCC  >  W3110S.gb/4183973‑4184102
                                                              |                                                                   
aGCGTGTCCTGCACAAGGCATGAGTAGTGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCATTTCAt                                                               <  1:1504436/69‑1 (MQ=255)
    tgtCCTGCACAAGGCATGAGTAGTGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCATTTCAt                                                               <  1:389504/65‑1 (MQ=255)
    tgtCCTGCACAAGGCATGAGTAGTGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCATTTCATGACGCa                                                         <  1:1798826/71‑1 (MQ=255)
           cacaAGGCATGAGTAGTGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCCTTTCATGACGCAGTCcatt                                                  >  1:2253851/1‑71 (MQ=255)
           cacaAGGCATGAGTAGTGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCATTTCATGACGCAGTCcatt                                                  <  1:948976/71‑1 (MQ=255)
              aaGGCATGAGTAGTGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCAttt                                                                  <  1:436775/52‑1 (MQ=255)
                ggCATGAGTAGTGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCATTTCa                                                                >  1:122500/1‑52 (MQ=255)
                    tGAGTAGTGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCATTTCATGACGCAGTCCATTCATTCCTg                                          >  1:2319639/1‑70 (MQ=255)
                      agtagtGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCCTTTCATGACGCAGTCCattcattc                                             >  1:1061858/1‑65 (MQ=255)
                      agtagtGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCATTTCATGACGCAGTCca                                                    >  1:846580/1‑58 (MQ=255)
                      agtagtGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCATTTCATGACGCAGTCca                                                    >  1:1899451/1‑58 (MQ=255)
                      agtagtGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCATTTCATGACGCAGTCCATTCATTCCTGt                                         >  1:125958/1‑69 (MQ=255)
                       gtagtGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGTAAg                                      <  1:1567945/71‑1 (MQ=255)
                          gtGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCATTTCATGACg                                                           <  1:2304045/47‑1 (MQ=255)
                           tgctgcCAGCAAGCAGCAGGTGATTTTATTGAGTCATTTCATGACGCAGTCCATTCATTCCTGTAAGtt                                    <  1:1638466/69‑1 (MQ=255)
                             ctgcCAGCAAGCAGCAGGTGATTTTATTGAGTCATTTCATGACGCAGTCCATTCATTCCTGTAAGTTcccc                                >  1:1148691/1‑71 (MQ=255)
                             ctgcCAGCAAGCAGCAGGTGATTTTATTGAGTCATTTCATGACGCAGTCCATTCATTCCTGTAAGTTcccc                                <  1:1165596/71‑1 (MQ=255)
                                  agcaagcaGCAGGTGATTTTATTGAGTCATTTCATGACGCAGTCCATTCATTCCTGTAAGtt                                    >  1:2039713/1‑62 (MQ=255)
                                  agcaagcaGCAGGTGATTTTATTGAGTCATTTCATGACGCAGTCCATTCATTCCTGTAAGTTcccc                                <  1:1455201/66‑1 (MQ=255)
                                     aagcaGCAGGTGATTTTATTGAGTCATTTCATGACGCAGTCCATTCATTCCTGTAAGTTCCCCGACGGaaa                        <  1:2334358/71‑1 (MQ=255)
                                     aagcaGCAGGTGATTTTATTGAGTCATTTCATGACGCAGTCCATTCATTCCTGTAAGTTCCCCGACGGaaa                        <  1:485712/71‑1 (MQ=255)
                                          gcagTTGATTTTATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGTAAGTTCCCCGACGGAAACAttt                   <  1:302223/71‑1 (MQ=255)
                                               tGATTTTATTGAGTCATTTCATGACGCAGTCcattc                                                 <  1:2512558/36‑1 (MQ=255)
                                               tGATTTTATTGAGTCATTTCATGACGCAGTCCATTCATTCCTGTAAg                                      <  1:772515/47‑1 (MQ=255)
                                                gATTTTATTGAGTCATTTCATGACGCAGTCcattca                                                >  1:2458113/1‑36 (MQ=255)
                                                gATTTTATTGAGTCATTTCATGACGCAGTCCattcattc                                             >  1:2050617/1‑39 (MQ=255)
                                                gATTTTATTGAGTCATTTCATGACGCAGTCCATTCATTcc                                            >  1:2307403/1‑40 (MQ=255)
                                                gATTTTATTGAGTCATTTCATGACGCAGTCCATTCATTCCTGTAAGTTCCCCGACGGAAACATTTa                  >  1:1529587/1‑66 (MQ=255)
                                                gATTTTATTGAGTCATTTCATGACGCAGTCCATTCATTCCTGTAAGTTCCCCGACGGAAACATTTAt                 >  1:2578533/1‑67 (MQ=255)
                                                gATTTTATTGAGTCATTTCATGACGCAGTCCATTCATTCCTGTAAGTTCCCCGACGGAAACATTTATCAAg             >  1:488543/1‑71 (MQ=255)
                                                gATTTTATTGAGTCATTTCATGACGCAGTCCATTCATTCCTGTAAGTTCCCCGACGGAAACATTTATCAAg             >  1:2411942/1‑71 (MQ=255)
                                                 aTTTTATTGAGTCATTTCATGACGCAGTCCATTCATTCCTGTAAGtt                                    >  1:533151/1‑47 (MQ=255)
                                                 aTTTTATTGAGTCATTTCATGACGCAGTCCATTCATTCCTGTAAGTTCCCCGACGGAAACATTTATCAAg             <  1:2776280/70‑1 (MQ=255)
                                                 aTTTTATTGAGTCATTTCATGACGCAGTCCATTCATTCCTGTAAGTTCCCCGACGGAAACATTTATCAAg             <  1:579464/70‑1 (MQ=255)
                                                 aTTTTATTGAGTCATTTCATGACGCAGTCCATTCATTCCTGTAAGTTCCCCGACGGAAACATTTATCAAg             <  1:23435/70‑1 (MQ=255)
                                                   tttATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGTAAGTTCCCCGACGGAAACATTTATCAAGCgg          >  1:698161/1‑71 (MQ=255)
                                                    ttATTGAGTCATTTCATGACGCAGTCCATTCATTCCTGTAAGTTCCCCGACGGAAACATTTATCAAGCggg         >  1:2227951/1‑71 (MQ=255)
                                                     tATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGTa                                        <  1:859716/39‑1 (MQ=255)
                                                           gTCATTTCATGACGCAGTCAATTCATTCCTGTAAGTTCCCCGACGGAAACATTTATCAAGCGGGGTACGcc  <  1:501436/71‑1 (MQ=255)
                                                              |                                                                   
AGCGTGTCCTGCACAAGGCATGAGTAGTGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGTAAGTTCCCCGACGGAAACATTTATCAAGCGGGGTACGCC  >  W3110S.gb/4183973‑4184102

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: