New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 442659134 (1.130)6 (0.230) 6/196 2.2 15.4% intergenic (+34/+37) rplI/yjfZ 50S ribosomal subunit protein L9/uncharacterized protein
?NC_000913 = 4426592 36 (1.360)intergenic (+35/+36) rplI/yjfZ 50S ribosomal subunit protein L9/uncharacterized protein
Rejected: Frequency below/above cutoff threshold.

GCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4426484‑4426591
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aaaacgccgaccaATGGTCGGCGTTTTTACGTCTCGTTGAATAACGAATTATTCAGCTACTACGTTTACGATCACTTTCGCGAATACTTCGCTGTGAACCTGGAAGCTCACTTCGTGTTCGCCA  <  NC_000913/4426592‑4426482
                                                                                                                                                                                                                           
GCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCAATGGTCGGCGTTT                                                                                                    >  5:89055/1‑121
GCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCAATGGTCGGCGTTT                                                                                                    <  6:89055/121‑1
     CACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCAATGGTCGGCGTTTTTACGTCTCGTTGAATAACGAATTATTCAGCTA                                                                   >  1:16502/1‑149
              AGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCAATGGTCGGCGTTTTTACGTCTCGTTGAATAACGAATTATTCAGCTACTACGTTTA                                                          <  2:16502/149‑1
                                                                          TTCGTTATTCAACGAGACGTAAAAACGCCGACCAATGGTCGGCGTTTTTACGTCTCGTTGAATAACGAATTATTCAGCTACTACGTTTACGAT                                                      >  1:57564/1‑93
                                                                          TTCGTTATTCAACGAGACGTAAAAACGCCGACCAATGGTCGGCGTTTTTACGTCTCGTTGAATAACGAATTATTCAGCTACTACGTTTACGAT                                                      <  2:57564/93‑1
                                                                                                       GACCAATGGTCGGCGTTTTTACGTCTCGTTGAATAACGAATTATTAAGCTACTACGTTTACGATCACTTTCGCGAATACTTCGCTGTGAACCTGGAAGCTCACTTCGTGTTCGCCA  >  6:171075/1‑116
                                                                                                                                                                                                                           
GCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4426484‑4426591
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aaaacgccgaccaATGGTCGGCGTTTTTACGTCTCGTTGAATAACGAATTATTCAGCTACTACGTTTACGATCACTTTCGCGAATACTTCGCTGTGAACCTGGAAGCTCACTTCGTGTTCGCCA  <  NC_000913/4426592‑4426482

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.