New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 553366 =10 (0.190)13 (0.250) 6/60 NT 44.1% coding (1526/2817 nt) ECB_00512 putative tail component of prophage
?REL606 4169224 = 23 (0.380)intergenic (+40/‑173) rpoC/htrC DNA‑directed RNA polymerase subunit beta'/heat shock protein

TAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATT  >  REL606/4169189‑4169258
                                   |                                  
tAATCCGCAAATAACTTAAAAACCCGCTTCGGCggg                                    <  1:6996673/36‑1 (MQ=255)
 gtcgggtcataaaccgcaagatgcggcgtggcgggg                                   >  1:686839‑M2/35‑35 (MQ=255)
 gtcggatcataaaccgcaagatgcggcgtggcggGt                                   >  1:4387317‑M2/35‑36 (MQ=255)
 aaTCCGCAAATAACGTAAAAACCCGCTTCGGCGGGt                                   >  1:3958501/1‑36 (MQ=255)
 aaTCCGCAAATAACATAAAAACCCGCTTCGGCGGGt                                   <  1:3956816/36‑1 (MQ=255)
  tcggatcataaaccgcaagatgcggcgtggcggggc                                  >  1:2730510‑M2/34‑34 (MQ=255)
  tcggatcataaaccgcaagatgcggcgtggcggGtt                                  >  1:1181699‑M2/34‑36 (MQ=255)
  aTCCGCAAATCACGTAAAAACCCGCTTCGGCGGGtt                                  <  1:5864990/36‑1 (MQ=255)
  aTCCGCAAATAACGTAAAAACCCGCTTCGGCGGGtt                                  <  1:7654807/36‑1 (MQ=255)
   cggatcataaaccgcaagatgcggcgtggcggGttt                                 >  1:6168679‑M2/33‑36 (MQ=255)
   cggatcataaaccgcaagatgcggcgtggcggGttt                                 >  1:5837054‑M2/33‑36 (MQ=255)
    ccGAAAATAACGTAAAAACCCGCTGCGGCGGGtttt                                >  1:7702982/1‑36 (MQ=17)
     gatcataaaccgcaagatgcggcgtggcgggtttct                               >  1:695112‑M2/31‑34 (MQ=255)
     gatcataaaccgcaagatgcggcgtggcggGttttt                               >  1:5480154‑M2/31‑36 (MQ=255)
     gatcataaaccgcaagatgcggcgtggcggGttttt                               >  1:5397397‑M2/31‑36 (MQ=255)
     gatcataaaccgcaagatgcggcgtggcggGttttt                               >  1:380547‑M2/31‑36 (MQ=255)
     cGCAAATAACGTAAAAACCCGCTTCGGCGGGttttt                               >  1:5840969/1‑36 (MQ=255)
      gCAAATAACGTAAAAACCCGCTTCGGCGGGtttttt                              >  1:6969209/1‑36 (MQ=255)
       tcataaaccgcaagatgcggcgtggcggGttttttt                             >  1:6241526‑M2/29‑36 (MQ=255)
        cataaaccgcaagatgcggcgtggcggTTTTTTTTa                            >  1:1259593‑M2/28‑36 (MQ=255)
        cataaaccgcaagatgcggcgtggcggGTTTTTTTa                            >  1:1411725‑M2/28‑36 (MQ=255)
        aaaTAACGTAAAAACCCGCTTCGGCGGGTTTTTTTa                            <  1:2651584/36‑1 (MQ=255)
         aaTAACGTAAAAACCCGCTTCGGCGGGttttttttt                           >  1:2140482/1‑34 (MQ=255)
          aTAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATg                          >  1:900176/1‑36 (MQ=255)
             aCGTAAAAACCCGCTTCGGCGGGTTTTTTTATgggg                       >  1:2537715/1‑36 (MQ=255)
                tAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGa                    <  1:2981891/36‑1 (MQ=255)
                 aaaaaCCCGCTTCGGCGGGTTTTTTTATGGGGGGAg                   >  1:7509454/1‑36 (MQ=255)
                   aaaCCCGCTTCGGCGGGTTTTTTTATGGGGGGAGtt                 <  1:4698938/36‑1 (MQ=255)
                    aaCCCGCTTCGGCGGGTTTTTTTTTGTGGGGTGttt                >  1:738488/1‑36 (MQ=14)
                    aaCACGCTTCGGCGGGTTTTTTTATGGGGGGAGttt                >  1:3156754/1‑36 (MQ=255)
                     aCCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTa               <  1:6825626/36‑1 (MQ=255)
                     aCCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTa               <  1:7768469/36‑1 (MQ=255)
                      cccGCTTCGGCGGGTTTTTTTATGGGGGGTGTTTAg              >  1:4172645/1‑36 (MQ=255)
                      cccGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAg              <  1:57905/36‑1 (MQ=255)
                          cTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGaa          <  1:2347359/36‑1 (MQ=255)
                            tCGGCGGGTTTTTTTATGGGTGGAGTTTAGGGAAag        <  1:1499383/36‑1 (MQ=255)
                              ggcggGTTTTTTTATGGGGGGAGTTTAGGGAAagag      <  1:4167616‑M2/31‑1 (MQ=255)
                               gcggGTTTTTTTATGTGGTGAGTTTAGGGACAGAGc     <  1:3825337/36‑1 (MQ=25)
                                 ggGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCAt   >  1:6943287‑M2/3‑36 (MQ=255)
                                  gGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCAtt  <  1:7329965‑M2/35‑1 (MQ=255)
                                  gGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCAtt  >  1:371049‑M2/2‑36 (MQ=255)
                                  gGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCAtt  >  1:2625228‑M2/2‑36 (MQ=255)
                                   |                                  
TAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATT  >  REL606/4169189‑4169258

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 21 ≤  ATCG  < 31 ≤  ATCG  < 37 ≤  ATCG