New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 267891387 (1.340)7 (0.120) 4/132 7.4 8.9% intergenic (+267/‑16) ECB_02509/ECB_02510 Fels‑2 prophage protein/Recombinase
?REL606 = 2679078 67 (1.180)coding (150/1557 nt) ECB_02510 Recombinase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TATTCTTTACATAATCAGTCCGAATATATCAAAGATTATGCTGAAAAGAACAATATGGAAATCGCTTATACCTACGATGATTCAGGTAAGAGCGGAGTCAGTATCGTAGGCAGGCATTCTTTGCAGCAGTTACTTAGCGATGTAGAACAA  >  REL606/2679004‑2679153
                                                                          |                                                                           
tatTCTTTACATAATCAGTCCGAATATATCAAAGATTATGCTGAAAAGAACAATATGGAAATCGCTTATACCTACg                                                                            <  1:480496/76‑1 (MQ=255)
tatTCTTTACATAATCAGTCCGAATATATCAAAGATTATGCTGAAAAGAACAATATGGAAATCGCTTATACCTACg                                                                            <  1:1765184/76‑1 (MQ=255)
tatTCTTTACATAATCAGTCCGAATATATCAAAGATTATGCTGAAAAGAACAATATGGAAATCGCTTATACCTACg                                                                            <  1:1765161/76‑1 (MQ=255)
  ttCTTTACATATTCAGTCCGAATATATCAAAGATTATGCTGAAAAGAACAATATGGAAATCGCTTATACCTACgat                                                                          >  1:1031319/1‑76 (MQ=255)
  ttCTTTACATAATCAGTCCGAATATATCAAAGATTATGCTGAAAAGAACAATATGGAAATCGCTTATACCTACgat                                                                          >  1:1414123/1‑76 (MQ=255)
  ttCTTTACATAATCAGTCCGAATATATCAAAGATTATGCTGAAAAGAACAATATGGAAATCGCTTATACCTACgat                                                                          >  1:1543733/1‑76 (MQ=255)
   tCTTTACATAATCAGTCCGAATATATCAAAGATTATGCTGAAAAGAACAATATGGAAATCGCTTATACCTACgatg                                                                         >  1:1058704/1‑76 (MQ=255)
   tCTTTACATAATCAGTCCGAATATATCAAAGATTATGCTGAAAAGAACAATATGGAAATCGCTTATACCTACgatg                                                                         >  2:1995100/1‑76 (MQ=255)
    cTTTACATAATCAGTCCGAATATATCAAAGATTATGCTGAAAAGAACAATATGGAAATCGCTTATACCTACgatga                                                                        >  2:1996499/1‑76 (MQ=255)
    cTTTACATAATCAGTCCGAATATATCAAAGATTATGCTGAAAAGAACAATATGGAAATCGCTTATACCTACgatga                                                                        >  1:455588/1‑76 (MQ=255)
         cATAATCAGTCCGAATATATCAAAGATTATGCTGAAAAGAACAATATGGAAATCGCTTATACCTACGATGATTCAg                                                                   >  1:2020175/1‑76 (MQ=255)
          aTAATCAGTCCGAATATATCAAAGATTATGCTGAAAAGAACAATATGGAAATCGCTTATACCTACGATGATTCAgg                                                                  >  2:1005416/1‑76 (MQ=255)
          aTAATCAGTCCGAATATATCAAAGATTATGCTGAAAAGAACAATATGGAAATCGCTTATACCTACGATGATTCAgg                                                                  <  1:1681542/76‑1 (MQ=255)
              tCAGTCCGAATATATCAAAGATTATGCTGAAAAGAACAATATGGAAATCGCTTATACCTACGATGATTCAGGTAag                                                              >  2:903206/1‑76 (MQ=255)
               cAGTCCGAATATATCAAAGATTATGCTGAAAAGAACAATATGGAAATCGCTTATACCTACGATGATTCAGGTAaga                                                             >  1:500129/1‑76 (MQ=255)
               cAGTCCGAATATATCAAAGATTATGCTGAAAAGAACAATATGGAAATCGCTTATACCTACGATGATTCAGGTAaga                                                             <  1:25711/76‑1 (MQ=255)
                   ccGAATATATCAAAGATTATGCTGAAAAGAACAATATGGAAATCGCTTATACCTACGATGATTCAGGTAAGAGCgg                                                         >  1:785291/1‑76 (MQ=255)
                   ccGAATATATCAAAGATTATGCTGAAAAGAACAATATGGAAATCGCTTATACCTACGATGATTCAGGTAAGAGCgg                                                         >  1:785272/1‑76 (MQ=255)
                    cGAATATATCAAAGATTATGCTGAAAAGAACAATATGGAAATCGCTTATACCTACGATGATTCAGGTAAGAGCGGa                                                        <  2:120456/76‑1 (MQ=255)
                    cGAATATATCAAAGATTATGCTGAAAAGAACAATATGGAAATCGCTTATACCTACGATGATTCAGGTAAGAGCGGa                                                        <  1:1742195/76‑1 (MQ=255)
                    cGAATATATCAAAGATTATGCTGAAAAGAACAATATGGAAATCGCTTATACCTACGATGATTCAGGTAAGAGCGGa                                                        <  2:229723/76‑1 (MQ=255)
                    cGAATATATCAAAGATTATGCTGAAAAGAACAATATGGAAATCGCTTATACCTACGATGATTCAGGTAAGAGCGGa                                                        >  1:483919/1‑76 (MQ=255)
                        tatatCAAAGATTATGCTGAAAAGAACAATATGGAAATCGCTTATACCTACGATGATTCAGGTAAGAGCGGAGTCa                                                    <  1:773461/76‑1 (MQ=255)
                        tatatCAAAGATTATGCTGAAAAGAACAATATGGAAATCGCTTATACCTACGATGATTCAGGTAAGAGCGGAGTCa                                                    >  2:516738/1‑76 (MQ=255)
                         atatCAAAGATTATGCTGAAAAGAACAATATGGAAATCGCTTATACCTACGATGATTCAGGTAAGAGCGGAGTCAg                                                   <  1:1827086/76‑1 (MQ=255)
                         atatCAAAGATTATGCTGAAAAGAACAATATGGAAATCGCTTATACCTACGATGATTCAGGTAAGAGCGGAGTCAg                                                   <  1:1827112/76‑1 (MQ=255)
                          tatCAAAGATTATGCTGAAAAGAACAATATGGAAATCGCTTATACCTACGATGATTCAGGTAAGAGCGGAGTCAGt                                                  >  1:111355/1‑76 (MQ=255)
                            tCAAAGACTATGCTGAAAAGAACAATATGGAAATCGCTTATACCTACGATGATTCAGGTAAGAGCGGAGTCAGTAt                                                <  2:1838667/76‑1 (MQ=255)
                             cAAAGATTATGCTGAAAAGAACAATATGGAAATCGCTTATACCTACGATGATTCAGGTAAGAGCGGAGTCAGTATc                                               >  2:1852947/1‑76 (MQ=255)
                             cAAAGATTATGCTGAAAAGAACAATATGGAAATCGCTTATACCTACGATGATTCAGGTAAGAGCGGAGTCAGTATc                                               >  1:675510/1‑76 (MQ=255)
                             cAAAGATTATGCTGAAAAGAACAATATGGAAATCGCTTATACCTACGATGATTCAGGTAAGAGCGGAGTCAGTATc                                               >  1:1647482/1‑76 (MQ=255)
                               aaGATTATGCTGAAAAGAACAATATGGAAATCGCTTATACCTACGATGATTCAGGTAAGAGCGGAGTCAGTATCGt                                             <  2:121862/76‑1 (MQ=255)
                                aGATTATGCTGAAAAGAACAATATGGAAATCGCTTATACCTACGATGATTCAGGTAAGAGCGGAGTCAGTATCGTa                                            >  2:1246668/1‑76 (MQ=255)
                                 gATTATGCTGAAAAGAACAATATGGAAATCGCTTATACCTACGATGATTCAGGTAAGAGCGGAGTCAGTATCGTag                                           >  1:1162469/1‑76 (MQ=255)
                                 gATTATGCTGAAAAGAACAATATGGAAATCGCTTATACCTACGATGATTCAGGTAAGAGCGGAGTCAGTATCGTag                                           >  1:2023207/1‑76 (MQ=255)
                                  aTTATGCTGAAAAGAACAATATGGAAATCGCTTATACCTACGATGATTCAGGTAAGAGCGGAGTCAGTATCGTagg                                          >  2:1719352/1‑76 (MQ=255)
                                       gCTGAAAAGAACAATATGGAAATCGCTTATACCTACGATGATTCAGGTAAGAGCGGAGTCAGTATCGTaggcaggc                                     <  1:1978988/76‑1 (MQ=255)
                                         tGAAAAGAACAATATGGAAATCGCTTATACCTACGATGATTCAGGTAAGAGCGGAGTCAGTATCGTAGGCAGGCAt                                   >  1:1016776/1‑76 (MQ=255)
                                         tGAAAAGAACAATATGGAAATCGCTTATACCTACGATGATTCAGGTAAGAGCGGAGTCAGTATCGTAGGCAGGCAt                                   >  1:1613416/1‑76 (MQ=255)
                                          gAAAAGAACAATATGGAAATCGCTTATACCTACGATGATTCAGGTAAGAGCGGAGTCAGTATCGTAGGCAGGCAtt                                  >  1:984480/1‑76 (MQ=255)
                                          gAAAAGAACAATATGGAAATCGCTTATACCTACGATGATTCAGGTAAGAGCGGAGTCAGTATCGTAGGCAGGCAtt                                  <  1:1135180/76‑1 (MQ=255)
                                          gAAAAGAACAATATGGAAATCGCTTATACCTACGATGATTCAGGTAAGAGCGGAGTCAGTATCGTAGGCAGGCAtt                                  >  2:458214/1‑76 (MQ=255)
                                           aaaaGAACAATATGGAAATCGCTTATACCTACGATGATTCAGGTAAGAGCGGAGTCAGTATCGTAGGCAGGCATTc                                 >  1:413899/1‑76 (MQ=255)
                                              aGAACAATATGGAAATCGCTTATACCTACGATGATTCAGGTAAGAGCGGAGTCAGTATCGTAGGCAGGCATTCttt                              <  2:1654206/76‑1 (MQ=255)
                                              aGAACAATATGGAAATCGCTTATACCTACGATGATTCAGGTAAGAGCGGAGTCAGTATCGTAGGCAGGCATTCttt                              <  2:158399/76‑1 (MQ=255)
                                               gAACAATATGGAAATCGCTTATACCTACGATGATTCAGGTAAGAGCGGAGTCAGTATCGTAGGCAGGCATTCTTTg                             >  2:612437/1‑76 (MQ=255)
                                                aaCAATATGGAAATCGCTTATACCTACGATGATTCAGGTAAGAGCGGAGTCAGTATCGTAGGCAGGCATTCTTTgc                            <  2:1071357/76‑1 (MQ=255)
                                                 aCAATATGGAAATCGCTTATACCTACGATGATTCAGGTAAGAGCGGAGTCAGTATCGTAGGCAGGCATTCTTTgca                           >  2:233459/1‑76 (MQ=255)
                                                  cAATATGGAAATCGCTTATACCTACGATGATTCAGGTAAGAGCGGAGTCAGTATCGTAGGCAGGCATTCTTTgcag                          >  2:550702/1‑76 (MQ=255)
                                                  cAATATGGAAATCGCTTATACCTACGATGATTCAGGTAAGAGCGGAGTCAGTATCGTAGGCAGGCATTCTTTgata                          <  2:111355/76‑4 (MQ=255)
                                                     tatGGAAATCGCTTATACCTACGATGATTCAGGTAAGAGCGGAGTCAGTATCGTAGGCAGGCATTCTTTGcagcag                       <  1:1481438/76‑1 (MQ=255)
                                                      atGGAAATCGCTTATACCTACGGTGATTCAGGTAAGAGCGGAGTCAGTATCGTAGGCAGGCATTCTTTGCAGCAGt                      >  2:378701/1‑76 (MQ=255)
                                                        ggAAATCGCTTATACCTACGATGATTCAGGTAAGAGCGGAGTCAGTATCGTAGGCAGGCATTCTTTGCAGCAGTTa                    <  2:496766/76‑1 (MQ=255)
                                                         gAAATCGCTTATACCTACGATGATTCAGGTAAGAGCGGAGTCAGTATCGTAGGCAGGCATTCTTTGCAGCAGTTAc                   >  2:429229/1‑76 (MQ=255)
                                                          aaaTCGCTTATACCTACGATGATTCAGGTAAGAGCGGAGTCAGTATCGTAGGCAGGCATTCTTTGCAGCAGTTACt                  <  1:158250/76‑1 (MQ=255)
                                                           aaTCGCTTATACCTACGATGATTCAGGTAAGAGCGGAGTCAGTATCGTAGGCAGGCATTCTTTGCAGCAGTTACtt                 <  2:578517/76‑1 (MQ=255)
                                                           aaTCGCTTATACCTACGATGATTCAGGTAAGAGCGGAGTCAGTATCGTAGGCAGGCATTCTTTGCAGCAGTTACtt                 <  2:365535/76‑1 (MQ=255)
                                                           aaTCGCTTATACCTACGATGATTCAGGTAAGAGCGGAGTCAGTATCGTAGGCAGGCATTCTTTGCAGCAGTTACtt                 <  1:477339/76‑1 (MQ=255)
                                                           aaTCGCTTATACCTACGATGATTCAGGTAAGAGCGGAGTCAGTATCGTAGGCAGGCATTCTTTGCAGCAGTTACtt                 <  1:1192135/76‑1 (MQ=255)
                                                             tCGCTTATACCTACGATGATTCAGGTAAGAGCGGAGTCAGTATCGTAGGCAGGCATTCTTTGCAGCAGTTACTTAg               >  2:1210895/1‑76 (MQ=255)
                                                             tCGCTTATACCTACGATGATTCAGGTAAGAGCGGAGTCAGTATCGTAGGCAGGCATTCTTTGCAGCAGTTACTTAg               >  1:1507343/1‑76 (MQ=255)
                                                             tCGCTTATACCTACGATGATTCAGGTAAGAGCGGAGTCAGTATCGTAGGCAGGCATTCTTTGCAGCAGTTACTTAg               >  1:1507347/1‑76 (MQ=255)
                                                                 ttATACCTACGATGATTCAGGTAAGAGCGGAGTCAGTATCGTAGGCAGGCATTCTTTGCAGCAGTTACTTAGCGAt           >  2:1030728/1‑76 (MQ=255)
                                                                  tataCCTACGATGATTCAGGTAAGAGCGGAGTCAGTATCGTAGGCAGGCATTCTTTGCAGCAGTTTCTTAGCGATg          >  2:1939297/1‑76 (MQ=255)
                                                                   ataCCTACGATGATTCAGGTAAGAGCGGAGTCAGTATCGTAGGCAGGCATTTTTTGCAGTAGTTACTTAGCGATGt         >  2:1562436/1‑76 (MQ=255)
                                                                   ataCCTACGATGATTCAGGTAAGAGCGGAGTCAGTATCGTAGGCAGGCATTCTTTGCAGCAGTTa                    >  1:107122/1‑65 (MQ=255)
                                                                   ataCCTACGATGATTCAGGTAAGAGCGGAGTCAGTATCGTAGGCAGGCATTCTTTGCAGCAGTTa                    <  2:107122/65‑1 (MQ=255)
                                                                   ataCCTACGATGATTCAGGTAAGAGCGGAGTCAGTATCGTAGGCAGGCATTCTTTGCAGCAGTTACTTAGCGATGt         >  1:1297430/1‑76 (MQ=255)
                                                                   ataCCTACGATGATTCAGGTAAGAGCGGAGTCAGTATCGTAGGCAGGCATTCTTTGCAGCAGTTACTTAGCGATGt         >  1:4176/1‑76 (MQ=255)
                                                                   ataCCTACGATGATTCAGGTAAGAGCGGAGTCAGTATCGTAGGCAGGCATTCTTTGCAGCAGTTACTTAGCGATGt         >  2:1562415/1‑76 (MQ=255)
                                                                     aCCTACGATGATTCAGGTAAGAGCGGAGTCAGTATCGTAGGCAGGCATTCTTTGCAGCAGTTACTTAGCGATGTAg       >  1:266275/1‑76 (MQ=255)
                                                                     aCCTACGATGATTCAGGTAAGAGCGGAGTCAGTATCGTAGGCAGGCATTCTTTGCAGCAGTTACTTAGCGATGTAg       >  1:1382589/1‑76 (MQ=255)
                                                                     aCCTACGATGATTCAGGTAAGAGCGGAGTCAGTATCGTAGGCAGGCATTCTTTGCAGCAGTTACTTAGCGATGTAg       <  2:1523924/76‑1 (MQ=255)
                                                                      ccTACGATGATTCAGGTAAGAGCGGAGTCAGTATCGTAGGCAGGCATTCTTTGCAGCAGTTACTTAGCGATGTAGa      <  2:1632794/76‑1 (MQ=255)
                                                                       cTACGATGATTCAGGTAAGAGCGGAGTCAGTATCGTAGGCAGGCATTCTTTGCAGCAGTTACTTAGCGATGTAGaa     >  2:1964198/1‑76 (MQ=255)
                                                                        tACGATGATTCAGGTAAGAGCGGAGTCAGTATCGTAGGCAGGCATTCTTTGCAGCAGTTACTTAGCGATGTAGAAc    >  2:1035751/1‑76 (MQ=255)
                                                                          cGATGATTCAGGTAAGAGCGGAGTCAGTATCGTAGGCAGGCATTCTTTGCAGCAGTTACTTAGCGATGTAGAACaa  <  2:1500517/76‑1 (MQ=255)
                                                                          |                                                                           
TATTCTTTACATAATCAGTCCGAATATATCAAAGATTATGCTGAAAAGAACAATATGGAAATCGCTTATACCTACGATGATTCAGGTAAGAGCGGAGTCAGTATCGTAGGCAGGCATTCTTTGCAGCAGTTACTTAGCGATGTAGAACAA  >  REL606/2679004‑2679153

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