New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_012967	= 16976	0 (0.000)	6 (0.680)	6/46	NT	100%	intergenic (‑75/‑512)	mokC/nhaA	type I toxin‑antitoxin system toxin MokC/Na+/H+ antiporter NhaA
	?	NC_012967	= 590469	NA (NA)					intergenic (+24/‑74)	ECB_RS24845/hokE	IS3 family transposase/type I toxin‑antitoxin system toxin HokE

TCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGATAGTAAAATATCACTCGGGGCTTTTCTCTATCTG  >  NC_012967/16941‑17011                                    |                                    tttcatcgtccaacaaaatgggtgcagtacaGGTTg                                     >  2:167270‑M1/32‑36 (MQ=255)  ttcatcgtccaacaaaatgggtgcagtacaGGTTGa                                    >  1:744184‑M1/31‑36 (MQ=255)    catcgtccaacaaaatgggtgcagtacaGGTTGATa                                  <  2:125881‑M1/8‑1 (MQ=255)     atcgtgcgacaaaatgggtgcagtacaGGTTGATAg                                 >  2:833501‑M1/28‑36 (MQ=255)     atcgtccaacaaaatgggtgcagtacaGGTTGATAg                                 >  2:735209‑M1/28‑36 (MQ=255)     atcgtccaacaaaatgggtgcagtacaGGTTGATAg                                 >  1:984453‑M1/28‑36 (MQ=255)       tgtgTTCAGCATATAGGGGGCCTCGGGTTGtactga                               <  2:464831‑M1/36‑7 (MQ=255)        gtgtTCAGCATATAGGGGGCCTCGGGTTGtactgac                              >  1:260324‑M1/1‑29 (MQ=255)            tCAGCATATAGGGGGCCTCGGGTTGtactgacccca                          <  2:594939‑M1/36‑12 (MQ=255)            aacaaaatgggtgcagtacaGGTTGATAGTAAAata                          >  1:125870‑M1/21‑36 (MQ=255)            aacaaaatgggtgcagtacaGGTTGATAGTAAAata                          <  1:73612‑M1/16‑1 (MQ=255)             cAGCATATAGGGGGCCTCGGGTTGtactgaccccaa                         >  1:97758‑M1/1‑24 (MQ=255)                agatgggtgcggtacaGGTTGATAGTAAAATATCAc                      >  1:558839‑M1/17‑36 (MQ=255)                 aatgggtgcagtacaGGTTGATAGTAAAATATCACt                     >  1:270010‑M1/16‑36 (MQ=255)                           gtacaGGTTGATAGTAAAATATCACTCGGGGCtttt           >  1:857310‑M1/6‑36 (MQ=255)                            cTCGGGTTGtactgaccccaaaaagttgaacagtta          >  2:814709‑M1/1‑9 (MQ=255)                             tCGGGTTGtactgaccccaaaaagttggacagttaa         >  1:812306‑M1/1‑8 (MQ=255)                             acaGGTTGATAGTAAAATATCACTCGGGGCTTTtct         >  1:922425‑M1/4‑36 (MQ=255)                              cGGGTTGtactgaccccaaaaagttggacagttaaa        <  1:466506‑M1/36‑30 (MQ=255)                               aGGTTGATAGTAAAATATCACTCGGGGCTTTTctct       >  1:50335‑M1/2‑36 (MQ=255)                               aGGTTGATAGTAAAATATCACTCGGGGCTTTTctct       >  2:208611‑M1/2‑36 (MQ=255)                                ggTTGtactgaccccaaaaagttggacagttaaaca      <  2:31727‑M1/36‑32 (MQ=255)                                   tGATAGTAAAATATCACTCGGGGCTTTTCTCTATCt   >  2:203178/1‑36 (MQ=255)                                    gtactgaccccaaaaagttggacagttaaatacgag  >  1:557980‑M1/1‑1 (MQ=255)                                    gtactgaccccaaaaagttggacagttaaacacgag  <  2:259282‑M1/36‑36 (MQ=255)                                    gtactgaccccaaaaagttggacagttaaacacgag  >  1:706402‑M1/1‑1 (MQ=255)                                    gtactgaccccaaaaagttggacagttaaacacgag  <  2:146077‑M1/36‑36 (MQ=255)                                    gtactgaccccaaaaagttggacagttaaacacgag  >  1:595800‑M1/1‑1 (MQ=255)                                    gtactgaccccaaaaagttggacagttaaacacgag  <  1:570019‑M1/36‑36 (MQ=255)                                    gtactgaccccaaaaagttggacagttaaacacgag  >  1:521239‑M1/1‑1 (MQ=255)                                    |                                    TCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGATAGTAAAATATCACTCGGGGCTTTTCTCTATCTG  >  NC_012967/16941‑17011

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