New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AP009048	2109295 =	40 (0.710)	4 (0.080)	3/90	NT	9.4%	coding (69/1167 nt)	wbbH	O‑antigen polymerase
	?	AP009048	2109315 =	42 (0.850)					coding (49/1167 nt)	wbbH	O‑antigen polymerase

ATTAATGATCTATCTTGTAATTAGTGTCTTTCTCATTACAGCATTTATCTGTTTATATCTTAAGAAGGATATA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AP009048/2109367‑2109295 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gatataAACAGATAAATGCTGTAATGAGAAAGACACTAATTACAAGATAGATCATTAATCCGTACTCATTATATT  >  AP009048/2109315‑2109383                                                                                                                                                ATTAATGATCTATCTTGTAATTAGTGTCTTTCTCATTACAGCATTTATCTGTTTATATCTTAAGAAGGATA                                                                         >  1:3664604/1‑71  TTAATGATCTATCTTGTAATTAGTGTCTTTCTCATTACAGCATTTATCTGTTTATATCTTAAGAAGGATAT                                                                        <  1:2506331/71‑1   TAATGATCTATCTTGTAATTAGTGTCTTTCTCATTACAGCATTTATCTGTTTATATCTTAAGAAGGATATA                                                                       >  1:4653668/1‑71          CTATCTTGTAATTAGTGTCTTTCTCATTACAGCATTTATCTGTTTATATCTTAAGAAGGATATA                                                                       <  1:306850/64‑1          CTATCTTGTAATTAGTGTCTTTCTCATTACAGCATTTATCTGTTTATATCTTAAGAAGGATA                                                                         >  1:4267201/1‑62                       AGTGTCTTTCTCATTACAGCATTTATCTGTTTATATCCTTCTTAAGATATAAACAGATAAATGCTGTAAT                                                    <  1:4172186/70‑1                                                     TTATATCTTAAGAAGGATATAAACAGATAAATGCTGTAATGAGAAAGACACTAAT                                     <  1:2436248/55‑1                                                         ATCTTAAGAAGTATATAAACAGATAAATGCTGTAAT                                                    <  1:11928/36‑1                                                         ATCTTAAGAAGGATATAAACAGATAAATGCTGTAAT                                                    <  1:1868489/36‑1                                                                     ATATAAACAGATAAATGCTGTAATGAGAAAGACACTAATTACAAGATAG                           <  1:671615/49‑1                                                                       ATAAACAGATAAATGCTGTAATGAGAAAGACACTAATTACAAGATA                            >  1:4421315/1‑46                                                                        TAAACAGATAAATGCTGTAATGAGAAAGACACTAATTACAAGATAGATCATTAATCCGTACTCATTATATT  >  1:1051049/1‑71                                                                        TAAACAGATAAATGCTGTAATGAGAAAGACACTAATTACAAGATAGATCATTAATCCGTACTCATTATA    <  1:1213527/69‑1                                                                         AAACAGATAAATGCTGTAATGAGAAAGACACTAATTACAAGA                              >  1:4971064/1‑42                                                                                                                                                ATTAATGATCTATCTTGTAATTAGTGTCTTTCTCATTACAGCATTTATCTGTTTATATCTTAAGAAGGATATA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AP009048/2109367‑2109295 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gatataAACAGATAAATGCTGTAATGAGAAAGACACTAATTACAAGATAGATCATTAATCCGTACTCATTATATT  >  AP009048/2109315‑2109383

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