New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 93900243 (0.710)3 (0.050) 3/210 5.7 7.4% coding (1009/1344 nt) rarA recombination intermediate processing DNA‑dependent ATPase
?NC_000913 = 939036 32 (0.530)coding (1043/1344 nt) rarA recombination intermediate processing DNA‑dependent ATPase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GCGCGTCGCTGTCTGGCGATTGCGTCTGAAGACGTCGGTAATGCCGATCCACGGGCGATGCAGGTGGCAATTGCGGCCTGGGATTGCTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/938867‑939002
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aGCGTTGCTTTTTGGCGCGCAGGCCAGGTAAACA  <  NC_000913/939036‑939004
                                                                                                                                                                         
GCGCGTCGCTGTCTTGCGATTGCGTCTGAAGACGTCGGTAATGCCGATCCACGGGCGATGCCGGTGGCAATTGCGGCGTGGGTTTGCTTTACTCGCGTTGGCCCGGCGGTAGGTGAACACGCCATTGCTCATGCGAGtg                                >  1:464588/1‑137
                                                                TGGCAATTGCGGCCTGGGATTGCTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGAGCGTTGCTTTTTGGCGCGCAGGCCAGGTAAACA  <  1:877486/105‑1
                                                                TGGCAATTGCGGCCTGGGATTGCTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGAGCGTTGCTTTTTGGCGCGCAGGCCAGGTAAACA  >  2:877486/1‑105
                                                                                                                                                                         
GCGCGTCGCTGTCTGGCGATTGCGTCTGAAGACGTCGGTAATGCCGATCCACGGGCGATGCAGGTGGCAATTGCGGCCTGGGATTGCTTTACTCGCGTTGGCCCGGCGGAAGGTGAACGCGCCATTGCTCAGGCGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/938867‑939002
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aGCGTTGCTTTTTGGCGCGCAGGCCAGGTAAACA  <  NC_000913/939036‑939004

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 40 ≤ 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.