New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1179565 =	15 (0.480)	4 (0.130)	4/254	3.5	21.4%	intergenic (+61/‑66)	nagK/cobB	N‑acetyl‑D‑glucosamine kinase/deacetylase of acs and cheY, chemotaxis regulator
	?	NC_000913	1179588 =	15 (0.510)					intergenic (+84/‑43)	nagK/cobB	N‑acetyl‑D‑glucosamine kinase/deacetylase of acs and cheY, chemotaxis regulator
Rejected: Coverage evenness skew score above cutoff.

ATCTAACCGATTAAACAACAGAGGTTGCTATGCTGTCGCGTCGGGGTCATCGGTTAAGTCGTTTTCGTAAAAATAAACGCCGCCTGCGCGAGCGTTTGCGTCAGCGTATTTTTTTCAGAGATAAAGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGT  >  NC_000913/1179491‑1179718                                                                                                  |                                                                                                                                   aTCTAACCGATTAAACAACAGAGGTTGCTATGCTGTCGCGTCGGGGTCATCGGTTAAGTCGTTTTCGTAAAAATAAACGCCGCCTGCGCGAGCGTTTGCGTCAGCGTATTTTTTTCAGAGATAAAGTGGTGCCGGAAGCAATGGaaaaa                                                                                 >  1:119386/1‑149 (MQ=255)                    agagGTTGCTATGCTGTCGCGTCGGGGTCATCGGTTAAGTCGTTTTCGTAAAAATAAACGCCGCCTGCGCGAGCGTTTGCGTCAGCGTATTTTTTTCAGAGATAAAGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGa                                                              <  2:119386/149‑1 (MQ=255)                             tATGCTGTCGCGTCGGGGTCATCGGTTAAGTCGTTTTCGTAAAAATAAACGCCGCCTGCGCGAGCGTTTGCGTCAGCGTATTTTTTTCAGAGATAAAGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAg                                                     >  6:35229/1‑149 (MQ=255)                                              gTCATCGGTTAAGTCGTTTTCGTAAAAATAAACGCCGCCTGCGCGAGCGTTTGCGTCAGCGTATTTTTTTCAGAGATAAAGTGGTGCCGGAAGCAATgg                                                                                      <  5:124295/99‑1 (MQ=255)                                              gTCATCGGTTAAGTCGTTTTCGTAAAAATAAACGCCGCCTGCGCGAGCGTTTGCGTCAGCGTATTTTTTTCAGAGATAAAGTGGTGCCGGAAGCAATgg                                                                                      >  6:124295/1‑99 (MQ=255)                                                     gTTAAGTCGTTTTCGTAAAAATAAACGCCGCCTGCGCGAGCGTTTGCGTCAGCGTATTTTTTTCAGAGATAAAGTGGTGCAGGAAGCAATGGAAAAACCAAGAGTACTCTTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTc                             >  1:136486/1‑149 (MQ=255)                                                          gTCGTTTTCGTAAAAATAAACGCCGCCTGCGCGAGCGTTTGCGTCAGCGTATTTTTTTCAGAGATAAAGTGGTGCCGGAAGCAATGGAAAAACCaa                                                                             <  7:75684/96‑1 (MQ=255)                                                          gTCGTTTTCGTAAAAATAAACGCCGCCTGCGCGAGCGTTTGCGTCAGCGTATTTTTTTCAGAGATAAAGTGGTGCCGGAAGCAATGGAAAAACCaa                                                                             >  8:75684/1‑96 (MQ=255)                                                           tCGTTTTCGTAAAAATAAACGCCGCCTGCGCGAGCGTTTGCGTCAGCGTATTTTTTTCAGAGATAAAGTGGTGCCGGAAGCAATGGAAAAACCaa                                                                             <  7:19044/95‑1 (MQ=255)                                                           tCGTTTTCGTAAAAATAAACGCCGCCTGCGCGAGCGTTTGCGTCAGCGTATTTTTTTCAGAGATAAAGTGGTGCCGGAAGCAATGGAAAAACCaa                                                                             >  8:19044/1‑95 (MQ=255)                                                                   gTAAAAATAAACGCCGCCTGCGCGAGCGTTTGCGTCAGCGTATTTTTTTCAGAGATAAAGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGcgcc               >  2:1127/1‑149 (MQ=255)                                                                    tAAAAATAAACGCCGCCTGCGCGAGCGTTTGCGTCAGCGTATTTTTTTCAGAGATAAAGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCg                                          >  5:110625/1‑121 (MQ=255)                                                                    tAAAAATAAACGCCGCCTGCGCGAGCGTTTGCGTCAGCGTATTTTTTTCAGAGATAAAGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCg                                          <  6:110625/121‑1 (MQ=255)                                                                               gccgccTGCGCGAGCGTTTGCGTCAGCGTATTTTTTTCAGAGATAAAGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCtg   >  2:89584/1‑149 (MQ=255)                                                                                ccgccTGCGCGAGCGTTTGCGTCAGCGTATTTTTTTCAGAGATAAAGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCtgt  <  1:1127/149‑1 (MQ=255)                                                                                                  |                                                                                                                                   ATCTAACCGATTAAACAACAGAGGTTGCTATGCTGTCGCGTCGGGGTCATCGGTTAAGTCGTTTTCGTAAAAATAAACGCCGCCTGCGCGAGCGTTTGCGTCAGCGTATTTTTTTCAGAGATAAAGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGT  >  NC_000913/1179491‑1179718

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