New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 268044734 (1.540)3 (0.150) 3/134 2.6 10.8% coding (1519/1557 nt) ECB_02510 Recombinase
?REL606 = 2680655 19 (0.960)coding (181/897 nt) ECB_02511 ParB family protein
Rejected: Frequency below/above cutoff threshold.

AAAGCATAAAAAAAAGAGGATTAAGCAAGCCTATAAGCGTAAGAGCTATTGATGAAGGCGATTTCAAATATGCTTTAATTTGTGGTCAAGGGAGAATAGAGGCTCTCGTTGCATTAGGTGAAACTATTATTCCAGCAATTATAAGAGA  >  REL606/2680581‑2680728
                                                                          |                                                                         
aaaGCATAAAAAAAAGAGGATTAAGCAAGCCTATAAGCGTAAGAGCTATTGATGAAGGCGATTTCAAATATGCTTt                                                                          <  2:427329‑M2/76‑2 (MQ=255)
           aaaaGAGGATTAAGCAAGCCTATAAGCGTAAGAGCTATTGATGAAGGCGATTTCAAATATGCTTTAATTTGTGGTc                                                               >  2:328131/1‑76 (MQ=255)
                   aTTAAGCAAGCCTATAAGCGTAAGAGCTATTGATGAAGGCGATTTCAAATATGCTTtaatttgtagaagtttatca                                                       >  2:237780‑M2/1‑56 (MQ=255)
                   aTTAAGCAAGCCTATAAGCGTAAGAGCTATTGATGAAGGCGATTTCAAATATGCTTTAATTTGTGGTCAAGGgaga                                                       >  1:480004/1‑76 (MQ=255)
                    ttAAGCAAGCCTATAAGCGTAAGAGCTATTGATGAAGGCGATTTCAAATATGCTTTAATTTGTGGTCAAGGGAGaa                                                      >  1:539329/1‑76 (MQ=255)
                     tAAGCAAGCCTATAAGCGTAAGAGCTATTGATGAAGGCGATTTCAAATATGCTTTAATTTGTGGTCAAGGGAGAAt                                                     <  2:419303/76‑1 (MQ=255)
                     tAAGCAAGCCTATAAGCGTAAGAGCTATTGATGAAGGCGATTTCAAATATGCTTTAATTTGTGGTCAAGGGAGAAt                                                     <  1:608510/76‑1 (MQ=255)
                       agcaagcCTATAAGCGTAAGAGCTATTGATGAAGGCGATTTCAAATATGCTTtaatttgtagaagtttatcaagat                                                   >  2:95478‑M2/1‑52 (MQ=255)
                         caagcCTATAAGCGTAAGAGCTATTGATGAAGGCGATTTCAAATATGCTTtaatttgtagaagtttatcaagatta                                                 >  2:64102‑M2/1‑50 (MQ=255)
                            gcCTATAAGCGTAAGAGCTATTGATGAAGGCGATTTCAAATATGCTTTAATTTGTGGTCAAGGGAGAATAGAGGct                                              >  1:335960/1‑76 (MQ=255)
                               tataAGCGTAAGAGCTATTGATGAAGGCGATTTCAAATATGCTTTAATTTGTGGTCAAGGGAGAATAGAGGCtctc                                           <  1:402698/76‑1 (MQ=255)
                                 taAGCGTAAGAGCTATTGATGAAGGCGATTTCAAATATGCTTTAATTTGTGGTCAAGGGAGAATAGAGGCTCTCGt                                         >  2:339750/1‑76 (MQ=255)
                                      gTAAGAGCTATTGATGAAGGCGATTTCAAATATGCTTTAATTTGTGGTCAAGGGAGAATAGAGGCTCTCGTTGCAt                                    >  2:408009/1‑76 (MQ=255)
                                      gTAAGAGCTATTGATGAAGGCGATTTCAAATATGCTTTAATTTGTGGTCAAGGGAGAATAGAGGCTCTCGTTGCAt                                    >  2:408012/1‑76 (MQ=255)
                                         agagCTATTGATGAAGGCGATTTCAAATATGCTTTAATTTGTGGTCAAGGGAGAATAGAGGCTCTCGTTGCATTAg                                 <  2:400013/76‑1 (MQ=255)
                                                  gatgaAGGCGATTTCAAATATGCTTTAATTTGTGGTCAAGGGAGAATAGAGGCTCTCGTTGCATTAGGTGAAACta                        >  2:657896/1‑76 (MQ=255)
                                                           gATTTCAAATATGCTTTAATTTGTGGTCAAGGGAGAATAGAGGCTCTCGTTGCATTAGGTGAAACTATTATTCCAg               <  1:443985/76‑1 (MQ=255)
                                                           gATTTCAAATATGCTTTAATTTGTGGTCAAGGGAGAATAGAGGCTCTCGTTGCATTAGGTGAAACTATTATTCCAg               <  2:603335/76‑1 (MQ=255)
                                                           gATTTCAAATATGCTTTAATTTGTGGTCAAGGGAGAATAGAGGCTCTCGTTGCATTAGGTGAAACTATTATTCCAg               <  1:147845/76‑1 (MQ=255)
                                                                 aaaTATGCTTTAATTTGTGGTCAAGGGAGAATAGAGGCTCTCGTTGCATTAGGTGAAACTATTATTCCAGCAATta         <  2:681966/76‑1 (MQ=255)
                                                                  aaTATGCTTTAATTTGTGGTCAAGGGAGAATAGAGGCTCTCGTTGCATTAGGTGAAACTATTATTCCAGCAATtat        >  1:437444/1‑76 (MQ=255)
                                                                       gCTTTAATTTGTGGTCAAGGGAGAATAGAGGCTCTCGTTGCATTAGGTGAAACTATTATTCCAGCAATTATAagag   <  1:275616/76‑1 (MQ=255)
                                                                        cTTTAATTTGTGGTCAAGGGAGAATAGAGGCTCTCGTTGCATTAGGTGAAACTATTATTCCAGCAATTATAAgaga  >  1:559770/1‑76 (MQ=255)
                                                                          |                                                                         
AAAGCATAAAAAAAAGAGGATTAAGCAAGCCTATAAGCGTAAGAGCTATTGATGAAGGCGATTTCAAATATGCTTTAATTTGTGGTCAAGGGAGAATAGAGGCTCTCGTTGCATTAGGTGAAACTATTATTCCAGCAATTATAAGAGA  >  REL606/2680581‑2680728

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