New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	2039235 =	26 (0.540)	6 (0.130)	6/144	4.6	16.4%	coding (460/747 nt)	wbbA	putative glycosyltransferase WbbA
	?	REL606	2039555 =	36 (0.790)					coding (140/747 nt)	wbbA	putative glycosyltransferase WbbA
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TCGATATTTTCAAAGTGTTGATATAACGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTTATATTTAGATTCTTTTATTGTATTTATAAGATTTATAAAACAATTTTCATGACTAAACCACT  >  REL606/2039480‑2039630                                                                            |                                                                            tCGATATTTTCAAAGTGTTGATATAACGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAATg                                                                             <  1:1127290/76‑1 (MQ=255)    atatTTTCAAAGTGTTGATATAACGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAATGCTa                                                                          >  1:1269872/1‑76 (MQ=255)         ttCAAAGTTTTGATATAACGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAATGCTATTTgg                                                                     <  1:768749/76‑1 (MQ=255)          tCAAAGTGTTGATATAACGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAATGCTATTTggn                                                                    <  1:1011726/76‑2 (MQ=255)          tCAAAGTGTTGATATAACGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGt                                                                    <  1:875504/76‑1 (MQ=255)          tCAAAGTGTTGATATAACGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGt                                                                    <  1:578807/76‑1 (MQ=255)             aaGTGTTGATATAACGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCtt                                                                 <  1:1081214/76‑1 (MQ=255)                 gtTGATATAACGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTtata                                                             <  1:1192106/76‑1 (MQ=255)                 gtTGATATAACGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTtata                                                             <  1:1192119/76‑1 (MQ=255)                 gtTGATATAACGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTtata                                                             >  1:359267/1‑76 (MQ=255)                  ttGATATAACGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTTatat                                                            <  1:873417/76‑1 (MQ=255)                    gATATAACGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTTATAttt                                                          <  1:373338/76‑1 (MQ=255)                    gATATAACGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTTATAttt                                                          <  1:1523052/76‑1 (MQ=255)                     atatAACGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTTATATTTa                                                         >  2:548784/1‑76 (MQ=255)                           cGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTTATATTTAGActct                                                   <  1:572869/76‑1 (MQ=255)                           cGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTTATATTTAGATTCt                                                   <  1:572895/76‑1 (MQ=255)                            ggATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTTATATTTAGATTCtt                                                  <  1:679926/76‑1 (MQ=255)                            ggATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTTATATTTAGATTCtt                                                  <  1:679934/76‑1 (MQ=255)                               tCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTTATATTTAGATTCTTTTa                                               <  1:343242/76‑1 (MQ=255)                                 tAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTTATATTTAGATTCTTTTAtt                                             <  1:708601/76‑1 (MQ=255)                                   aGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTTATATTTAGATTCTTTTATTGt                                           <  1:830279/76‑1 (MQ=255)                                       aGCAACGTTTAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTTATATTTAGATTCTTTTATTGTAttt                                       <  1:927905/76‑1 (MQ=255)                                          aaCGTTTAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTTATATTTAGATTCTTTTATTGTATTtata                                    <  1:539238/76‑1 (MQ=255)                                              tttAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTTATATTTAGATTCTTTTATTGTATTTATAAGAt                                <  1:525484/76‑1 (MQ=255)                                                tAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTTATATTTAGATTCTTTTATTGTATTTATAAGAttt                              <  1:1089545/76‑1 (MQ=255)                                                tAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTTATATTTAGATTCTTTTATTGTATTTATAAGAttt                              <  1:1089558/76‑1 (MQ=255)                                                 aaTGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTTATATTTAGATTCTTTTATTGTATTTATAAGATTta                             <  1:1476057/76‑1 (MQ=255)                                                  aTGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTTATATTTAGATTCTTTTATTGTATTTATAAGATTtat                            <  2:1102673/76‑1 (MQ=255)                                                     aTCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTTATATTTAGATTCTTTTATTGTATTTATAAGATTTATaaa                         <  1:280419/76‑1 (MQ=255)                                                     aTCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTTATATTTAGATTCGTTTATTGTATTTATAAGATTTATaaa                         <  1:280416/76‑1 (MQ=255)                                                      tCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTTATATTTAGATTCTTTTATTGTATTTATAAGATTTATaaaa                        >  1:706458/1‑76 (MQ=255)                                                        cATCAAAAATAAAAGTAAATGCTATTTGGTCTTTATATTTAGATTCTTTTATTGTATTTATAAGATTTATAAAACa                      <  1:266597/76‑1 (MQ=255)                                                                 tAAAAGTAAATGCTATTTGGTCTTTATATTTAGATTCTTTTATTGTATTTATAAGATTTATAAAACAATTTTCATg             <  1:581630/76‑1 (MQ=255)                                                                      gTAAATGCTATTTGGTCTTTATATTTAGATTCTTTTATTGTATTTATAAGATTTATAAAACAATTTTCATGACTaa        <  1:107862/76‑1 (MQ=255)                                                                       tAAATGCTATTTGGTCTTTATATTTAGATTCTTTTATTGTATTTATAAGATTTATAAAACAATTTTCATGACTaaa       >  1:1434846/1‑76 (MQ=255)                                                                        aaaTGCTATTTGGTCTTTATATTTAGATTCTTTTATTGTATTTATAAGATTTATAAAACAATTTTCATGACTAAAc      >  1:88203/1‑76 (MQ=255)                                                                            gCTATTTGGTCTTTATATTTAGATTCTTTTATTGTATTTATAAGATTTATAAAACAATTTTCATGACTAAACCACt  >  1:505887/1‑76 (MQ=255)                                                                            |                                                                            TCGATATTTTCAAAGTGTTGATATAACGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAATGCTATTTGGTCTTTATATTTAGATTCTTTTATTGTATTTATAAGATTTATAAAACAATTTTCATGACTAAACCACT  >  REL606/2039480‑2039630

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