Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,624,5140AT38.5% 28.3 / 12.6 26I491F (ATT→TTT) ppkpolyphosphate kinase
Reads supporting (aligned to +/- strand):  ref base A (3/13);  new base T (10/0);  total (13/13)
Fisher's exact test for biased strand distribution p-value = 1.08e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGCCAACATTCACTGGGCGAAGCGCCTGACCGAAGCAGGCGTGCACGTTATCTTCTCTGCGCCGGGGCTGAAAATTCACGCCAAACTGTTCCTGATTTCACGTAAAGAAAACGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGG‑AACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTG‑AAAACCCATACCGTCCGGTGAC‑ATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGT‑ATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGC‑AAGGGCTGCCCAGTGGTATCACCCTGAAGCTAAATAACCTTGTCGATAAA‑GGCCTGGTTGATCGTCTGTATGCGGCCTCCAGCTCCGGCGTACCGGTTAATCTGCTGGTTCGCGGAATGTGTTCGCTGATCCC  >  NC_000913/2624270‑2624755
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      
agCCAACATTCACTGGGCGAAGCGCCTGACCGAAGCAGGCGTGCACGTTATCTTCTCTGCGCCGGGGCTGAAAATTCACGCCAAACTGTTCCTGATTTCACGTAAAGAAAACGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGG‑AACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCTATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTG‑aaa                                                                                                                                                                                                                                                 <  1:410197/251‑1 (MQ=255)
                               gAAGCAGGCGTGCACGTTATCTTCTCTGCGCCGGGGCTGAAAATTCACGCCAAACTGTTCCTGATTTCACGTAAAGAAAACGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGG‑AACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTTACTTTTTTT‑AAAACCCATACCCTCCGGGGAC‑aatttatttttt                                                                                                                                                                                                                 >  1:30019/1‑240 (MQ=255)
                                aaGCAGGCGTGCACGTTATCTTCTCTGCGCCGGGGCTGAAAATTCACGCCAAACTGTTCCTGATTTCACGTAAAGAAAACGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGG‑AACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTTACTTTTTTT‑TAAAACCATACCCGCCGGGGGC‑ATTTGATTAtttt                                                                                                                                                                                                                >  1:132795/1‑250 (MQ=255)
                                                      ctctGCGCCGGGGCTGAAAATTCACGCCAAACTGTTCCTGATTTCACGTAAAGAAAACGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGG‑AACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGGTTTAAAATTTTTT‑AAAAACCATACCCCCCGGGGAG‑ATTTGTTTATTTTATGGTaaagcccccaaaatccc                                                                                                                                                                                          >  1:415884/1‑235 (MQ=255)
                                                        ctGCGCCGGGGCTGAAAATTCACGCCAAACTGTTCCTGATTTCACGTAAAGAAAACGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGG‑AACTTTAACGAAAAAACCGCGCGTCTTTATACTGTCTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTG‑AAAACCCATACCGTCCGGTGAC‑ATTTGATTATTTAATGGTATCGCCGCAAAACTccc                                                                                                                                                                                          <  2:38508/249‑1 (MQ=255)
                                                                       aaaTTCACGCCAAACTGTACCTGATTTCACGTAAAGAAAACGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGG‑AACTTTAACGAAAAAACCGCGCTGCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTG‑AAAACCCATACCGTCCGGTGTC‑ATTTGGGTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGT‑ATGaa                                                                                                                                                                        <  1:390997/251‑1 (MQ=255)
                                                                             cagCCAAACTGTTCCTGATTTCACGTAAAAAAAACGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGG‑AACTTTAACGAAAAAACCGCGCGTCTTTATACTGCCTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTG‑AAAACCCATACCGTCCGGTGAC‑ATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGT‑ATGAAAtggtg                                                                                                                                                                  <  2:280556/249‑1 (MQ=255)
                                                                                                  cACGTAAAGAAAACGGTGAAGTGGTGCGTTATGCACACAGCGGGACCGGG‑AACTTTAACGAAAAACCCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTG‑AAAACCCATACCGTCCGGTGAC‑ATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGT‑ATGAAATGGTGGACCGCGAGATCGCCAAcgcg                                                                                                                                             <  2:195720/251‑1 (MQ=255)
                                                                                                    cGTAAAGAAAACGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGG‑AACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTTACTTTATTT‑AAAACCCATACCGTCCGGTGAC‑ATTTTATTATTTTATGGTATCGCCGCAAAACTCCCCCCGCCTATTGT‑ATGAAATGGTGGGCCGCGAGAACGcccacccgcc                                                                                                                                           >  1:309669/1‑243 (MQ=255)
                                                                                                       aaaGAAAACGGTGAAGGGGGGCGTTAGGCACACAGCGGGACCGGG‑AACTTGAACGAAAAAACCGCGCGGCTTTATACGGACTATTCGGTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTG‑AAAACCCATACCGTCCGGTGAC‑ATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGT‑ATGAAATGGTGGACCGCGAGATCGCCAACGCgcagc‑a                                                                                                                                       <  2:200704/251‑1 (MQ=255)
                                                                                                           aaaaCGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGG‑AACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTAAATTTTTTG‑AAAACCCCTACCGGCCCGTGGC‑ATTTGGTTATTTTATGGTATTGCCCCAAAACTACCCCCCCCCCTTGA‑TTGAAATTGGGGAACGCGCGATCGGCCACGCCCAGC‑AAggg                                                                                                                                   >  1:388980/1‑251 (MQ=255)
                                                                                                            aaaCGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGG‑AACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTTACTTTTTTT‑AAAACCCATACCGTCCGGAGAC‑ATTTGGTTATTTTATTGTTTTGCCCCCAAACTCCCGCCCCCTATTTT‑ATTAAAAGGTGGGCCGCGAGATCGcccacaccccgcaagggg                                                                                                                                   >  1:286005/1‑235 (MQ=255)
                                                                                                                      gtggggCGTTATGCACACTTCGGG‑CCGGGCACCTTTACCGAAAAAACCGCGCTTTTTTATACGGACTATCCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCTGGTATTTAACTTTATTG‑AAAACCCATACCGGCTGGTGCC‑ATTTGATTATTTAATGGTATCGCCGCAAACCTCCCGCCGCCTATTGT‑ATGAAATGGTGGACCGCGAGATCGCCAACGCTCAGC‑ATGGGCTGCCCAGTgg                                                                                                                        <  2:227279/251‑1 (MQ=255)
                                                                                                                          tgCGTTATGCACACATCGGGACCGGG‑AACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTTACTTTTTTG‑AAAAACCATACCCTCCGGGGAC‑ATTTGATTATTTTATGGTATCGCCGCAAAAATCCCCCCGCCTATTGT‑ATGAAATGGGGGGCCGCGAGAACGCCAACCCCCAGCCAAGGG‑TGCCCAGTggggtt                                                                                                                    >  1:375737/1‑247 (MQ=255)
                                                                                                                                   cacacaTCGGGACCGGG‑AACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTAACTTTTTTT‑TAAAACCATAACGTCCGGGGAC‑ATTTGTTTATTTTAAGGTATCGCCGCCAAAATCCCCCCGCCCATTGT‑ATGAAAAGGGGGACCGCGAGATCGCCAACGCGCAGC‑AAGGGGTGCCCACGGGGATCAccccga                                                                                                             >  1:318768/1‑246 (MQ=255)
                                                                                                                                                 ggg‑ACCTTTAACGAAAAACCCGCGCGTCTTTATACTGACTGTTCGTTGCGGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTG‑AAAACCCATACCGTCCGGTGAC‑ATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGT‑ATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGC‑AAGGGCTGCCCAGTGGTATCACCCTGAAGCTAAATAACCTTGt                                                                                             <  1:12861/251‑1 (MQ=255)
                                                                                                                                                        tttAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAAATTTATTT‑AAAAACCCTACCGTCCGGGGAC‑AATTGATTTTTTTATGGGATCCCCGCCAAACACCCGCCCCCTCTTTT‑TTGAAAAGGGGGGCCGCGCGATCGCCAACGCGCAGC‑AAGGGGTGCCCCGCGGGATCACCCCTAAGCCAAAAAACCTCTTCGATa                                                                                        >  1:375518/1‑250 (MQ=255)
                                                                                                                                                                           gTCTTTAAACTGACTATTCGTTGCTGCCCGCCGATGCGCGCATCCCCAACGAAGTACGGCGGGTATTTAACTTTATGG‑AAAACCCATACCGTCCGGTGAC‑ATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGT‑ATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGC‑AAGGGCTGCCCAGTGGTATCACCCTGAAGCTAAATAACCGTGTCGATAAA‑GGCCTGGTTGATCGTCTg                                                                   <  2:124017/251‑1 (MQ=255)
                                                                                                                                                                                aaaccttacctttCGTGGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTG‑AAAACCCATACCGTCCGGTGAC‑ATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGT‑ATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGC‑AAGGGCTGCCCAGTGGTATCACCCTGAAGCTAAATAACCTTGTCGATAAA‑GGCCTGGTTGATCGTCTGTATGc                                                              <  2:250367/240‑1 (MQ=255)
                                                                                                                                                                                  taCTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTTACCTTTTTGTAAAACCC‑TACCCGCCCGTGGA‑ATTTTATTTTTTTATGGGGTCGCCGCAAAAAACCCCCCGCCTCTTGT‑ATGAAAAGGGGGAACGCGCGATCGCCACCCCGCCGC‑AAAGGCTGGCCACTGGAGTCCACCTGAAGCCAAATAACCTTGTCGGAAAA‑GGCCCGGTTTATCGGCTGTATGCgg                                                            >  1:89900/1‑251 (MQ=255)
                                                                                                                                                                                             cGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTG‑AAAACCCATACCGTCCGGTGAC‑ATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGT‑ATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGC‑AAGGGCTGCCCAGTGGTATCACCCTGAAGCTAAATAACCTTGGCGATAAA‑GGCCTGGTTGATCGTCTGTATGCGGCCTCCAGCTcc                                                 <  1:121114/251‑1 (MQ=255)
                                                                                                                                                                                                 gCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGTTTTAACTTTTTTG‑GAAAACCATAACGCCCGGGGACCATTT‑TTTATTTTATTGTATCGCCGCAAAAACCCCCCCCCCCATTGTTAT‑AAAAGGTGGGGCGCCAGAACGCCAACGCCCCGC‑CAGGGCTGGCCAGGGGGATCACCCCCGAGCTAAATAAACTTGTTGGTAAA‑AGCCTGGGTGATCGTCTGTTTTTGGGCTCCAGCTCCcgcg                                             >  1:286236/1‑251 (MQ=255)
                                                                                                                                                                                                    aaCCGCCGATGCGCGCATCACCAACGAAGTACGGCGTGTATTTAACTTTATTG‑AAAACCCATACCGTCCGGTGAC‑ATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGT‑ATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGC‑AAGGGCTGCCCAGTGGTAGCACCCTGAAGCTAAATAACCTTGTCGATAAA‑GGCCTGGTTGATCGGCTGTATGCGGCCTCCAGCTCCGGCGt                                            <  2:413591/248‑1 (MQ=255)
                                                                                                                                                                                                     aCCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTTACTTTTTTT‑AAAACCCCTACCGCCCCGGGGC‑CTTTTATTTTTTTATAGGATCGCCCCAAAAATCCCGCCGCCCTTTTT‑TTGAAATGGGGGGCCGCGAGATAGCCCACACGCAGC‑CAGGGCTGCCCCGTGGGGTCAACCCGAAGCTAAATAAACCTGTCGATAAAAGGCC‑GGTTGGTAGTCTGTATGCGGGCTcccgcccccgcggacc                                         >  1:278932/1‑237 (MQ=255)
                                                                                                                                                                                                                    aTCACCAACGAAGTACGGCGGGTATTTAACTTTATTG‑AAAACCCATACCGTCCGGTGAC‑ATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGT‑ATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGC‑AAGGGCTGCCCAGTGGTATCACCCTGAAGCTAAATAACCTTGTCGATAAA‑GGCCTGGTTGATCGTCTGTATGCGGCCTCCAGCTCCGGCGTACCGGTTAATCTGCTGGt                          <  2:132795/251‑1 (MQ=255)
                                                                                                                                                                                                                                            tttAACTTTATTG‑AACACCCATACCGTCCGGTGAC‑ATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGT‑ATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGC‑AAGGGCTGCCCAGTGGTATCACCCTGAAGCTAAATAACCTTGTCGATAAA‑GGCCTGGTTGATCGTCTGTATGCGGCCTCCAGCTCCGGCGTACCGGTTAATCTGCTGGTTCGCGGAATGTGTTCGCTGATccc  >  1:313390/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      
AGCCAACATTCACTGGGCGAAGCGCCTGACCGAAGCAGGCGTGCACGTTATCTTCTCTGCGCCGGGGCTGAAAATTCACGCCAAACTGTTCCTGATTTCACGTAAAGAAAACGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGG‑AACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTG‑AAAACCCATACCGTCCGGTGAC‑ATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGT‑ATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGC‑AAGGGCTGCCCAGTGGTATCACCCTGAAGCTAAATAACCTTGTCGATAAA‑GGCCTGGTTGATCGTCTGTATGCGGCCTCCAGCTCCGGCGTACCGGTTAATCTGCTGGTTCGCGGAATGTGTTCGCTGATCCC  >  NC_000913/2624270‑2624755

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