Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,151,6860GT51.5% ‑5.0 / 22.3 33intergenic (+40/+25)yegJ/yegKuncharacterized protein/ser/thr phosphatase‑related protein
Reads supporting (aligned to +/- strand):  ref base G (12/4);  new base T (9/8);  total (21/12)
Fisher's exact test for biased strand distribution p-value = 2.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.49e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGC  >  NC_000913/2151579‑2151906
                                                                                                             |                                                                                                                                                                                                                                      
ttGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTTAATTTTTTTCTTTAACCGCCCTTAGAAACTCAATCTGCCCACAGCGCCCGTGCCACAGTTTTTTCATCATCCGTACCCCCCCTTATCCCCCGGCTGGACAAAAAATGCTTTAACAAATTAGGGAAATAAATCAAgt                                                                                               >  2:1280289/1‑239 (MQ=255)
 tGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAACGCCAGTGCCCACGtttttt                                                                                                                                                                          >  1:1236309/1‑164 (MQ=255)
  gATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTTTTTTCTTTTACCGCCCTGAGAAac                                                                                                                                                                                                             >  1:2335928/1‑130 (MQ=255)
  gATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTTTTTTCTTCCAACGCCCAGAGAAAATCACTCTTCCCCCACCGCCCGTGGCCACGGTTTTTCAACATCCCTACGCTCGTTTACCCCCCGGCTGGACAAAAAATGCTTTTAAACATTAGGTAATAAATCAAgttgt                                                                                            >  2:1180809/1‑239 (MQ=255)
  gATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCCGGCCGGACAAAAACTGCTTTAACAATTCAGGTAATAAAACCAgttgt                                                                                            >  1:1642386/1‑239 (MQ=255)
         tATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTTGATTGTTGTTTTTAACCCCCATGAGAAAATCACTCTTCCCACCGAGCCCGTGGCCACGTTTTTTTATTATCCGTACGCTCGTTTTCTCCCCGGCCGGACCACAAATCCTTTAACAATTCACGGTATAAAATCAGTTGTTCCcgcg                                                                                     >  2:1375740/1‑239 (MQ=255)
              gTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGCTTGAATGTTTTTTTTAAACCCCAAGAGAAACTCAACTTCCCCCCCACGCCCGTGGCCACGTGCTTTTAATAACCCCACCCCCGTTTTCCCCCCGGCTCGTCCAAAACCGCCTTTACAATTTCCGTAATAAATCTAATTTTTCTTGCGTcgct                                                                                >  2:1215874/1‑239 (MQ=255)
              gTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTTAATTTTTTTTTTCAACCGCCTTTAGAAACTCATCCTTCCCCCAACTACCGTGCCTAACCTTTATCATCAACCTTACCCCTCCTTCCCGCCGGGCTTGACAAAAACTGCTTTAACCAATCCGGGAATAAAGACAGTTGTTTCTg                                                                                       >  2:104845/1‑232 (MQ=255)
              gTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTTTTTAACCGCCAAGAGaaa                                                                                                                                                                                                              >  1:848699/1‑117 (MQ=255)
              gTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGGTGAATTGTTTTCTTCCAACCCCAGTTGAAAATCAATCTGCCCACCACGCCAGGCGCCACCTTTTATCATCATCCGTACGCTCGCTTACAGCCGGGCGGGACAAAAAAGGCTTTAACAATTAACGTACTAAATTCAgtcgtttccgcat                                                                                    >  2:1936973/1‑228 (MQ=255)
                                            cTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTTTTTTTTTCAACCCCCATTGGAAAATCAATCTGCCCCCATCGCCAGTGCCAACGTTTTATACTCCTCCTGTCGCTCGTTTACCCCCCGGCCGGGCAAAAAACGCTTTAACAATTCAGGGAATAAATCACATTGTGTCTGTGGTTTTTCCCCCCACCCATTAGAAAAACGGGGTaaa                                                 >  2:458491/1‑239 (MQ=255)
                                            atAAGAGAAAAATTATTTCGAAAATCGGCCTTTATATTTATTGTTAAACCTGGGGTTTGATTTTTTGTTTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGc                                                                                                                                           <  1:970500/152‑1 (MQ=255)
                                              aaGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTTTTTAACCGCCATGAGAAAATCAATCTGCCCACAACCCCAGTGCCAAAGGTTTATTATTATTCCCACGGTCGTTTCTCGCCTGGGCGGGAAAAAAATGCCTTAAACAATCAAGGtataa                                                                                                       >  2:507178/1‑182 (MQ=255)
                                                     aaaTTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTTTTTTCTTCAACCCCCATTAGAAAATCACTCTGCCCCCAGCGCCAGGGGCAACGTTTTTTTATCATCCCTACGCTCGTTTACCGCCCGGCGGGGCAAAAACTGCTTTAACAAATTAGGGAATAAATCGAGTTGTTCCTGAGTTGCTGCCGCCACTCCAATGAAAAAAGGGGTAAAGAAAAGcaca                                      >  2:367181/1‑239 (MQ=255)
                                                      aaTTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTTTTTTTAACCGCCCTGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGGTTTTTCATCAACCCTACGCCCGTTTTCCGCCGGGGGGGACaaaaa                                                                                                                                >  2:970500/1‑153 (MQ=255)
                                                           gTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTTTTTTCTTCAAACGCCCTGAGAAAATCACTCTTCCCACAGCGCCAGTGCCAAAGATTCATTATCTACCCTACACCCCTTTTCCACCCGGCGGGGGAAAAAATGCGTTTACAAATTAGGTAATAAAATGAGTGTTTCTTGATGGTCTGGCGCCAGCCCACTTAAAAACGGGGTAGAAGAAACGACATgagga                                >  2:1871307/1‑239 (MQ=255)
                                                                 tAATCAGCCTTTATTTTTTTTTTTAAAACGGGGGTTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGtt                                                                                                                                                                              <  1:246099/98‑1 (MQ=255)
                                                                   atccccctttATTTTTTTTTTTTAAACGGGGGTTAAAATTTTTTTCTTCACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATc                        <  1:2108269/231‑1 (MQ=255)
                                                                   aTCACCCTTTATTTTTTTTGTTAAAACGGGCGTTGAAATTGTTGTCTTAACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAg                                                                                                <  1:606614/170‑1 (MQ=255)
                                                                        cttttATTTTTTTTGTTAAACCGGGGGTTTGATTTTTTGTCTCCACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGc                    <  1:1375740/235‑1 (MQ=255)
                                                                        ccGTTATAATTATTGTTTAAACTGGCGGTTGAAttgtt                                                                                                                                                                                                                                      >  2:2254363/1‑38 (MQ=255)
                                                                         ttttATATTTTTGTTTAAAGTGGGGTGTTGATTTTTTTTCTTCACCCGCCATGAGAAACTCACGCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGc                    <  1:1419041/235‑1 (MQ=255)
                                                                           ttATTTTTTTTTTTAAACCGGGGGGTTAATTTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCcac                                                                                                                                                                                                <  2:659808/70‑1 (MQ=255)
                                                                           ttATATTTATTGTTTAAACTGGCGGTTGAATTTTTTTTTTCAACCGCCCTTAGGAAATCACTCTGCCCCCAGCGCCAGTTCCCACCTTTTttcatcat                                                                                                                                                                    >  2:246099/1‑98 (MQ=255)
                                                                             atatTTATTGTTTAAACTGGCGGTTGAATTGTTGTTTTTCAACCCCCTGAGAAAATCACTCTTCCCACAGCGCCAGTGCCAACGTTTTTTCATCATCCGTACGCCCGTTTTCCGCCGGGCTGGACAAAAACTGCTTTTACAATTCAGGGAATAAAACGAGTTGTTCCTGc                                                                                      >  2:606614/1‑170 (MQ=255)
                                                                               tttttttttTAAAACGGGGGGTTGATTTTTTTCTTACAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACg                                                         <  1:1488634/195‑1 (MQ=255)
                                                                                ttttttttttAACCGGGGGTTTGATTTGTTGTCTCCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTatcatc                                                                                                                                                                      <  1:1362926/91‑1 (MQ=255)
                                                                                     tGTTTAAACTGGCGGTTTAATTTTTGTCTTCAACCGCCCTGAGAAAATCACTCTGCCCCCAGCCCCCGTGCCAACGTTTTTTTATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAAATGCTTTAACAAtt                                                                                                                 >  1:43929/1‑136 (MQ=255)
                                                                                     tGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCCCAGCGCCCGTg                                                                                                                                                                                      >  1:659808/1‑70 (MQ=255)
                                                                                      ttttAAACGGGGGGTTTAATTTTTTTTCTCAACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAg      <  1:1603026/238‑1 (MQ=255)
                                                                                      tgttAAAAGGGGGGTTGAATTTGTTTCTTCAATCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAg      <  1:536202/237‑1 (MQ=255)
                                                                                      gTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCCCCCTGAGAAAATCCATCTCCCCACAGCGCCCGTGCCAACGCTTTTTCCTCATCCGTACGCGCGTTTTCCGCCGGGCTGGACACAAAATTCCTTTAAAAATTAGGTGATAAATGCAGTTGTTCCTGCGTTGATGCCGCCACTCACTTTAAAAAAGGGGGAAAGAAAGGACACAGAGGGGAATTTTCCAGCATATTTAActgcaac     >  2:2010500/1‑233 (MQ=255)
                                                                                          aaaCTGGCGGTTGAATTGTTGTCTTTAACCGCCCTGAGAAAATCACTCTTCCCCCAGCGCCCGTGCCAACGTTTTATCATCATCCGTACGc                                                                                                                                                            >  2:1362926/1‑91 (MQ=255)
                                                                                          aaaCCGGGGGTTGAATTTTTTTTCTCAACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGc  <  1:833681/239‑1 (MQ=255)
                                                                                                             |                                                                                                                                                                                                                                      
TTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGC  >  NC_000913/2151579‑2151906

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: