Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0009131,528,8820G→A100.0% 1.0 2pseudogene (981/2037 nt)rhsEpseudogene, Rhs family
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (1/0):  ref base (1/0):  total (2/0)

TATCCAGACGGTATACGAGCCGGGAAGCTTCACGCCGCTCATCCGGGTCGAGACAGAGAACGGCGAGCGGGAAAAAGCGCAGCGGCGCAGCCTGGCAGAGACGCTCCAGCAGGAAGGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTT  >  NC_000913/1528876‑1529176
      |                                                                                                                                                                                                                                                                                                      
tATCCAGACGGTATACGAGCCGGGAAGCTTCACGCCGCTCATCCGGGTCGAGACAGAGAACGGCGAGCGGGAAAAAGCGCAGCGGCGCAGCCTGGCAGAGACGCTCCAGCAGGAAGGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGCCCGGTGGCGCCGGACTCCACCCCGGCCCGACAAGtg  >  1:399342/1‑300 (MQ=2)
tATCCAAACGGTATACGAGCCGGGAAGCTTCACGCCGCTCATCCGGGTCGAGACAGAGAACGGCGAGCGGGAAAAAGCGCAGCGGCGCAGCCTGGCAGAGACGCTCCAGCAGGATGGGAGTGCGAACGGCCCCGGAGTGGTGTTCCCTGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGCAAGCCGGGGGTGGCTTGCGCAGTGCGGGCTGAGGGTGTAGCCACTGGCCCGACCGGTGGCGCCGGAAAAACGACCGGCGCGACAAGtt  >  2:117448/1‑301 (MQ=11)
      |                                                                                                                                                                                                                                                                                                      
TATCCAGACGGTATACGAGCCGGGAAGCTTCACGCCGCTCATCCGGGTCGAGACAGAGAACGGCGAGCGGGAAAAAGCGCAGCGGCGCAGCCTGGCAGAGACGCTCCAGCAGGAAGGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTT  >  NC_000913/1528876‑1529176

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 22 ≤ 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.