New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4053426 =184 (1.420)4 (0.050)
+CAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATT
4/180 7.9 3.5% intergenic (+8/+221) hemN/yshB coproporphyrinogen III oxidase, SAM and NAD(P)H dependent, oxygen‑independent/uncharacterized protein
?NC_000913 = 4103596 NA (NA)intergenic (+100/+6) yiiM/cpxA 6‑N‑hydroxylaminopurine resistance protein/sensory histidine kinase in two‑component regulatory system with CpxR
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GCACATTGATGGCGCTTTCGTTTTCTACTTTGTAAACGAAGCGCCAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4053472‑4053426
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GCATTTGCAGGCAGGATAAGGCGTTTACGCCGCATCCGGCATTTGAGCAGGATGC  <  NC_000913/4103596‑4103542
                                               |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||                                                       
GCACATTGATGGCGCTTTCGTTTTCTACTTTGTAAACGAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGAAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATct             >  2:265125/1‑148
        ATGGCGCTTTCGTTTTCTACTTTGTAAACGAAGCGCCATCAATGTGCTGAGATCAAGATTGCCTGATGCGCTGGGCAACAAAAACTGTTGAAATTATTGCTATTGTAGGACGGATAAGGCGTTTTCGCCGCATCAGGCATCTGAGCAtta     >  1:258303/1‑147
          GGCGCTTTCGTTTTCTACTTTGTAAACGAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATG   >  1:75097/1‑150
           GCGCTTTCGTTTTCTACTTTGTAAACGAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGC  >  1:725853/1‑150
                                               |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||                                                       
GCACATTGATGGCGCTTTCGTTTTCTACTTTGTAAACGAAGCGCCAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4053472‑4053426
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GCATTTGCAGGCAGGATAAGGCGTTTACGCCGCATCCGGCATTTGAGCAGGATGC  <  NC_000913/4103596‑4103542

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