New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4624213 =138 (0.900)5 (0.030) 4/168 NT 3.9% coding (577/645 nt) ytjB SMP_2 family putative membrane‑anchored periplasmic protein
?NC_000913 4624245 = 119 (0.820)coding (545/645 nt) ytjB SMP_2 family putative membrane‑anchored periplasmic protein
Rejected: Frequency below cutoff threshold.

GATAACACCACTAACATTTTACGCCTGATGTTGCTGCTCTCACTGGCAATCGGTGTAGTGCTGACCCGCACGCTGCTACAGGGTAAACGCACCCGCTGGCAGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4624315‑4624213
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcagcGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGT  >  NC_000913/4624245‑4624342
                                                                                                                                                                                                         
GATAACACCACTAACATTTTACGCCTGATGTTGCTGCTCTCACTGGCAATCGGTGTAGTGCTGACCCGCACGCTGCTACAGGGTAAACGCACCCGCTGGC                                                                                                       >  1:1101659/1‑100
GATAACACCACTAACATTTTACGCCTGATGTTGCTGCTCTCACTGGCAATCGGTGTAGTGCTGACCCGCACGCTGCTACAGGGTAAACGCACCCGCTGGC                                                                                                       <  1:559128/100‑1
  TAACACCACTAACATTTTACGCCTGATGTTGCTGCTCTCACTGGCAATCGGTGTAGTGCTGACCCGCACGCTGCTACAGGGTAAACGCACCCGCTGGCAG                                                                                                     >  2:884197/1‑100
  TAACACCACTAACATTTTACGCCTGATGTTGCTGCTCTCACTGGCAATCGGTGTAGTGCTGACCCGCACGCTGCTACAGGGTAAACGCACCCGCTGGCAG                                                                                                     >  2:2075372/1‑100
  TAACACCACTAACATTTTACGCCTGATGTTGCTGCTCTCACTGGCAATCGGTGTAGTGCTGACCCGCACGCTGCTACAGGGTAAACGCACCCGCTGGCAG                                                                                                     <  2:1757474/100‑1
    ACACCACTAACATTTTACGCCTGATGTTGCTGCTCTCACTGGCAATCGGTGTAGTGCTGACCCGCACGCTGCTACAGGGTAAACGCACCCGCTGGCAGC                                                                                                    <  1:4246136/99‑1
                                                                  CGCACGCTGCTACAGGGTAAACGCACCCGCTGGCAGCGTGCGGGTCAGCACTACACCGCTTGCCAGTGAGTGCAGCAACATCAGGCGTAAAATGTTAGTG                                     >  1:3946215/1‑100
                                                                  CGCACGCTGCTACAGGGTAAACGCACCCGCTGGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGTAGCAACATCAGGCGTAAAATGTTAGTG                                     >  1:2033557/1‑100
                                                                   GCACGCTGCTACAGGGTAAACGCACCCGCTGGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGG                                    >  1:1377094/1‑100
                                                                    CACGCTGCTACAGGGTAAACGCACCCGCTGGCAGCGTGCGGGTCAGCACAACACCGATTGACAGTGAGAGCAGCAACATCAGGCGCAAAATGTTAGTGG                                    >  1:1077942/1‑99
                                                                        CTGCTACAGGGTAAACGCACACGCTGGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTA                               >  1:230338/1‑100
                                                                                                  GCAGTGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAG     >  2:461234/1‑100
                                                                                                  GCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAG     >  2:179091/1‑100
                                                                                                  GCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAG     <  2:3038487/100‑1
                                                                                                   CAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGC    >  1:3937999/1‑100
                                                                                                   CAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGC    <  2:3119333/100‑1
                                                                                                    AGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCG   <  2:3553339/100‑1
                                                                                                     GCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGT  >  2:2304526/1‑100
                                                                                                     GCGGGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGT  <  1:1004737/100‑1
                                                                                                                                                                                                         
GATAACACCACTAACATTTTACGCCTGATGTTGCTGCTCTCACTGGCAATCGGTGTAGTGCTGACCCGCACGCTGCTACAGGGTAAACGCACCCGCTGGCAGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4624315‑4624213
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcagcGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGT  >  NC_000913/4624245‑4624342

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