Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,177,230 C→T 8.2% D89N (GAT→AAT)  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,2300CT 8.2% 181.3 / 4.8 61D89N (GAT→AAT) gspKgeneral secretory pathway component, cryptic
Reads supporting (aligned to +/- strand):  ref base C (29/27);  new base T (2/3);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

ATTGAGATTAAAACAGTGCTGGGCATCGCGTAACTGCCACTTCACTGTATTGCCATCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAGTACTGCTGCAACGTGGTCTGTCTGTCGTTATCTTTAACG  >  W3110S.gb/4177174‑4177300
                                                        |                                                                      
aTTGAGATTAAAACAGTGCTGGGCATCGCGTAACTGCCACTTCACTGTATTGCCATCTTCCAGCTGCa                                                             <  1:3993453/68‑1 (MQ=255)
 ttGAGATTAAAACAGTGCTGGGCATCGCGTAACTGCCACTTCACTGTATTGCCATCTTCCAGc                                                                 <  1:1230924/63‑1 (MQ=255)
   gagaTTAAAACAGTGCTGGGCATCGCGTAACTGCCACTTCACTGTATTGCCATCTTCCAGCTGCAACtgtt                                                       >  1:560367/1‑71 (MQ=255)
   gagaTTAAAACAGTGCTGGGCATCGCGTAACTGCCACTTCACTGTATTGCCATCTTCCAGCTGCAACtgt                                                        <  1:1891067/70‑1 (MQ=255)
     gaTTAAAACAGTGCTGGGCATCGCGTAACTGCCACTTCACTGTATTGCCATCTTCCAGCTGCAACtgttgt                                                     <  1:3131860/71‑1 (MQ=255)
     gaTTAAAACAGTGCTGGGCATCGCGTAACTGCCACTTCACTGTATTGCCATCTTCCAGCTGCAACtgttgt                                                     <  1:1340085/71‑1 (MQ=255)
      aTTAAAACAGTGCTGGGCATCGCGTAACTGCCACTTCACTGTATTGCCATCTTCCAGCTGCAACt                                                          <  1:820513/65‑1 (MQ=255)
       ttAAAACAGTGCTGGGCATCGCGTAACTGCCACTTCACTGTATTGCCATCTTCCAGCTGCAACTGTTGTgg                                                   <  1:764987/71‑1 (MQ=255)
       ttAAAACAGTGCTGGGCATCGCGTAACTGCCACTTCACTGTATTGCCATCTTCCAGCTGCAACTGTTGTgg                                                   <  1:513414/71‑1 (MQ=255)
        tAAAACAGTGCTGGGCATCGCGTAACTGCCACTTCACTGTATTGCCATCTTCCAGCTGCAACTGTTGTGGc                                                  >  1:274791/1‑71 (MQ=255)
         aaaaCAGTGCTGGGCATCGCGTAACTGCCACTTCACTGTATTGCCATCTTCCAGCTGCAACTGTTGTGGCt                                                 >  1:2521245/1‑71 (MQ=255)
         aaaaCAGTGCTGGGCATCGCGTAACTGCCACTTCACTGTATTGCCATCTTCCAGCTGCAACTGTTGTGGCt                                                 >  1:362480/1‑71 (MQ=255)
          aaaCAGTGCTGGGCATCGCGTAACTGCCACTTCACTGTATTGCCATCTTCCAGCTGCAACTGTTGTGGCt                                                 <  1:817325/70‑1 (MQ=255)
           aaCAGTGCTGGGCATCGCGTAACTGCCACTTCACTGTATTGCCATCTTACAGCTGCAACTGTTGTGGCtgt                                               <  1:1778612/71‑1 (MQ=255)
            aCAGTGCTGGGCATCGCGTAACTGCCACTTCACTGTATTGCCATCTTCCAGCTGCAACTGTTGTgg                                                   <  1:3260007/66‑1 (MQ=255)
              aGTGCTGGGCATCGCGTAACTGCCACTTCACTGTATTGCCATCTTCCAGCTGCAACTGTTTTGGCTGTGcc                                            <  1:588914/71‑1 (MQ=255)
               gTGCTGGGCATCGCGTAACTGCCACTTCACTGTATTGCCATCTTCCAGCTGCAACTGTTGTGGCTGTGccc                                           <  1:435361/71‑1 (MQ=255)
                  cTGGGCATCGCGTAACTGCCACTTCACTGTATTGCCATCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAg                                         >  1:4001299/1‑70 (MQ=255)
                  cTGGGCATCGCGTAACTGCCACTTCACTGTATTGCCATCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAg                                         >  1:3171283/1‑70 (MQ=255)
                  cTGGGCATCGCGTAACTGCCACTTCACTGTATTGCCATCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAg                                         >  1:944627/1‑70 (MQ=255)
                  cTGGGCATCGCGTAACTGCCACTTCACTGTATTGCCATCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAg                                         >  1:2201822/1‑70 (MQ=255)
                    gggCATCGCGTAACTGCCACTTCACTGTATTGCCATCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAGTAc                                      >  1:650065/1‑71 (MQ=255)
                     ggCATCGCGTAACTGCCACTTCACTGTATTGCCATCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAGTAct                                     <  1:221556/71‑1 (MQ=255)
                     ggCATCGCGTAACTGCCACTTCACTGTATTGCCATCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAGTAct                                     >  1:3199629/1‑71 (MQ=255)
                         tCGCGTAACTGCCACTTCACTGTATTGCCATCTTCCAGCTGc                                                              <  1:2906205/42‑1 (MQ=255)
                           gcgTAACTGCCACTTCACTGTATTGCCATCTTCCAGCTGCAACtgt                                                        >  1:4253867/1‑46 (MQ=255)
                           gcgTAACTGCCACTTCACTGTATTGCCATCTTCCAGCTGCAACtgt                                                        >  1:3576519/1‑46 (MQ=255)
                           gcgTAACTGCCACTTCACTGTATTGCCATCTTCCAGCTGCAACtgt                                                        >  1:1764002/1‑46 (MQ=255)
                           gcgTAACTGCCACTTCACTGTATTGCCATCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAGTACtgctgc                                >  1:2561590/1‑70 (MQ=255)
                           gcgTAACTGCCACTTCACTGTATTGCCATCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAGTACtgctgc                                >  1:4297783/1‑70 (MQ=255)
                              tAACTGCCACTTCACTGTATTGCCATCTTCCAGCTGCAACtgt                                                        <  1:3297955/43‑1 (MQ=255)
                                   gCCACTTCACTGTATTGCCATCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAGTACTGCTGCAACGTGgtc                       >  1:627478/1‑71 (MQ=255)
                                    ccACTTCACTGTATTGCCATTTTCCAGCTGCAACTGTTGCGGTTGTGCCCAGTACTGCTGCAACGTGgttt                      <  1:725319/71‑3 (MQ=255)
                                    ccACTTCACTGTATTGCCATTTTCCAGCTGCAACTGTTGCGGTTGTGCCCAGTACTGCTGCAACGTGgttt                      >  1:2431399/1‑69 (MQ=255)
                                    ccACTTCACTGTATTGCCATTTTCCAGCTGCAACTGTTGCGGTTGTGCCCAGTACTGCTGCAACGTGgttt                      >  1:4013523/1‑69 (MQ=255)
                                    ccACTTCACTGTATTGCCATTTTCCAGCTGCAACTGTTGCGGTTGTGCCCAGTACTGCTGCAACGTGgttt                      <  1:4030693/71‑3 (MQ=255)
                                    ccACTTCACTGTATTGCCATTTTCCAGCTGCAACTGTTGCGGTTGTGCCCAGTACTGCTGCAACGTGgttt                      <  1:882044/71‑3 (MQ=255)
                                       cctCACTGTATTGCCATCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAGTACTGCTGCAACGTGgtctgt                    >  1:4341851/3‑70 (MQ=255)
                                       cTTCACTGTATTGCCATCTTCCATCTGCAACTGTTGTTGCTGTTCCCAg                                         <  1:521503/49‑1 (MQ=255)
                                       cTTCACTGTATTGCCATCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAGTACTGCTGCa                               <  1:1482204/59‑1 (MQ=255)
                                       cTTCACTGTATTGCCATCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAGTACTGCTGCAACGTGgtctgtc                   >  1:451569/1‑71 (MQ=255)
                                        ttCACTGTATTGCCATCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAGTACTGCTGCAACGTGgtctg                     <  1:3548315/68‑1 (MQ=255)
                                         tCACTGTATTGCCATCTTCCAGTTGCAACTGTTGTGGCTGTGCCCAGTACTGCTGCAACGTGGtctgtctg                 >  1:1836292/1‑71 (MQ=255)
                                         tCACTGTATTGACATCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAGTa                                       >  1:1174001/1‑49 (MQ=255)
                                          cACTGTATTGCCATCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAGTa                                       >  1:2609270/1‑48 (MQ=255)
                                          cACTGTATTGCCATCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAGTACTGCTGCAACGTGGTctgtctgt                >  1:2846190/1‑71 (MQ=255)
                                          cACTGTATTGCCATCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAGTACTGCTGCAACGTGGTctgtctgt                >  1:3987804/1‑71 (MQ=255)
                                           aCTGTATTGCCATCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAGTACTGCTGCAACGTg                          >  1:4201991/1‑60 (MQ=255)
                                               tATTGCCATCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAGTACTGCTGCAACGTGGTCTGTCTGTCGtt            <  1:4174818/70‑1 (MQ=255)
                                               tATTGCCATCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAGTACTGCTGCAACGTGGTCTGTCTGTCGtt            <  1:284345/70‑1 (MQ=255)
                                                aTTGCCATCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAGTACTGCTGCAACGTGGTCTGTCTGTCGTTa           >  1:3087324/1‑70 (MQ=255)
                                                 ttGCCATCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAGTACTGCTGCAACGTGGTCTGTCTGTCGTTAt          <  1:3944301/70‑1 (MQ=255)
                                                 ttGCCATCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAGTACTGCTGCAACGTGGTCTGTCTGTCGTTATc         >  1:2471506/1‑71 (MQ=255)
                                                   gCCATCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAGTACTGCTGCAACGTGGTCTGTCTGTCGtt            <  1:2736415/66‑1 (MQ=255)
                                                   gCCATCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAGTACTGCTGCAACGTGGTCTGTCTGTCGTTATCt        <  1:4108917/70‑1 (MQ=255)
                                                     cATCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAGTACTGCTGCAACGTGGTCTGTCTGTCGTTAt          >  1:1347052/1‑66 (MQ=255)
                                                     cATCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAGTACTGCTGCAACGTGGTCTGTCTGTCGTTATCtt       >  1:3921157/1‑69 (MQ=255)
                                                      aTCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAGTACTGCTGCAACGTGGTCTGTCTGTCGTTAt          <  1:3921078/65‑1 (MQ=255)
                                                       tCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAGTACTGCTGCa                               <  1:4236001/43‑1 (MQ=255)
                                                       tCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAGTACTGCTGCAACGTg                          <  1:837752/48‑1 (MQ=255)
                                                       tCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAGTACTGCTGCAACGTGgtctgt                    >  1:775225/1‑54 (MQ=255)
                                                        cTTCCAGCTGCAACTGTTGTGGCTGTGCCCAGTACTGCTGCAACGTGGTCTGTCTGTCGTTATCTTTAACg  >  1:2192045/1‑71 (MQ=255)
                                                        |                                                                      
ATTGAGATTAAAACAGTGCTGGGCATCGCGTAACTGCCACTTCACTGTATTGCCATCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAGTACTGCTGCAACGTGGTCTGTCTGTCGTTATCTTTAACG  >  W3110S.gb/4177174‑4177300

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: