New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 306111726 (0.500)3 (0.060) 3/60 NT 8.7% coding (752/1173 nt) ECB_02854 conserved hypothetical protein
?REL606 = 4141016 37 (0.620)coding (89/705 nt) yijC DNA‑binding transcriptional repressor

CGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAAT  >  REL606/4140981‑4141051
                                   |                                   
cGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAc                                     <  1:2696614/36‑1 (MQ=255)
 gTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACc                                    <  1:3830225‑M2/36‑2 (MQ=255)
 gTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACc                                    <  1:6153664‑M2/36‑2 (MQ=255)
 gTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACc                                    >  1:301733‑M2/1‑35 (MQ=255)
  tGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACcc                                   <  1:3051194‑M2/36‑3 (MQ=255)
  tGATAGGCGTAAGAGCGCAACAAAAAGAAAAAACcc                                   >  1:7280637‑M2/1‑34 (MQ=255)
   gATGGGCGTAAGAGCGCAACAAAAAGAAAAAACccg                                  <  1:6895984‑M2/36‑4 (MQ=255)
    aTGGGCGTAAGAGCGCAACAAAAAGAAAAAAAccga                                 <  1:6208368‑M2/36‑5 (MQ=255)
      gggCGTAAGAGCGCAACAAAAAGAAAAAACccgccg                               <  1:5990074‑M2/36‑7 (MQ=255)
            aaGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCt                         >  1:3325303/1‑36 (MQ=255)
            aaGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCt                         >  1:6119711/1‑36 (MQ=255)
             agagCGCAACAAAAAGAAAAAACCCGCCGTTCGCtg                        >  1:1912004/1‑36 (MQ=255)
              gagCGCAACAAAAAGAAAAAACCCGCCGTTCGCtgg                       >  1:171745/1‑36 (MQ=255)
               agCGCAACAAAAAGAAAAAACCCGCCGTTCGCtggt                      >  1:3239305/1‑36 (MQ=255)
                gcgcAACAAAAAGAAAAAACCCGCCGTTCGCcggtg                     >  1:5140247/1‑36 (MQ=255)
                     acaaAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGc                >  1:816878/1‑36 (MQ=255)
                     acaaAAAGAAAAAACCCGCAGTTCGCTGGTGGAAGc                >  1:666657/1‑36 (MQ=255)
                      aaaaaaGAAAAACCCCGCGGTTCGCTGGTGGAAGcc               <  1:4902354/35‑1 (MQ=25)
                       aaaaaGAAAAAACCCGCCGTTCGCTGGTGGAAGCCg              >  1:7280252/1‑36 (MQ=255)
                       aaaaaAAAAAACCCCGCCGTTCGCTGGTGGAAGCCg              <  1:6014068/36‑1 (MQ=37)
                        aaaaGAAAAACCCCGCCGTTCGCTGGTGGAAGCCGc             <  1:4991079/36‑1 (MQ=255)
                        aaaaGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGc             >  1:4042973/1‑36 (MQ=255)
                        aaaaAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGc             <  1:4902753/36‑1 (MQ=37)
                        aaaaAAAAAAACccgccaatcgaggccagactttta             >  1:4749566‑M2/1‑12 (MQ=255)
                        aaaaAAAAAAACccgccaatattttacataacgcca             >  1:4587770‑M2/1‑12 (MQ=255)
                         aaaGAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCa            <  1:7539899/36‑1 (MQ=255)
                         aaaGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCa            >  1:97323/1‑36 (MQ=255)
                         aaaGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCa            >  1:286534/1‑36 (MQ=255)
                         aaaGAAAAAACCCGACGCTCGCTGGTGGAATCCGCa            >  1:7603257/1‑36 (MQ=25)
                         aaaAAAAAAACccgccgatggcccgcgcaagcggga            >  1:3839781‑M2/1‑11 (MQ=255)
                         aaaAAAAAAACccgccaatattttacataacgccac            >  1:5391464‑M2/1‑11 (MQ=255)
                          aaGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCAt           >  1:3653493/1‑36 (MQ=255)
                          aaAAAAAAACccgccaatcgaggccagacttttact           >  1:5372163‑M2/1‑10 (MQ=255)
                           aGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCAtt          <  1:1353354/36‑1 (MQ=255)
                           aGAAAAAACCCGCCGTTCGCTGGGGGAAGCCGCAtt          >  1:6939911/1‑36 (MQ=255)
                            gAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCAttt         >  1:5027081/1‑36 (MQ=255)
                            gAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCAttt         >  1:245752/1‑36 (MQ=255)
                             aacaaaCCCGCCGTTCGCTGGTGGAAGCCGCATTTa        >  1:3552447/4‑36 (MQ=255)
                             aaaaaaCCCGCCGTTCGCTGGTGGAAGCCGCATTTa        <  1:4297595/36‑1 (MQ=255)
                               caaacccgccgatggcccgcgcaagcgggatcaggt      >  1:7487252‑M2/2‑5 (MQ=255)
                               acaacccgccgatggcccgcgcaagcgggatcaggt      >  1:3739923‑M2/3‑5 (MQ=255)
                               aaacCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGc      <  1:2481368/36‑1 (MQ=255)
                               aaacCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGc      <  1:4557850/36‑1 (MQ=255)
                               aaaaCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGc      >  1:371666/1‑36 (MQ=255)
                               aaaaCCCGCCGTTCGCTGGTGGAAGCCGAATTTAGc      >  1:3475725/1‑36 (MQ=255)
                                aaaCccgccaatcgaggccagacttttactgaaggt     >  1:6439791‑M2/1‑4 (MQ=255)
                                aaaCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGcc     >  1:2994291/1‑36 (MQ=255)
                                aaaCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGcc     >  1:1910805/1‑36 (MQ=255)
                                 aaCccgccgatggcccgcgcaagcgggatcaggtaa    >  1:2608735‑M2/1‑3 (MQ=255)
                                 aaCccgccaatattttacataacgccaccgagtaca    >  1:5249996‑M2/1‑3 (MQ=255)
                                 aaCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCa    <  1:7628911/36‑1 (MQ=255)
                                 aaCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCa    >  1:7875305/1‑36 (MQ=255)
                                  aCccgccgatggcccgcgcaagcgggatcaggtaag   >  1:137695‑M2/1‑2 (MQ=255)
                                  aCccgccaatattttacataacgccaccgagtacag   >  1:7163836‑M2/1‑2 (MQ=255)
                                  aCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCaa   >  1:6610445/1‑36 (MQ=255)
                                  aCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCaa   >  1:5472244/1‑36 (MQ=255)
                                  aCCCGCCGTTCGCTGGCGGAAGCCGCATTTAGCCaa   >  1:7256204/1‑36 (MQ=255)
                                   cccgccgatggcccgcgcaagcgggatcaggtaagg  >  1:6106475‑M2/1‑1 (MQ=255)
                                   cccGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAAt  <  1:6286732/36‑1 (MQ=255)
                                   |                                   
CGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAAT  >  REL606/4140981‑4141051

Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 19 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG