New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4477260 =50 (1.380)4 (0.120) 4/220 2.8 8.5% intergenic (+9/+47) yjgL/argI SopA‑central‑domain‑like hexapeptide repeat protein/ornithine carbamoyltransferase 1
?NC_000913 4477291 = 40 (1.200)intergenic (+40/+16) yjgL/argI SopA‑central‑domain‑like hexapeptide repeat protein/ornithine carbamoyltransferase 1
Rejected: Frequency below/above cutoff threshold.

CTATAAAGATGCCGGGAAATACAAACCATCCCGGCATC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4477297‑4477260
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tcccggcatcTTTATAGCGATAGCAATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTA  >  NC_000913/4477291‑4477428
                                                                                                                                                                                
CTATAAAGATGCCGGGAAATACAAACCATCCCGGCATCTTTATAGCGATAGCAATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGAT                                       >  1:5137/1‑139
CTATAAAGATGCCGGGAAATACAAACCATCCCGGCATCTTTATAGCGATAGCAATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGAT                                       <  2:5137/139‑1
 tcTAAAGAGGCCGGGAAATACACACCCTCCCGGCATCTTTATAGCGAGAGCAAGTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGGATGCATACGGTTTTCCGCCTGATCAAACACAA                                                     <  4:113079/122‑1
 TATAAAGATGCCGGGAAATACAAACCATCCCGGCATCTTTATAGCGATAGCAATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGGTTTCCGCCTGATCAAAAACAA                                                     >  3:113079/1‑124
                            TCCCGGCATCTTTATAGCGATAGCAATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTA  <  1:758/148‑1
                                 GCATCTTTATAGCGATAGCAATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATG    <  5:66787/141‑1
                                 GCATCTTTATAGCGATAGCAATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATG    >  6:66787/1‑141
                                                                                                                                                                                
CTATAAAGATGCCGGGAAATACAAACCATCCCGGCATC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4477297‑4477260
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tcccggcatcTTTATAGCGATAGCAATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGTGACTTCCATACCGCCATGTA  >  NC_000913/4477291‑4477428

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.