New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 3941883NA (NA)7 (0.120) 5/280 8.0 29.5% noncoding (76/1542 nt) rrsC 16S ribosomal RNA of rrnC operon
?NC_000913 = 4208238 17 (0.280)noncoding (92/1542 nt) rrsE 16S ribosomal RNA of rrnE operon
Rejected: Coverage evenness skew score above cutoff.

GATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTG‑ATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCT  >  NC_000913/3941822‑3942023
                                                             |                                                                                                                                             
gagcATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTG‑ATGGAGGGGGATAactact                                                             <  2:233127/140‑1 (MQ=1)
                                                 caaaggaaaCAGTAAACAGCTGGTTGTTTCGCTGACGAGTGGTGGACGGGTGTGTAATGTCTGGGAAACTTCCTG‑ATGTAGGGGGATAACTACTGGAAACGGGAGCTACCACCGCATAACGTCGCAAGACCAAAGAGGGGGAc            <  2:934510/136‑1 (MQ=1)
                                                      ggaaCATGACACAGCTTGCTGTTGCGCTGACGAGTGGCGGACAGGTGAGTAATGTCTGGGAAACTGCCTG‑ATGGAGGGGGATAACTACTAGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCg       <  2:294051/141‑1 (MQ=11)
                                                      gTAACAGTAGACGGCCTGCGGTTGCGCTGACGAGTGGCGGACGGGTGAGGAATGTCTGGGAAACTGCCTG‑ATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCg       <  1:799/143‑1 (MQ=11)
                                                       ttaCATGACACCGCTTGCTGTCTCGATGACGAGTCGCGGAGGGGTGAGTTATGTCGGGGAAACTGCCTG‑ATGTAGGGGGATAACTACTGGAAACGGGAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCgg      <  2:820486/141‑1 (MQ=11)
                                                         aCAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGAGTGAGTAATCTGCAGGAAA‑TGCCCGCATGGCACGGGAAAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACATGCGGGc    >  2:486053/1‑143 (MQ=11)
                                                          agcgAAACAGCAAGCTGTTTC‑CTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTG‑ATGGAGGGGGATAACTACTGAATGTAAATGCTAGGACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCct  >  1:797388/4‑143 (MQ=11)
                                                             |                                                                                                                                             
GATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTG‑ATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCT  >  NC_000913/3941822‑3942023

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