Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,184,035 C→A 26.9% 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,0350CA26.9% 119.0 / 26.5 52G3* (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 (29/9);  new base A (4/10);  total (33/19)
Fisher's exact test for biased strand distribution p-value = 2.85e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.71e-01

CTCAGCGTGTCCTGCACAAGGCATGAGTAGTGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGTAAGTTCCCCGACGGAAACATTTATCAAGCGGGGTACGCCG  >  W3110S.gb/4183970‑4184103
                                                                 |                                                                    
ctcAGCGTGTCCTGCACAAGGCATGAGTAGTGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCCttt                                                                   <  1:965411/69‑1 (MQ=255)
ctcAGCGTGTCCTGCACAAGGCATGAGTAGTGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCCttt                                                                   <  1:3687648/69‑1 (MQ=255)
     cGTGTCCTGCACAAGGCATGAGTAGTGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCCTTTCAt                                                                <  1:2608385/67‑1 (MQ=255)
     cGTGTCCTGCACAAGGCATGAGTAGTGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCCTTTCAt                                                                <  1:1462070/67‑1 (MQ=255)
              cacaAGGCATGAGTAGTGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCCTTTCATGACGCAGTcc                                                      >  1:453771/1‑68 (MQ=255)
                   ggCATGAGTAGTGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCATTTCATGACGCAGTCcattcat                                                <  1:2830605/69‑1 (MQ=255)
                         agtagtGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCCTTTCATGACg                                                            >  1:850683/1‑51 (MQ=255)
                         agtagtGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCCTTTCATGACg                                                            >  1:2849266/1‑51 (MQ=255)
                         agtagtGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCCTTTCATGACg                                                            >  1:1614154/1‑51 (MQ=255)
                         agtagtGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCCTTTCATGACg                                                            >  1:3254874/1‑51 (MQ=255)
                         agtagtGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCCTTTCATGACGCAGTCCATTCATTCCt                                            >  1:1822434/1‑67 (MQ=255)
                         agtagtGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCCTTTCATGACGCAGTCCATTCATTCCt                                            >  1:3400690/1‑67 (MQ=255)
                         agtagtGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCCTTTCATGACGCAGTCCATTCATTCCt                                            >  1:829166/1‑67 (MQ=255)
                           tagtGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGTa                                         <  1:1297727/68‑1 (MQ=255)
                           tagtGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGTa                                         <  1:847745/68‑1 (MQ=255)
                            agtGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCATTTCATGACGCAGTCCATTCATTCCTGTAAg                                       <  1:3581296/69‑1 (MQ=255)
                            agtGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCATTTCATGACGCAGTCCATTCATTCCTGTAAg                                       <  1:2214830/69‑1 (MQ=255)
                            agtGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCATTTCATGACGCAGTCCATTCATTCCTGTAAg                                       <  1:1290904/69‑1 (MQ=255)
                            agtGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCATTTCATGAAGCAGTCCATTCATTCCTGTAAg                                       <  1:919480/69‑1 (MQ=255)
                             gtGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGTAAGt                                      >  1:3386099/1‑69 (MQ=255)
                             gtGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGTAAGt                                      >  1:1159293/1‑69 (MQ=255)
                                ctgcCAGCAAGCAGCAGGTGATTTTATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGTAAGTTcc                                   >  1:1969366/1‑69 (MQ=255)
                                ctgcCAGCAAGCAGCAGGTGATTTTATTGAGTCATTTCATGACGCAGTCCATTCATTCCTGTAAGTTcc                                   >  1:2829416/1‑69 (MQ=255)
                                  gcCAGCAAGCAGCAGGTGATTTTATTGAGTCATTTCATGACGCAGTCCATTCATTCCTGTAAGTTcccc                                 <  1:3641956/69‑1 (MQ=255)
                                   ccAGCAAGCAGCAGGTGATTTTATTGAGTCATTTCATGACGCAGTCCATTCATTCCTGTAAGTTc                                    >  1:401013/1‑65 (MQ=255)
                                   ccAGCAAGCAGCAGGTGATTTTATTGAGTCATTTCATGACGCAGTCCATTCATTCCTGTAAGTTc                                    >  1:2173883/1‑65 (MQ=255)
                                   ccAGCAAGCAGCAGGTGATTTTATTGAGTCATTTCATGACGCAGTCCATTCATTACTGTAAGTTc                                    >  1:2405577/1‑65 (MQ=255)
                                     agcaagcaGCAGGTGATTTTATTGAGTCATTTCATGACGCAGTCcat                                                    <  1:1391815/47‑1 (MQ=255)
                                     agcaagcaGCAGGTGATTTTATTGAGTCATTTCATGACGCAGTCcat                                                    <  1:2488586/47‑1 (MQ=255)
                                     agcaagcaGCAGGTGATTTTATTGAGTCATTTCATGACGCAGTCcat                                                    <  1:1733136/47‑1 (MQ=255)
                                        aagcaGCAGGTGATTTTATTGAGTCCTTTCATGACGCAGTCcatt                                                   >  1:3234853/1‑45 (MQ=255)
                                        aagcaGCAGGTGATTTTATTGAGTCCTTTCATGACGCAGTCcatt                                                   >  1:47489/1‑45 (MQ=255)
                                        aagcaGCAGGTGATTTTATTGAGTCCTTTCATGACGCAGTCcat                                                    >  1:2929182/1‑44 (MQ=255)
                                        aagcaGCAGGTGATTTTATTGAGACCTTTCATGACGCAGTCcatt                                                   >  1:699369/1‑45 (MQ=255)
                                         agcagcAGGTGATTTTATTGAGTCATTTCATGACGCAGTCcatt                                                   <  1:162891/44‑1 (MQ=255)
                                                ggTGATTTTATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGTAAGTTCCCCGACGGAAACAttt                    >  1:2660146/1‑68 (MQ=255)
                                                  tGATTTTATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGTAAGTTCCCCGACGGAAACATTTa                   >  1:2220109/1‑67 (MQ=255)
                                                   gATTTTATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGTAAGTTCCCCGACGGAAACATTTATCa                <  1:2076236/69‑1 (MQ=255)
                                                    aTTTTATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGTAAGTTCCCCGACGGAAACATTTATCaa               >  1:3225805/1‑69 (MQ=255)
                                                       ttATTGAGTCCTTTCATGACGCAGTCCattcattc                                              >  1:3567849/1‑35 (MQ=255)
                                                       ttATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGTAAGTTCCCCGACGGAATCAtt                     >  1:1235275/1‑60 (MQ=255)
                                                       ttATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGTAAGTTCCCCGACGGAAACAtt                     >  1:3550336/1‑60 (MQ=255)
                                                       ttATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGTAAGTTCCCCGACGGAAACAtt                     >  1:1704599/1‑60 (MQ=255)
                                                       ttATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGTAAGTTCCCCGACGGAAACAtt                     >  1:3248867/1‑60 (MQ=255)
                                                       ttATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGTAAGTTCCCCGACGGAAACAtt                     <  1:30319/60‑1 (MQ=255)
                                                       ttATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGTAAGTTCCCCGACGGAAACAtt                     >  1:2913415/1‑60 (MQ=255)
                                                       ttATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGTAAGTTCCCCGACGGAAACAtt                     <  1:1547812/60‑1 (MQ=255)
                                                       ttATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGTAAGTTCCCCGACGGAAACATTTATCAAGCg            >  1:3549220/1‑69 (MQ=255)
                                                       ttATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGTAAGTTCCCCGACGGAAACATTTATCAAGCg            >  1:2287039/1‑69 (MQ=255)
                                                        tATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGTAAGTTCCCCGACGGAAACATTTATCAAGCgg           >  1:3709874/1‑69 (MQ=255)
                                                        tATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGTAAGTTCCCCGACGGAAACATTTATCAAGCgg           >  1:894710/1‑69 (MQ=255)
                                                        tATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGTAAGTTCCCCGACGGAAACATTTATCAAGCgg           >  1:1146247/1‑69 (MQ=255)
                                                                 cTTTCATGACGCAGTCCATTCATTCCTGTAAGTTCCCCGACGGAAACATTTATCAAGCGGGGTACGCCg  >  1:935288/1‑69 (MQ=255)
                                                                 |                                                                    
CTCAGCGTGTCCTGCACAAGGCATGAGTAGTGCTGCCAGCAAGCAGCAGGTGATTTTATTGAGTCCTTTCATGACGCAGTCCATTCATTCCTGTAAGTTCCCCGACGGAAACATTTATCAAGCGGGGTACGCCG  >  W3110S.gb/4183970‑4184103

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: