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.

TGTGTACTTTAATAAAATAAATACTAAAAGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAACACTAACCAAAACAACCATCTCATATTCATTGCATATAGGTCATTCGCAATAATTTTAGATAG  >  REL606/2036985‑2037125                                                                            |                                                                  tgtgtACTTTAATAAAATAAATACTAAAAGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATaaa                                                                   >  1:1164235/1‑76 (MQ=255) tgtgtACTTTAATAAAATAAATACTAAAAGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATaaa                                                                   >  1:238999/1‑76 (MQ=255)  gtgtACTTTAATAAAATAAATACTAAAAGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATaaaa                                                                  <  1:1115729/76‑1 (MQ=255)    gtACTTTAATAAACTAAATACTAAAAGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAac                                                                >  1:415472/1‑76 (MQ=255)      aCTTTAATAAAATAAATACTAAAAGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAacac                                                              <  1:782574/76‑1 (MQ=255)       cTTTAATAAAATAAATACTAAAAGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAACACt                                                             >  2:302763/1‑76 (MQ=255)               aaataaatACTAAAAGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAACACTAACCaaaa                                                     <  1:351316/76‑1 (MQ=255)                 ccaaatcCTAAAAGCGAAAATCTAACTGCTGTGAATTGGAATACCCACCCCGGCACTACCCACACTAACCCAaacc                                                   >  2:1211350/3‑75 (MQ=255)                 ataaataCTAAAAGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAACACTAACCAAaaca                                                   <  2:685550/76‑1 (MQ=255)                  taaataCTAAAAGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAACACTAACCAAaacaa                                                  <  1:345225/76‑1 (MQ=255)                  taaataCTAAAAGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAACACTAACCAAaacaa                                                  <  1:914582/76‑1 (MQ=255)                   aaataCTAAAAGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAACACTAACCAAaacaac                                                 <  1:218199/76‑1 (MQ=255)                     ataCTAAAAGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAACACTAACCAAAACAACCa                                               >  2:494001/1‑76 (MQ=255)                      taCTAAAAGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAACACTAACCAAAACAACCAt                                              <  1:922596/76‑1 (MQ=255)                       aCTAAAAGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAACACTAACCAAAACAACCAtc                                             >  2:598854/1‑76 (MQ=255)                        cTAAAAGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAACACTAACCAAAACAACCAtct                                            <  1:966842/76‑1 (MQ=255)                        cTAAAAGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAACACTAACCAAAACAACCAtct                                            <  1:966837/76‑1 (MQ=255)                          aaaaGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAACACTAACCAAAACAACCATCTCa                                          >  1:1152361/1‑76 (MQ=255)                                  aaaTATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAACACTAACCAAAACAACCATCTCATATTCAtt                                  <  2:107769/76‑1 (MQ=255)                                    atatAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAACACTAACCAAAACAACCATCTCATATTCATTGc                                >  1:1104893/1‑76 (MQ=255)                                      ataACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAACACTAACCAAAACAACCATCTCATATTCATTGCat                              >  1:60273/1‑76 (MQ=255)                                      ataACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAACACTAACCAAAACAACCATCTCATATTCATTGCat                              <  2:275413/76‑1 (MQ=255)                                       taACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAACACTAACCAAAACAACCATCTCATATTCATTGCata                             >  2:1004486/1‑76 (MQ=255)                                        aaCTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAACACTAACCAAAACAACCATCTCATATTCATTGCatat                            <  2:761325/76‑1 (MQ=255)                                        aaCTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAACACTAACCAAAACAACCATCTCATATTCATTGCatat                            <  2:387816/76‑1 (MQ=255)                                         aCTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAACACTAACCAAAACAACCATCTCATATTCATTGCAtata                           <  2:1063024/76‑1 (MQ=255)                                            gctgTGAATTGGAATATCCAACCCGGAAATAAAAACACTAACCAAAACAACCATCTCATATTCATTGCATATAGGt                        <  2:135013/76‑1 (MQ=255)                                             ctgTGAATTGGAATATCCAACCCGGAAATAAAAACACTAACCAAAACAACCATCTCATATTCATTGCATATAGGTc                       <  1:24413/76‑1 (MQ=255)                                              tgtgAATTGGAATATCCAACCCGGAAATAAAAACACTAACCAAAACAACCATCTCATATTCATTGCATATAGGTCa                      >  2:600939/1‑76 (MQ=255)                                                     tGGAATATCCAACCCGGAAATAAAAACACTAACCAAAACAACCATCTCATATTCATTGCATATAGGTCATTCGCaa               >  2:62183/1‑76 (MQ=255)                                                         atatCCAACCCGGAAATAAAAACACTAACCAAAACAACCATCTCATATTCATTGCATATAGGTCATTCGCaataat           >  2:657284/1‑76 (MQ=255)                                                            tCCAACCCGGAAATAAAAACACTAACCAAAACAACCATCTCATATTCATTGCATATAGGTCATTCGCAATAAtttt        <  2:612286/76‑1 (MQ=255)                                                             ccAACCCGGAAATNAAAACACTAACCAAAACAACCATCTCATATTCATTGCATATAGGTCATTCGCAATAATTTTa       >  1:858108/1‑76 (MQ=255)                                                              cAACCCGGAAATAAAAACACTAACCAAAACAACCATCTCATATTCATTGCATATAGGTCATTCGCAATAATTTTAg      <  2:50912/76‑1 (MQ=255)                                                                aCCCGGAAATAAAAACACTAACCAAAACAACCATCTCATATTCATTGCATATAGGTCATTCGCAATAATTTTAGAt    <  1:1142108/76‑1 (MQ=255)                                                                  ccGGAAATAAAAACACTAACCAAAACAACCATCTCATATTCATTGCATATAGGTCATTCGCAATAATTTTAGATAg  <  1:672520/76‑1 (MQ=255)                                                                            |                                                                  TGTGTACTTTAATAAAATAAATACTAAAAGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAACACTAACCAAAACAACCATCTCATATTCATTGCATATAGGTCATTCGCAATAATTTTAGATAG  >  REL606/2036985‑2037125

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