New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 5224500 =1 (0.020)8 (0.190)
+CTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCA
8/174 4.8 29.3% intergenic (‑204/+175) PP_4603/PP_4604 transposase/DMT superfamily permease
?NC_002947 5545626 = 59 (0.890)intergenic (‑766/+111) rpsF/rlmB 30S ribosomal protein S6/23S rRNA (guanosine(2251)‑2'‑O)‑methyltransferase RlmB
Rejected: Coverage evenness skew score above cutoff.

ACAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGCGGT  >  NC_002947/5224493‑5224639
       |                                                                                                                                           
aCAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGCGGt  <  1:368718/147‑1 (MQ=32)
       |                                                                                                                                           
ACAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGCGGT  >  NC_002947/5224493‑5224639

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