New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_012967	2772857 =	NA (NA)	13 (1.540)	7/44	NT	93.7%	intergenic (+378/‑54)	iap/ECB_RS13760	alkaline phosphatase isozyme conversion aminopeptidase/IS4 family transposase
	?	NC_012967	= 4524530	1 (0.100)					coding (502/549 nt)	fimA	type 1 fimbrial major subunit FimA

GGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGT  >  NC_012967/2772827‑2772882                               |                          ggAACACCCTGAACGGTGGGGAAAAACCCCATAAgc                      >  1:469934/1‑36 (MQ=255)         cagAACAGTGAGGAAAAACCCCATAAGCGCTAACtt              <  1:722775/34‑1 (MQ=17)         cTGAACGGTGGGTAAAAACCCCATAAGCGCTAACtt              <  2:789720/36‑1 (MQ=21)              cGGTGGGGAAAAACCCCATAAGCGCTAACTTAAggg         >  2:868949/1‑36 (MQ=255)                gTGGAGAAAAACCCCATAAGCGCTAACTTAAGGGtt       <  1:637595/36‑1 (MQ=14)                     gAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGt  <  2:254619/36‑1 (MQ=40)                               |                          GGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGT  >  NC_012967/2772827‑2772882

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