New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3854888	14 (0.540)	3 (0.130)	3/114	2.2	21.4%	coding (967/1185 nt)	emrD	multidrug efflux system protein
	?	NC_000913	= 3854930	10 (0.450)					coding (1009/1185 nt)	emrD	multidrug efflux system protein

TGCTCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAA  >  NC_000913/3854826‑3854951                                                               |                                                                tgctCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGcgcg                                                            <  2:132459/68‑1 (MQ=255) tgctCGGTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGcgcg                                                            <  2:3598/68‑1 (MQ=255)    tCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATg                                                         >  2:150992/1‑68 (MQ=255)       ttCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAg                                                      >  1:592886/1‑68 (MQ=255)         accgccgCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGAGATGGAGca                                                    >  1:299054/2‑67 (MQ=255)           cgccgcGCTGCTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGt                                                  <  2:202688/68‑1 (MQ=255)             ccgcGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGcgccagccgtgcccgcc                                                <  1:708056‑M1/68‑18 (MQ=255)                           ttCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCggg                                  >  1:25701/1‑68 (MQ=255)                           ttCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCggg                                  >  2:804431/1‑68 (MQ=255)                           ttCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCggg                                  >  2:726771/1‑68 (MQ=255)                           ttCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCggg                                  >  1:922455/1‑68 (MQ=255)                            tCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGcgccagccgtgcccgccaggaaggggaacggc                                 <  1:267775‑M1/68‑33 (MQ=255)                                           tGTTTCCGCTGGCGACCAGCGcgccagccgtgcccgccaggaaggggaacggctccatcgcgccgctg                  >  1:561130‑M1/1‑21 (MQ=255)                                                 cGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTcggc            >  1:301818/1‑68 (MQ=255)                                                      gCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCTCGGCTGGCGCGCTGGGCGGCGGTCt       <  2:299054/68‑1 (MQ=255)                                                        gACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGc     <  1:85394/68‑1 (MQ=255)                                                        gACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGc     <  1:150992/68‑1 (MQ=255)                                                           cAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGcgctggtcgccagcggaaaca  >  1:561130‑M2/1‑47 (MQ=255)                                                               |                                                                TGCTCGTTCCCGCCGCGCTGTTCTTTTTCGGTGCCGGGATGCTGTTTCCGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAA  >  NC_000913/3854826‑3854951

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