New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 5224500 =5 (0.070)10 (0.220)
+CTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCA
9/178 4.9 34.3% intergenic (‑204/+175) PP_4603/PP_4604 transposase/DMT superfamily permease
?NC_002947 5545626 = 54 (0.770)intergenic (‑766/+111) rpsF/rlmB 30S ribosomal protein S6/23S rRNA (guanosine(2251)‑2'‑O)‑methyltransferase RlmB
Rejected: Coverage evenness skew score above cutoff.

AACAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGCGGTTTACC  >  NC_002947/5224492‑5224644
        |                                                                                                                                                
aaCAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGCGGtt      <  5:378861/149‑1 (MQ=32)
   aaaGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGAACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTCGGTGATCTCCTACAGGTCGGTTTAc   >  2:31315/1‑149 (MQ=16)
    aaagaaaaGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCCGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGCGGTTTAcc  >  8:366070/5‑149 (MQ=11)
    aaGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGCGGTTTAcc  >  2:112401/1‑149 (MQ=35)
    aaGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGCGGTTTAcc  <  7:203321/149‑1 (MQ=35)
    aaGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGGCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGCGGTTTAcc  <  2:282635/149‑1 (MQ=25)
        |                                                                                                                                                
AACAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGCGGTTTACC  >  NC_002947/5224492‑5224644

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.