New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_012967	16972 =	0 (0.000)	9 (0.940)	3/50	NT	100%	intergenic (‑71/‑516)	mokC/nhaA	type I toxin‑antitoxin system toxin MokC/Na+/H+ antiporter NhaA
	?	NC_012967	588495 =	NA (NA)					intergenic (‑222/‑47)	ybdK/ECB_RS24840	YbdK family carboxylate‑amine ligase/IS3‑like element IS150 family transposase

ACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGATAGTAAAATATCACTCGGGGCTTTTCTCTA  >  NC_012967/16939‑17007                                  |                                    aCTCTACATGTGTTCAGCATATAGGGGGCCTCGGGt                                   >  1:59974/1‑36 (MQ=255)  attacatcgtccaacaaaatgggtgcagtacaGGtt                                  >  2:925445‑M1/33‑36 (MQ=255)   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                             >  1:984453‑M1/28‑36 (MQ=255)       atcgtccaacaaaatgggtgcagtacaGGTTGATAg                             >  2:735209‑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   >  2:208611‑M1/2‑36 (MQ=255)                                 aGGTTGATAGTAAAATATCACTCGGGGCTTTTctct   >  1:50335‑M1/2‑36 (MQ=255)                                  ggTTGtactgaccccaaaaagttggacagttaaaca  <  2:31727‑M1/36‑32 (MQ=255)                                  |                                    ACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGATAGTAAAATATCACTCGGGGCTTTTCTCTA  >  NC_012967/16939‑17007

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