New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 302154745 (0.670)4 (0.070) 3/124 8.4 9.9% coding (1138/1140 nt) ECB_02821 hypothetical protein
?REL606 = 3021727 35 (0.630)coding (958/1140 nt) ECB_02821 hypothetical protein
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

AACAACAACAAGAGATATAAATTCTTAATCGCATTCATAAAAAATAATCATCATCCTGACAACAAGAATATTTTAATATATTGACTCTAGGAGTTTTGTCTTTATTTTCTTATCTTCATACACATAATTACATAGCCATTGTATTCTCT  >  REL606/3021473‑3021621
                                                                          |                                                                          
aacaacaacAAGAGATATAAATTCTTAATCGCATTCATAAAAAATAATCATCATCCTGACAACAAGAATATTTTa                                                                            >  2:1067512/1‑75 (MQ=255)
 acaacaacaAGAGATATAAATTCTTAATCGCATTCATAAAAAATAATCATCATCCTGACAACAAGAATATTTTaa                                                                           >  1:2054920/1‑75 (MQ=255)
     caacaaGAGATATAAATTCTTAATCGCATTCATAAAAAATAATCATCATCCTGACAACAAGAATATTTTAatata                                                                       <  2:1767975/75‑1 (MQ=255)
     caacaaGAGATATAAATTCTTAATCGCATTCATAAAAAATAATCATCATCCTGACAACAAGAATATTTTAatata                                                                       <  1:1875861/75‑1 (MQ=255)
       acaaGAGATATAAATTCTTAATCGCATTCATAAAAAATAATCATCATCCTGACAACAAGAATATTTTAATATAtt                                                                     <  1:695355/75‑1 (MQ=255)
         aaGAGATATAAATTCTTAATCGCATTCATAAAAAATAATCATCATCCTGACAACAAGAATATTTTAATATATTGa                                                                   >  1:1898001/1‑75 (MQ=255)
            agaTATAAATTCTTAATCGCATTCATAAAAAATAATCATCATCCTGACAACAAGAATATTTTAATATATTGActc                                                                <  1:1052968/75‑1 (MQ=255)
             gaTATAAATTCTTAATCGCATTCATAAAAAATAATCATCATCCTGACAACAAGAATATTTTAATATATTGActct                                                               <  2:1334208/75‑1 (MQ=255)
             gaTATAAATTCTTAATCGCATTCATAAAAAATAATCATCATCCTGACAACAAGAATATTTTAATATATTGActct                                                               >  1:18456/1‑75 (MQ=255)
                            tCGCATTCATAAAAAATAATCATCATCCTGACAACAAGAATATTTTAATATATTGACTCTAGGAGTTTTGTCttt                                                >  2:1068088/1‑75 (MQ=255)
                               cATTCATAAAAAATAATCATCATCCTGACAACAAGAATATTTTAATATATTGACTCTAGGAGTTTTGTCTTTAtt                                             <  2:1835442/75‑1 (MQ=255)
                                 ttCATAAAAAATAATCATCATCCTGCCAACAAGAATATTTTAATATCTTGACTCTAGGAGTTTTGTCTTTAtttt                                           >  1:2039468/1‑75 (MQ=255)
                                   cATAAAAAATAATCATCATCCTGACAACAAGAATATTTTAATATATTGACTCTAGGAGTTTTGTCTTTATTTTCt                                         <  1:2113753/75‑1 (MQ=255)
                                   cATAAAAAATAATCATCATCCTGACAACAAGAATATTTTAATATATTGACTCTAGGAGTTTTGTCTTTATTTTCt                                         <  1:2113742/75‑1 (MQ=255)
                                    aTAAAAAATAATCATCATCCTGACAACAAGAATATTTTAATATATTGACTCTAGGAGTTTTGTCTTTATTTTCtt                                        >  1:2110203/1‑75 (MQ=255)
                                         aaaTAATCATCATCCTGACAACAAGAATATTTTAATATATTGACTCTAGGAGTTTTGTCTTTATTTTCTTATCtt                                   >  1:1937872/1‑75 (MQ=255)
                                          aataatCATCATCCTGACAACAAGAATATTTTAATATATTGACTCTAGGAGTTTTGTCTTTATTTTCTTATCTTc                                  <  1:1254226/75‑1 (MQ=255)
                                             aatCATCATCCTGACAACAAGAATATTTTAATATATTGACTCTAGGAGTTTTGTCTTTATTTTCTTATCTTCATa                               >  2:173536/1‑75 (MQ=255)
                                                 atcatcCTGACAACAAGAATATTTTAATATATTGACTCTAGGAGTTTTGTCTTTATTTTCTTATCTTCATAcaca                           >  2:340428/1‑75 (MQ=255)
                                                  tcatcCTGACAACAAGAATATTTTAATATATTGACTCTAGGAGTTTTGTCTTTATTTTCTTATCTTCATACACAt                          <  1:1296920/75‑1 (MQ=255)
                                                  tcatcCTGACAACAAGAATATTTTAATATATTGACTCTAGGAGTTTTGTCTTTATTTTCTTATCTTCATACACAt                          <  1:1177433/75‑1 (MQ=255)
                                                   catcCTGACAACAAGAATATTTTAATATATTGACTCTAGGAGTTTTGTCTTTATTTTCTTATCTTCATACACATa                         <  1:471238/75‑1 (MQ=255)
                                                      ccTGACAACAAGAATATTTTAATATATTGACTCTAGGAGTTTTGTCTTTATTTTCTTATCTTCATACACATAAtt                      <  2:2129039/75‑1 (MQ=255)
                                                      ccTGACAACAAGAATATTTTAATATATTGACTCTAGGAGTTTTGTCTTTATTTTCTTATCTTCATACACATAAtt                      <  2:2129028/75‑1 (MQ=255)
                                                      ccTGACAACAAGAATAGTTTAATATATTGACTCTAGGAGTTTTGTCTTTATTTTCTTATCTTCATACACATAAtt                      <  2:2129063/75‑1 (MQ=255)
                                                        tGACAACAAGAATATTTTAATATATTGACTCTAGGAGTTTTGTCTTTATTTTCTTATCTTCATACACATAATTAc                    <  1:544359/75‑1 (MQ=255)
                                                          acaacaAGAATATTTTAATATATTGACTCTAGGAGTTTTGTCTTTATTTTCTTATCTTCATACACATAATTACAt                  <  1:261892/75‑1 (MQ=255)
                                                          acaacaAGAATATTTTAATATATTGACTCTAGGAGTTTTGTCTTTATTTTCTTATCTTCATACACATAATTACAt                  <  1:552812/75‑1 (MQ=255)
                                                          acaacaAGAATATTTTAATATATTGACTCTAGGAGTTTTGTCTTTATTTTCTTATCTTCATACACATAATTACAt                  <  1:996583/75‑1 (MQ=255)
                                                           caacaaGAATATTTTAATATATTGACTCTAGGAGTTTTGTCTTTATTTTCTTATCTTCATACACATAATTACATa                 >  2:1160369/1‑75 (MQ=255)
                                                            aacaaGAATATTTTAATATATTGACTCTAGGAGTTTTGTCTTTATTTTCTTATCTTCATACACATAATTACATAg                >  2:347699/1‑75 (MQ=255)
                                                            aacaaGAATATTTTAATATATTGACTCTAGGAGTTTTGTCTTTATTTTCTTATCTTCATACACATAATTACATAg                >  1:1369091/1‑75 (MQ=255)
                                                            aacaaGAATATTTTAATATATTGACTCTAGGAGTTTTGTCTTTATTTTCTTATCTTCATACACATAATTACATAg                <  1:529671/75‑1 (MQ=255)
                                                            aacaaGAATATTTTAATATATTGACTCTAGGAGTTTTGTCTTTATTTTCTTATCTTCATACACATAATTACATAg                <  1:2016617/75‑1 (MQ=255)
                                                            aacaaGAATATTTTAATATATTGACTCTAGGAGTTTTGTCTTTATTTTCTTATCTTCATACACATAATTACATAg                >  1:1369086/1‑75 (MQ=255)
                                                              caaGAATATTTTAATATATTGACTCTAGGAGTTTTGTCTTTATTTTCTTATCTTCATACACATAATTACATAGcc              <  1:1602614/75‑1 (MQ=255)
                                                              caaGAATATTTTAATATATTGACTCTAGGAGTTTTGTCTTTATTTTCTTATCTTCATACACATAATTACATAGcc              <  1:1602610/75‑1 (MQ=255)
                                                                 gAATATTTTAATATATTGACTCTAGGAGTTTTGTCTTTATTTTCTTATCTTCATACACATAATTACATAGCCAtt           >  2:902235/1‑75 (MQ=255)
                                                                 gAATATTTTAATATATTGACTCTAGGAGTTTTGTCTTTATTTTCTTATCTTCATACACATAATTACATAGCCAtt           <  1:1231531/75‑1 (MQ=255)
                                                                 gAATATTTTAATATATTGACTCTAGGAGTTTTGTCTTTATTTTCTTATCTTCATACACATAATTACATAGCCAtt           >  1:1089343/1‑75 (MQ=255)
                                                                   atatTTTAATATATTGACTCTAGGAGTTTTGTCTTTATTTTCTTATCTTCATACACATAATTACATAGCCATTGt         <  1:1067512/75‑1 (MQ=255)
                                                                        ttAATATATTGACTCTAGGAGTTTTGTCTTTATTTTCTTATCTTCATACACATAATTACATAGCCATTGTATtct    >  1:704924/1‑75 (MQ=255)
                                                                         tAATATATTGACTCTAGGAGTTTTGTCTTTATTTTCTTATCTTCATACACATAATTACATAGCCATTGTATtctc   <  1:917378/75‑1 (MQ=255)
                                                                         tAATATATTGACTCTAGGAGTTTTGTCTTTATTTTCTTATCTTCATACACATAATTACATAGCCATTGTATtctc   <  1:917370/75‑1 (MQ=255)
                                                                          aaTATATTGACTCTAGGAGTTTTGTCTTTATTTTCTTATCTTCATACACATAATTACATAGCCATTGTATtctct  <  1:655642/75‑1 (MQ=255)
                                                                          aaTATATTGACTCTAGGAGTTTTGTCTTTATTTTCTTATCTTCATACACATAATTACATAGCCATTGTATtctct  >  2:635602/1‑75 (MQ=255)
                                                                          |                                                                          
AACAACAACAAGAGATATAAATTCTTAATCGCATTCATAAAAAATAATCATCATCCTGACAACAAGAATATTTTAATATATTGACTCTAGGAGTTTTGTCTTTATTTTCTTATCTTCATACACATAATTACATAGCCATTGTATTCTCT  >  REL606/3021473‑3021621

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