Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,177,469 G→A 100% A9V (GCG→GTG)  gspK ← general secretory pathway component, cryptic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,177,4690GA85.0% 95.3 / 10.9 40A9V (GCG→GTG) gspKgeneral secretory pathway component, cryptic
Reads supporting (aligned to +/- strand):  ref base G (1/5);  new base A (22/12);  total (23/17)
Fisher's exact test for biased strand distribution p-value = 6.66e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.28e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ATGTCTGCTGCCAGAGCAGCCATGAGCGCCAGTAACATGAGTACAATGAGCAGCGCGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGCGTCTGAGCTGGCGACTCT  >  W3110S.gb/4177415‑4177532
                                                      |                                                               
aTGTCTGCTGCCAGAGCAGCCATGAGCGCCAGTAACATGAGTACAATGAGCAGCACGACGCCGCgttgtt                                                  <  1:2644954/70‑1 (MQ=255)
   tCTGCTGCCAGAGCAGCCATGAGCGCCAGTAACATGAGTACAATGAGCAGCACGACGCCGCGTTGTTCat                                               <  1:2697740/70‑1 (MQ=255)
   tCTGCTGCCAGAGCAGCCATGAGCGCCAGTAACATGAGTACAATGAGCAGCACGACGCCGCGTTGTTCat                                               <  1:1825800/70‑1 (MQ=255)
       ctgcCAGAGCAGCCATGAGCGCCAGTAACATGAGTACAATGAGCAGCACGACGCCGCGTTGTTCattatt                                           <  1:1171641/70‑1 (MQ=255)
       ctgcCAGAGCAGCCATGAGCGCCAGTAACATGAGTACAATGAGCAGCACGACGCCGCGTTGTTCattatt                                           <  1:3298341/70‑1 (MQ=255)
                  gcCATGAGCGCCAGTAACATGAGTACAATGAGCAGCACGACGCCGCGTTGTTCATTATTCATGGTtattat                               <  1:845982/71‑1 (MQ=255)
                   ccATGAGCGCCAGTAACATGAGTACAATGAGCAGCACGACGCCGCGTTGTTCattatt                                           <  1:1001653/58‑1 (MQ=255)
                   ccATGAGCGCCAGTAACATGAGTACAATGAGCAGCACGACGCCGCGTTGTTCattatt                                           <  1:262659/58‑1 (MQ=255)
                   ccATGAGCGCCAGTAACATGAGTACAATGAGCAGCACGACGCCGCGTTGTTCattatt                                           <  1:1329343/58‑1 (MQ=255)
                        aGCGCCAGTAACATGAGTACAATGAGCAGCACGACGCCGCGTTGTTCATTATTCATGGttat                                  >  1:2780493/1‑62 (MQ=255)
                        aGCGCCAGTAACATGAGTACAATGAGCAGCACGACGCCGCGTTGTTCATTATTCATGGttat                                  >  1:2708545/1‑62 (MQ=255)
                        aGCGCCAGTAACATGAGTACAATGAGCAGCACGACGCCGCGTTGTTCATTATTCATGGttat                                  >  1:3647239/1‑62 (MQ=255)
                        aGCGCCAGTAACATGAGTACAATGAGCAGCACGACGCCGCGTTGTTCATTATTCATGGttat                                  >  1:1249790/1‑62 (MQ=255)
                        aGCGCCAGTAACATGAGTACAATGAGCAGCACGACGCCGCGTTGTTCATTATTCATGGttat                                  >  1:2364125/1‑62 (MQ=255)
                               gTAACATGAGTACAATGAGCAGCACGACGCCGCgttgtt                                                  <  1:2272637/39‑1 (MQ=255)
                                 aaCATGAGTACAATGAGCAGCACGACGCCGCgttgtt                                                  >  1:3599384/1‑37 (MQ=255)
                                 aaCATGAGTACAATGAGCAGCACGACGCCGCgttgtt                                                  >  1:879277/1‑37 (MQ=255)
                                 aaCATGAGTACAATGAGCAGCACGACGCCGCgttgtt                                                  >  1:1068183/1‑37 (MQ=255)
                                 aaCATGAGTACAATGAGCAGCACGACGCCGCgttgtt                                                  >  1:1743801/1‑37 (MQ=255)
                                 aaCATGAGTACAATGAGCAGCACGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGCGTCt                >  1:2480123/1‑71 (MQ=255)
                                 aaCATGAGTACAATGAGCAGCACGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGCGTCt                >  1:214628/1‑71 (MQ=255)
                                 aaCATGAGTACAATGAGCAGCACGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGCGTCt                >  1:1965189/1‑71 (MQ=255)
                                 aaCATGAGTACAATGAGCAGCACGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGCGTCt                >  1:1162461/1‑71 (MQ=255)
                                  ctgtGAGTACAATGAGCAGCGCGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGCGTCTg               <  1:2315897/68‑1 (MQ=255)
                                  ctgtGAGTACAATGAGCAGCGCGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGCGTCTg               <  1:1298804/68‑1 (MQ=255)
                                  ctgtGAGTACAATGAGCAGCGCGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGCGTCTg               <  1:2974220/68‑1 (MQ=255)
                                  aCATGAGTACAATGAGCAGCACGACGCCGCGTTGTTCa                                                >  1:1357602/1‑38 (MQ=255)
                                  aCATGAGTACAATGAGCAGCACGACGCCGCGTTGTTCa                                                >  1:1270198/1‑38 (MQ=255)
                                  aCATGAGTACAATGAGCAGCACGACGCCGCGTTGTTCa                                                >  1:3318116/1‑38 (MQ=255)
                                  aCATGAGTACAATGAGCAGCACGACGCCGCGTTGTTCa                                                >  1:3413629/1‑38 (MQ=255)
                                   tgtGAGTACAATGAGCAGCGCGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGCGTCTg               <  1:3090465/68‑1 (MQ=255)
                                      gAGTACAATGAGCAGCGCGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGCGTCTGAGCt           >  1:1349522/1‑71 (MQ=255)
                                       aGTACAATGAGCAGCACGACGCCGCGTTGTTCatta                                             >  1:1289343/1‑36 (MQ=255)
                                         tACAATGAGCAGCACGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGCGTCTGAGCTGGc        <  1:2830/71‑1 (MQ=255)
                                         tACAATGAGCAGCACGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGCGGCTGAGCTGGc        <  1:192607/71‑1 (MQ=255)
                                           cAATGAGCAGCGCGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGCGTCTGAGCTGGCGa      <  1:327281/71‑1 (MQ=255)
                                           cAATGAGCAGCACGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGCGTCTGAGCTGGCGa      >  1:3625495/1‑71 (MQ=255)
                                              tGAGCAGCACGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGCGTCTGAGCTGGCGActc   >  1:2367404/1‑71 (MQ=255)
                                               gAGCAGCACGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGCGTCTGAGCTGGCGActct  >  1:3022665/1‑71 (MQ=255)
                                               gAGCAGCACGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGCGTCTGAGCTGGCGActct  >  1:2703765/1‑71 (MQ=255)
                                                      |                                                               
ATGTCTGCTGCCAGAGCAGCCATGAGCGCCAGTAACATGAGTACAATGAGCAGCGCGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGCGTCTGAGCTGGCGACTCT  >  W3110S.gb/4177415‑4177532

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: