New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	= 3553762	87 (1.880)	6 (0.140)	5/138	4.3	6.9%	coding (3883/4236 nt)	rhsB	rhsB element core protein RshB
	?	REL606	= 3553969	83 (1.950)					coding (4090/4236 nt)	rhsB	rhsB element core protein RshB
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GTATCAGCTTATCATACTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGTCATTTTATTTAGGTACGCCCAGTGGTAATTTTAAAGCAATAG  >  REL606/3553894‑3554044                                                                            |                                                                            gTATCAGCTTATCATACTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAGGATgaga                                                                             <  1:934473/76‑1 (MQ=255)  tATCAGCTTATCATACTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGaa                                                                            <  2:961185/76‑1 (MQ=255)  tATCAGCTTATCATACTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGaa                                                                            <  1:974042/76‑1 (MQ=255)  tATCAGCTTATCATACTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGaa                                                                            >  1:1218937/1‑76 (MQ=255)  tATCAGCTTATCATACTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGaa                                                                            >  2:519636/1‑76 (MQ=255)  tATCAGCTTATCATACTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGaa                                                                            <  2:340118/76‑1 (MQ=255)   aTCAGCTTATCATACTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGaaa                                                                           <  1:307870/76‑1 (MQ=255)    tCAGCTTATCATACTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGaaaa                                                                          <  2:665929/76‑1 (MQ=255)     cAGCTTATCATACTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAAt                                                                         <  1:151374/76‑1 (MQ=255)     cAGCTCATCATACTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAAt                                                                         <  2:1478476/76‑1 (MQ=255)      aGCTTATCATACTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATa                                                                        <  2:449840/76‑1 (MQ=255)       gCTTATCATACTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAg                                                                       >  1:242632/1‑76 (MQ=255)       gCTTATCATACTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAg                                                                       >  1:242620/1‑76 (MQ=255)        cTTATCATACTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGt                                                                      >  2:1084230/1‑76 (MQ=255)         ttATCATACTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGtt                                                                     <  2:1246307/76‑1 (MQ=255)         ttATCATACTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGtt                                                                     >  2:1161306/1‑76 (MQ=255)         ttATCATACTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGtt                                                                     >  2:1286126/1‑76 (MQ=255)         ttATCATACTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGtt                                                                     <  2:1388091/76‑1 (MQ=255)          tATCATACTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTa                                                                    >  1:653715/1‑76 (MQ=255)            tCATACTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTaaa                                                                  >  1:52837/1‑76 (MQ=255)             cATACTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTaaaa                                                                 <  1:48162/76‑1 (MQ=255)              aTACTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAg                                                                >  2:254072/1‑76 (MQ=255)              aTACTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAg                                                                >  1:948582/1‑76 (MQ=255)              aTACTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAg                                                                <  2:864517/76‑1 (MQ=255)               tACTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAgn                                                               <  1:718741/76‑2 (MQ=255)               tACTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGt                                                               >  2:715088/1‑76 (MQ=255)                aCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTa                                                              >  1:173079/1‑76 (MQ=255)                aCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTa                                                              <  2:1060300/76‑1 (MQ=255)                       gcgCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATaa                                                       >  2:395416/3‑76 (MQ=255)                       gTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATaa                                                       <  2:1118204/76‑1 (MQ=255)                       gTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATaa                                                       <  2:1210466/76‑1 (MQ=255)                       gTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATaa                                                       <  2:1024145/76‑1 (MQ=255)                       gTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATaa                                                       <  1:433707/76‑1 (MQ=255)                           aGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAat                                                   <  2:228154/76‑1 (MQ=255)                           aGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAat                                                   <  2:526690/76‑1 (MQ=255)                            gATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAata                                                  >  2:272585/1‑76 (MQ=255)                             aTAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAatat                                                 <  1:1462195/76‑1 (MQ=255)                             aTAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAatat                                                 >  1:357852/1‑76 (MQ=255)                               aGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCa                                               <  1:1092791/76‑1 (MQ=255)                               aGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCa                                               <  1:1279688/76‑1 (MQ=255)                               aGTCACGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCa                                               <  2:1460937/76‑1 (MQ=255)                                gTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCaa                                              >  2:572515/1‑76 (MQ=255)                                gTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCaa                                              <  2:652348/76‑1 (MQ=255)                                gTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCaa                                              >  2:921375/1‑76 (MQ=255)                                  cATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGt                                            >  1:462628/1‑76 (MQ=255)                                    tGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGTCa                                          <  1:504628/76‑1 (MQ=255)                                    tGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGTCa                                          >  2:655741/1‑76 (MQ=255)                                    tGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGTCa                                          <  2:114473/76‑1 (MQ=255)                                     ggAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATANCAATATCAAGTCag                                         <  1:812124/76‑2 (MQ=255)                                     ggAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGTCAt                                         <  2:1108936/76‑1 (MQ=255)                                       agaATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGTCAttt                                       >  2:1001497/1‑76 (MQ=255)                                        gaATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGTCAtttt                                      <  1:1402731/76‑1 (MQ=255)                                         aaTATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGTCATttta                                     >  1:333787/1‑76 (MQ=255)                                         aaTATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGTCATttta                                     >  2:202828/1‑76 (MQ=255)                                          atatTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGTCATtttat                                    <  2:1384947/76‑1 (MQ=255)                                          atatTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGGTAAAAGTAAAGATAACAATATCAAGTCATtttat                                    <  1:9/76‑1 (MQ=255)                                             ttGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGTCATtttattta                                 >  2:808322/1‑76 (MQ=255)                                              tGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGTCATTTTATTTAg                                <  1:320038/76‑1 (MQ=255)                                               gggACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGTCATTTTATTTAgg                               <  2:1392600/76‑1 (MQ=255)                                                ggACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGTCATTTTATTTAGGt                              <  2:1214883/76‑1 (MQ=255)                                                 gACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGTCATTTTATTTAGGTa                             <  2:342834/76‑1 (MQ=255)                                                  aCGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGTCATTTTATTTAGGTAc                            <  2:578277/76‑1 (MQ=255)                                                  aCGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGTCATTTTATTTAGGTAc                            >  1:1097402/1‑76 (MQ=255)                                                   cGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGTCATTTTATTTAGGTACg                           <  2:1403217/76‑1 (MQ=255)                                                   cGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGTCATTTTATTTAGGTACg                           <  2:459941/76‑1 (MQ=255)                                                     aaaTATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGTCATTTTATTTAGGTACGcc                         >  1:328865/1‑76 (MQ=255)                                                     aaaTATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGTCATTTTATTTAGGTACGcc                         <  2:75555/76‑1 (MQ=255)                                                      aaTATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGTCATTTTATTTAGGTACGccc                        <  2:1042428/76‑1 (MQ=255)                                                       atatTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGTCATTTTATTTAGGTACGCCCa                       <  1:570230/76‑1 (MQ=255)                                                       atatTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGTCATTTTATTTAGGTACGCCCa                       <  2:366206/76‑1 (MQ=255)                                                        tatTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGTCATTTTATTTAGGTACGCCCAg                      >  2:797869/1‑76 (MQ=255)                                                        tatTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGTCATTTTATTTAGGTACGCCCAg                      <  2:487670/76‑1 (MQ=255)                                                        tatTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGTCATTTTATTTAGGTACGCCCAg                      <  1:906700/76‑1 (MQ=255)                                                         atTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGTCATTTTATTTAGGTACGCCCAGt                     >  1:1211909/1‑76 (MQ=255)                                                         atTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGTCATTTTATTTAGGTACGCCCAGt                     >  2:788226/1‑76 (MQ=255)                                                           tttCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGTCATTTTATTTAGGTACGCCCAGTgg                   <  2:266069/76‑1 (MQ=255)                                                           tttCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGTCATTTTATTTAGGTACGCCCAGTgg                   >  1:980044/1‑76 (MQ=255)                                                            ttCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGTCATTTTATTTAGGTACGCCCAGTGGt                  >  2:85840/1‑76 (MQ=255)                                                             tCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGTCATTTTATTTAGGTACGCCCAGTGGTa                 >  1:231056/1‑76 (MQ=255)                                                             tCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGTCATTTTATTTAGGTACGCCCAGTGGTa                 <  1:61219/76‑1 (MQ=255)                                                              cAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGTCATTTTATTTAGGTACGCCCAGTGGTaa                <  2:555439/76‑1 (MQ=255)                                                              cAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGTCATTTTATTTAGGTACGCCCAGTGGTaa                >  2:815911/1‑76 (MQ=255)                                                              cAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGTCATTTTATTTAGGTACGCCCAGTGGTaa                <  2:81237/76‑1 (MQ=255)                                                                 gTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGTCATTTTATTTAGGTACGCCCAGTGGTAAttt             >  2:695354/1‑76 (MQ=255)                                                                  tAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGTCATTTTATTTAGGTACGCCCAGTGGTAAtttt            >  1:607902/1‑76 (MQ=255)                                                                  tAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGTCATTTTATTTAGGTACGCCCAGTGGTAAtttt            >  2:1284921/1‑76 (MQ=255)                                                                   aaaGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGTCATTTTATTTAGGTACGCCCAGTGGTAATTTTa           <  1:1443943/76‑1 (MQ=255)                                                                   aaaGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGTCATTTTATTTAGGTACGCCCAGTGGTAATTTTa           >  1:1118260/1‑76 (MQ=255)                                                                   aaaGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGTCATTTTATTTAGGTACGCCCAGTGGTAATTTTa           >  1:1221016/1‑76 (MQ=255)                                                                     aGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGTCATTTTATTTAGGTACGCCCAGTGGTAATTTTaaa         <  1:1325443/76‑1 (MQ=255)                                                                     aGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGTCATTTTATTTAGGTACGCCCAGTGGTAATTTTaaa         <  2:418211/76‑1 (MQ=255)                                                                        tGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGTCATTTTATTTAGGTACGCCCAGTGGTAATTTTAAAGCa      <  1:1015198/76‑1 (MQ=255)                                                                        tGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGTCATTTTATTTAGGTACGCCCAGTGGTAATTTTAAAGCa      <  1:720907/76‑1 (MQ=255)                                                                         gagaAAATAGTTAAAAGTAAAGATAACAATATCAAGTCATTTTATTTAGGTACGCCCAGTGGTAATTTTAAAGCaa     >  1:1379212/1‑76 (MQ=255)                                                                          agaAAATAGTTAAAAGTAAAGATAACAATATCAAGTCATTTTATTTAGGTACGCCCAGTGGTAATTTTAAAGCAAt    >  2:721688/1‑76 (MQ=255)                                                                            aaaaTAGTTAAAAGTAAAGATAACAATATCAAGTCATTTTATTTAGGTACGCCCAGTGGTAATTTTAAAGCAATAg  <  1:160996/76‑1 (MQ=255)                                                                            |                                                                            GTATCAGCTTATCATACTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATCAAGTCATTTTATTTAGGTACGCCCAGTGGTAATTTTAAAGCAATAG  >  REL606/3553894‑3554044

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