New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	2036832 =	26 (0.660)	4 (0.110)	3/136	4.5	11.8%	coding (127/1113 nt)	wbbC	putative glycosyltransferase WbbC
	?	REL606	2037060 =	36 (1.010)					coding (1071/1176 nt)	wzy	O antigen polymerase (Wzy)
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

AGGTAATATTATTGTTTTCTGGAAAATTGAACAATTCTTTAACAGGGATATTAGCCTGGTGAACAACTATTATTTCACAATCAATCCCATCGAATATTTCATTAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCAA  >  REL606/2036763‑2036903                                                                      |                                                                        aGGTAATATTATTGTTTTCTGGAAAATTGAACAATTCTTTAACAGGGATATTAGCCTGGTGAACAACTATTATTTc                                                                   >  2:583278/1‑76 (MQ=255)     aaTATTATTGTTTTCTGGAAAATTGAACAATTCTTTAACAGGGATATTAGCCTGGTGAACAACTATTATTTCAcaa                                                               >  2:186876/1‑76 (MQ=255)       tattatTGTTTTCTGGAAAATTGAACAATTCTTTAACAGGGATATTAGCCTGGTGAACAACTATTATTTCAcaatc                                                             >  2:1115729/1‑76 (MQ=255)        attattGTTTTCTGGAAAATTGAACAATTCTTTAACAGGGATATTAGCCTGGTGAACAACTATTATTTCAcaatca                                                            <  2:1090312/76‑1 (MQ=255)        attattGTTTTCTGGAAAATTGAACAATTCTTTAACAGGGATATTAGCCTGGTGAACAACTATTATTTCAcaatca                                                            <  2:1268131/76‑1 (MQ=255)         ttattGTTTTCTGGAAAATTGAACAATTCTTTAACAGGGATATTAGCCTGGTGAACAACTATTATTTCAcaatcaa                                                           <  1:1146021/76‑1 (MQ=255)          tattGTTTTCTGGAAAATTGAACAATTCTTTAACAGGGATATTAGCCTGGTGAACAACTATTATTTCAcaatcaat                                                          <  1:316237/76‑1 (MQ=255)          tattGTTTTCTGGAAAATTGAACAATTCTTTAACAGGGATATTAGCCTGGTGAACAACTATTATTTCAcaatcaat                                                          >  2:732237/1‑76 (MQ=255)           attGTTTTCTGGAAAATTGAACAATTCTTTAACAGGGATATTAGCCTGGTGAACAACTATTATTTCACaatcaatc                                                         <  2:1052593/76‑1 (MQ=255)            ttGTTTTCTGGAAAATTGAACAATTCTTTAACAGGGATATTAGCCTGGTGAACAACTATTATTTCACAATCAATcc                                                        <  2:375419/76‑1 (MQ=255)              gTTTTCTGGAAAATTGAACAATTCTTTAACAGGGATATTAGCCTGGTGAACAACTATTATTTCACAATCAATCCCa                                                      >  1:15836/1‑76 (MQ=255)                 ttCTGGAAAATTGAACAATTCTTTAACAGGGGTATTAGCCTGGTGAACAACTATTATTTCACAATCAATCCCATCg                                                   >  2:1015843/1‑76 (MQ=255)                 ttCTGGAAAATTGAACAATTCTTTAACAGGGATATTAGCCTGGTGAACAACTATTATTTCACAATCAATCCCATCg                                                   >  1:761325/1‑76 (MQ=255)                     ggAAAATTGAACAGTTCTTTAACAGGGATATTAGCCTGGTGAACAACTATTATTTCACAATCAATCCCATCGAata                                               >  1:844264/1‑76 (MQ=255)                        aaaTTGAACAATTCTTTAACAGGGATATTAGCCTGGTGAACAACTATTATTTCACAATCAATCCCATCGAATAttt                                            >  2:1202231/1‑76 (MQ=255)                           ttGAACAATTCTTTAACAGGGATATTAGCCTGGTGAACAACTATTATTTCACAATCAATCCCATCGAATATTTCAt                                         <  2:767645/76‑1 (MQ=255)                                cAATTCTTTAACAGGGATATTAGCCTGGTGAACAACTATTATTTCACAATCAATCCCATCGAATATTTCATTAata                                    >  2:464804/1‑76 (MQ=255)                                   ttCTTTAACAGGGATATTAGCCTGGAGAACAACTATTATTTCACAATCAATCCCATCGCATATTTCATTAATAttt                                 >  1:808623/1‑76 (MQ=255)                                           cAGGGATATTAGCCTGGTGAACAACTATTATTTCACAATCAATCCCATCGAATATTTCATTAATATTTTTTGCAAt                         <  2:43237/76‑1 (MQ=255)                                              ggATATTAGCCTGGTGAACAACTATTATTTCACAATCAATCCCATCGAATATTTCATTAATATTTTTTGCAATTaa                      <  1:138082/76‑1 (MQ=255)                                              ggATATTAGCCTGGTGAACAACTATTATTTCACAATCAATCCCATCGAATATTTCATTAATATTTTTTGCAATTaa                      >  1:275413/1‑76 (MQ=255)                                              ggATATTAGCCTGGTGAACAACTATTATTTCACAATCAATCCCATCGAATATTTCATTAATATTTTTTGCAATTaa                      >  1:776385/1‑76 (MQ=255)                                                 tatTAGCCTGGTGAACAACTATTATTTCACAATCAATCCCATCGAATATTTCATTAATATTTTTTGCAATTAAAtc                   <  2:1180774/76‑1 (MQ=255)                                                             gAACAACTATTATTTCACAATCAATCCCATCGAATATTTCATTAATACTTTTTGCAATTAAATCTCTAAAGTATAc       <  2:1194484/76‑1 (MQ=255)                                                                  acTATTATTTCACAATCAATCCCATCGAATATTTCATTAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCaa  >  2:345225/1‑76 (MQ=255)                                                                  acTATTATTTCACAATCAATCCCATCGAATATTTCATTAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCaa  >  2:914582/1‑76 (MQ=255)                                                                      |                                                                        AGGTAATATTATTGTTTTCTGGAAAATTGAACAATTCTTTAACAGGGATATTAGCCTGGTGAACAACTATTATTTCACAATCAATCCCATCGAATATTTCATTAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCAA  >  REL606/2036763‑2036903

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