New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	2041245 =	35 (0.610)	5 (0.090)	5/140	6.3	12.6%	coding (139/1113 nt)	vioA	VioA, involved in dTDP‑N‑acetylviosamine synthesis
	?	REL606	2041495 =	37 (0.690)					coding (1326/1428 nt)	wzx	O antigen translocase (Wzx)
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GGCTTGTAATGCTGTGATCAATGCAATTGTTGCATTGTTAAAAAGCGAGATATGCTGCACACCCAAAAATTCACATAATTTTTCCTCGAGTTCTTGATGAAATGGTCCATTATTTGTCAACCATTTATTTTTCCATATTTTTTCCAGGTA  >  REL606/2041171‑2041320                                                                           |                                                                            ggCTTGTAATGCTGTGATCAATGCAATTGTTGCATTGTTAAAAAGCGAGATATGCTGCACACCCAAAAATTCACAt                                                                            <  1:749812/76‑1 (MQ=255)  gCTTGTAATGCTGTGATCAATGCAATTGTTGCATTGTTAAAAAGCGAGATATGCTGCACACCCAAAAATTCACATa                                                                           >  1:1398896/1‑76 (MQ=255)     tGTAATGCTGTGATCAATGCAATTGTTGCATTGTTAAAAAGCGAGATATGCTGCACACCCAAAAATTCACATAAtt                                                                        <  1:844587/76‑1 (MQ=255)       tAATGCTGTGATCAATGCAATTGTTGCATTGTTAAAAAGCGAGATATGCTGCACACCCAAAAATTCACATAAtttt                                                                      <  1:236605/76‑1 (MQ=255)       tAATGCTGTGATCAATGCAATTGTTGCATTGTTAAAAAGCGAGATATGCTGCACACCCAAAAATTCACATAAtttt                                                                      <  1:236578/76‑1 (MQ=255)            cTGTGATCAATGCAATTGTTGCATTGTTAAAAAGCGAGATATGCTGCACACCCAAAAATTCACATAATTTTTCCTc                                                                 >  2:495667/1‑76 (MQ=255)               tgATCAATGCAATTGTTGCATTGTTAAAAAGCGAGATATGCTGCACACCCAAAAATTCACATAATTTTTCCTCGAg                                                              <  1:1452600/76‑1 (MQ=255)               tgATCAATGCAATTGTTGCATTGTTAAAAAGCGAGATATGCTGCACACCCAAAAATTCACATAATTTTTCCTCGAg                                                              >  1:263932/1‑76 (MQ=255)                gATCAATGCAATTGTTGCATTGTTAAAAAGCGAGATATGCTGCACACCCAAAAATTCACATAATTTTTCCTCGAGt                                                             <  1:1690305/76‑1 (MQ=255)                 aTCAATGCAATTGTTGCATTGTTAAAAAGCGAGATATGCTGCACACCCAAAAATTCACATAATTTTTCCTCGAGtt                                                            >  2:1464879/1‑76 (MQ=255)                 aTCAATGCAATTGTTGCATTGTTAAAAAGCGAGATATGCTGCACACCCAAAAATTCACATAATTTTTCCTCGAGtt                                                            >  2:209848/1‑76 (MQ=255)                  tCAATGCAATTGTTGCATTGTTAAAAAGCGAGATATGCTGCACACCCAAAAATTCACATAATTTTTCCTCGAGTTc                                                           <  1:641358/76‑1 (MQ=255)                   cAATGCAATTGTTGCATTGTTAAAAAGCGAGATATGCTGCACACCCAAAAATTCACATAATTTTTCCTCGAGTTCt                                                          <  1:303086/76‑1 (MQ=255)                    aaTGCAATTGTTGCATTGTTAAAAAGCGAGATATGCTGCACACCCAAAAATTCACATAATTTTTCCTCGAGTTCtt                                                         <  1:214668/76‑1 (MQ=255)                    aaTGCAATTGTTGCATTGTTAAAAAGCGAGATATGCTGCACACCCAAAAATTCACATAATTTTTCCTCGAGTTCtt                                                         <  1:334029/76‑1 (MQ=255)                      tGCAATTGTTGCATTGTTAAAAAGCGAGATATGCTGCACACCCAAAAATTCACATAATTTTTCCTCGAGTTCTtga                                                       <  1:1586927/76‑1 (MQ=255)                       gCAATTGTTGCATTGTTAAAAAGCGAGATATGCTGCACACCCAAAAATTCACATAATTTTTCCTCGAGTTCTtgat                                                      <  2:926362/76‑1 (MQ=255)                       gCAATTGTTGCATTGTTAAAAAGCGAGATATGCTGCACACCCAAAAATTCACATAATTTTTCCTCGAGTTCTtgat                                                      <  1:720425/76‑1 (MQ=255)                         aaTTGTTGCATTGTTAAAAAGCGAGATATGCTGCACACCCAAAAATTCACATAATTTTTCCTCGAGTTCTtgatga                                                    <  2:1642644/76‑1 (MQ=255)                           ttgttgCATTGTTAAAAAGCGAGATATGCTGCACACCCAAAAATTCACATAATTTTTCCTCGAGTTCTTGATGaaa                                                  <  2:1695638/76‑1 (MQ=255)                                   ttGTTAAAAAGCGAGATATGCTGCACACCCAAAAATTCACATAATTTTTCCTCGAGTTCTTGATGAAATGGTCCat                                          <  2:1356257/76‑1 (MQ=255)                                        aaaaaGCGAGATATGCTGCACACCCAAAAATTCACATAATTTTTCCTCGAGTTCTTGATGAAATGGTCCATTAttt                                     >  1:232783/1‑76 (MQ=255)                                           aaGCGAGATATGCTGCACACCCAAAAATTCACATAATTTTTCCTCGAGTTCTTGATGAAATGGTCCATTATTTGTc                                  <  1:652222/76‑1 (MQ=255)                                              cGAGATATGCTGCACACCCAAAAATTCACATAATTTTTCCTCGAGTTCTTGATGAAATGGTCCATTATTTGTCAAc                               >  2:250318/1‑76 (MQ=255)                                               gagaTATGCTGCACACCCAAAAATTCACATAATTTTTCCTCGAGTTCTTGATGAAATGGTCCATTATTTGTCAAcc                              >  2:1348007/1‑76 (MQ=255)                                                agaTATGCTGCACACCCAAAAATTCACATAATTTTTCCTCGAGTTCTTGATGAAATGGTCCATTATTTGTCAACCa                             >  2:26041/1‑76 (MQ=255)                                                 gaTATGCTGCACACCCAAAAATTCACATAATTTTTCCTCGAGTTCTTGATGAAATGGTCCATTATTTGTCAACCat                            <  2:258448/76‑1 (MQ=255)                                                   tatGCTGCACACCCAAAAATTCACATAATTTTTCCTCGAGTTCTTGATGAAATGGTCCATTATTTGTCAACCattt                          >  2:1480783/1‑76 (MQ=255)                                                    atGCTGCACACCCAAAAATTCACATAATTTTTCCTCGAGTTCTTGATGAAATGGTCCATTATTTGTCAACCattta                         <  2:557480/76‑1 (MQ=255)                                                     tgctgcACACCCAAAAATTCACATAATTTTTCCTCGAGTTCTTGATGAAATGGTCCATTATTTGTCAACCatttat                        <  1:875096/76‑1 (MQ=255)                                                           acacCCAAAAATTCACATAATTTTTCCTCGAGTTCTTGATGAAATGGTCCATTATTTGTCAACCATTTATTTTTcc                  <  2:1313391/76‑1 (MQ=255)                                                               ccAAAAATTCACATAATTTTTCCTCGAGTTCTTGATGAAATGGTCCATTATTTGTCAACCATTTATTTTTCCatat              <  2:177034/76‑1 (MQ=255)                                                                cAAAAATTCACATAATTTTTCCTCGAGTTCTTGATGAAATGGTCCATTATTTGTCAACCATTTATTTTTCCATAtt             <  1:1446587/76‑1 (MQ=255)                                                                   aaaTTCACATAATTTTTCCTCGAGTTCTTGATGAAATGGTCCATTATTTGTCAACCATTTATTTTTCCATAttttt          <  1:1348007/76‑1 (MQ=255)                                                                      ttCACATAATTTTTCCTCGAGTTCTTGATGAAATGGTCCATTATTTGTCAACCATTTATTTTTCCATATTTTTTcc       <  1:1238986/76‑1 (MQ=255)                                                                           aTAATTTTTCCTCGAGTTCTTGATGAAATGGTCCATTATTTGTCAACCATTTATTTTTCCATATTTTTTCCAGGTa  <  1:1274356/76‑1 (MQ=255)                                                                           |                                                                            GGCTTGTAATGCTGTGATCAATGCAATTGTTGCATTGTTAAAAAGCGAGATATGCTGCACACCCAAAAATTCACATAATTTTTCCTCGAGTTCTTGATGAAATGGTCCATTATTTGTCAACCATTTATTTTTCCATATTTTTTCCAGGTA  >  REL606/2041171‑2041320

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