New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	= 2042003	34 (1.540)	4 (0.190)	3/144	2.7	11.2%	coding (818/1428 nt)	wzx	O antigen translocase (Wzx)
	?	REL606	= 2042213	31 (1.460)					coding (608/1428 nt)	wzx	O antigen translocase (Wzx)
Rejected: Frequency below/above cutoff threshold.

ATATATTTGAATAAAAAGCTTCAAGTACTCCAGCAATTAGTAAGTTACTGCTAAATGAAAATAGCTTTTTGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGTACAATACAATTTATGACACTATAGG  >  REL606/2042139‑2042288                                                                           |                                                                            atatatTTGAATAAAAAGCTTCAAGTACTCCAGCAATTAGTAAGTTACTGCTAAATGAAAATAGCTTTTTGAATGt                                                                            >  1:566264‑M2/1‑75 (MQ=255)        tGAATAAAAAGCTTCAAGTACTCCAGCAATTAGTAAGTTACTGCTAAATGAAAATAGCTTTTTGAATGTTTCATAg                                                                     <  2:287501/76‑1 (MQ=255)            tAAAAAGCTTCAAGTACTCCAGCAATTAGTAAGTTACTGCTAAATGAAAATAGCTTTTTGAATGTTTCATAGCaaa                                                                 <  1:536860/76‑1 (MQ=255)                  gCTTCAAGTACTCCAGCAATTAGTAAGTTACTGCTAAATGAAAATAGCTTTTTGAATGTTTCATAGCAAAATTtct                                                           <  1:665666/76‑1 (MQ=255)                   cTTCAAGTACTCCAGCAATTAGTAAGTTACTGCTAAATGAAAATAGCTTTTTGAATGtttataatggcaaagggaa                                                          >  2:375743‑M2/1‑57 (MQ=255)                   cTTCAAGTACTCCAGCAATTAGTAAGTTACTGCTAAATGAAAATAGCTTTTTGAATGTTTCATAGCAAAATTtctc                                                          >  1:505492/1‑76 (MQ=255)                      cAAGTACTCCAGCAATTAGTAAGTTACTGCTAAATGAAAATAGCTTTTTGAATGtttataatggcaaagggaaaat                                                       >  2:391182‑M2/1‑54 (MQ=255)                      cAAGTACTCCAGCAATTAGTAAGTTACTGCTAAATGAAAATAGCTTTTTGAATGtttataatggcaaagggaaaat                                                       >  2:391175‑M2/1‑54 (MQ=255)                        aGTACTCCAGCAATTAGTAAGTTACTGCTAAATGAAAATAGCTTTTTGAATGTTTCATAGCAAAATTTCTCTTTTg                                                     <  2:156583/76‑1 (MQ=255)                           aCTCCAGCAATTAGTAAGTTACTGCTAAATGAAAATAGCTTTTTGAATGTTTCATAGCAAAATTTCTCTTTTGGAt                                                  >  1:174703/1‑76 (MQ=255)                           aCTCCAGCAATTAGTAAGTTACTGCTAAATGAAAATAGCTTTTTGAATGTTTCATAGCAAAATTTCTCTTTTGGAt                                                  >  1:174711/1‑76 (MQ=255)                           aCTCCAGCAATTAGTAAGTTACTGCTAAATGAAAATAGCTTTTTGAATGTTTCATAGCAAAATTTCTCTTTTGGAt                                                  <  2:370293/76‑1 (MQ=255)                           aCTCCAGCAATTAGTAAGTTACTGCTAAATGAAAATAGCTTTTTGAATGTTTCATAGCAAAATTTCTCTTTTGGAt                                                  <  2:370291/76‑1 (MQ=255)                               cAGCAATTAGTAAGTTACTGCTAAATGAAAATAGCTTTTTGAATGtttcatagcaaagggaaaatagatgatgctt                                              >  2:89325‑M2/1‑45 (MQ=255)                                aGCAATTAGTAAGTTACTGCTAAATGAAAATAGCTTTTTGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCAt                                             <  2:14804/76‑1 (MQ=255)                                    aTTAGTAAGTTACTGCTAAATGAAAATAGCTTTTTGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGAt                                         <  1:150589/76‑1 (MQ=255)                                       aGTAAGTTACTGCTAAATGAAAATAGCTTTTTGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATtat                                      <  1:559040/76‑1 (MQ=255)                                             tttcTGCTTAATGAAAATAGCTTTTTGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAatat                                <  1:105193/73‑1 (MQ=255)                                                   cTAAATGAAAATAGCTTTTTGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGta                          >  2:18250/1‑76 (MQ=255)                                                    tAAATGAAAATAGCTTTTTGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGtac                         >  1:705357/1‑76 (MQ=255)                                                     aaaTGAAAATAGCTTTTTGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGtaca                        >  2:563410/1‑76 (MQ=255)                                                       aTGAAAATAGCTTTTTGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGtacaat                      <  1:395943/76‑1 (MQ=255)                                                       aTGAAAATAGCTTTTTGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGtacaat                      <  2:525220/76‑1 (MQ=255)                                                        tGAAAATAGCTTTTTGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGtacaata                     <  1:6634/76‑1 (MQ=255)                                                        tGAAAATAGCTTTTTGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGtacaata                     <  1:6628/76‑1 (MQ=255)                                                             aTAGCTTTTTGAAGGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGTACAATACAAtt                <  1:290876/76‑1 (MQ=255)                                                              tAGCTTTTTGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGTACAATACAAttt               <  2:357960/76‑1 (MQ=255)                                                               aGCTTTTTGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGTACAATACAATTTa              <  1:79257/76‑1 (MQ=255)                                                               aGCTTTTTGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGTACAATACAATTTa              >  1:565612/1‑76 (MQ=255)                                                                gCTTTTTGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGTACAATACAATTTAt             <  1:566191/76‑1 (MQ=255)                                                                gCTTTTTGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGTACAATACAATTTAt             >  1:359049/1‑76 (MQ=255)                                                                gCTTTTTGAATGTTTCATAGCAAAATTTCTCGTTTGGATACCATGGATTATAAATATTCAGTACAATACAATTTAt             <  1:566190/76‑1 (MQ=255)                                                                  tttttGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGTACAATACAATTTATGa           >  2:400437/1‑76 (MQ=255)                                                                   ttttGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGTACAATACAATTTATGac          >  1:380199/1‑76 (MQ=255)                                                                     ttGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGTACAATACAATTTATGacac        >  2:230714/1‑76 (MQ=255)                                                                           gTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGTACAATACAATTTATGACACTATAgg  <  1:32698/76‑1 (MQ=255)                                                                           |                                                                            ATATATTTGAATAAAAAGCTTCAAGTACTCCAGCAATTAGTAAGTTACTGCTAAATGAAAATAGCTTTTTGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGTACAATACAATTTATGACACTATAGG  >  REL606/2042139‑2042288

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