New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 150634120 (0.420)8 (0.180) 6/140 4.5 30.0% intergenic (+203/‑37) ECB_01417/ydcC hypothetical protein/hypothetical protein
?REL606 = 1506614 NA (NA)coding (237/792 nt) ydcC hypothetical protein
Rejected: Coverage evenness skew score above cutoff.

AGCTGTTTTTTATTGCTTATACATGATCAAATACTCCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAAATGATGGAACATATTTCTATTAT  >  REL606/1506319‑1506415
                      |                                                                          
aGCTGTTTTTTATTGCTTATACATGATCAAATACTCCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAAat                       >  1:472934/1‑76 (MQ=32)
 gCTGTTTTTTATTGCTTATACATGATCAAATACTCCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAAatg                      <  2:88346/76‑1 (MQ=34)
  cTGTTTTTTATTGCTTATACATGATCAAATACTCCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAAatga                     <  1:329913/76‑1 (MQ=34)
    gTTTTTTATTGCTTATACATGATCAAATACTCCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAAatgatg                   <  1:1088978/76‑1 (MQ=34)
    gTTTTTTATTGCTTATACATGATCAAATACTCCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAAatgatg                   >  1:1339216/1‑76 (MQ=34)
    gTTTTTTATTGCTTATACATGATCAAATACTCCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAAatgatg                   >  1:1044480/1‑76 (MQ=34)
      tttttATTGCTTATACATGATCAAATACTCCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAAATGATGGa                 >  1:609314/1‑76 (MQ=34)
       ttttATTGCTTATACATGATCAAATACTCCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAAATGATGGaa                <  1:178769/76‑1 (MQ=34)
       ttttATTGCTTATACATGATCAAATACTCCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAAATGATGGaa                <  1:684728/76‑1 (MQ=34)
         ttATTGCTTATACATGATCAAATACTCCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAAATGATGGAACa              <  1:839833/76‑1 (MQ=34)
           aTTGCTTATACATGATCAAATACTCCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAAATGATGGAACata            <  1:534471/76‑1 (MQ=34)
           aTTGCTTATACATGATCAAATACTCCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAAATGATGGAACata            <  1:998966/76‑1 (MQ=34)
            ttGCTTATACATGATCAAATACTCCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAAATGATGGAACatat           <  2:524414/76‑1 (MQ=34)
            ttGCTTATACATGATCAAACACTCCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAAATGATGGAACatat           <  2:610117/76‑1 (MQ=21)
             tGCTTATACATGATCAAATACTCCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAAATGATGTAACATAtt          >  1:259079/1‑76 (MQ=21)
             tGCTTATACATGATCAAATACTCCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAAATGATGGAACATAtt          >  1:259087/1‑76 (MQ=34)
             tGCTTATACATGATCAAATACTCCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAAATGATGGAACATAtt          >  1:584597/1‑76 (MQ=34)
                ttATACATGATCAAATACTCCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAAATGATGGAACATATTTCt       >  2:562304/1‑76 (MQ=34)
                   taCATGATCAAATACTCCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAAATGATGGAACATATTTCtatt    >  2:1332178/1‑76 (MQ=34)
                     cATGATCAAATACTCCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAAATGATGGAACATATTTCtattat  >  2:1376582/1‑76 (MQ=34)
                      |                                                                          
AGCTGTTTTTTATTGCTTATACATGATCAAATACTCCTTACATAATTAAGGAGAACAAAATGGAACTTAAAAAAATGATGGAACATATTTCTATTAT  >  REL606/1506319‑1506415

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