New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	= 2042596	31 (1.400)	4 (0.200)	4/136	2.2	13.5%	coding (225/1428 nt)	wzx	O antigen translocase (Wzx)
	?	REL606	= 2042873	23 (1.150)					coding (494/546 nt)	rmlC	dTDP‑4‑dehydrorhamnose 3,5‑epimerase (RmlC)
Rejected: Frequency below/above cutoff threshold.

ATATAAGAAGATAAACAATAAAGTACAATATATATTAGTATTGATAGACAAATATTAAATATAAATACTGTCGAATAATCT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  REL606/2042516‑2042596 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ataatctCTGAGAAAGACATGAAACAACCATCATTCAGAGAACTTTTTGAGAAAGTTTAATTGATGAATAACAAACTAGCCA  <  REL606/2042873‑2042799                                                                                                                                                              ATATAAGAAGATAAACAATAAAGTACAATATATATTAGTATTGATAGACAAATATTAAATATAAATACTGTCGAAT                                                                                  >  2:45800/1‑76   ATAAGAAGATAAACAATAAAGTACAATATATATTAGTATTGATAGACAAATATTAAATATAAATACTGTCGAATAA                                                                                >  2:441189/1‑76                                                     TATTAAATATAAATACTGTCGAATAATCTCTGAGAAAGACATGAAACAACCATCATTCAGAGAACTTTTTGAGAAA                              >  2:529311/1‑76                                                      ATTAAATATAAATACTGTCGAATAATCTCTGAGAAAGACATGAAACAACCATCATTCAGAGAACTTTTTGAGAAAG                             >  2:463976/1‑76                                                          AATATAAATACTGTCGAATAATCTCTGAGAAAGACATGAAACAACCATCATTCAGAGAACTTTTTGAGAAAGTTTA                         >  2:661906/1‑76                                                                 ATACTGTCGAATAATCTCTGAGAAAGACATGAAACAACCATCATTCAGAGAACTTTTTGAGAAAGTTTAATTGATG                  >  2:537426/1‑76                                                                              ATCTCTGAGAAAGACATGAAACAACCATCATTCAGAGAACTTTTTGAGAAAGTTTAATTGATGAATAACAAACTAG     >  1:647005/1‑76                                                                                CTCTGAGAAAGACATGAAACAACCATAATTCAGAGAACTTTTTGAGAAAGTTTAATTGATGAATAACAAACTAGCC   >  2:429265/1‑76                                                                                 TCTGAGAAAGACATGAAACAACCATCATTCAGAGAACTTTTTGAGAAAGTTTAATTGATGAATAACAAACTAGCCA  <  2:146219/76‑1                                                                                                                                                              ATATAAGAAGATAAACAATAAAGTACAATATATATTAGTATTGATAGACAAATATTAAATATAAATACTGTCGAATAATCT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  REL606/2042516‑2042596 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ataatctCTGAGAAAGACATGAAACAACCATCATTCAGAGAACTTTTTGAGAAAGTTTAATTGATGAATAACAAACTAGCCA  <  REL606/2042873‑2042799

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