Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,426,6740TC72.4% 93.3 / 41.1 58D393D (GAT→GACtrkGpotassium transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (0/16);  new base C (15/27);  total (15/43)
Fisher's exact test for biased strand distribution p-value = 5.58e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ATAAATCAGCTTTCTCATCCCAGAGCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGATCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATACTCTCTTCACGGTGTTTTTTA  >  W3110S.gb/1426609‑1426735
                                                                 |                                                             
aTAAATCAGCTTTCTCATCCCAGAGCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCgtgt                                                          <  1:2605806/71‑1 (MQ=255)
aTAAATCAGCTTTCTCATCCCAGAGCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCgtgt                                                          <  1:1279014/71‑1 (MQ=255)
aTAAATCAGCTTTCTCATCCCAGAGCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCgtgt                                                          <  1:2520064/71‑1 (MQ=255)
 tAAATCAGCTTTCTCATCCCAGAGCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCgtgt                                                          <  1:1826529/70‑1 (MQ=255)
     tCAGCTTTCTCATCCCAGAGCGTTGTTGAGTGTAAATGTAGGAGGGAAGATTGTTACAGACCgtgt                                                          <  1:2953661/66‑1 (MQ=255)
     tCAGCTTTCTCATCCCAGAGCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCgtgt                                                          <  1:1843720/66‑1 (MQ=255)
            tctcATCCCAGAGCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGATCGTGTAATGAGGTCTg                                               <  1:25617/70‑1 (MQ=255)
            tctcATCCCAGAGCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGATCGTGTAATGAGGTCTg                                               <  1:1742582/70‑1 (MQ=255)
            tctcATCCCAGAGCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGATCGTGTAATGAGGTCTg                                               <  1:1984547/70‑1 (MQ=255)
                       agCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCgtgt                                                          <  1:1637059/48‑1 (MQ=255)
                        gCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGtttctt                                  <  1:86385/71‑1 (MQ=255)
                        gCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGtttctt                                  <  1:769656/71‑1 (MQ=255)
                        gCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGtttctt                                  <  1:795627/71‑1 (MQ=255)
                        gCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGtttctt                                  <  1:2181878/71‑1 (MQ=255)
                        gCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGtttctt                                  <  1:3019370/71‑1 (MQ=255)
                        gCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGtttctt                                  <  1:3033391/71‑1 (MQ=255)
                        gCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGtttctt                                  <  1:3175607/71‑1 (MQ=255)
                        gCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGtttctt                                  <  1:3395262/71‑1 (MQ=255)
                          gttgttGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGtttctttt                                >  1:307037/1‑71 (MQ=255)
                          gttgttGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGtttctttt                                >  1:1267524/1‑71 (MQ=255)
                          gttgttGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGtttctttt                                >  1:3700994/1‑71 (MQ=255)
                          gttgttGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGtttctttt                                >  1:3985214/1‑71 (MQ=255)
                          gttgttGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGtttctttt                                >  1:258515/1‑71 (MQ=255)
                          gttgttGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGtttctttt                                >  1:2196670/1‑71 (MQ=255)
                          gttgttGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGtttctttt                                >  1:2039064/1‑71 (MQ=255)
                                 aGTGTAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTta                         <  1:727401/71‑1 (MQ=255)
                                 aGTGTAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTta                         <  1:671454/71‑1 (MQ=255)
                                 aGTGTAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTta                         <  1:3768883/71‑1 (MQ=255)
                                   tgtAAATGTAGGAGGGAAGATAGTTACAGATCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTtata                       <  1:1229433/71‑1 (MQ=255)
                                   tgtAAATGTAGGAGGGAAGATAGTTACAGATCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTtata                       <  1:4076120/71‑1 (MQ=255)
                                   tgtAAATGTAGGAGGGAAGATAGTTACAGATCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTtata                       <  1:4009763/71‑1 (MQ=255)
                                   tgtAAATGTAGGAGGGAAGATAGTTACAGATCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTtata                       <  1:3336915/71‑1 (MQ=255)
                                   tgtAAATGTAGGAGGGAAGATAGTTACAGATCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTtata                       <  1:2731004/71‑1 (MQ=255)
                                   tgtAAATGTAGGAGGGAAGATAGTTACAGATCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTtata                       <  1:3062510/71‑1 (MQ=255)
                                   tgtAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTgatc                       >  1:2104109/1‑70 (MQ=255)
                                   tgtAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTgatc                       >  1:825719/1‑70 (MQ=255)
                                   tgtAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTgatc                       >  1:2903592/1‑70 (MQ=255)
                                   tgtAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTgatc                       >  1:444801/1‑70 (MQ=255)
                                   tgtAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTgatc                       >  1:1814236/1‑70 (MQ=255)
                                   tgtAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTgatc                       >  1:1635458/1‑70 (MQ=255)
                                   tgtAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTATTCTTgatc                       >  1:3431345/1‑70 (MQ=255)
                                    gtAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGtttctt                                  <  1:2787728/59‑1 (MQ=255)
                                     tAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATAc                      <  1:2034428/70‑1 (MQ=255)
                                     tAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATAc                      <  1:3502673/70‑1 (MQ=255)
                                      aaaTGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATAc                      <  1:611682/69‑1 (MQ=255)
                                                 ggAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGtttctt                                  <  1:1812448/46‑1 (MQ=255)
                                                      aTAGTTACAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATACTCTCTTCACGGTGttttt    <  1:3451986/71‑1 (MQ=255)
                                                      aTAGTTACAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATACTCTCTTCACGGTGttttt    <  1:601629/71‑1 (MQ=255)
                                                      aTAGTTACAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATACTCTCTTCACGGTGttttt    <  1:743156/71‑1 (MQ=255)
                                                      aTAGTTACAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATACTCTCTTCACGGTGttttt    <  1:1590553/71‑1 (MQ=255)
                                                        aGTTACAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATACTCTCTTCACGGTGTTTTTta  >  1:3467264/1‑71 (MQ=255)
                                                          ttACAGATCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATACTCTCTTCACGGTGTTTTTta  <  1:469454/69‑1 (MQ=255)
                                                          ttACAGATCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATACTCTCTTCACGGTGTTTTTta  <  1:1607119/69‑1 (MQ=255)
                                                          ttACAGATCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATACTCTCTTCACGGTGTTTTTta  <  1:794334/69‑1 (MQ=255)
                                                          ttACAGATCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATACTCTCTTCACGGTGTTTTTta  <  1:1584075/69‑1 (MQ=255)
                                                          ttACAGATCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATACTCTCTTCACGGTGTTTTTta  <  1:3821906/69‑1 (MQ=255)
                                                          ttACAGATCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATACTCTCTTCACGGTGTTTTTta  <  1:3791860/69‑1 (MQ=255)
                                                          ttACAGATCGTGTAAAGAGGTCTGTATGGAGTTTCTTTTTTCTTTATACTCTCTTCACGGTGTTTTTta  <  1:3440438/69‑1 (MQ=255)
                                                                 |                                                             
ATAAATCAGCTTTCTCATCCCAGAGCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGATCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATACTCTCTTCACGGTGTTTTTTA  >  W3110S.gb/1426609‑1426735

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: