New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 54095278 (1.440)7 (0.150) 3/132 6.1 8.6% coding (161/1527 nt) ybcK predicted recombinase
?REL606 = 541418 80 (1.680)coding (627/1527 nt) ybcK predicted recombinase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TTTGTTCCAGACCCCGACAGGGTGAAAACTATTGAGCTAATTTTTAAACTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGACC  >  REL606/541350‑541492
                                                                    |                                                                          
tttGTTCCAGACCCCGACAGGGTGAAAACTATTGAGCTAATTTTTAAACTCAGGATGGAAAGGCGCTCATTGAATg                                                                     <  2:32458/76‑1 (MQ=255)
 ttGTTCCAGACCCCGACAGGGTGAAAACTATTGAGCTAATTTTTAAACTCAGGATGGAAAGGCGCTCATTGAATGc                                                                    <  1:864515/76‑1 (MQ=255)
   gTTCCAGACCCCGACAGGGTGAAAACTATTGAGCTAATTTTTAAACTCAGGATGGAAAGGCGCTCATTGAATGCaa                                                                  <  2:709618/76‑1 (MQ=255)
     tCCAGACCCCGACAGGGTGAAAACTATTGAGCTAATTTTTAAACTCAGGATGGAAAGGCGCTCATTGAATGCAATa                                                                >  2:1437215/1‑76 (MQ=255)
       cAGCCCCCGACAGGGTGAAAACTATTGAGCTAATTTTTTAACTCAGGATGGAAAGGCGCTCATTGAATGCAATAGc                                                              <  1:467268/76‑1 (MQ=255)
       cAGACCCCGACAGGGTGAAAACTATTGAGCTAATTTTTAAACTCAGGATGGAAAGGCGCTCATTGAATGCAATAGc                                                              >  1:534730/1‑76 (MQ=255)
       cAGACCCCGACAGGGTGAAAACTATTGAGCTAATTTTTAAACTCAGGATGGAAAGGCGCTCATTGAATGCAATAGc                                                              >  2:23552/1‑76 (MQ=255)
        aGACCCCGACAGTGTGAAAACTATTGAGCTAATTTTTAAACTCAGGATGGAAAGGCGCTCATTGAATGCAATAGcc                                                             >  1:1505054/1‑76 (MQ=255)
        aGACCCCGACAGTGTGAAAACTATTGAGCTAATTTTTAAACTCAGGATGGAAAGGCGCTCATTGAATGCAATAGcc                                                             >  1:1505042/1‑76 (MQ=255)
         gACCCCGACAGGGTGAAAACTATTGAGCTAATTTTTAAACTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCa                                                            <  2:1146502/76‑1 (MQ=255)
          aCCCCGACAGGGTGAAAACTATTGAGCTAATTTTTAAACTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCaa                                                           <  1:1547061/76‑1 (MQ=255)
           ccccGACAGGGTGAAAACTATTGAGCTAATTTTTAAACTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAg                                                          <  2:688649/76‑1 (MQ=255)
           ccccGACAGGGTGAAAACTATTGAGCTAATTTTTAAACTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAg                                                          <  2:1447397/76‑1 (MQ=255)
            cccGACAGGGTGAAAACTATTGAGCTAATTTTTAAACTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGt                                                         >  2:328902/1‑76 (MQ=255)
            cccGACAGGGTGAAAACTATTGAGCTAATTTTTAAACTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGt                                                         >  2:328916/1‑76 (MQ=255)
            cccGACAGGGTGAAAACTATTGAGCTAATTTTTAAACTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGt                                                         <  2:701182/76‑1 (MQ=255)
            cccGACAGGGTGAAAACTATTGAGCTAATTTTTAAACTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGt                                                         <  2:668067/76‑1 (MQ=255)
                 cAGGGTGAAAACTATTGAGCTAATTTTTAGACTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTa                                                    <  2:293791/76‑1 (MQ=255)
                 cAGGGTGAAAACTATTGAGCTAATTTTTAAACTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTa                                                    <  2:121616/76‑1 (MQ=255)
                 cAGGGTGAAAACTATTGAGCTAATTTTTAAACTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTa                                                    <  2:442135/76‑1 (MQ=255)
                 cAGGGTGAAAACTATTGAGCTAATTTTTAAACTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTa                                                    >  1:574542/1‑76 (MQ=255)
                 cAGGGTGAAAACTATTGAGCTAATTTTTAAACTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTa                                                    <  1:577391/76‑1 (MQ=255)
                  aGGGTGAAAACTATTGAGCTAATTTTTAAACTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTaa                                                   >  2:625333/1‑76 (MQ=255)
                  aGGGTGAAAACTATTGAGCTAATTTTTAAACTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTaa                                                   >  2:834149/1‑76 (MQ=255)
                  aGGGTGAAAACTATTGAGCTAATTTTTAAACTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTaa                                                   >  1:1625359/1‑76 (MQ=255)
                  aGGGTGAAAACTATTGAGCTAATTTTTAAACTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTaa                                                   <  1:1684289/76‑1 (MQ=255)
                  aGGGTGAAAACTATTGAGCTAATTTTTAAACTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTaa                                                   <  2:25419/76‑1 (MQ=255)
                        aaaaCTATTGAGCTAATTTTTAAACTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATc                                             <  1:1412751/76‑1 (MQ=255)
                         aaaCTATTGAGCTAATTTTTAAACTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCa                                            <  1:321181/76‑1 (MQ=255)
                            cTATTGAGCTAATTTTTAAACTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGc                                         <  2:156268/76‑1 (MQ=255)
                            cTATTGAGCTAATTTTTAAACTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGc                                         <  2:1155235/76‑1 (MQ=255)
                             tATTGAGCTAATTTTTAAACTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCt                                        >  1:1711771/1‑76 (MQ=255)
                              aTTGAGCTAATTTTTAAACTCAGGATGGAAAGGCGCTGATTGAATGCAATAGCCNAGTATTTAAATGATCATGCTg                                       >  2:273039/1‑76 (MQ=255)
                              aTTGAGCTAATTTTTAAACTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTg                                       >  2:1651595/1‑76 (MQ=255)
                              aTTGAGCTAATTTTTAAACTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTg                                       <  2:721156/76‑1 (MQ=255)
                               ttGAGCTAATTTTTAAACTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGt                                      <  1:1228866/76‑1 (MQ=255)
                               ttGAGCTAATTTTTAAACTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGt                                      >  1:16025/1‑76 (MQ=255)
                                tGAGCTAATTTTTAAACTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTa                                     <  2:217883/76‑1 (MQ=255)
                                 gAGCTAATTTTTAAACTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTaa                                    >  2:61533/1‑76 (MQ=255)
                                 gAGCTAATTTTTAAACTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTaa                                    >  2:74853/1‑76 (MQ=255)
                                   gCTAATTTTTAAACTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAg                                  <  1:1152702/76‑1 (MQ=255)
                                     tAATTTTTAAACTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGaa                                <  1:398670/76‑1 (MQ=255)
                                      aaTTTTTAAACTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAAt                               >  2:1413530/1‑76 (MQ=255)
                                      aaTTTTTAAACTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAAt                               >  1:184097/1‑76 (MQ=255)
                                       aTTTTTAAACTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAAtt                              >  2:343254/1‑76 (MQ=255)
                                        tttttAAACTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAAttt                             >  1:1539667/1‑76 (MQ=255)
                                        tttttAAACTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAAttt                             >  1:1539654/1‑76 (MQ=255)
                                           ttAAACTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTtctc                          >  2:1627044/1‑76 (MQ=255)
                                            tAAACTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCa                         <  1:1030062/76‑1 (MQ=255)
                                            tAAACTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCa                         >  2:206637/1‑76 (MQ=255)
                                             aaaCTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAg                        >  2:532948/1‑76 (MQ=255)
                                             aaaCTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAg                        >  2:1068546/1‑76 (MQ=255)
                                             aaaCTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAg                        >  2:1520308/1‑76 (MQ=255)
                                             aaaCTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAg                        >  2:1520291/1‑76 (MQ=255)
                                             aaaCTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAg                        >  2:1503372/1‑76 (MQ=255)
                                                cTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGaa                     <  1:911170/76‑1 (MQ=255)
                                                cTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGaa                     >  1:884433/1‑76 (MQ=255)
                                                cTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGaa                     >  1:739517/1‑76 (MQ=255)
                                                  cAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGaaaa                   >  2:743234/1‑76 (MQ=255)
                                                  cAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGaaaa                   >  2:100013/1‑76 (MQ=255)
                                                   aGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAaaag                  <  2:1544655/76‑1 (MQ=255)
                                                     gATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAaaagaa                <  1:435038/76‑1 (MQ=255)
                                                      aTGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAaaagaaa               >  1:1234692/1‑76 (MQ=255)
                                                       tGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAaaagaaag              >  1:356903/1‑76 (MQ=255)
                                                       tGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAaaagaaag              >  1:36264/1‑76 (MQ=255)
                                                       tGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAaaagaaag              >  1:365729/1‑76 (MQ=255)
                                                       tGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAaaagaaag              >  2:875176/1‑76 (MQ=255)
                                                       tGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAaaagaaag              >  2:508754/1‑76 (MQ=255)
                                                       tGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAaaagaaag              >  1:920905/1‑76 (MQ=255)
                                                       tGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAaaagaaag              >  1:921528/1‑76 (MQ=255)
                                                         gAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTg            <  1:87686/76‑1 (MQ=255)
                                                         gAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTg            >  2:17050/1‑76 (MQ=255)
                                                          aaaGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGc           >  2:605381/1‑76 (MQ=255)
                                                          aaaGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGc           <  2:1331626/76‑1 (MQ=255)
                                                           aaGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCa          >  2:697119/1‑76 (MQ=255)
                                                           aaGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCa          >  2:697122/1‑76 (MQ=255)
                                                           aaGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCa          >  2:1493153/1‑76 (MQ=255)
                                                             ggCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATg        <  2:632269/76‑1 (MQ=255)
                                                               cgcTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATGTCTCAGGAAAAGAAAGTGCATggg      <  2:1211514/76‑1 (MQ=255)
                                                                  tCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGAc   >  2:513252/1‑76 (MQ=255)
                                                                  tCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGAc   >  1:921367/1‑76 (MQ=255)
                                                                   cATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGAcc  >  1:81364/1‑76 (MQ=255)
                                                                    |                                                                          
TTTGTTCCAGACCCCGACAGGGTGAAAACTATTGAGCTAATTTTTAAACTCAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGACC  >  REL606/541350‑541492

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