New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4638146	87 (0.870)	6 (0.070)	5/156	NT	6.7%	intergenic (+26/‑32)	creC/creD	sensory histidine kinase in two‑component regulatory system with CreB or PhoB/inner membrane protein
	?	NC_000913	= 4638170	92 (1.050)					intergenic (+50/‑8)	creC/creD	sensory histidine kinase in two‑component regulatory system with CreB or PhoB/inner membrane protein
Rejected: Frequency below cutoff threshold.

TTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4638037‑4638146 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑acatagtcttcTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGACGGTGAAGTCGAAGCGAGGCCAGCACGCCACCTTCC  <  NC_000913/4638170‑4638072                                                                                                                                                                                                                   TTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCA                                                                                                               <  1:714987/100‑1  TGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCAC                                                                                                              <  2:1643506/100‑1   GTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCAATTCACATAGCTTCAAATTCTTCCCACA                                                                                                             >  2:2394828/1‑100      TAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAG                                                                                                          <  2:465030/100‑1      TAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAG                                                                                                          <  2:1762300/100‑1        ACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTC                                                                                                        >  2:2409881/1‑100         CGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCT                                                                                                       <  2:2188206/100‑1          GGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTT                                                                                                      >  1:2020343/1‑100                                                GCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCTCCTTTGCAATGGCAGCAGGATATGAAGACTATGTGG                                                                <  1:1870156/100‑1                                                GCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGG                                                                >  1:388922/1‑100                                                GCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCGTCAAATTCTTCCCACATAGGCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGG                                                                <  2:1499279/100‑1                                                 CTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGG                                                               <  2:42467/100‑1                                                  TGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGA                                                              <  2:388921/100‑1                                                   GGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAA                                                             <  1:79271/100‑1                                                                                                    ACATAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGACGGTGAAGTCGAAGCGAGGCCAGCAC            <  1:2608266/100‑1                                                                                                    ACATAGTCTTCTCCGTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGACGGTGAAGTCGAAGCGAGGCCAGCA             <  2:1790776/99‑1                                                                                                       TAGTCTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGACGGTGAAGTCGAAGCGAGGCCAGCACGCC         <  2:1539576/100‑1                                                                                                           CTTCTCCTTTGCAATGGCAGCAGGATACGAAGACTATGCGGGAAGAATTTGAAGCTATGTGAAGTGACGGTGAAGTCGAAGCGAGGCCAGCACGCCACCT     >  2:1334116/1‑100                                                                                                              CTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGACGGTGAAGTCGAAGCGAGGCCAGCACGCCACCTTCC  <  2:501347/100‑1                                                                                                              CTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGACGGTGAAGTCGAAGCGAGGCCAGCACGCCACCTTCC  <  2:1705136/100‑1                                                                                                                                                                                                                   TTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4638037‑4638146 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑acatagtcttcTCCTTTGCAATGGCAGCAGGATACGAAGACTATGTGGGAAGAATTTGAAGCTATGTGAAGTGACGGTGAAGTCGAAGCGAGGCCAGCACGCCACCTTCC  <  NC_000913/4638170‑4638072

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