New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3854880 =	125 (1.440)	3 (0.040)	3/246	NT	2.6%	coding (959/1185 nt)	emrD	multidrug efflux system protein
	?	NC_000913	3854940 =	107 (1.330)					coding (1019/1185 nt)	emrD	multidrug efflux system protein
Rejected: Frequency below/above cutoff threshold.

ATTAAGGTCATCAACAACCCCAGGCTGCCCTGACCGGTTTGCGGCAACATCGCAGAGAGCGACGCCAGCACGCCGGAACCAATGTTTTGCAGACCGCCGACCAGCGCGCCAGCCGTGCCCGCCAGGAAGGGGAACGGCTCCATCGCGCCGCTGGTCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3855036‑3854880 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgctggtcgGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCA  >  NC_000913/3854940‑3855087                                                                                                                                                                                                                                                                                                                   ATTAAGGTCATCAACAACCCCAGGCTGCCCTGACCGGTTTGCGGCAACATCGCAGAGAGCGACGCCAGCACGCCGGAACCAATGTTTTGCAGACCGCCGACCAGCGCGCCAGCCGTGCCCGCCAGGAAGGGGAACGGCTCCATCGCGCC                                                                                                                                                              >  3:35887/1‑149  TTAAGGTCATCAACAACCCCAGGCTGCCCTGACCGGTTTGCGGCAACATCGCAGAGAGCGACGCCAGCACGCCGGAACCAATGTTTTGCAGACCGCCGACCAGCGCGCCAGCCGTGCCCGCCAGGAAGGGGAACGGCTCCATCGCGCCG                                                                                                                                                             >  2:131211/1‑149    AAGGTCATCAACAACCCCAGGCTGCCCTGACCGGTTTGCGGCAACATCGCAGAGAGCGACGCCAGCACGCCGGAACCAATGTTTTGCAGACCGCCGACCAGCGCGCCAGCCGTGCCCGCCAGGAAGGGGAACGGCTCCATCGCGCCGC                                                                                                                                                            <  8:114298/148‑1     AGGTCATCAACAACCCCAGGCTGCCCTGACCGGTTTGCGGCAACATCGCAGAGAGCGACGCCAGCACGCCGGAACCAATGTTTTGCAGACCGCCGACCAGCGCGCCAGCCGTGCCCGCCAGGAAGGGGAATGGCTCCATCGCGCCGCTG                                                                                                                                                          <  3:217879/149‑1        TCATCAACAACCCCAGGCTGCCCTGACCGGTTTGCGGCAACATCGCAGAGAGCGACGCCAGCACGCCGGAACCAATGTTTTGCAGACCGCCGACCAGCGCGCCAGCCGTGCCCGCCAGGAAGGGGAACGGCTCCATCGCGCCGCTGGTC                                                                                                                                                       <  1:339787/149‑1        TCATCAACAACCCCAGGCTGCCCTGACCGGTTTGCGGCAACATCGCAGAGAGCGACGCCAGCACGCCGGAACCAATGTTTTGCAGACCGCCGACCAGCGCGCCAGCCGTGCCCGCCAGGAAGGGGAACGGCTCCATCGCGCCGCTGGTC                                                                                                                                                       >  8:3830/1‑149             AACAACCCCAGGCTGCCCTGACCGGTTTGCGGCAACATCGCAGAGAGCGACGCCAGCACGCCGGAACCAATGTTTTGCAGACCGCCGACCAGCGCGCCAGCCGTGCCCGCCAGGAAGGGGAACGGCTCCATCGCGCCGCTGGTC                                                                                                                                                       <  8:129165/144‑1             AACAACCCCAGGCTGCCCTGACCGGTTTGCGGCAACATCGCAGAGAGCGACGCCAGCACGCCGGAACCAATGTTTTGCAGACCGCCGACCAGCGCGCCAGCCGTGCCCGCCAGGAAGGGGAACGGCTCCATCGCGCCGCTGGTC                                                                                                                                                       >  7:129165/1‑144                                                                                               CGCCGACCAGCGCGCCAGCCGTGCCCGCCAGGAAGGGGAACGGCTCCATCGCGCCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTT                                                                 >  8:386169/1‑148                                                                                                 CCGACCAGCGCGCCAGCCGTGCCCGCCAGGAAGGGGAACGGCTCCATCGCGCCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGAT                                                              >  5:139657/1‑149                                                                                                                  CGTGCCCGGCAGGCAGGGGAACGGCTCCATCCCGCCGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGGGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGGTGTTGATGACCTTAATGGGATTGT                                             <  6:139657/149‑1                                                                                                                                                     CGCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGG          <  5:86237/149‑1                                                                                                                                                      GCTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGA         <  8:2971/149‑1                                                                                                                                                       CTGGTCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGAT        >  3:321992/1‑149                                                                                                                                                           TCGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCG    <  2:186718/149‑1                                                                                                                                                            CGGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGC   <  3:96361/149‑1                                                                                                                                                             GGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCA  >  6:297714/1‑149                                                                                                                                                                                                                                                                                                                   ATTAAGGTCATCAACAACCCCAGGCTGCCCTGACCGGTTTGCGGCAACATCGCAGAGAGCGACGCCAGCACGCCGGAACCAATGTTTTGCAGACCGCCGACCAGCGCGCCAGCCGTGCCCGCCAGGAAGGGGAACGGCTCCATCGCGCCGCTGGTCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3855036‑3854880 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgctggtcgGCGGTCTGCAAAACATTGGTTCCGGCGTGCTGGCGTCGCTCTCTGCGATGTTGCCGCAAACCGGTCAGGGCAGCCTGGGGTTGTTGATGACCTTAATGGGATTGTTGATCGTGCTGTGCTGGCTGCCGCTGGCGACGCGGATGTCGCA  >  NC_000913/3854940‑3855087

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