New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	2037389 =	30 (0.760)	3 (0.080)	3/142	4.7	7.9%	coding (742/1176 nt)	wzy	O antigen polymerase (Wzy)
	?	REL606	2037953 =	42 (1.120)					coding (178/1176 nt)	wzy	O antigen polymerase (Wzy)
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TTATATAAATAAAATTTTCTTCAAAAAAATCAGTAGATTCTCCATTATTTTTTAAATTCTGAAGTTGATTAATATTATCTGTTATAAATGCAACGTGTGAAACTGTTTGAACTCTGCCAACAACCTTCATAAATTGTTGATTAAAAGAATC  >  REL606/2037314‑2037464                                                                            |                                                                            ttATATAAATAAAATTTTCTTCAAAAAAATCAGTAGATTCTCCATTATTTTTTAAATTCTGAAGTTGATTAAtatt                                                                             <  1:1032254/76‑1 (MQ=255)   atataAATAAAATTTTCTTCAAAAAAATCAGTAGATTCTCCATTATTTTTTAAATTCTGAAGTTGATTAAtattat                                                                           >  2:39891/1‑76 (MQ=255)            aaaTTTTCTTCAAAAAAATCAGTAGATTCTCCATTATTTTTTAAATTCTGAAGTTGATTAATATTATCTGTTATaa                                                                  >  1:753520/1‑76 (MQ=255)               ttttCTTCAAAAAAATCAGTAGATTCTCCATTATTTTTTAAATTCTGAAGTTGATTAATATTATCTGTTATAAATg                                                               <  1:1026795/76‑1 (MQ=255)               ttttCTTCAAAAAAATCAGTAGATTCTCCATTATTTTTTAAATTCTGAAGTTGATTAATATTATCTGTTATAAATg                                                               <  1:1038414/76‑1 (MQ=255)               ttttCTTCAAAAAAATCAGTAGATTCTCCATTATTTTTTAAATTCTGAAGTTGATTAATATTATCTGTTATAAATg                                                               >  1:1134884/1‑76 (MQ=255)                  tcttcAAAAAAATCAGTAGATTCTCCATTATTTTTTAAATTCTGAAGTTGATTAATATTATCTGTTATAAATGCaa                                                            <  1:46084/76‑1 (MQ=255)                      cAAAAAAATCAGTAGATTCTCCATTATTTTTTAAATTCTGAAGTTGATTAATATTATCTGTTATAAATGCAACgtg                                                        >  1:212172/1‑76 (MQ=255)                       aaaaaaaTCAGTAGATTCTCCATTATTTTTTAAATTCTGAAGTTGATTAATATTATCTGTTATAAATGCAACgtgt                                                       >  2:771876/1‑76 (MQ=255)                           aaaTCAGTAGATTCTCCATTATTTTTTAAATTCTGAAGTTGATTAATATTATCTGTTATAAATGCAACGTGTGaaa                                                   >  1:1002972/1‑76 (MQ=255)                           aaaTCAGTAGATTCTCCATTATTTTTTAAATTCTGAAGTTGATTAATATTATCTGTTATAAATGCAACGTGTGaaa                                                   <  2:212172/76‑1 (MQ=255)                           aaaTCAGTAGATTCTCCATTATTTTTTAAATTCTGAAGTTGATTAATATTATCTGTTATAAATGCAACGTGTGaaa                                                   >  1:577094/1‑76 (MQ=255)                            aaTCAGTAGATTCTCCATTATTTTTTAAATTCTGAAGTTGATTAATATTATCTGTTATAAATGCAACGTGTGAAAc                                                  <  2:997493/76‑1 (MQ=255)                                  tAGATTCTCCATTATTTTTTAAATTCTGAAGTTGATTAATATTATCTGTTATAAATGCAACGTGTGAAACTGTTTg                                            >  1:244411/1‑76 (MQ=255)                                  aaGATTCTCCATTATTTTTTAAATTCTGAAGTTGATTAATATTATCTGTTATAAATGCAACGTGTGAAACTGTTTg                                            >  1:699561/2‑76 (MQ=255)                                   aGATTCTCCATTATTTTTTAAATTCTGAAGTTGATTAATATTATCTGTTATAAATGCAACGTGTGAAACTGTTTGa                                           >  1:992494/1‑76 (MQ=255)                                   aGATTCTCCATTATTTTTTAAATTCTGAAGTTGATTAATATTATCTGTTATAAATGCAACGTGTGAAACTGTTTGa                                           >  1:77013/1‑76 (MQ=255)                                   aGATTCTCCATTANTTTTTAAATTCTGAAGTTGATTAATATTATCTGTTATAAATGCAACGTGTGAAACTGTTTGa                                           >  1:850264/1‑76 (MQ=255)                                     aTTCTCCATTATTTTTTAAATTCTGAAGTTGATTAATATTATCTATTATAAATGCAACGTGTGAAACTGTTTGAAc                                         >  1:238591/1‑76 (MQ=255)                                              tattTTTTAAATTCTGAAGTTGATTAATATTATCTGTTATAAATGCAACGTGTGAAACTGTTTGAACTCTGCcaac                                >  2:693636/1‑76 (MQ=255)                                                   tttAAATTCTGAAGTTGATTAATATTATCTGTTATAAATGCAACGTGTGAAACTGTTTGAACTCTGCCAACAACCt                           >  1:501376/1‑76 (MQ=255)                                                           cTGAAGTTGATTAATATTATCTGTTATAAATGCAACGTGTGAAACTGTTTGAACTCTGCCAACAACCTTCATAAAt                   >  1:1266041/1‑76 (MQ=255)                                                            tGAAGTTGATTAATATTATCTGTTATAAATGCAACGTGTGAAACTGTTTGAACTCTGCCAACAACCTTCATAAAtt                  >  1:957907/1‑76 (MQ=255)                                                               aGTTGATTAATATTATCTGTTATAAATGCAACGTGTGAAACTGTTTGAACTCTGCCAACAACCTTCATAAAttgtt               <  1:886122/76‑1 (MQ=255)                                                                gTTGATTAATATTATCTGTTATAAATGCAACGTGTGAAACTGTTTGAACTCTGCCAACAACCTTCATAAAttgttg              >  1:933170/1‑76 (MQ=255)                                                                gTTGATTAATATTATCTGTTATAAATGCAACGTGTGAAACTGTTTGAACTCTGCCAACAACCTTCATAAAttgttg              <  2:536772/76‑1 (MQ=255)                                                                gTTGATTAATATTATCTGTTATAAATGCAACGTGTGAAACTGTTTGAACTCTGCCAACAACCTTCATAAAttgttg              <  2:564001/76‑1 (MQ=255)                                                                 ttGATTAATATTATCTGTTATAAATGCAACGTGTGAAACTGTTTGAACTCTGCCAACAACCTTCATAAATTGTTGa             >  1:1207517/1‑76 (MQ=255)                                                                     ttAATATTATCTGTTATAAATGCAACGTGTGAAACTGTTTGAACTCTGCCAACAACCTTCATAAATTGTTGATTaa         <  1:1238815/76‑1 (MQ=255)                                                                         tattatCTGTTATAAATGCAACGTGTGAAACTGTTTGAACTCTGCCAACAACCTTCATAAATTGTTGATTAAAAGa     >  1:1141457/1‑76 (MQ=255)                                                                            tatCTGTTATAAATGCAACGTGTGAAACTGTTTGAACTCTGCCAACAACCTTCATAAATTGTTGATTAAAAGAATc  >  2:743583/1‑76 (MQ=255)                                                                            |                                                                            TTATATAAATAAAATTTTCTTCAAAAAAATCAGTAGATTCTCCATTATTTTTTAAATTCTGAAGTTGATTAATATTATCTGTTATAAATGCAACGTGTGAAACTGTTTGAACTCTGCCAACAACCTTCATAAATTGTTGATTAAAAGAATC  >  REL606/2037314‑2037464

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