Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,426,674 T→C 71.2% D393D (GAT→GAC trkG → potassium transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,426,6740TC71.2% 93.7 / 43.0 59D393D (GAT→GACtrkGpotassium transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (8/9);  new base C (23/19);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 7.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

AAATCAGCTTTCTCATCCCAGAGCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGATCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATACTCTCTTCACGGTGTTTTTTATA  >  W3110S.gb/1426611‑1426737
                                                               |                                                               
aaaTCAGCTTTCTCATCCCAGAGCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCgtgt                                                            <  1:3453314/69‑1 (MQ=255)
 aaTCAGCTTTCTCATCCCAGAGCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCgtgt                                                            <  1:734552/68‑1 (MQ=255)
    cAGCTTTCTCATCCCAGAGCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGATCGTGTAATg                                                        <  1:2915352/69‑1 (MQ=255)
    cAGCTTTCTCATCCCAGAGCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGATCGTGTAATg                                                        <  1:3655926/69‑1 (MQ=255)
    cAGCTTTCTCATCCCAGAGCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATg                                                        <  1:2258361/69‑1 (MQ=255)
    cAGCTTTCTCATCCCAGAGCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATg                                                        <  1:3594594/69‑1 (MQ=255)
    cAGCTTTCTCATCCCAGAGCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATg                                                        <  1:399438/69‑1 (MQ=255)
    cAGCTTTCTCATCCCAGAGCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCGGGTAATg                                                        <  1:1291166/69‑1 (MQ=255)
         ttCTCATCCCAGAGCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCgtgt                                                            >  1:2510371/1‑60 (MQ=255)
         ttCTCATCCCAGAGCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCgtgt                                                            >  1:2328821/1‑60 (MQ=255)
         ttCTCATCCCAGAGCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCgtgt                                                            >  1:501147/1‑60 (MQ=255)
         ttCTCATCCCAGAGCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCgtgt                                                            >  1:1397418/1‑60 (MQ=255)
         ttCTCATCCCAGAGCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCgtgt                                                            >  1:630647/1‑60 (MQ=255)
         ttCTCATCCCAGAGCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCgtgt                                                            >  1:782831/1‑60 (MQ=255)
            tcATCCCAGAGCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGATCGTGTAATGAGGTCTg                                                 <  1:468666/68‑1 (MQ=255)
                cccAGAGCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATgg                                            <  1:2668316/69‑1 (MQ=255)
                cccAGAGCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATgg                                            <  1:2228652/69‑1 (MQ=255)
                 ccAGAGCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGa                                           >  1:2565422/1‑69 (MQ=255)
                          tgttgAGTGTAAATGTAGGAGGGAAGATAGTTACAGATCGTGTAATGAGGTCTGTATGGAGtttctttt                                  >  1:513786/1‑69 (MQ=255)
                          tgttgAGTGTAAATGTAGGAGGGAAGATAGTTACAGATCGTGTAATGAGGTCTGTATGGAGtttcttt                                   >  1:1465541/1‑68 (MQ=255)
                          tgttgAGTGTAAATGTAGGAGGGAAGATAGTTACAGATCGTGTAATGAGGTCTGTATGGAGtttcttt                                   >  1:1374717/1‑68 (MQ=255)
                          tgttgAGTGTAAATGTAGGAGGGAAGATAGTTACAGATCGTGTAATGAGGTCTGTATGGAGtttcttt                                   >  1:1856900/1‑68 (MQ=255)
                              gAGTGTAAATGTAGGAGGGAAGATAGTTACAGATCgtgt                                                            >  1:195154/1‑39 (MQ=255)
                              gAGTGTAAATGTAGGAGGGAAGATAGTTACAGATCgtgt                                                            >  1:748226/1‑39 (MQ=255)
                                gtgtAAATGTAGGAGGGAAGATAGTTACAGATCGTGTAATGAGGTCTGTATGGAGtttctttt                                  >  1:463744/1‑63 (MQ=255)
                                gtgtAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGtttcttttttctt                             <  1:2021141/68‑1 (MQ=255)
                                gtgtAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGtttcttttttctt                             <  1:900888/68‑1 (MQ=255)
                                 tgtAAATGTAGGAGGGAAGATAGTTACAGATCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTta                           <  1:3551467/69‑1 (MQ=255)
                                 tgtAAATGTAGGAGGGAAGATAGTTACAGATCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTta                           <  1:573028/69‑1 (MQ=255)
                                 tgtAAATGTAGGAGGGAAGATAGTTACAGATCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTta                           <  1:1172118/69‑1 (MQ=255)
                                 tgtAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTta                           <  1:2094186/69‑1 (MQ=255)
                                   tAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGtttcttttttcttt                            >  1:1299534/1‑66 (MQ=255)
                                   tAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTta                           >  1:1212901/1‑67 (MQ=255)
                                     aaTGTAGGAGGGAAGATAGTTACAGATCGTGTAAt                                                         >  1:1318512/1‑35 (MQ=255)
                                     aaTGTAGGAGGGAAGATAGTTACAGATCGTGTAAt                                                         <  1:2890702/35‑1 (MQ=255)
                                     aaTGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATAct                       <  1:1486919/69‑1 (MQ=255)
                                     aaTGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATAct                       <  1:3391548/69‑1 (MQ=255)
                                     aaTGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATAct                       <  1:3197936/69‑1 (MQ=255)
                                     aaTGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATAct                       >  1:3064383/1‑69 (MQ=255)
                                     aaTGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATAc                        >  1:2368093/1‑68 (MQ=255)
                                         tAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATActctc                    >  1:640842/1‑68 (MQ=255)
                                         tAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATActctc                    >  1:636121/1‑68 (MQ=255)
                                          aggaggGAAGATAGTTACAGATCGTGTAATGAGGTCTg                                                 <  1:1744121/38‑1 (MQ=255)
                                               ggAAGATAGTTACAGATCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATACTCTCTTCACgg             <  1:2649119/69‑1 (MQ=255)
                                                 aaGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGtttcttttttctt                             >  1:2047257/1‑51 (MQ=255)
                                                 aaGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGtttcttttttctt                             >  1:1742961/1‑51 (MQ=255)
                                                 aaGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATACTCTCTTCACg              >  1:103701/1‑66 (MQ=255)
                                                 aaGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATACTCTCTTCACGgtg           >  1:1582786/1‑69 (MQ=255)
                                                 aaGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATACTCTCTTCACGgtg           >  1:1571146/1‑69 (MQ=255)
                                                 aaGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATACTCTCTTCACGgtg           >  1:3323062/1‑69 (MQ=255)
                                                 aaGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATACTCTCTTCACGgtg           >  1:2614036/1‑69 (MQ=255)
                                                     tAGTTACAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATACTCTCTTCACGGTGttt        >  1:1717234/1‑68 (MQ=255)
                                                      aGTTACAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATACTCTCTTCACGGTGttttt      <  1:1329056/69‑1 (MQ=255)
                                                      aGTTACAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATACTCTCTTCACGGTGttttt      >  1:2021281/1‑69 (MQ=255)
                                                      aGTTACAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATACTCTCTTCACGGTGttttt      <  1:1298521/69‑1 (MQ=255)
                                                      aGTTACAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATACTCTCTTCACGGTGttttt      <  1:2420284/69‑1 (MQ=255)
                                                         tACAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATACTCTCTTCACGGTGTTTTTta    <  1:1554168/68‑1 (MQ=255)
                                                         tACAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATACTCTCTTCACGGTGTTTTTta    <  1:2484533/68‑1 (MQ=255)
                                                          aCAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATACTCTCTTCACGGTGTTTTTtata  >  1:3325738/1‑69 (MQ=255)
                                                               |                                                               
AAATCAGCTTTCTCATCCCAGAGCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGATCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATACTCTCTTCACGGTGTTTTTTATA  >  W3110S.gb/1426611‑1426737

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: