Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,184,035 C→A 25.5% 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,0350CA25.5% 104.7 / 20.8 47G3* (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 (18/17);  new base A (5/7);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 7.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.83e-01

AGCGTGTCCTGCACAAGGCATGAGTAGTGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGTAAGTTCCCCGACGGAAACATTTATCAAGCGGGGTACGCC  >  W3110S.gb/4183973‑4184102
                                                              |                                                                   
aGCGTGTCCTGCACAAGGCATGAGTAGTGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCATTTCAt                                                               <  1:1211612/69‑1 (MQ=255)
aGCGTGTCCTGCACAAGGCATGAGTAGTGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCATTTCAt                                                               <  1:1782978/69‑1 (MQ=255)
      tCCTGCACAAGGCATGAGTAGTGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCCTTTCATGACGCAGt                                                       >  1:2279252/1‑71 (MQ=255)
      tCCTGCACAAGGCATGAGTAGTGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCCTTTCATGACGCAGt                                                       >  1:1944551/1‑71 (MQ=255)
      tCCTGCACAAGGCATGAGTAGTGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCCTTTCATGACGCAGt                                                       >  1:1738147/1‑71 (MQ=255)
           cacaAGGCATGAGTAGTGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCCTTTCATGa                                                             >  1:1156466/1‑60 (MQ=255)
           cacaAGGCATGAGTAGTGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCCTTTCATGa                                                             >  1:300221/1‑60 (MQ=255)
           cacaAGGCATGAGTAGTGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCCTTTCATGa                                                             >  1:1827548/1‑60 (MQ=255)
           cacaAGGCATGAGTAGTGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCCTTTCATGACGCAGTCcatt                                                  <  1:3976947/71‑1 (MQ=255)
           cacaAGGCATGAGTAGTGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCATTTCATGACGCAGTCcatt                                                  <  1:2391140/71‑1 (MQ=255)
              aaGGCATGAGTAGTGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCCTTTCATGACGCAGTCcattcat                                               <  1:2292263/71‑1 (MQ=255)
              aaGGCATGAATAGTGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCCTTTCATGACGCAGTCcattcat                                               <  1:3897413/71‑1 (MQ=255)
                   aTGAGTAGTGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTg                                          >  1:3805107/1‑71 (MQ=255)
                   aTGAGTAGTGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTg                                          >  1:3100809/1‑71 (MQ=255)
                    tGAGTAGTGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGt                                         <  1:1536296/71‑1 (MQ=255)
                    tGAGTAGTGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCATTTCAt                                                               <  1:549811/49‑1 (MQ=255)
                      agtagtGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCATTTCATg                                                              >  1:2536109/1‑48 (MQ=255)
                      agtagtGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCATTTCATGa                                                             >  1:1195419/1‑49 (MQ=255)
                      agtagtGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCATTTCATGa                                                             >  1:3045206/1‑49 (MQ=255)
                      agtagtGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCATTTCATGa                                                             >  1:2629465/1‑49 (MQ=255)
                      agtagtGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCATTTCATGa                                                             >  1:2440149/1‑49 (MQ=255)
                      agtagtGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCATTTCATGACGCAGTCCATTCATTCCTGTa                                        <  1:3039497/70‑1 (MQ=255)
                             ctgcCAGCAAGCAGCAGGTGATTTTATTGAGTCATTTCATGACGCAGTCCATTCATTCCTGTAAGTTcccc                                <  1:2031073/71‑1 (MQ=255)
                                ccAGCAAGCAGCAGGTGATTTTATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGTAAGTTCCCCg                               <  1:3342109/69‑1 (MQ=255)
                                ccAGCAAGCAGCAGGTGATTTTATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGTAAGTTCCCCg                               <  1:2110433/69‑1 (MQ=255)
                                 cAGCAAGCAGCAGGTGATTTTATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGTaa                                       <  1:1600538/60‑1 (MQ=255)
                                  agcaagcaGCAGGTGATTTTATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGTaa                                       <  1:106593/59‑1 (MQ=255)
                                         agcagGTGATTTTATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGTAAGTTccc                                 >  1:3178075/1‑58 (MQ=255)
                                         agcagGTGATTTTATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGTAAGTTCCCCGACGGAAACAtt                    <  1:2232397/71‑1 (MQ=255)
                                           cagGTGATTTTATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGTAAGTTCCCCGACGGAAACATTTa                  <  1:85100/71‑1 (MQ=255)
                                              gTGATTTTATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGTAAGTTCCCCGACGGAAACATTTATCa               <  1:1828918/71‑1 (MQ=255)
                                               tGATTTTATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGTAAGTTCCCCGACGGAAACATTTATc                >  1:1362470/1‑69 (MQ=255)
                                               tGATTTTATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGTAAGTTCCCCGACGGAAACATTTATc                >  1:2186412/1‑69 (MQ=255)
                                                gATTTTATTGAGTCCTTTCATGACGCAGTCcattca                                                >  1:2020238/1‑36 (MQ=255)
                                                gATTTTATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGTAAg                                      >  1:2815439/1‑46 (MQ=255)
                                                gATTTTATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGTAAg                                      >  1:3852148/1‑46 (MQ=255)
                                                gATTTTATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGTAAg                                      >  1:1455929/1‑46 (MQ=255)
                                                 aTTTTATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGt                                         <  1:1313347/42‑1 (MQ=255)
                                                 aTTTTATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGt                                         <  1:384628/42‑1 (MQ=255)
                                                 aTTTTATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGt                                         <  1:1692355/42‑1 (MQ=255)
                                                 aTTTTATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGt                                         <  1:1173363/42‑1 (MQ=255)
                                                 aTTTTATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGTAAGTTCCCCGACGGAAACATTTATCAAGc            >  1:2174016/1‑71 (MQ=255)
                                                 aTTTTATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGTAAGTTCCCCGACGGAAACATTTATCAAGc            >  1:2918607/1‑71 (MQ=255)
                                                    ttATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGTAAGTTCACCGAc                             >  1:2339804/1‑51 (MQ=255)
                                                      aTTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGTAAGTTCCCCGACGGAAACATTTATCAAGCgggg        <  1:3587646/70‑1 (MQ=255)
                                                           gTCCTTTCATGACGCAGTCCATTCATTCCTGTAAGTTCCCCGACGGAAACATTTATCAAGCGGGGTACGcc  <  1:2470577/71‑1 (MQ=255)
                                                           gTCATTTCATGACGCAGTCCATTCATTCCTGTAAGTTCCCCGACGGAAACATTTATCAAGCGGGGTACGcc  <  1:2095031/71‑1 (MQ=255)
                                                              |                                                                   
AGCGTGTCCTGCACAAGGCATGAGTAGTGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGTAAGTTCCCCGACGGAAACATTTATCAAGCGGGGTACGCC  >  W3110S.gb/4183973‑4184102

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: