New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4559341 =	51 (1.110)	4 (0.090)	3/248	NT	7.9%	coding (186/1173 nt)	iadA	isoaspartyl dipeptidase
	?	NC_000913	4559360 =	46 (1.070)					coding (167/1173 nt)	iadA	isoaspartyl dipeptidase

CGCATTTGTATGCGCCGGAAGATCGGGGAATTTGCGATGTCCTCGTCGCTAACGGCAAAATTATCGCCGTTGCCAGCAATATCCCTTCTGACATTGTACCGAACTGCACGGTTGTCGATCTCAGTGGGCAGATCCTCTGCCCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4559483‑4559341 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tctgcccaCTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAGG  >  NC_000913/4559360‑4559497                                                                                                                                                                                                                                                                                           CGCATTTGTATGCGCCGGAAGATCGGGGAATTTGCGATGTCCTCGTCGCTAACGGCAAACTTATCGCCGTTGCCAGCAATATCCCTTCTGACATTGTACCGAACTGCACGGTTGTCGATCTCAGTGGGCAGATCCTCTGCCCA                                                                                                                                            >  2:150999/1‑143                                                                          CAGCAATATCCCTTCTGACATTGTACCGAACTGCACGGTTGTCGATCTCAGTGGGCAGATCCTCTGCCCACTGAGATCGACAACCGTGCAGTTCGGT                                                                                                                 >  1:60957/1‑97                                                                          CAGCAATATCCCTTCTGACATTGTACCGAACTGCACGGTTGTCGATCTCAGTGGGCAGATCCTCTGCCCACTGAGATCGACAACCGTGCAGTTCGGT                                                                                                                 <  2:60957/97‑1                                                                                                                             TGGGCAGATCCTCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCG                >  1:463198/1‑143                                                                                                                             TGGGCAGATCCTCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATAGTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCG                >  2:14349/1‑143                                                                                                                                       agCTGCCCACTGAGAGCGACAGCCGTGCAGTTCGGTACAAGGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAG      <  1:228614/141‑1                                                                                                                                        TCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGC     >  1:554021/1‑143                                                                                                                                           GCCCTCTGAGATCGACAACCGTGCAGTTAGATAAGGTGGCATAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAGG  <  2:185905/143‑1                                                                                                                                           GCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAGG  >  1:862837/1‑143                                                                                                                                                                                                                                                                                           CGCATTTGTATGCGCCGGAAGATCGGGGAATTTGCGATGTCCTCGTCGCTAACGGCAAAATTATCGCCGTTGCCAGCAATATCCCTTCTGACATTGTACCGAACTGCACGGTTGTCGATCTCAGTGGGCAGATCCTCTGCCCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4559483‑4559341 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tctgcccaCTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAGG  >  NC_000913/4559360‑4559497

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