Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,177,249 T→C 10.9% P82P (CCA→CCG 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,2490TC10.9% 160.8 / 11.1 64P82P (CCA→CCGgspKgeneral secretory pathway component, cryptic
Reads supporting (aligned to +/- strand):  ref base T (32/25);  new base C (4/3);  total (36/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.29e-01

GATTAAAACAGTGCTGGGCATCGCGTAACTGCCACTTCACTGTATTGCCATCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAGTACTGCTGCAACGTGGTCTGTCTGTCGTTATCTTTAACGTCCTGGGTGAGGCTGGCAA  >  W3110S.gb/4177179‑4177319
                                                                      |                                                                      
gaTTAAAACAGTGCTGGGCATCGCGTAACTGCCACTTCACTGTATTGCCATCTTCCAGCTGCAACtgttgt                                                                        <  1:3131860/71‑1 (MQ=255)
gaTTAAAACAGTGCTGGGCATCGCGTAACTGCCACTTCACTGTATTGCCATCTTCCAGCTGCAACtgttgt                                                                        <  1:1340085/71‑1 (MQ=255)
  ttAAAACAGTGCTGGGCATCGCGTAACTGCCACTTCACTGTATTGCCATCTTCCAGCTGCAACTGTTGTgg                                                                      <  1:513414/71‑1 (MQ=255)
  ttAAAACAGTGCTGGGCATCGCGTAACTGCCACTTCACTGTATTGCCATCTTCCAGCTGCAACTGTTGTgg                                                                      <  1:764987/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:3171283/1‑70 (MQ=255)
             cTGGGCATCGCGTAACTGCCACTTCACTGTATTGCCATCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAg                                                            >  1:2201822/1‑70 (MQ=255)
             cTGGGCATCGCGTAACTGCCACTTCACTGTATTGCCATCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAg                                                            >  1:4001299/1‑70 (MQ=255)
             cTGGGCATCGCGTAACTGCCACTTCACTGTATTGCCATCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAg                                                            >  1:944627/1‑70 (MQ=255)
               gggCATCGCGTAACTGCCACTTCACTGTATTGCCATCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAGTAc                                                         >  1:650065/1‑71 (MQ=255)
                ggCATCGCGTAACTGCCACTTCACTGTATTGCCATCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAGTAct                                                        >  1:3199629/1‑71 (MQ=255)
                ggCATCGCGTAACTGCCACTTCACTGTATTGCCATCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAGTAct                                                        <  1:221556/71‑1 (MQ=255)
                      gcgTAACTGCCACTTCACTGTATTGCCATCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAGTACtgctgc                                                   >  1:2561590/1‑70 (MQ=255)
                      gcgTAACTGCCACTTCACTGTATTGCCATCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAGTACtgctgc                                                   >  1:4297783/1‑70 (MQ=255)
                              gCCACTTCACTGTATTGCCATCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAGTACTGCTGCAACGTGgtc                                          >  1:627478/1‑71 (MQ=255)
                               ccACTTCACTGTATTGCCATTTTCCAGCTGCAACTGTTGCGGTTGTGCCCAGTACTGCTGCAACGTGgttt                                         <  1:4030693/71‑3 (MQ=255)
                               ccACTTCACTGTATTGCCATTTTCCAGCTGCAACTGTTGCGGTTGTGCCCAGTACTGCTGCAACGTGgttt                                         >  1:2431399/1‑69 (MQ=255)
                               ccACTTCACTGTATTGCCATTTTCCAGCTGCAACTGTTGCGGTTGTGCCCAGTACTGCTGCAACGTGgttt                                         >  1:4013523/1‑69 (MQ=255)
                               ccACTTCACTGTATTGCCATTTTCCAGCTGCAACTGTTGCGGTTGTGCCCAGTACTGCTGCAACGTGgttt                                         <  1:725319/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:3987804/1‑71 (MQ=255)
                                     cACTGTATTGCCATCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAGTACTGCTGCAACGTGGTctgtctgt                                   >  1:2846190/1‑71 (MQ=255)
                                      aCTGTATTGCCATCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAGTACTGCTGCAACGTg                                             >  1:4201991/1‑60 (MQ=255)
                                          tATTGCCATCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAGTACTGCTGCAACGTGGTCTGTCTGTCGtt                               <  1:284345/70‑1 (MQ=255)
                                          tATTGCCATCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAGTACTGCTGCAACGTGGTCTGTCTGTCGtt                               <  1:4174818/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)
                                                    ttCCAGCTGCAACTGTTGTGGCTGTGCCCAGTACTGCTGCAACGTGGTCTGTCTGTCGTTATCttt                         <  1:1428335/66‑1 (MQ=255)
                                                    ttCCAGCTGCAACTGTTGTGGCTGTGCCCAGTACTGCTGCAACGTGGTCTGTCTGTCGTTATCTTTAACGt                    >  1:3722340/1‑71 (MQ=255)
                                                    ttCAAGCTGCAACTGTTGTGGCTGTGCCCAGTACTGCTGCAACGTGGTCtgtctgtc                                  >  1:436747/1‑57 (MQ=255)
                                                      ccAGCTGCAACTGTTGTGGCTGTGCCCAGTACTGCTGCAACGTGGTCTGTCTGTCGTTATCTTTAAc                      <  1:3291414/67‑1 (MQ=255)
                                                      ccAGCTGCAACTGTTGTGGCTGTGCCCAGTACTGCTGCAACGTGGTCTGTCTGTCGTTATCTTTAACGt                    <  1:3335257/69‑1 (MQ=255)
                                                      ccAGCTGCAACTGTTGTGGCTGTGCCCAGTACTGCTGCAACGTGGTCTGTCTGTCGTTATCTTTAACGTcc                  >  1:1032016/1‑71 (MQ=255)
                                                      ccAGCTGCAACTGTTGTGGCTGTGCCCAGTACTGCTGCAACGTGGTCTGTCTGTCGTTATCTTTAACGTcc                  >  1:1536568/1‑71 (MQ=255)
                                                       cAGCTGCAACTGTTGTGGCTGTGCCCAGTACTGCTGCAACGTGGTCTGTCTGTCGTTATCTTTAACGTCCt                 <  1:440742/71‑1 (MQ=255)
                                                         gCTGCAACTGTTGTGGCTGTGCCCAGTACTGCTGCAACGTGGTCTGTCTGTCGTTATCTTTAACGTCCt                 >  1:2444103/1‑69 (MQ=255)
                                                          cTGCAACTGTTGCGGTTGTGCCCAGTACTGCTGCAACGTGGTTTGTCTGTCGTTATCTTTAACGTCCTggg              >  1:1050781/1‑71 (MQ=255)
                                                          cTGCAACTGTTGCGGTTGTGCCCAGTACTGCTGCAACGTGGTTTGTCTGTCGTTATCTTTAACGTCCTggg              >  1:1878187/1‑71 (MQ=255)
                                                             cAACTGTTGTGGCTGTGCCCAGTACTGCTGCAACGTGGTCTGTCTGTCGTTATCTTTAACGTCCTGGGt             <  1:1109840/69‑1 (MQ=255)
                                                             cAACTGTTGTGGCTGTGCCCAGTACTGCTGCAACGTGGTCTGTCTGTCGTTATCTTTAACGTCCTGGGTGa           <  1:2269015/71‑1 (MQ=255)
                                                                cTGTTGTGGCTGTGCCCAGTACTGCTGCAACGTGGTCTGTCTGTCGTTATCTTTaa                       >  1:1709015/1‑56 (MQ=255)
                                                                 tgttgtGGCTGTGCCCAGTACTGCTGCAACGTGGTCTGTCTGTCGTTATCttt                         <  1:3404429/53‑1 (MQ=255)
                                                                   ttgtGGCTGTGCCCAGTACTGCTGCAACGTGGTCTGTCTGTCGTTATCTTTAACGTCCTGGGTGAGGCTgg     <  1:3648563/71‑1 (MQ=255)
                                                                   ttgtGGCTGTGCCCAGTACTGCTGCAACGTGGTCTGTCTGTCGTTATCTTTAAAGTCCTGGGTGAGGCTgg     >  1:3518953/1‑71 (MQ=255)
                                                                     gtgGCTGTGCCCAGTACTGCTGCAACGTGGTCTGTCTGTCGTTATCTTTAACGTCCTGGGTGAGGCTGGc    <  1:4037224/70‑1 (MQ=255)
                                                                     gtgGCTGTGCCCAGTACTGCTGCAACGTGGTCTGTCTGTCGTTATCTTTAACGTCCTGGGTGAGGCTGGc    <  1:3850624/70‑1 (MQ=255)
                                                                     gtgGCTGTGCCCAGTACTGCTGCAACGTGGTCTGTCTGTCGTTATCTTTAACGTCCTGGGTGAGGCTGGc    <  1:3275284/70‑1 (MQ=255)
                                                                     gtgGCTGTGCCCAGTACTGCTGCAACGTGGTCTGTCTGTCGTTATCTTTAACGTCCTGGGTGAGGCTGGc    >  1:1981047/1‑70 (MQ=255)
                                                                      tgGCTGTGCCCAGTACTGCTGCAACGTGGTCTGTCTGTCGTTATCttt                         <  1:3480868/48‑1 (MQ=255)
                                                                      tgGCTGTGCCCAGTACTGCTGCAACGTGGTCTGTCTGTCGTTATCTTTAACGTCCTGGGTGAGGCTGGCaa  <  1:4256991/71‑1 (MQ=255)
                                                                      |                                                                      
GATTAAAACAGTGCTGGGCATCGCGTAACTGCCACTTCACTGTATTGCCATCTTCCAGCTGCAACTGTTGTGGCTGTGCCCAGTACTGCTGCAACGTGGTCTGTCTGTCGTTATCTTTAACGTCCTGGGTGAGGCTGGCAA  >  W3110S.gb/4177179‑4177319

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: