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

TGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGG‑TTGAATTGTTGTCTTCAA‑CCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATC‑CGTACG‑CTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAAC‑AATTCAGGTAATAAAT‑CGAGTTGTTCCTG‑CGTTGCTG‑CCGCCAG‑TCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCAT‑ATTTAACGCCAGT  >  NC_000913/2151580‑2151903
                                                                                                           |                                                                                                                                                                                                                                  
tGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGG‑GTGAATTGTTGTTTTTAA‑CCCCCATgcga                                                                                                                                                                                                              >  1:1918815/1‑128 (MQ=255)
    ccTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGG‑TTGAATTGTTTTTTTTAA‑ACGCCCTGAGAAAACCACTCTTCCCACAGCGCCCGTGCCAACGGTTTATTATAATC‑CGttag                                                                                                                                                          >  1:75726/1‑172 (MQ=255)
    ccTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGG‑TTGAATTGTTGTCTTCAA‑CCGCCATGAGAAAATCACTCTGCCCACAGCGCCAGTGCCAACGGTTTATCATCATC‑CGTACG‑CTCGTTTACCGCCGGGCTGGACAAAAAATGCTTTTAC‑AATTCAGGTAAAAAAA‑CGAgttgtt                                                                                         >  1:1529489/1‑237 (MQ=255)
       ttATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGG‑TTTAATTTTTTTCTTCTA‑CCCCCCTTGGAAAATCACCCTGCCCACCGCCCCAGTGcccaac                                                                                                                                                                              >  1:2222521/1‑149 (MQ=255)
            ggTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGG‑TTTAATTTTTTTCTTCCA‑ACGCCCTTAGAAAATCACTCTTCCCCCCACGCCCAGGCCCACGTTTTATCATCATC‑CCTACC‑CTCGTTTACAGCCCGGGTGGGCACAAAATCCCTTAAC‑ACTTCAGGTAATAAAT‑CTATGTGTTCCTgag                                                                                   >  2:1832768/1‑233 (MQ=255)
             gTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTTGAATT‑TTGTCTTTTA‑CCCCCCTGAGAAACTCACACTGcccccccc                                                                                                                                                                                           >  1:2196733/1‑129 (MQ=255)
                aCCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGG‑TTGAATTGTTGTCTTCAA‑CCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTTTCATCATC‑CGTACG‑CTCGTTTACCCCCGGGCTGGAc                                                                                                                                   >  1:268943/1‑185 (MQ=255)
                     ccaaaaaaaAGAAGCCGAGAAGATAAGAGAAAAATTATTTCCCTAATCACCCTTTATTTTTTTTTTTAAAACTGGGGG‑TTAATTTTTTGTCTTCAA‑CCGCCATGAGAAactc                                                                                                                                                                                                         <  1:1062124/109‑1 (MQ=255)
                            gaagaagCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGG‑TTGAATTGTTTTCTTCAA‑CCGCCCTTAGAAACTCACTCTGCCCCCAACGCCAGTGCCCACGTTTTTTCATCATC‑CCTACG‑CTCGTTTACCGCCGGGCCGGGCCAAAAATCCTTTTAA‑AATTCAGGGAATAAAA‑CGAGTTGTTTCTT‑CGTTGCTt‑ccgcccggtccatg                                                             >  2:735320/1‑230 (MQ=255)
                               gaagCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGG‑TTGAATTGTTGTCTTCAA‑CCGCCCTGAGAAACTCACTCTGCCca                                                                                                                                                                                               >  2:1062124/1‑112 (MQ=255)
                                  gCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGG‑TTGAATTGTTTTCTTCAA‑CCCCCATTAGAAAATCACTCTGCCCCCAGCGCCAGTGCCAACGCTTTATCCTCATC‑CCTACG‑CTCGTTTAcccc                                                                                                                                             >  1:1121087/1‑155 (MQ=255)
                                             aaGAGCAAAATTAGTTCGCTAATAAGCCGTTATATTTATTGTTTAAACTGGCGG‑TTTGATTGTTTTCTTTAA‑TCCCCACGTTAAAAATACTCTCTCCCAAGGCCCAGGGCCAAAGATTTTTTAT‑ATCACCTACGCCTCG‑TTACCACCCGGCTGGGCATAAAATACTTTAAC‑CATTTAAGGAAGTAAT‑ACAGTTGTTTCTGACGTGGCTC‑CGGCCGT‑GCCTTTAAAAAATGGGGTaacga                                            >  2:429147/1‑239 (MQ=255)
                                                    aaaTTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGG‑TTGAATTTTTTTCTTCAA‑ACCCCCTGAGAAAATCACTCTTCCCCCAGCGCCCGGGGCCACGTTTTTTCATCATC‑CCTACC‑CTCGTTTACCGCCGGGCTGGACAAAAAATGCTTTTAC‑AATTTCGGGAATAAAT‑AGAGTTGTTCCTG‑CGTTTCTG‑CCGCCCC‑TCCCTTTAAAAAACGGGTGAAGAAAAGcaca                                    >  2:2147639/1‑239 (MQ=255)
                                                          gTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGG‑TTTGATTTTTTTCTTTAA‑CCGCCCTGAGAAAATCACCTTTCCCAAAGCGCCCGAGTCAAAGTTTTTTTATCTAC‑CATACG‑ACCGGTTACCCCCCGGGTGGACAAAAAAATCTTTAACAAATTATGGTAA‑AAAT‑CGAGTTGTTCCTT‑CCGTGTCT‑CCGCCAG‑TTCACTTAAAAATGGGGTCAATAAATGcacaagtgta                              >  2:2255945/1‑233 (MQ=255)
                                                          gTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGG‑TTTAATTTTTTTCTTCCA‑ACCCCCTTTGAAAATCACTCTGCCCCCAGCGCCAGGGCCCCCGTTTTTTCCTCATC‑CCGGCG‑CTCCTTTCCCCCCCGGGTGGGCAAAAACTTCTTTTAC‑AAATTAGGGAATAAAT‑CCCGTTGTTCCTG‑CCGTTCTG‑CCGCCAG‑ACCCTTTAAAAAAGGGGTAAAGAAAGGAACACtagga                              >  2:1660698/1‑239 (MQ=255)
                                                          gTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGG‑TTGAAATGTTGTCTTCAA‑CCGCCCAGGGAAACTCACTCTGCCCACAACGCCCGTGCCAAc                                                                                                                                                                               >  1:2295704/1‑101 (MQ=255)
                                                             cGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGG‑TTTTATTGTTGTCTTCCA‑CACCCATGCGGAACATACTCTTCCca                                                                                                                                                                                               >  1:1056369/1‑82 (MQ=255)
                                                               cTAATCAGCCGTTATATTTATTGTTTAAACTGGCGG‑TTGTATTTTTTTCTTCAA‑CCGCCATGAGAAAATCACTCTGCCCCCAGCCCCAGGGCCAAAGTTTTATCATCATC‑CGTACC‑CTCGTTTTCCGCCGGGCTCGGCACAAACTGCTTTTAC‑AATTTCCGTAATAAAT‑CGAGCTGTTCCTC‑CGTTTCTT‑CCCCCGA‑CCTATTTATAAACTGGGGAAAAAAAGAC‑CATtaggaggaaat                        >  2:1483821/1‑234 (MQ=255)
                                                                       tcattATTTTTGTTGTTAAAGTTGGGGG‑TTGTGTT‑TTTTCCTCAACCCGCCATGAGAAACTCCCTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATC‑CGTACG‑CTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAAC‑AATTCAGGTAATAAAT‑CGAGTTGTTCCTG‑CGTTGCTG‑CCGCCAG‑TCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGc                  <  1:1936503/235‑1 (MQ=255)
                                                                            atatTTTTTGTTTAAACGGGCGG‑TTGAATTGTTTTCTCCAA‑CCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATC‑CGTACG‑CTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAAC‑AATTCAg                                                                                                            <  1:1581486/147‑1 (MQ=255)
                                                                            atatTTATTGTTTAAACTGGCGG‑TTGAATTTTTTTTTTCAA‑CCGCCCTTAGAAAATTACTCTGCCCCCAGCGCCAGTTCCAACCTTTTTTTCTCATC‑CCTACG‑CTCGTTTACCACCGGGGTGGGCAAAAACACCTTTTAC‑AATTCATGGAATTAAA‑CGCGTTTGTCCTG‑CGTTGTTC‑ACGCCAG‑GTCATTTAAAAACGGGGGAAAGAAAAGCACAATGGGGGAAATAACCACCAT‑Attt           >  2:960917/1‑239 (MQ=255)
                                                                                 tattGTTTAAACTGGCGG‑TTGAATTGTTGTTTTCAA‑CCCCCATGAGAAAATCACTCTGCCCACAGCGCCCGTGCCAACGTTTTTTTATCATT‑CCTACG‑CCCGTTTACCCCCGGGGCGGGCAAAAAACTCTTTAAC‑AATTCAGGGAATAAAA‑CGAATTGTTTCCG‑CGTTTCTG‑TCGCCAGATCATTTGAAAAACGGGTAAAAAAAAGGCA‑ATTAGGGGAATTATCCAGCAG‑ATTTAAcc       >  2:1936503/1‑237 (MQ=255)
                                                                                    tGTTTAAACTGGCGG‑TTGAATTTTTTTTTTTAA‑CCCCCCTGAGAAACTCACTCTTCCCCCAGCGCCCGTGCCAACg                                                                                                                                                                              >  1:1238404/1‑76 (MQ=255)
                                                                                    tGTTTAAACTGGCGG‑TTGAATTGTTTTCTTCAA‑CCCCCATTAAAAAATTACTTTTCCCACAGCGCCAGTGGCCACGGTTTTTCATCATC‑CGTACG‑CTCCTTTAACACCCGGGCGGACCAAAACTGCTTTTAC‑AATTTAAGGAATAAAA‑CCAGTTGTTTCTG‑CGTTTCTG‑CCGGCAA‑TTCATTTTAAAAAGGGGGAAAAAAAAGCACAATAGGAGAAGTATACATCAC‑AATTaaaagaca   >  2:1106475/1‑233 (MQ=255)
                                                                                    tGTTAAATCTGGTGG‑TGTATTTTTTTTTTTCAC‑CCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATC‑CGTACG‑CTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAAC‑AATTCAGGTAATAAAT‑CGAGTTGTTCCTG‑CGTTGCTG‑CCGCCAG‑TCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCAT‑ATTTAACGCCAg   <  1:884160/239‑1 (MQ=255)
                                                                                      tttAAGCTGGCGG‑TTGGATTTTTGTTTTTAA‑CCGCCCTGAGAAAATCACCCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATC‑CGTACG‑CCCGTTTACCGCCGGGCCGGACAAAAACTGCTTTAAC‑AAATCAGGTAATAAAT‑c                                                                                                 >  2:1581486/1‑147 (MQ=255)
                                                                                      ttaaatcgggggg‑TTGATTTTTTTTCTCCAC‑CCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATC‑CGTACG‑CTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAAC‑AATTCAGGTAATAAAT‑CGAGTTGTTCCTG‑CGTTGCTG‑CCGCCAG‑TCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCAT‑ATTTAACGCCAGt  <  1:299007/227‑1 (MQ=255)
                                                                                       ttAAACTGGCGG‑TTGAATTGTTGTCTTTAA‑CCGCCCTGAGAAAATCACTCTCCCCCCCACCCCAGGTCCCACGGTTTATCATTCTC‑CCGACG‑CCCGTTTTACCCCGGGGTGGACAAAAAATGCTTTTAC‑AATTTAGGGAATAAATCCGAG‑TGTTCCCT‑CGTTTCTGCCCCCCCG‑CCCTTTAAAAACCGGGGAAAAAAAAGGCCCTTG‑GGAGAAATTTTCAACATAATTTAAcccagac  >  2:2208600/1‑232 (MQ=255)
                                                                                       ttAAACTGGCGG‑TTGAATTGTTGTCTTCAA‑CCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCCAAGGTTTTTCATCATC‑CGTACG‑CTCGTTTACCCCCGGGCTGGAAAAAAACTGCttttaaa                                                                                                                   >  1:1689256/1‑126 (MQ=255)
                                                                                                           |                                                                                                                                                                                                                                  
TGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGG‑TTGAATTGTTGTCTTCAA‑CCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATC‑CGTACG‑CTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAAC‑AATTCAGGTAATAAAT‑CGAGTTGTTCCTG‑CGTTGCTG‑CCGCCAG‑TCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCAT‑ATTTAACGCCAGT  >  NC_000913/2151580‑2151903

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