Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 47,828 A→G 36.8% V20V (GTA→GTG kefC → potassium:proton antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb47,8280AG36.8% 43.6 / 57.9 57V20V (GTA→GTGkefCpotassium:proton antiporter
Reads supporting (aligned to +/- strand):  ref base A (19/17);  new base G (14/7);  total (33/24)
Fisher's exact test for biased strand distribution p-value = 4.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.21e-01

GGCCCATCATGGATAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGTACGTCTTGGTCTGGGATCGGTACTTGGCTACCTGATCGCCGGCTGCATTATTGGCC  >  W3110S.gb/47761‑47895
                                                                   |                                                                   
ggCCCATCATGGATAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTAcc                                                                   >  1:2945001/1‑70 (MQ=255)
ggCCCATCATGGATAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTAcc                                                                   >  1:1094658/1‑70 (MQ=255)
      tcatGGATAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGCCCATTg                                                              <  1:965919/69‑1 (MQ=255)
      tcatGGATAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGCCCATTg                                                              <  1:2531447/69‑1 (MQ=255)
      tcatGGATAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGCCCATTg                                                              <  1:1790265/69‑1 (MQ=255)
      tcatGGATAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGCCCATTg                                                              <  1:172518/69‑1 (MQ=255)
      tcatGGATAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGCCCATTg                                                              <  1:3256834/69‑1 (MQ=255)
      tcatGGATAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTTACCCATTGCgg                                                           <  1:3040908/71‑1 (MQ=255)
      tcatGGATAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTGCGGCAGCGCTGATTGTGCCCATTg                                                              <  1:503364/69‑1 (MQ=255)
      tcatGGATAGCCATACGCTGATTCAGGCGCTGATATATCTCGGTTCGGCAGCGCTGATTGTGCCCATTg                                                              <  1:3272307/69‑1 (MQ=255)
        atGGATAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCgg                                                           <  1:3291115/70‑1 (MQ=255)
        atGGATAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCgg                                                           <  1:487892/70‑1 (MQ=255)
        atGGATAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCgg                                                           <  1:2054535/70‑1 (MQ=255)
        atGGATAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCgg                                                           <  1:246273/70‑1 (MQ=255)
        atGGATAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGGTCGGCAGCGCTGATTGTACCCATTGCgg                                                           <  1:3215925/70‑1 (MQ=255)
         tGGATAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCgg                                                           <  1:1641517/69‑1 (MQ=255)
             tAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGTACGTc                                                     <  1:455093/71‑1 (MQ=255)
             tAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGTACGTc                                                     <  1:459038/71‑1 (MQ=255)
             tAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGTACGTc                                                     <  1:2706996/71‑1 (MQ=255)
             tAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGTACGTc                                                     <  1:2546234/71‑1 (MQ=255)
             tAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGTACGTc                                                     <  1:2542347/71‑1 (MQ=255)
             tAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGTACGTc                                                     <  1:1343769/71‑1 (MQ=255)
              aGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGCCCa                                                                 >  1:3009611/1‑58 (MQ=255)
              aGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGCCCa                                                                 >  1:492641/1‑58 (MQ=255)
              aGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGCCCa                                                                 >  1:674905/1‑58 (MQ=255)
              aGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGCCCa                                                                 >  1:377586/1‑58 (MQ=255)
              aGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGCCCa                                                                 >  1:752498/1‑58 (MQ=255)
              aGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGCCCa                                                                 >  1:2296830/1‑58 (MQ=255)
              aGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGCCCa                                                                 >  1:1743163/1‑58 (MQ=255)
              aGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGCCCa                                                                 >  1:3141944/1‑58 (MQ=255)
              aGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGCCCa                                                                 >  1:2078267/1‑58 (MQ=255)
              aGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTGCCCa                                                                 >  1:2057594/1‑58 (MQ=255)
                     cGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGTACGTCTTGGTCTg                                             >  1:1906999/1‑71 (MQ=255)
                     cGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGTACGTCTTGGTCTg                                             >  1:2375115/1‑71 (MQ=255)
                     cGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGTACGTCTTGGTCTg                                             >  1:2436276/1‑71 (MQ=255)
                     cGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGTACGTCTTGGTCTg                                             >  1:3357832/1‑71 (MQ=255)
                     cGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGTACGTCTTGGTCTg                                             >  1:1880669/1‑71 (MQ=255)
                     cGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGTACGTCTTGGTCTg                                             >  1:2990748/1‑71 (MQ=255)
                     cGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGTACGTCTTAGTCTg                                             >  1:76227/1‑71 (MQ=255)
                                        ttATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGTACGTCTTGGTCTGGGATCGGTACTTGGCTAcc                          >  1:873994/1‑71 (MQ=255)
                                        ttATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGTACGTCTTGGTCTGGGATCGGTACTTGGCTAcc                          >  1:2549674/1‑71 (MQ=255)
                                        ttATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGTACGTCTTGGTCTGGGATCGGTACTTGGCTAcc                          >  1:2362762/1‑71 (MQ=255)
                                         tATCTCGGTTCGGCAGCGCTGATTGTGCCCATTGCGGTACGTCTTGGTCTGGGATTGGTACTTGGCTACCt                         >  1:1688498/1‑71 (MQ=255)
                                         tATCTCGGTTCGGCAGCGCTGATTGTGCCCATTGCGGTACGTCTTGGTCTGGGATCGGTACTTGGCTAcc                          >  1:770962/1‑70 (MQ=255)
                                         tATCTCGGTTCGGCAGCGCTGATTGTGCCCATTGCGGTACGTCTTGGTCTGGGATCGGTACTTGGCTACCt                         >  1:1021653/1‑71 (MQ=255)
                                         tATCTCGGTTCGGCAGCGCTGATTGTGCCCATTGCGGTACGTCTTGGTCTGGGATCGGTACTTGGCTACCt                         >  1:2268078/1‑71 (MQ=255)
                                         tATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGTACGTCTTGGTCTGGGATCGGTACTTGGCTACCt                         >  1:1693531/1‑71 (MQ=255)
                                         tATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGTACGTCTTGGTCTGGGATCGGTACTTGGCTACCt                         >  1:1035111/1‑71 (MQ=255)
                                         tATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGTACGTCTTGGTCTGGGATCGGTACTTGGCTACCt                         >  1:1113158/1‑71 (MQ=255)
                                         tATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGTACGTCTTGGTCTGGGATCGGTACTTGGCTACCt                         >  1:1524554/1‑71 (MQ=255)
                                                  tCGGCAGCGCTGATTGTACCCATTGCGGTACGTCTTGGTCTGGGATCGGTACTTGGCTAcc                          >  1:1667819/1‑61 (MQ=255)
                                                  tCGGCAGCGCTGATTGTACCCATTGCGGTACGTCTTGGTCTGGGATCGGTACTTGGCTACCTGATCGCCgg                >  1:938306/1‑71 (MQ=255)
                                                  tCGGCAGCGCTGATTGTACCCATTGCGGTACGTCTTGGTCTGGGATCGGTACTTGGCTACCTGATCGCCg                 >  1:1032406/1‑70 (MQ=255)
                                                         cgcTGATTGTACCCATTGCGGTACGTCTTGGTCTGGGATCGGTACTTGGCTACCTGATCGCCGGCTGCatt         <  1:1891390/71‑1 (MQ=255)
                                                         cgcTGATTGTACCCATTGCGGTACGTCTTGGTCTGGGATCGGTACTTGGCTACCTGATCGCCGGCTGCatt         <  1:2421807/71‑1 (MQ=255)
                                                          gcTGATTGTACGCATTGCGGTACGTCTTGGTCTGGGATCGGTACTTGGCTACCTGATCGCCGGCTGCatt         <  1:302580/70‑1 (MQ=255)
                                                                 gTACCCATTGCGGTACGTCTTGGTCTGGGATCGGTACTTGGCTACCTGATCGCCGGCTGCATTATTGGcc  <  1:553039/70‑1 (MQ=255)
                                                                   |                                                                   
GGCCCATCATGGATAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGTACGTCTTGGTCTGGGATCGGTACTTGGCTACCTGATCGCCGGCTGCATTATTGGCC  >  W3110S.gb/47761‑47895

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