Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0009134,166,5790A→C100.0% 1.0 4intergenic (+294/‑80)murI/rrsBglutamate racemase/16S ribosomal RNA of rrnB operon
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (3/0):  ref base (0/1):  total (3/1)

ATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAGGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCGAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGC  >  NC_000913/4166295‑4166872
                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                     
aTGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCACAAATAAATGCTTGACTCTGTCGCGGGAAGGCGTATTATGCACACCCCGCGCCGCTGAGACAAAGCTAAGCGGCACTGCTCTTTAACAATTTATCCGACAATCTGTGTTGGCACTCGAAGATACGGATTATTAACGTCGACAGACGAAAAATGCATa                                                                                                                                                                                                                                                                                       >  1:258485/1‑301 (MQ=255)
aTGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAGGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCGAAGCGGCCCTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCCAGACGAAAAATGAATa                                                                                                                                                                                                                                                                                       >  1:98641/1‑301 (MQ=14)
aTGTTGCGCGGTCAGAAAATTATGTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCATGAGAACTACGGCAGAGGAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAGGCGTATTATGCAGAACACGCGCCGATGAGAAAAAGCGAAGCGGAACTGCTCTTAAACAATTGATCAGACACTCAGATTGGGCACTACAAGATACGGACTCTTAAGGTCGACAGACGAAAAATGAATa                                                                                                                                                                                                                                                                                       >  1:338790/1‑301 (MQ=11)
                                                                                                                                                                                                                                                                                     aaggcGCAAGAAGAAAAATGAATAACAAGTCTCAAGAGTGAACACGTAATTAATTAAGAAGTTTAATTCTGTGAGGGTCAAACGTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGc  <  1:632632/297‑1 (MQ=2)
                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                     
ATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAGGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCGAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGC  >  NC_000913/4166295‑4166872

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.