Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,177,469 G→A 57.3% 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,4690GA57.3% 29.4 / 97.5 75A9V (GCG→GTG) gspKgeneral secretory pathway component, cryptic
Reads supporting (aligned to +/- strand):  ref base G (22/10);  new base A (24/19);  total (46/29)
Fisher's exact test for biased strand distribution p-value = 3.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.00e-01

AACTGAGCGTCATGTCTGCTGCCAGAGCAGCCATGAGCGCCAGTAACATGAGTACAATGAGCAGCGCGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGCGTCTGAGCTGGCGACTCTTC  >  W3110S.gb/4177404‑4177534
                                                                 |                                                                 
atcTGAGCGTCATGTCTGCTGCCAGAGCAGCCATGAGCGCCAGTAACATGAGTACAATGAGCAGcacg                                                                 <  1:2724898/66‑1 (MQ=255)
          cATGTCTGCTGCCAGAGCAGCCATGAGCGCCAGTAACATGAGTACAATGAGCAGCACGACGCCGCgttg                                                      >  1:733744/1‑69 (MQ=255)
          cATGTCTGCTGCCAGAGCAGCCATGAGCGCCAGTAACATGAGTACAATGAGCAGCACGACGCCGCgtt                                                       >  1:610172/1‑68 (MQ=255)
           aTGTCTGCTGCCAGAGCAGCCATGAGCGCCAGTAACATGAGTACAATGAGCAGCACGACGCCGCgttg                                                      >  1:1639303/1‑68 (MQ=255)
            tGTCTGCTGCCAGAGCAGCCATGAGCGCCAGTAACATGAGTACAATGAGCAGCGCGACGCCGCgttgtt                                                    >  1:2312486/1‑69 (MQ=255)
             gTCTGCTGCCAGAGCAGCCATGAGCGCCAGTAACATGAGTACAATGAGCAGCGCGACGCCGCgttgtt                                                    <  1:3716712/68‑1 (MQ=255)
             gTCTGCTGCCAGAGCAGCCATGAGCGCCAGTAACATGAGTACAATGAGCAGCGCGACGCCGCgttgtt                                                    <  1:1222016/68‑1 (MQ=255)
             gTCTGCTGCCAGAGCAGCCATGAGCGCCAGTAACATGAGTACAATGAGCAGCGCGACGCCGCgttgtt                                                    <  1:3584922/68‑1 (MQ=255)
             gTCTGCTGCCAGAGCAGCCATGAGCGCCAGTAACATGAGTACAATGAGCAGCGCGACGCCGCgttgtt                                                    <  1:3580547/68‑1 (MQ=255)
             gTCTGCTGCCAGAGCAGCCATGAGCGCCAGTAACATGAGTACAATGAGCAGCGCGACGCCGCgttgtt                                                    <  1:1581322/68‑1 (MQ=255)
                tgctgcCAGAGCAGCCATGAGCGCCAGTAACATGAGTACAATGAGCAGCGCGACGCCGCGTTGTTCat                                                 <  1:750100/68‑1 (MQ=255)
                    gcCAGAGCAGCCATGAGCGCCAGTAACATGAGTAGAATGAGCAGCGCGACGCCGCGTTGTTCattat                                              <  1:67938/67‑1 (MQ=255)
                    gcCAGAGCAGCCATGAGCGCCAGTAACATGAGTACAATGAGCAGCGCGACGCCGCGTTGTTCattat                                              <  1:1032423/67‑1 (MQ=255)
                     ccAGAGCAGCCATGAGCGCCAGTAACATGAGTACAATGAGCAGCACGACGCCGCGTTGTTCATTATTCa                                           >  1:860956/1‑69 (MQ=255)
                     ccAGAGCAGCCATGAGCGCCAGTAACATGAGTACAATGAGCAGCACGACGCCGCGTTGTTCATTATTCa                                           >  1:3716011/1‑69 (MQ=255)
                        gagCAGACATGAGGGCCAGTAACATGAGTACAATGAGCAGCACGACGCCGCgttgtt                                                    >  1:3632002/1‑57 (MQ=255)
                           cagcCATGAGCGCCAGTAACATGAGTACAATGAGCAGCGCGACGCCGCGTTGTTCATTATTCATGGtt                                      >  1:3477494/1‑68 (MQ=255)
                            agcCATGAGCGCCAGTAACATGAGTACAATGAGCAGCACGACGCCGCGTTGTTCATTATTCATGGttat                                    <  1:1877073/69‑1 (MQ=255)
                            agcCATGAGCGCCAGTAACATGAGTACAATGAGCAGCACGACGCCGCGTTGTTCATTATTCATGGttat                                    <  1:3260861/69‑1 (MQ=255)
                             gcCATGAGCGCCAGTAACATGAGTACAATGAGCAGCACGACGCCGCGTTGTTCATTATTCATGGttatt                                   >  1:151539/1‑69 (MQ=255)
                             gcCATGAGCGCCAGTAACATGAGTACAATGAGCAGCACGACGCCGCGTTGTTCATTATTCATGGttat                                    >  1:1809946/1‑68 (MQ=255)
                               cATGAGCGCCAGTAACATGAGTACAATGAGCAGCACGACGCCGCGTTGTTCATTATTCATGGTtattat                                 >  1:2361453/1‑69 (MQ=255)
                                aTGAGCGCCAGTAACATGAGTACAATGAGCAGCACGACGCCGCGTTGTTCATTATTCATGGTTattatt                                >  1:3098139/1‑69 (MQ=255)
                                 tGAGCGCCAGTAACATGAGTACAATGAGCAGCACGACGCCGCGTTGTTCat                                                 <  1:686665/51‑1 (MQ=255)
                                 tGAGCGCCAGTAACATGAGTACAATGAGCAGCACGACGCCGCGTTGTTCat                                                 <  1:1845713/51‑1 (MQ=255)
                                 tGAGCGCCAGTAACATGAGTACAATGAGCAGCACGACGCCGCGTTGTTCat                                                 <  1:1901520/51‑1 (MQ=255)
                                 tGAGCGCCAGTAACATGAGTACAATGAGCAGCACGACGCCGCGTTGTTCat                                                 <  1:2107952/51‑1 (MQ=255)
                                 tGAGCGCCAGTAACATGAGTACAATGAGCAGCACGACGCCGCGTTGTTCat                                                 <  1:339885/51‑1 (MQ=255)
                                 tGAGCGCCAGTAACATGAGTACAATGAGCAGCACGACGCCGCGTTGTTCATTATTCATGGTTATTATTc                               >  1:3599972/1‑69 (MQ=255)
                                           tAACATGAGTACAATGAGCAGCACGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGCg                     <  1:568976/69‑1 (MQ=255)
                                            aaCATGAGTACAATGAGCAGCGCGACGCCGCGTTGTTCat                                                 >  1:1639215/1‑40 (MQ=255)
                                            aaCATGAGTACAATGAGCAGCGCGACGCCGCGTTGTTCat                                                 >  1:3244300/1‑40 (MQ=255)
                                            aaCATGAGTACAATGAGCAGCGCGACGCCGCGTTGTTCat                                                 >  1:380008/1‑40 (MQ=255)
                                            aaCATGAGTACAATGAGCAGCGCGACGCCGCGTTGTTCat                                                 >  1:3213045/1‑40 (MQ=255)
                                            aaCATGAGTACAATGAGCAGCGCGACGCCGCGTTGTTCat                                                 >  1:2720682/1‑40 (MQ=255)
                                            aaCATGAGTACAATGAGCAGCACGACGCCGCGTTGTTCat                                                 >  1:2719407/1‑40 (MQ=255)
                                            aaCATGAGTACAATGAGCAGCACGACGCCGCGTTGTTCat                                                 >  1:1031831/1‑40 (MQ=255)
                                            aaCATGAGTACAATGAGCAGCACGACGCCGCGTTGTTCat                                                 >  1:2303377/1‑40 (MQ=255)
                                            aaCATGAGTACAATGAGCAGCACGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGcc                         <  1:541739/64‑1 (MQ=255)
                                             aCATGAGTACAATGAGCAGCGCGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGCGt                    >  1:2030080/1‑68 (MQ=255)
                                             aCATGAGTACAATGAGCAGCGCGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGCGt                    >  1:2935235/1‑68 (MQ=255)
                                             aCATGAGTACAATGAGCAGCGCGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGCGt                    >  1:141014/1‑68 (MQ=255)
                                             aCATGAGTACAATGAGCAGCACGACGCCGCGTTGTTCATTATTCATGGt                                       >  1:3591235/1‑49 (MQ=255)
                                             aCATGAGTACAATGAGCAGCACGACGCAGCGTTGTTCATTATTCATGGt                                       >  1:1654298/1‑49 (MQ=255)
                                               aTGAGTACAATGAGCAGCGCGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGCGTCt                  >  1:1030642/1‑68 (MQ=255)
                                               aTGAGTACAATGAGCAGCGCGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGCGTCTg                 >  1:543254/1‑69 (MQ=255)
                                               aTGAGTACAATGAGCAGCGCGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGCGTCTg                 <  1:3069696/69‑1 (MQ=255)
                                                tGAGTACAATGAGCAGCGCGACGCCGCGTTGTTCATTATTCAt                                          >  1:1243105/1‑43 (MQ=255)
                                                tGAGTACAATGAGCAGCACGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGCGTCTGa                <  1:2319799/69‑1 (MQ=255)
                                                tGAGTACAATGAGCAGCACGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGCGTCTGa                <  1:3405778/69‑1 (MQ=255)
                                                  aGTACAATGAGCAGCACGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGcc                         <  1:2462816/58‑1 (MQ=255)
                                                  aGTACAATGAGCAGCACGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGcc                         <  1:1472018/58‑1 (MQ=255)
                                                      cAATGAGCAGCGCGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGc                      >  1:1044974/1‑57 (MQ=255)
                                                      cAATGAGCAGCGCGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGc                      >  1:1948928/1‑57 (MQ=255)
                                                      cAATGAGCAGCGCGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGc                      >  1:2883506/1‑57 (MQ=255)
                                                      cAATGAGCAGCGCGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGc                      >  1:1170389/1‑57 (MQ=255)
                                                      cAATGAGCAGCACGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGCGTCTGAGCTGGc          >  1:454244/1‑69 (MQ=255)
                                                      cAATGAGCAGCACGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGCGTCTGAGCTGGc          >  1:2258835/1‑69 (MQ=255)
                                                      cAATGAGCAGCACGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGCGTCTGAGCTGGc          >  1:449902/1‑69 (MQ=255)
                                                        aTGAGCAGCACGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGCGTCTGAGCTGGCg         <  1:3334681/68‑1 (MQ=255)
                                                         tGAGCAGCACGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGCGTCTGAGCTGGCGAc       >  1:1153469/1‑69 (MQ=255)
                                                         tGAGCAGCACGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGCGTCTGAGCTGGCGAc       >  1:334169/1‑69 (MQ=255)
                                                           agcagcGCGAGGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGCGTCTGAGCTGGCGActc     >  1:3306821/1‑69 (MQ=255)
                                                           agcagcGCGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGCGTCTGAGCTGGCGActc     >  1:1610499/1‑69 (MQ=255)
                                                           agcagcGCGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGCGTCTGAGCTGGCGActc     >  1:3466252/1‑69 (MQ=255)
                                                           agcagcGCGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGCGTCTGAGCTGGCGAct      >  1:2735344/1‑68 (MQ=255)
                                                           agcagcACGACGCCGCGTTGTTCATTATTCATGGTTattatt                                >  1:2813278/1‑42 (MQ=255)
                                                           agcagcACGACGCCGCGTTGTTCATTATTCATGGTTattatt                                >  1:2617691/1‑42 (MQ=255)
                                                           agcagcACGACGCCGCGTTGTTCATTATTCATGGTTattatt                                >  1:2613292/1‑42 (MQ=255)
                                                           agcagcACGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGCGTCTGAGCTGGCGActc     <  1:202096/69‑1 (MQ=255)
                                                           agcagcACGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGCGTCTGAGCTGGCGActc     <  1:3678460/69‑1 (MQ=255)
                                                           agcagcACGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGCGTCTGAGCTGGCGActc     <  1:1415813/69‑1 (MQ=255)
                                                           agcagcACGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGCGTCTGAGCTGGCGActc     <  1:3495667/69‑1 (MQ=255)
                                                             cagcGCGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGCGTCTGAGCTGGCGActct    <  1:1090808/68‑1 (MQ=255)
                                                             cagcACGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGCGGCTGAGCTGGCGACtctt   <  1:1461804/69‑1 (MQ=255)
                                                              agcGCGACGCCGCGTTGTTCATTATTCATGGTTATTACTCCCGGCCTGCGTCTGAGCTGGCGACtcttc  >  1:2652708/1‑69 (MQ=255)
                                                                 |                                                                 
AACTGAGCGTCATGTCTGCTGCCAGAGCAGCCATGAGCGCCAGTAACATGAGTACAATGAGCAGCGCGACGCCGCGTTGTTCATTATTCATGGTTATTATTCCCGGCCTGCGTCTGAGCTGGCGACTCTTC  >  W3110S.gb/4177404‑4177534

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: