New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 3941882NA (NA)8 (0.130) 5/280 8.2 31.1% noncoding (76/1542 nt) rrsC 16S ribosomal RNA of rrnC operon
?NC_000913 = 4208237 18 (0.290)noncoding (92/1542 nt) rrsE 16S ribosomal RNA of rrnE operon
Rejected: Coverage evenness skew score above cutoff.

GAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCT  >  NC_000913/3941792‑3942022
                                                                                          |                                                                                                                                            
gAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGCAGCTTGCTGCTTCGCTGACGAGTGGCGGACGCGCGCCGAATGTCTgcg                                                                                          >  1:317324/1‑141 (MQ=11)
                             gagcATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAactact                                                             <  1:233127/140‑1 (MQ=1)
                                            tGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGCCAGCTACGGGCGACGGTGGcgga                                              >  2:492873/1‑140 (MQ=1)
                                                                         aagttgaggggTAACAGGAAACAGCTGGCTGTTGCGCGGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAgg                 <  1:398740/134‑1 (MQ=1)
                                                                              caaaggaaaCAGTAAACAGCTGGTTGTTTCGCTGACGAGTGGTGGACGGGTGTGTAATGTCTGGGAAACTTCCTGATGTAGGGGGATAACTACTGGAAACGGGAGCTACCACCGCATAACGTCGCAAGACCAAAGAGGGGGAc            <  1:934510/136‑1 (MQ=1)
                                                                                   ggaaCATGACACAGCTTGCTGTTGCGCTGACGAGTGGCGGACAGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTAGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCg       <  1:294051/141‑1 (MQ=11)
                                                                                   gTAACAGTAGACGGCCTGCGGTTGCGCTGACGAGTGGCGGACGGGTGAGGAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCg       <  2:799/143‑1 (MQ=11)
                                                                                    ttaCATGACACCGCTTGCTGTCTCGATGACGAGTCGCGGAGGGGTGAGTTATGTCGGGGAAACTGCCTGATGTAGGGGGATAACTACTGGAAACGGGAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCgg      <  1:820486/141‑1 (MQ=11)
                                                                                       agcgAAACAGCAAGCTGTTTC‑CTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGAATGTAAATGCTAGGACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCct  >  2:797388/4‑143 (MQ=11)
                                                                                          |                                                                                                                                            
GAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCT  >  NC_000913/3941792‑3942022

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.