New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_012967	1733645 =	0 (0.000)	14 (1.580)	13/46	NT	100%	coding (681/1413 nt)	pykF	pyruvate kinase PykF
	?	NC_012967	= 2775875	NA (NA)					intergenic (‑45/+204)	ECB_RS25125/cysH	IS3 family transposase/phosphoadenosine phosphosulfate reductase

CACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGA  >  NC_012967/1733610‑1733677                                    |                                 attccatttcatcgtccaacaaaatgggtgcagtac                                  >  2:258912‑M1/36‑36 (MQ=255) attccatttcatcgtccaacaaaatgggtgcagtac                                  >  2:527704‑M1/36‑36 (MQ=255) attccatttcatcgtccaacaaaatgggtgcagtac                                  <  1:730833‑M1/1‑1 (MQ=255) actccatttcatcgtccaacaaaatgggtgcagtac                                  <  1:168861‑M1/1‑1 (MQ=255)  ttccatttcattgtccaacaaaatgggtgcagtaca                                 >  1:961481‑M1/35‑36 (MQ=255)  ttccatttcatcgtccaacaaaatgggtgctgtaca                                 <  2:408006‑M1/2‑1 (MQ=255)  ttccatttcatcgtccaacaaaatgggtgcagtaca                                 >  1:798407‑M1/35‑36 (MQ=255)  ttccatttcatcgtccaacaaaatgggtgcagtaca                                 >  2:562385‑M1/35‑36 (MQ=255)  ttccatttcatcgtccaacaaaatgggtgcagtaca                                 <  1:725966‑M1/2‑1 (MQ=255)  acaTCATCTCCAAAATCGAAAACCAGGAAGGCCTca                                 >  2:655813/1‑36 (MQ=255)    ccatttcatcgtccaacaaaatgggtgcagtacaac                               >  1:518720‑M1/33‑36 (MQ=255)          tcCAAAATCGAAAACCAGGAAGGCCTCAACAtgtac                         >  1:802586‑M1/1‑31 (MQ=255)          catcgtccaacaaaatgggtgcagtaCAACAACTTc                         >  2:181002‑M1/27‑36 (MQ=255)           ccAAAATCGAAAACCAGGAAGGCCTCAACAtgtact                        >  2:988009‑M1/1‑30 (MQ=255)           atcgtccaacaaaatgggtgcagtaCAACAACTTcg                        >  1:669229‑M1/26‑36 (MQ=255)            tcgtccaacaaaatgggtgcagtaCAACAACTTcga                       >  1:703351‑M1/25‑36 (MQ=255)             cgtccaacaaaatgggtgtagtaCAACAACTTcgac                      >  1:428739‑M1/24‑36 (MQ=255)               tccaacaaaatgggtgcagtaCAACAACTTcgacga                    >  1:237377‑M1/22‑36 (MQ=255)               aaTCGAAAACCAGGAAGGCCTCAACAtgtactgacc                    >  2:429035‑M1/1‑26 (MQ=255)                ccaacaaaatgggtgcagtaCAACAACTTCGACGaa                   <  2:14550‑M1/16‑1 (MQ=255)                ccaacaaaatgggtgcagtaCAACAACTTCGACGaa                   <  2:649258‑M1/16‑1 (MQ=255)                aTCGAAAACCAGGAAGGCCTCAACAtgtactgaccc                   <  1:691462‑M1/36‑12 (MQ=255)                    aaaaCCAGGAAGGCCTCAACAtgtactgaccccaaa               <  2:323199‑M1/36‑16 (MQ=255)                     aaagtgggtgcagtaCAACAACTTCGACGAAATCCt              >  2:762145‑M1/16‑36 (MQ=255)                     aaaatgggtgcagtaCAACAACTTCGACGAAATCCt              <  1:817429‑M1/21‑1 (MQ=255)                      aaatgggtgcagtaCAACAACTTCGACGAAATCCTc             >  2:589228‑M1/15‑36 (MQ=255)                       aatgggtgcagtaCAACAACTTCGACGAAATCCTCg            <  2:968700‑M1/23‑1 (MQ=255)                         tgggtgcagtaCAACAACTTCGACGAAATCCTCGaa          >  2:544944‑M1/12‑36 (MQ=255)                          aggaggGCCTCAACAtgtactgaccccaaaaagttg         <  1:589780‑M1/36‑22 (MQ=255)                          aGGAAGGTCTCAACAtgtactgacaccaaaaagttg         >  2:279697‑M1/1‑15 (MQ=255)                           ggtgcagtaCAACAACTTCGACGAAATCCTCGAAGc        <  2:781527‑M1/27‑1 (MQ=255)                           ggAAGGCCTCAACAtgtactgaccccaaaaagttgg        <  1:651316‑M1/36‑23 (MQ=255)                            gtgcagtaCAACAACTTCGACGAAATCCTCGAAGcc       <  2:314777‑M1/28‑1 (MQ=255)                               cagtaCAACAACTTCGACGAAATCCTCGAAGCctct    <  2:284243‑M1/31‑1 (MQ=255)                                agtaCAACAACTTCGACGAAATCCTCGCAGCCTCTg   <  1:381863‑M1/32‑1 (MQ=255)                                 gaaCAACAACTTCGACGAAATCCTCGAAGCCTCTGa  >  2:556875‑M1/4‑36 (MQ=255)                                    |                                 CACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGA  >  NC_012967/1733610‑1733677

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