New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? CP009273 = 3935260NA (NA)10 (0.020) 6/106 NT NA noncoding (93/1542 nt) rrsC 16S ribosomal RNA of rrnC operon
?CP009273 = 4028965 NA (NA)noncoding (75/1542 nt) rrsA 16S ribosomal RNA of rrnA operon

CCTAACACATGCAAGTCGAACGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGG  >  CP009273/3935214‑3935309
                                              |                                                 
ccaaaCACATGCAAGTCGAACGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACggg                               <  1:37317708/64‑1 (MQ=1)
                              aatcAGCTTGCTGTTTCGCTGACGAGTGTCGTACGGGTGAGTAATGTCTGGGAAACTGCCTGATgg  <  1:9475722/63‑1 (MQ=1)
                              aaccAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGAt    <  1:4967493/61‑1 (MQ=1)
                                              |                                                 
CCTAACACATGCAAGTCGAACGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGG  >  CP009273/3935214‑3935309

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