New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 556412 =34 (0.640)7 (0.130) 5/60 NT 18.7% coding (1697/2346 nt) ECB_00513 conserved hypothetical protein
?REL606 1539626 = 27 (0.440)coding (46/2400 nt) dos cAMP phosphodiesterase, heme‑regulated

GACGGCAGAAAGAGTCGGTGAATGGTGGTTGTTTACCGCCACTCTTTTTGCAGATGCGGAAATGATAATCA  >  REL606/556377‑556447
                                   |                                   
cgatgcggaaaatgccgccgatggcatttttttccc                                     <  1:967809‑M1/1‑1 (MQ=255)
cgatgcggaaaatgccgccgatggcatttttttccc                                     >  1:27282‑M1/36‑36 (MQ=255)
 gatgcggaaaatgccgccgatggcatttttttcccc                                    >  1:5982825‑M1/35‑36 (MQ=255)
 gatgcggaaaatgccgccgatggcatttttttcccc                                    >  1:4652702‑M1/35‑36 (MQ=255)
 acggCAGAAAGAGTCGGTGAATGGTGGTTGTTTAcc                                    <  1:5822412/36‑1 (MQ=255)
 acggCAGAAAGAGTCGGTGAATGGTGGTTGTTTAcc                                    <  1:6672952/36‑1 (MQ=255)
 acggCAGAAAGAGTCGGTGAATGGTGGTTGTTTAcc                                    <  1:4219774/36‑1 (MQ=255)
  cggCAGAAAGAGTCGGTGAATGGTGGTTGTTTACCg                                   <  1:8311220/36‑1 (MQ=255)
   tgcggaaaatgccgccgatggcatttttttccCCGc                                  >  1:8013653‑M1/33‑36 (MQ=255)
   ggCAGAAAGAGTCGGTGAATGGTGGTTGTTTACCGc                                  <  1:937814/36‑1 (MQ=255)
   ggCAGAAAGAGTCGGTGAATGGTGGTTGTTTACCGc                                  >  1:7842908/1‑36 (MQ=255)
    gcggaaaatgccgccgatggcatttttttccCCGcc                                 <  1:1529973‑M1/5‑1 (MQ=255)
    gCAGAAAGAGTCGGTGAATGGTGGTTGTTTACCGcc                                 <  1:4961612/36‑1 (MQ=255)
     cAGAAAGAGTCGGTGAATGGTGGGTGTTTACCGCCa                                >  1:6789882/1‑36 (MQ=255)
     aaGAAAGACTCGGTAAATGGGGTTTGTTTACCGCCa                                <  1:8387314/35‑1 (MQ=16)
      aGAAAGAGTCGGTGAATGGTGGTTGTTTACCGcccc                               >  1:1926440/1‑34 (MQ=255)
       gAAAGAGTCGGTGAATGGTGGTTGTTTACCGCCAct                              >  1:3750178/1‑36 (MQ=255)
          agagTCGGTGAATGGTGGTTGTTTACCGCCACTCtt                           >  1:6977950/1‑36 (MQ=255)
          agagTCGGTGAATGGTGGTTGTTTACCGCCACTCtt                           >  1:3978864/1‑36 (MQ=255)
            agTCGGGGAATGGTGGTTGTTTACCGCCAATAtttt                         >  1:7150077/1‑36 (MQ=21)
             gTCGGTGAATGGTGGTTGTTTACCGCCACTCttttt                        >  1:5259225/1‑36 (MQ=255)
               cGGTGAATGGTGGTTGTTTACCGCCACTCTTTTTGc                      >  1:1395103/1‑36 (MQ=255)
               cGGTGAATGGTGGTTGTTTACCGCCACTCTTTTTGc                      >  1:1550569/1‑36 (MQ=255)
               cGGTGAATGGTGGTTGTTTACCGCCACTCTTTTTGc                      >  1:7719761/1‑36 (MQ=255)
                ggTGAATGGTGGTTGTTTACCGCCACTCTTTTTGCa                     >  1:4943835/1‑36 (MQ=255)
                ggTGAATGGTGGTTGTTTACCGCCACTCTTTTTGCa                     >  1:4939417/1‑36 (MQ=255)
                ggTGAATGGTGGTTGTTTACCGCCACTCTTTTTGCa                     >  1:5714503/1‑36 (MQ=255)
                 gTGAATGGTGGTTGTTTACCGCCACTCTTTTTGCAg                    >  1:4965798/1‑36 (MQ=255)
                     aTGGTGGTTGTTTTCCGCCACTCTTTTTGCAGATGc                >  1:6825833/1‑36 (MQ=255)
                      tggtggTTGTTTACCGCCACTCTTTTTGCAGATGCg               >  1:7295084/1‑36 (MQ=255)
                       ggtggtTGTTTACCGCCACTCTTTTTGCAGATGCgg              >  1:6681138/1‑36 (MQ=255)
                        gtggtTGTTCCCCGCCACTCTTTTTGCAGATGCGGa             <  1:569001/36‑1 (MQ=37)
                         tgttttttccCCGCCACTCTTTTTGCAGATGCGGaa            >  1:801073‑M1/11‑36 (MQ=255)
                          tttttttccCCGCCACTCTTTTTGCAGATGCGGaaa           >  1:460327‑M1/10‑36 (MQ=255)
                          ggtTGTTTACCGCCACTCTTTTTGCAGATGCGGaaa           >  1:5620918/1‑36 (MQ=255)
                          ggtTGTTTACCGCCACTCTTTTTGCAGATGCGGaaa           >  1:4396984/1‑36 (MQ=255)
                           ttttttccCCGCCACTCTTTTTGCAGATGCGGAAAt          >  1:5792133‑M1/9‑36 (MQ=255)
                           gttgttTACCGCCACTCTTTTTGCAGATGCGGAAAt          >  1:2537739/1‑36 (MQ=255)
                            ttgttTACCGCCACTCTTTTTGCAGATGCGGAAATg         >  1:1287012/1‑36 (MQ=255)
                            ttgttTACCGCCACTCTTTTTGCAGATGCGGAAATg         >  1:462885/1‑36 (MQ=255)
                             ttttccCCGCCACTCTTTTTGCAGATGCGGAAATGa        >  1:7205936‑M1/7‑36 (MQ=255)
                              gttTACCGCCACTCTTTTTGCAGATGCGGAAATGAt       >  1:6139882/1‑36 (MQ=255)
                              gttTACCGCCACTCTTTTTGCAGATGCGGAAATGAt       >  1:8087555/1‑36 (MQ=255)
                               tttACCGCCACTCTTTTTGCAGATGCGGAAATGATa      >  1:4345209/1‑36 (MQ=255)
                                ttACCGCCACTCTTTTTGCAGATGCGGAAATGATaa     >  1:789278/1‑36 (MQ=255)
                                tccCCGCCACTCTTTTTGCAGATGCGGAAATGATaa     >  1:1045595‑M1/4‑36 (MQ=255)
                                tcaCCTCCACTCTTTTTGCATATTCGGAAATGATac     <  1:7369265/34‑2 (MQ=19)
                                 tACCGCCACTCTTTTTGCAGATGCGGAAATGATAAt    >  1:7201773/1‑36 (MQ=255)
                                 ccCCGCCACTCTTTTTGCAGATGCGGAAATGATAAt    >  1:8158424‑M1/3‑36 (MQ=255)
                                 ccCCGCCACTCTTTTTGCAGATGCGGAAATGATAAt    >  1:4747847‑M1/3‑36 (MQ=255)
                                  aCCGCCACTCTTTTTGCAGATGCGGAAATGATAATc   <  1:3047753/36‑1 (MQ=255)
                                  aCCGCCACTCTTTTTGCAGATGCGGAAATGATAATc   >  1:743450/1‑36 (MQ=255)
                                  aCCGCCACTCTTTTTGCAGATGCGGAAATGATAATc   >  1:2640777/1‑36 (MQ=255)
                                  aCCGCCACTCTTTTTGCAGATGCGGAAATGATAATc   >  1:6133408/1‑36 (MQ=255)
                                  aCCGCCACTCTTTTTGCAGATGCGGAAATGAGAATc   >  1:5365921/1‑36 (MQ=255)
                                   CCGCCACTCTTTTTGCAGATGCGGAAATGATAATCa  >  1:4264117‑M1/1‑36 (MQ=255)
                                   CCGCCACTCTTTTTGCAGATGCGGAAATGATAATCa  >  1:1468857‑M1/1‑36 (MQ=255)
                                   |                                   
GACGGCAGAAAGAGTCGGTGAATGGTGGTTGTTTACCGCCACTCTTTTTGCAGATGCGGAAATGATAATCA  >  REL606/556377‑556447

Base quality scores:  ATCG  < 3 ≤  ATCG  < 9 ≤  ATCG  < 20 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG