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

CGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGC  >  NC_000913/3854874‑3854993                                                         |                                                                cGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTcggc                                                      >  1:301818/1‑68 (MQ=255)      gCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCTCGGCTGGCGCGCTGGGCGGCGGTCt                                                 <  2:299054/68‑1 (MQ=255)        gACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGc                                               <  1:150992/68‑1 (MQ=255)        gACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGc                                               <  1:85394/68‑1 (MQ=255)           cAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGcgctggtcgccagcggaaaca                                            >  1:561130‑M2/1‑47 (MQ=255)                          gCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGcgctggtcgccagcggaaacagcatcccggcaccga                             <  1:267775‑M2/68‑37 (MQ=255)                                    ttCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTc                   >  1:766600/1‑68 (MQ=255)                                         ggCGGGCACGGCTGGCGcgctggtcgccagcggaaacagcatcccggcaccgaaaaagaacagcgcgg              <  1:708056‑M2/68‑52 (MQ=255)                                           cGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGctctct            <  1:23765/68‑1 (MQ=255)                                           cGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGctctct            >  1:8753/1‑68 (MQ=255)                                                     tGGCTCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGc  <  1:793587/68‑1 (MQ=255)                                                     tGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGc  <  2:835785/68‑1 (MQ=255)                                                     tGGCGCGCTGGGCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGc  <  2:592886/68‑1 (MQ=255)                                                         |                                                                CGCTGGCGACCAGCGGCGCGATGGAGCCGTTCCCCTTCCTGGCGGGCACGGCTGGCGCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGC  >  NC_000913/3854874‑3854993

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.