New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 2103888 =16 (0.250)11 (0.290)
+GCTGGCTGGCTGGCTGGCTGGCTGGCTGGCTGGCTGGCTGG
9/118 3.2 47.0% coding (155/216 nt) ECB_01992 hypothetical protein
?REL606 = 2103918 26 (0.400)coding (185/216 nt) ECB_01992 hypothetical protein
Rejected: Coverage evenness skew score above cutoff.

CTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCATTAGATTTCACAGCCAGCCAGCCAGCCAGCCAGCCAGCCAGGTGTGAATGTGGATTTTGTACTACGAAATAATTGCCCTTATTATACTATAAAAGGCAACACTAAAAAAGAAAACATTTCATCTTCTGCTGTTTTAACAAC  >  REL606/2103821‑2104018
                                                                                                 |                                                                                                    
cTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCATTAGATTTCACAGccagccagccagccagccagccagccagcca                                                                                                   <  1:179788/101‑4 (MQ=255)
cTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCATTAGATTTCACAGccagccagccagccagccagccagccagcca                                                                                                   >  1:1268249/1‑98 (MQ=255)
  ataCAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCATTAGATTTCACAGccagccagccagccagccagccagccagccagc                                                                                                 >  1:618914/1‑96 (MQ=255)
    aCAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCATTAGATTTCACAGccagccagccagccagccagccagccagccagcca                                                                                               <  1:1215902/101‑8 (MQ=255)
                                                            agccagccagccagccagccagccagccagccagcCAGGTGTGAATGGGGATTTTGTACTACGAAATAATTGCCCTTATTATACTATAAAAGGCAACACTa                                       >  1:451939/8‑101 (MQ=255)
                                                               cagccagccagccagccagccagccagccagcCAGGTGTGAATGTGGATTTTGTACTACGAAATAATTGCCCTTATTATACTATAAAAGGCAACACTaaaa                                    >  2:1175639/5‑101 (MQ=255)
                                                               cagccagccagccagccagccagccagccagcCAGGTGTGAATGTGGATTTTGTACTACGAAATAATTGCCCTTATTATACTATAAAAGGCAACACTaaaa                                    >  2:1175638/5‑101 (MQ=255)
                                                               cagccagccagccagccagccagccagccagcCAGGTGTGAATGTGGATTTTGTACTACGAAATAATTGCCCTTATTATACTATAAAAGGCAACACTaaaa                                    >  1:136185/5‑101 (MQ=255)
                                                                 gccggccagccagccagccagccagccagccaGGTGTGAATGTGGATTTTGTACTACGAAATAATTGCCCTTATTATACTATAAAAGGCAACACTaaaaaa                                  <  2:1568376/97‑1 (MQ=255)
                                                                     gtaaatcagccggccagccagccagccagGTGTGAATGTGGATTTTGTACTACGAAATAATTGCCCTTATTATACTATAAAAGGCAACACTAAAAAAGaaa                              <  2:1114518/95‑1 (MQ=255)
                                                                      ccagccagccagccagccagccagccagGTGTGAATGTGGATTTTGTACTACGAAATAATTGCCCTTATTATACTATAAAAGGCAACACTAAAAAAGaaaa                             <  2:10613/101‑1 (MQ=255)
                                                                       cagccagccagccagccagccagccagGTGTGAATGTGGATTTTGTACTACGAAATAATTGCCCTTATTATACTATAAAAGGCAACACTAAAAAAGAAAAc                            <  1:221267/101‑1 (MQ=255)
                                                                         gccagccagccagccagccagccagGTGTGAATGTGGATTTTGTACTACGAAATAATTTCCCTTATTATACTATAAAAGGCAACACTAAAAAAGAAAACAt                          >  1:811608/1‑101 (MQ=255)
                                                                         gccagccagccagccagccagccagGTGTGAATGTGGATTTTGTACTACGAAATAATTGCCCTTATTATACTATAAAAGGCAACACTAAAAAAGAAAACAt                          <  2:1609487/101‑1 (MQ=255)
                                                                         gccagccagccagccagccagccagGTGTGAATGTGGATTTTGTACTACGAAATAATTGCCCTTATTATACTATAAAAGGCAACACTAAAAAAGAAAACAt                          <  2:549587/101‑1 (MQ=255)
                                                                            agccagccagccagccagccagGTGTGAATGTGGATTTTGTACTACGAAATAATTGCCCTTATTATACTATAAAAGGCAACACTAAAAAAGAAAACATTTc                       >  2:798231/1‑101 (MQ=255)
                                                                              ccagccagccagccagccagGTGTGAATGTGGATTTTGTACTACGAAATAATTGCCCTTATTATACTATAAAAGGCAACACTAAAAAAGAAAACATTTCAt                     <  2:1268244/101‑1 (MQ=255)
                                                                               cagccagccagccagccagGTGTGAATGTGGATTTTGTACTACGAAATAATTGCCCTTATTATACTATAAAAGGCAACACTAAAAAAGAAAACATTTCAtc                    >  2:823730/1‑101 (MQ=255)
                                                                                  ccagccagccagccagGTGTGAATGTGGATTTTGTACTACGAAATAATTGCCCTTATTATACTATAAAAGGCAACACTAAAAAAGAAAACATTTCAtcttc                 >  2:1359627/1‑101 (MQ=255)
                                                                                     gccagccagccagGTGTGAATGTGGATTTTGTACTACGAAATAATTGCCCTTATTATACTATAAAAGGCAACACTAAAAAAGAAAACATTTCATCTTctgc              <  1:113636/101‑1 (MQ=255)
                                                                                      ccagccagccagGTGTGAATGTGGATTTTGTACTACGAAATAATTGCCCTTATTATACTATAAAAGGCAACACTAAAAAAGAAAACATTTCATCTTctgct             >  2:846567/1‑101 (MQ=255)
                                                                                      ccagccagccagGTGTGAATGTGGATTTTGTACTACGAAATAATTGCCCTTATTATACTATAAAAGGCAACACTAAAAAAGAAAACATTTCATCTTctgct             >  1:1702590/1‑101 (MQ=255)
                                                                                        agccagccagGTGTGAATGTGGATTTTGTACTACGAAATAATTGCCCTTATTATACTATAAAAGGCAACACTAAAAAAGAAAACATTTCATCTTCTGCTGt           <  1:953655/101‑1 (MQ=255)
                                                                                         gccagccagGTGTGAATGTGGATTTTGTACTACGAAATAATTGCCCTTATTATACTATAAAAGGCAACACTAAAAAAGAAAACATTTCATCTTCTGCTGtt          >  1:353703/1‑101 (MQ=255)
                                                                                          ccagccagGTGTGAATGTGGATTTTGTACTACGAAATAATTGCCCTTATTATACTATAAAAGGCAACACTAAAAAAGAAAACATTTCATCTTCTGCTGttt         >  1:598652/1‑101 (MQ=255)
                                                                                               cagGTGTGAATGTGGATTTTGTACTACGAAATAATTGCCCTTATTATACTATAAAAGGCAACACTAAAAAAGAAAACATTTCATCTTCTGCTGTTTTaaca    <  1:802209/101‑1 (MQ=255)
                                                                                                 ggTGTGAATGTGGATTTTGTACTACGAAATAATTGCCCTTATTATACTATAAAAGGCAACACTAAAAAAGAAAACATTTCATCTTCTGCTGTTTTaacaac  >  2:1242837/1‑101 (MQ=255)
                                                                                                 ggTGTGAATGTGGATTTTGTACTACGAAATAATTGCCCTTATTATACTATAAAAGGCAACACTAAAAAAGAAAACATTTCATCTTCTGCTGTTTTaacaac  >  2:261791/1‑101 (MQ=255)
                                                                                                 ggTGTGAATGTGGATTTTGTACTACGAAATAATTGCCCTTATTATACTATAAAAGGCAACACTAAAAAAGAAAACATTTCATCTTCTGCTGTTTTaacaac  >  2:1516846/1‑101 (MQ=255)
                                                                                                 ggTGTGAATGTGGATTTTGTACTACGAAATAATTGCCCTTATTATACTATAAAAGGCAACACTAAAAAAGAAAACATTTCATCTTCTGCTGTTTTaacaac  >  2:881831/1‑101 (MQ=255)
                                                                                                 |                                                                                                    
CTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCATTAGATTTCACAGCCAGCCAGCCAGCCAGCCAGCCAGCCAGGTGTGAATGTGGATTTTGTACTACGAAATAATTGCCCTTATTATACTATAAAAGGCAACACTAAAAAAGAAAACATTTCATCTTCTGCTGTTTTAACAAC  >  REL606/2103821‑2104018

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