Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913807,3430CG56.2% ‑3.0 / 22.9 16intergenic (‑62/+90)ybhC/ybhBouter membrane lipoprotein YbhC/putative kinase inhibitor
Reads supporting (aligned to +/- strand):  ref base C (0/7);  new base G (6/3);  total (6/10)
Fisher's exact test for biased strand distribution p-value = 1.14e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.17e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGATGGCATTGCGCCATCTGGCAGAGTGATTAACTAAACATCGCAGTAATCGAGGCGCTTGCCAGAGAGTGGAAATGAACGTTAAACCCGACCATCGCGCCGCTGGCACCTTCATCGACATCAATACGTTCTATATCCAGCGCGTGAACGGTAAAAATGTAGCGATGAGTTTCGCCTTTCGGCGGTGCTGCGCCATCGTACCCGGTTTTACCAAAGTCGGTACGCGTCTGCAA  >  NC_000913/807050‑807636
                                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                                     
catgcaggGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAAGTTggca                                                                                                                                                                                                                                                                                                 <  2:831844/297‑4 (MQ=255)
  ggcaggGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAAGTTggcatt                                                                                                                                                                                                                                                                                               >  2:487613/1‑295 (MQ=255)
        gTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCGGTTGAAGCCTGCTTTTTTATACTAAGTTggcattataaaa                                                                                                                                                                                                                                                                                         >  2:721760/1‑289 (MQ=255)
          aaTCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGTGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGCCGTTGATTGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGCAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATCCTAAGTTggaattataaaaaa                                                                                                                                                                                                                                                                                       >  2:329655/1‑287 (MQ=255)
               gggAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCAAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAAGTTggcattataaaaaagcatt                                                                                                                                                                                                                                                                                  >  1:328535/1‑282 (MQ=255)
                 gAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAAGTTggcattataaaaaagcattgc                                                                                                                                                                                                                                                                                <  1:308835/300‑21 (MQ=255)
                    ggATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACNGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAAGTTggcattataaaaaagcattgctta                                                                                                                                                                                                                                                                             >  1:597414/1‑277 (MQ=255)
                      aTTCCACGCTGCAGCACATGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAAGTTggcattataaaaaagcattgcttatc                                                                                                                                                                                                                                                                           >  2:264474/1‑275 (MQ=255)
                        tCCACGCTGCAACACCTGGTGTCAGGGGTGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAATGACGTTGATCGGGCGGGGTTGAGCGACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAAGTTggcattataaaaaagcattgcttatcaa                                                                                                                                                                                                                                                                         <  2:377535/300‑28 (MQ=255)
                                                                                                                                                                                                                                                            ttatatcattttacgtttctcgttcagCTTTTTTATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGATGGCATTGCGCCATCTGGCAGAGTGATTAACTAAACATCGCAGTAATCGAGGCGCTTGCCAGAGAGTGGAAATGAACGTTAAACCCGACCATCGCGCCGCTGGCACCTTCATCGACATCAATACGTTCTATATCCAGCGCGTGAACGGTAAAAATGTAGCGATGAGTTTCGCCTTTCGTCGGTGCTGCGCCATCg                                     <  1:551969/274‑1 (MQ=255)
                                                                                                                                                                                                                                                                 tcattttacgtttctcgttcagCTTTTTTATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGATGGCATTGCGCCATCTGGCAGAGTGATTAACTAAACATCGCAGTAATCGAGGCGCTTGCCAGAGAGTGGAAATGAACGTTAAACCCGACCATCGCGCCGCTGGCACCTTCATCGACATCAATACGTTCTATATCCAGCGCGTGAACGGTAAAAATGTAGCGATGAGTTTCGCCTTTCGGCGGTGCTGCGCCATCGTAccc                                <  1:517240/279‑1 (MQ=255)
                                                                                                                                                                                                                                                                    ttttacgtttctcgttcagCTTTTTTATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGATGGCATTGCGCCATCTGGCAGAGTGATTAACTAAACATCGCAGTAATCGAGGCGCTTGCCAGAGAGTGGAAATGAACGTTAAACCCGACCATCGCGCCGCTGGCACCTTCATCGACATCAATACGTTCTATATCCAGCGCGTGAACGGTAAAAATGTAGCGATGAGTTTCGCCTTTCGGCGGTGCTGCGCCATCGTACCCGGt                             <  1:315429/282‑1 (MQ=255)
                                                                                                                                                                                                                                                                          gtttctcgttcagCTTTTTTATACTAACTTGAGCTAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGATGGCATTGCGCCATCTGGCAGAGTGATTAACTAAACATCGCAGTAATCGAGGCGCTTGCCAGAGAGTGGAAATGAACGTTAAACCCGACCATCGCGCCGCTGGCACCTTCATCGACATCAATACGTTCTATATCCAGCGCGTGAACGGTAAAAATGTAGCGATGAGTTTCGCCTTTCGGCGGTGCTGCGCCATCGTACCCGGTTTTAcc                       <  2:287847/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                                      gCTTTTTTATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGATGGCATTGCGCCATCTGGCAGAGTGATTAACTAAACATCGCAGTAATCGAGGCGCTTGCCAGAGAGTGGAAATGAACGTTAAACCCGACCATCGCGCCGCTGGCACCTTCATCGACATCAATACGTTCTATATCCAGCGCGTGAACGGTAAAAATGTAGCGATGAGTTTCGCCTTTCGGCGGTGCTGCGCCATCGTACCCGGTTTTACCAAAGTCGGTAcg           <  1:793115/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                        ttttttATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAATTGATACCAGATGGCATTGCGCCATCTGGCAGAGTGATTAACTAAACATCGCAGTAATCGAGGCGCTTGCCAGAGAGTGGAAATGAACGTTAAACCCGACCATCGCGCCGCTGGCACCTTCATCGACATCAATACGTTCTATATCCAGCGCGTGAACGGTAAAAATGTAGCGATGAGTTTCGCCTTTCGGCGGTGCTGCGCCATCGTACCCGGTTTTACCAAAGTCGGTAcgcg         <  2:142209/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                               taCTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGATGGCATTGCGCCATCTGGCAGAGTGATTAACTAAACATCGCAGTAATCGAGGCGCTTGCCAGAGAGTGGAAATGAACGTTAAACCCGACCATCGCGCCGCTGGCACCTTCATCGACATCAATACGTTCTATATCCAGCGCGTGAACGGTAAAAATGTAGCGATGAGTTTCGCCTTTCGGCGGTGCTGCGCCATCGTACCCGGTTTTACCAAAGTCGGTACGCGTCTGCaa  <  1:193469/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                                     
CAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAGGACGTTGATCGGGCGGGGTTGAGCTACAGGCGGTCAGCGTCACGCCAAAAGCCAATGCCAGCGCCAGACGGGAAACTGAAAATGTGTTCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGATGGCATTGCGCCATCTGGCAGAGTGATTAACTAAACATCGCAGTAATCGAGGCGCTTGCCAGAGAGTGGAAATGAACGTTAAACCCGACCATCGCGCCGCTGGCACCTTCATCGACATCAATACGTTCTATATCCAGCGCGTGAACGGTAAAAATGTAGCGATGAGTTTCGCCTTTCGGCGGTGCTGCGCCATCGTACCCGGTTTTACCAAAGTCGGTACGCGTCTGCAA  >  NC_000913/807050‑807636

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