New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 151901656 (2.530)6 (0.300) 4/138 2.3 11.5% coding (850/957 nt) yddK hypothetical protein
?REL606 = 1519475 41 (2.020)coding (391/957 nt) yddK hypothetical protein
Rejected: Frequency below/above cutoff threshold.

TGTAATACTTACATTTCTGCCTTGAGGAAAATTTATTGACTCAAGAGCATTACTATTCAGACTTAGCTTTACGAGTCTATCGTAGTTGGTTGCATCAAAATGTGTGAGTTGGTTATTATCAAGATTTATAGATGTGAAGTTTTGCAG  >  REL606/1519401‑1519547
                                                                          |                                                                        
tgtAATACTTACATTTCTGCCTTGAGGAAAATTTATTGACTCAAGAGCATTACTATTCAGACTTAGCTTTACGAGt                                                                         <  1:576013‑M2/76‑2 (MQ=255)
tgtAATACTTACATTTCTGCCTTGAGGAAAATTTATTGACTCAAGAGCATTACTATTCAGACTTAGCTTTACGAGt                                                                         <  1:296774‑M2/76‑2 (MQ=255)
  tAATACTTACATTTCTGCCTTGAGGAAAATTTATTGACTCAAGAGCATTACTATTCAGACTTAGCTTTACGAGtct                                                                       >  2:200590‑M2/1‑73 (MQ=255)
        ttACATTTCTGCCTTGAGGCAAATTTATTGACTCAAGAGCATTACTATTCAGACTTAGCTTTACGAGTCTATCGTa                                                                 >  1:255120/1‑76 (MQ=255)
         tACATTTCTGCCTTGAGGAAAATTTATTGACTCAAGAGCATTACTATTCAGACTTAGCTTTACGAGTCTATCGTAg                                                                >  1:610549/1‑76 (MQ=255)
         tACATTTCTGCCTTGAGGAAAATTTATTGACTCAAGAGCATTACTATTCAGACTTAGCTTTACGAGTCTATCGTAg                                                                >  2:260831/1‑76 (MQ=255)
          aCATTTCTGCCTTGAGGAAAATTTATTGACTCAAGAGCATTACTATTCAGACTTAGCTTTACGAGTCTATCGTAgt                                                               >  1:256445/1‑76 (MQ=255)
          aCATTTCTGCCTTGAGGAAAATTTATTGACTCAAGAGCATTACTATTCAGACTTAGCTTTACGAGTCTATCGTAgt                                                               >  1:235754/1‑76 (MQ=255)
            ctgtCTGCCTTGAGGAAAATTTATTGACTCAAGAGCATTCCTATCCAGACTTAGCTTTACGAGTCTATCGTAgttg                                                             <  2:403860/73‑1 (MQ=255)
            aTTTCTGCCTTGAGGAAAATTTATTGACTCAAGAGCATTACTATTCAGACTTAGCTTTACGAGTCTATCGTAgttg                                                             <  1:677944/76‑1 (MQ=255)
               tCTGCCTTGAGGAAAATTTATTGACTCAAGAGCATTACTATTCAGACTTAGCTTTACGAGTCTATCGTAgttggtt                                                          <  2:419792/76‑1 (MQ=255)
                cTGCCTTGAGGAAAATTTATTGACTCAAGAGCATTACTATTCAGACTTAGCTTTACGAGTCTATCGTAgttggttg                                                         <  2:370003/76‑1 (MQ=255)
                 tGCCTTGAGGAAAATTTATTGACTCAAGAGCATTACTATTCAGACTTAGCTTTACGAGTCTATCGTAGTTGGTTGc                                                        >  2:695191/1‑76 (MQ=255)
                 tGCCTTGAGGAAAATTTATTGACTCAAGAGCATTACTATTCAGACTTAGCTTTACGAGTCTATCGTAGTTGGTTGc                                                        <  2:44470/76‑1 (MQ=255)
                  gCCTTGAGGAAAATTTATTGACTCAAGAGCATTACTATTCAGACTTAGCTTTACGAGTCTATCGTAGTTGGTTGCa                                                       >  2:197495/1‑76 (MQ=255)
                   ccTTGAGGAAAATTTATTGACTCAAGAGCATTACTATTCAGACTTAGCTTTACGAGtctatcaaaagtttaagact                                                      >  2:247218‑M2/1‑56 (MQ=255)
                   ccTTGAGGAAAATTTATTGACTCAAGAGCATTACTATTCAGACTTAGCTTTACGAGtctatcaaaagtttaagact                                                      >  2:247213‑M2/1‑56 (MQ=255)
                     ttGAGGAAAATTTATTGACTCAAGAGCATTACTATTCAGACTTAGCTTTACGAGTCTATCGTAGTTGGTTGCATCa                                                    >  1:280374/1‑76 (MQ=255)
                      tGAGGAAAATTTATTGACTCAAGAGCATTACTATTCAGACTTAGCTTTACGAGTCTATCGTAGTTGGTTGCATCaa                                                   <  1:144588/76‑1 (MQ=255)
                       gAGGAAAATTTATTGACTCAAGAGCATTACTATTCAGACTTAGCTTTACGAGTCTATCGTAGTTGGTTGCATCaaa                                                  <  2:30105/76‑1 (MQ=255)
                        aGGAAAATTTATTGACTCAAGAGCATTACTATTCAGACTTAGCTTTACGAGTCTATCGTAGTTGGTTGCATCaaaa                                                 <  1:134608/76‑1 (MQ=255)
                         ggAAAATTTATTGACTCAAGAGCATTACTATTCAGACTTAGCTTTACGAGTCTATCGTAGTTGGTTGCATCAAAAt                                                <  2:661404/76‑1 (MQ=255)
                              aTTTATTGACTCAAGAGCATTACTATTCAGACTTAGCTTTACGAGTCTATCGTAGTTGGTTGCATCAAAAtgtgtg                                           <  2:558438/76‑1 (MQ=255)
                                 tATTGACTCAAGAGCATTACTATTCAGACTTAGCTTTACGAGTCTATCGTAGTTGGTTGCATCAAAATGTGTGAGt                                        >  1:683430/1‑76 (MQ=255)
                                   ttGACTGAATAGCATTACTATTCAGACTTAGCTTTACGAGtctatcaaaagtttaagactggacaataataaaata                                      >  2:387281‑M2/1‑40 (MQ=255)
                                   ttGACTCAAGAGCATTACTATTCAGATTTAGCTTTACGAGtctatcaaaagtttaagactggacaataataaaata                                      >  2:387294‑M2/1‑40 (MQ=255)
                                     gACTCAAGAGCATTACTATTCAGACTTAGCTTTACGAGtctatcaaaagtttaagactggacaataataaaataac                                    >  2:50845‑M2/1‑38 (MQ=255)
                                     gACTCAAGAGCATTACTATTCAGACTTAGCTTTACGAGTCTATCGTAGTTGGTTGCATCAAAATGTGTGAGTTGGt                                    >  2:703734/1‑76 (MQ=255)
                                      aCTCAAGAGCATTACTATTCAGACTTAGCTTTACGAGTCTATCGTAGTTGGTTGCATCAAAATGTGTGAGTTGGtt                                   <  2:232366/76‑1 (MQ=255)
                                       cgcAAGAGCGTTACTATTCAGGCTTAGCTTTACGAGTCTATCGTAGTTGGTTGCATCAAAATGTGTGAGTTGGtta                                  <  1:711086/74‑1 (MQ=255)
                                         cAAGAGCATTACTATTCAGACTTAGCTTTACGAGTCTATCGTAGTTGGTTGCATCAAAATGTGTGAGTTGGttatt                                >  1:489332/1‑76 (MQ=255)
                                             agCATTACTATTCAGACTTAGCTTTACGAGTCTATCGTAGTTGGTTGCATCAAAATGTGTGAGTTGGTTATTATCa                            >  1:554944/1‑76 (MQ=255)
                                             agCATTACTATTCAGACTTAGCTTTACGAGTCTATCGTAGTTGGTTGCATCAAAATGTGTGAGTTGGTTATTATCa                            <  2:288307/76‑1 (MQ=255)
                                              gCATTACTATTCAGACTTAGCTTTACGAGTCTATCGTAGTTGGTTGCATCAAAATGTGTGAGTTGGTTATTATCaa                           >  1:391662/1‑76 (MQ=255)
                                               cATTACTATTCAGACTTAGCTTTACGAGTCTATCGTAGTTGGTTGCATCAAAATGTGTGAGTTGGTTATTATCAAg                          >  2:667561/1‑76 (MQ=255)
                                                     tATTCAGACTTAGCTTTACGAGTCTATCGTAGTTGGTTGCATCAAAATGTGTGAGTTGGTTATTATCAAGATTtat                    >  1:209629/1‑76 (MQ=255)
                                                     caTTCCGGCTTAGCGTTTTCAGTCTATCCTAGTTGGTTGCATCGAAATGTGTGAATTGGTTATTATCGATATTtag                    <  2:33536/75‑2 (MQ=255)
                                                      aTTCAGACTTAGCTTTACGAGTCTATCGTAGTTGGTTGCATCAAAATGTGTGAGTTGGTTATTATCAAGATTtata                   <  1:703734/76‑1 (MQ=255)
                                                      aTTCAGACTTAGCTTTACGAGTCTATCGTAGTTGGTTGCATCAAAATGTGTGAGTTGGTTATTATCAAGATTtata                   <  2:256445/76‑1 (MQ=255)
                                                         cAGACTTAGCTTTACGAGTCTATCGTAGTTGGTTGCATCAAAATGTGTGAGTTGGTTATTATCAAGATTTATAGAt                >  1:317543/1‑76 (MQ=255)
                                                          aGACTTAGCTTTACGAGTCTATCGTAGTTGGTTGCATCAAAATGTGTGAGTTGGTTATTATCAAGATTTATAGAtg               >  2:447970/1‑76 (MQ=255)
                                                            ctttagTTTTAGGAGtctatcaaaattttaagcctggacaataataaaataacccacattgatactaataatacat             >  2:271243‑M2/3‑15 (MQ=255)
                                                            aCTTAGCTTTACGAGTCTATCGTAGTTGGTTGCATCAAAATGTGTGAGTTGGTTATTATCAAGATTTATAGAtgtg             >  1:71633/1‑76 (MQ=255)
                                                              ttAGCTTTACGAGTCTATCGTAGTTGGTTGCATCAAAATGTGTGAGTTGGTTATTATCAAGATTTATAGATGTGaa           >  2:441400/1‑76 (MQ=255)
                                                                aGCTTTACGAGTCTATCGTAGTTGGTTGCATCAAAATGTGTGAGTTGGTTATTATCAAGATTTATAGATGTGAAGt         >  1:662142/1‑76 (MQ=255)
                                                                 gCTTTACGAGTCTATCGTAGTTGGTTGCATCAAAATGTGTGAGTTGGTTATTATCAAGATTTATAGATGTGATGtt        <  1:607218/76‑1 (MQ=255)
                                                                   tttACGAGTCTATCGTAGTTGGTTGCATCAAAATGTGTGAGTTGGTTATTATCAAGATTTATAGATGTGAAGtttt      <  1:356191/76‑1 (MQ=255)
                                                                   tttACGAGTCTATCGTAGTTGGTTGCATCAAAATGTGTGAGTTGGTTATATCAAGATTTATAGATGTGAAGTTTTg     >  2:199997/1‑76 (MQ=255)
                                                                     tACGAGTCTATCGTAGTTGGTTGCATCAAAATGTGTGAGTTGGTTATTATCAAGATTTATAGATGTGAAGTTTTGc    <  1:494425/76‑1 (MQ=255)
                                                                       cGAGTCTATCGTAGTTGGTTGCATCAAAATGTGTGAGTTGGTTATTATCAAGATTTATAGATGTGAAGTTTTGCAg  >  1:602122/1‑76 (MQ=255)
                                                                          |                                                                        
TGTAATACTTACATTTCTGCCTTGAGGAAAATTTATTGACTCAAGAGCATTACTATTCAGACTTAGCTTTACGAGTCTATCGTAGTTGGTTGCATCAAAATGTGTGAGTTGGTTATTATCAAGATTTATAGATGTGAAGTTTTGCAG  >  REL606/1519401‑1519547

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