New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 2312280 =75 (0.730)5 (0.050) 3/44 NT 9.1% intergenic (‑147/‑143) ais/yfbE hypothetical protein/uridine 5'‑(beta‑1‑threo‑pentapyranosyl‑4‑ulose diphosphate) aminotransferase, PLP‑dependent
?REL606 = 4042610 31 (0.330)intergenic (‑57/+10) yshA/yihO predicted outer membrane porin L/predited transporter

GGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTT  >  REL606/4042576‑4042644
                                  |                                  
gggAGGCCGCCACCCACCAGGATGGAAATAGCGGc                                    >  1:17137763/1‑35 (MQ=255)
gggAGGCCGCCACACACCAGGATGGAAATAGCGGc                                    <  1:15088418/35‑1 (MQ=255)
 ggAGGCCGCCACCCACCAGGATGGAAATAGCGGCt                                   <  1:7491178‑M2/35‑2 (MQ=255)
 ggAGGCCGCCACCCACCAGGATGGAAATAGCGGCt                                   >  1:6664507‑M2/1‑34 (MQ=255)
 ggAGGCCGCCACCCACCAGGATGGAAATAGCGGCt                                   <  1:19414125‑M2/35‑2 (MQ=255)
 ggAGGCCGCCACCCACCAGGATGGAAATAGCGGCt                                   >  1:14129374‑M2/1‑34 (MQ=255)
  gATGCCGCCACCCACCAGGATGGAAATAGCGGCta                                  <  1:3986566‑M2/35‑3 (MQ=255)
  gAGGCCGCCACCCACCAGTATGGAAATAGCGGCta                                  <  1:4657701‑M2/35‑3 (MQ=255)
  gAGGCCGCCACCCACCAGGATGGAAATAGCGGCta                                  >  1:8122389‑M2/1‑33 (MQ=255)
  gAGGCCGCCACCCACCAGGATGGAAATAGCGGCta                                  <  1:772473‑M2/35‑3 (MQ=255)
  gAGGCCGCCACCCACCAGGATGGAAATAGCGGCta                                  <  1:17172086‑M2/35‑3 (MQ=255)
  gAGCCCGCCACCCACCAGGATGGAAATAGCGGCta                                  >  1:3580627‑M2/1‑33 (MQ=255)
   aGGCCGCCACCCACCAGGATGGAAATAGCGGCtaa                                 <  1:9857231‑M2/35‑4 (MQ=255)
   aGGCCGCCACCCACCAGGATGGAAATAGCGGCtaa                                 <  1:15372953‑M2/35‑4 (MQ=255)
   aGGCCGCCACCCACCAGGATGGAAATAGCGGCtaa                                 <  1:13918176‑M2/35‑4 (MQ=255)
   aGGCCGCCACCCACCAGGATGGAAATAGCGGCTAt                                 <  1:5575313/35‑1 (MQ=255)
    ggCCGCCACCCACCAGGATGGAAATAGCGGCtaaa                                <  1:1902604‑M2/35‑5 (MQ=255)
    ggCCGCCACCCACCAGGATGGAAATAGCGGCtaaa                                <  1:14195482‑M2/35‑5 (MQ=255)
     gccgccACCCACCAGGATGGAAATAGCGGCTATTg                               <  1:9801432/35‑1 (MQ=255)
       cgccCCCCACCAGGATGGAAATAGCGGCTATTGAc                             <  1:8736454/35‑1 (MQ=255)
       cgccCCCCACCAGGATGGAAATAGCGGCTATTGAc                             <  1:10880030/35‑1 (MQ=255)
       cgccACCCACCAGGATGGAAATAGCGGCTATTGAc                             <  1:5640102/35‑1 (MQ=255)
        gccCCCCACCAGGATGGAAATAGCGGCTATTGACa                            <  1:1016252/35‑1 (MQ=255)
         ccacccacCAGGATGGAAATAGCGGCtaaaaaaaa                           <  1:13869675‑M2/35‑10 (MQ=255)
          ccccccccAGGAGGGAAATAGCGGCTATTGACAtt                          <  1:15616415/33‑1 (MQ=37)
          caccccccAGGAGGGAAATAGCGGCTATTGACAtt                          <  1:13783359/35‑1 (MQ=38)
          cacccaccAGGATGGAAATAGCGGCtaaaaaaaaa                          <  1:17227563‑M2/35‑11 (MQ=255)
          cacccaccAGGATGGAAATAGCGGCtaaaaaaaaa                          <  1:7749188‑M2/35‑11 (MQ=255)
          cacccaccAGGATGGAAATAGCGGCTATTGACAtt                          <  1:20306310/35‑1 (MQ=255)
           acccaccAGGATGGAAATAGCGGCtaaaaaaaaaa                         <  1:8886485‑M2/35‑12 (MQ=255)
           acccaccAGGATGGAAATAGCGGCTATTGACAttt                         <  1:12609566/35‑1 (MQ=255)
           acccaccAGGATGGAAATAGCGGCTATTGACAttt                         <  1:2712505/35‑1 (MQ=255)
           acacaccaGGATGGAAATAGCGGCTATTGACAttt                         <  1:13387022/32‑1 (MQ=255)
             ccaccaGGATGGAAATAGCGGCTATTGACATTtat                       <  1:10016120/35‑1 (MQ=255)
             ccaccaGGATGGAAATAGCGGCTATTGACATTtat                       <  1:16082600/35‑1 (MQ=255)
             ccaccaGGATGGAAATAGCGGCTATTGACATTtat                       <  1:1743319/35‑1 (MQ=255)
             ccaccaGGATGGAAATAGCGGCTATTGACATTtat                       <  1:12157127/35‑1 (MQ=255)
              caccaGGAGGGAAATAGCGGCTATTGACATTtata                      <  1:1138031/35‑1 (MQ=255)
               accaGGATGGAAATAGCGGCTATTGACATTTATaa                     <  1:3672691/35‑1 (MQ=255)
                  aGGATGGAAATAGCGGCtaaaaaaaaatataaaaa                  <  1:2680590‑M2/35‑19 (MQ=255)
                  aGGATGGAAATAGCGGCTATTGACATTTATAACGt                  <  1:18964838/35‑1 (MQ=255)
                  aGGATGGAAATAGCGGCTATTGACATTTATAACGt                  <  1:11906705/35‑1 (MQ=255)
                     ttGGAATTAGCGGCTATTGACATTTATAACGTTAc               <  1:4413014/34‑1 (MQ=255)
                     aTGGAAATAGCGTCTATTGACATTTATAACGTTAc               >  1:19838540/1‑35 (MQ=255)
                     aTGGAAATAGCGGCTATTGACATTTATAACGTTAc               <  1:17488342/35‑1 (MQ=255)
                      tGGAAATAGCGGCTATTGACATTTATAACGTTACg              <  1:6904458/35‑1 (MQ=255)
                        gAAATAGCGGCTATTGACATTTATAACGTTACGGa            <  1:2830170/35‑1 (MQ=255)
                        gAAATAGCGGCTATTGACATTTATAACGTTACGGa            <  1:15291800/35‑1 (MQ=255)
                        gAAATAGCGGCTATTGACATTTATAACGTTACGGa            <  1:5089339/35‑1 (MQ=255)
                         aaaTAGCGGCTATTGACATTTATAACGTTACGGaa           <  1:9115151/35‑1 (MQ=255)
                           aTAGCTGCTATTGACATTTATAACGTTACGGAAgg         >  1:4079124/1‑34 (MQ=255)
                           aTAGCGGCTATTGACATTTATAACGTTACGGAAGc         <  1:18956296/35‑1 (MQ=255)
                           aTAGCGGCTATTGACATTTATAACGTTACGGAAGc         <  1:3800737/35‑1 (MQ=255)
                           aTAGCGGCTATTGACATTTATAACGTTACGGAAGc         <  1:7686251/35‑1 (MQ=255)
                            tAGCGGCTATTGACATTTATAACGTTACGGAAGcc        <  1:13629523/35‑1 (MQ=255)
                            tAGCGGCTATTGACATTTATAACGTTACGGAAGcc        <  1:17842310/35‑1 (MQ=255)
                            tAGCGGCTATTGACATTTATAACGTTACGGAAGcc        <  1:9628691/35‑1 (MQ=255)
                             aGCGGCTATTGACATTTATAACGTTACGGAAGCCg       <  1:17293062/35‑1 (MQ=255)
                               cgtcTATTGACATTTATAACGTTACGGAAGGCGtt     >  1:9846369/4‑35 (MQ=37)
                                ggCTATTGACATTTATAAAGTTACGGAAGCCGttt    <  1:7473350/35‑1 (MQ=255)
                                 gCTATTGACATTTATAACGTTACGGAAGGCGcttt   >  1:15054750/1‑35 (MQ=37)
                                 gCTATTGACATTTATAACGTTACGGAAGCCGtttt   >  1:3036839/1‑35 (MQ=255)
                                 gCTATTGACATTTATAACGTTACGGAAGCCGgttt   >  1:12393060/1‑35 (MQ=255)
                                  cTATTGACATTTATAACGTTACGGAAGCCGttttt  >  1:18707881/1‑35 (MQ=255)
                                  |                                  
GGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCGTTTTT  >  REL606/4042576‑4042644

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