New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 1286707 =35 (0.670)6 (0.120) 3/60 NT 16.0% intergenic (+31/‑152) narI/ECB_01206 nitrate reductase 1, gamma (cytochrome b(NR)) subunit/hypothetical protein
?REL606 = 3571661 28 (0.460)coding (110/1470 nt) yhiP predicted transporter

GGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTC  >  REL606/3571626‑3571694
                                   |                                 
ggAGCGATTCGGCTACTACGGCGTGCAGGGCGTACt                                   <  1:602176/36‑1 (MQ=255)
  aGCGATTCGGCTACTACGGCGTGCAGGGCGTACTgg                                 <  1:5971273‑M2/36‑3 (MQ=255)
  aGCGATTCGGCTACTACGGCGTGCAGGGCGTACTgg                                 <  1:2070640‑M2/36‑3 (MQ=255)
  aGCGATTCGGCTACTACGGCGTGCAGGGCGTACTgg                                 <  1:2295332‑M2/36‑3 (MQ=255)
  aGCGATTCGGCTACTACGGCGTGCAGGGCGTACTgg                                 <  1:6138851‑M2/36‑3 (MQ=255)
   tcGATTCGGCTACTACGGCTTGCAGGGCGTACTGGc                                <  1:2168409/35‑1 (MQ=37)
   gCGATTCGGCTACTGCGGCGTGCAGGGCGTACTggc                                >  1:2648474‑M2/1‑33 (MQ=255)
   gCGATTCGGCTACTACGGCGTGCAGGGCGTACTggc                                >  1:2656143‑M2/1‑33 (MQ=255)
   gCGATTCGGCTACTACGGCGTGCAGGGCGTAATggc                                <  1:846079‑M2/36‑4 (MQ=255)
    cGATTCGGCTACTACGGCGTGCAGGGCGTACTggcg                               <  1:1268154‑M2/36‑5 (MQ=255)
     gATTCGGCTACTACGGCGTGCAGGGCGTACTggcgg                              <  1:2309910‑M2/36‑6 (MQ=255)
     gATTCGGCTACTACGGCGTGCAGGGCGTACTggcgg                              >  1:5452231‑M2/1‑31 (MQ=255)
     gATTCGGCTACTACGGCGTGCAGGGCGTACTggcgg                              >  1:4719901‑M2/1‑31 (MQ=255)
      aTTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGt                             <  1:5657538/36‑1 (MQ=255)
        tCGGCTACTACGGCGTGCAGGGCGTACTGGCGGttt                           <  1:5973031/36‑1 (MQ=255)
         cGGCTACTACGGCGTGCAGGGCGTACTGGCGGtttt                          <  1:2306141/36‑1 (MQ=255)
          ggCTACTACTGCGTGCAGGGCGTACTGGCGGTTttc                         <  1:7133567/36‑1 (MQ=255)
          ggCTACTACGGCGTGCAGGGCGTACTggcgggtttt                         <  1:2254335‑M2/36‑11 (MQ=255)
          ggCTACTACGGCGTGCAGGGCGTACTggcggggttt                         <  1:2927501‑M2/36‑11 (MQ=255)
           gCTACTACGGCGTGCAGGGCGTACTGGCGGTTttct                        <  1:5074108/36‑1 (MQ=255)
            cttctacGGCGTGCAGGGCGTACTGGCGGTTttctt                       <  1:897331/33‑1 (MQ=255)
             tactacGGCGTGCAGGGCGTACTggcggggtttttt                      <  1:6188374‑M2/36‑14 (MQ=255)
             tactacGGCGTGCAGGGCGTACTGGCGGTTttcttc                      <  1:537724/36‑1 (MQ=255)
              actacGGCGTGCAGGGCGTACTggcggggttttttt                     <  1:4072029‑M2/36‑15 (MQ=255)
               ctacGGCGTGCAGGGCGTACTggcgggttttttttt                    <  1:2116651‑M2/36‑16 (MQ=255)
               ctacGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGt                    <  1:4907518/36‑1 (MQ=255)
               caacGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGt                    <  1:2621753/34‑1 (MQ=255)
                tacGGCGTGCAGGGCGTACTggcggggttttttttt                   <  1:1880153‑M2/36‑17 (MQ=255)
                tacGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGtt                   <  1:2456522/36‑1 (MQ=255)
                tacGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGtt                   <  1:2797956/36‑1 (MQ=255)
                gacGGCGTGCGGGGCGTACTGGCGGTTTTCTTCGtt                   <  1:3333493/35‑1 (MQ=255)
                 acGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTa                  <  1:2347354/36‑1 (MQ=255)
                 acGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTa                  <  1:2663141/36‑1 (MQ=255)
                 acGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTa                  <  1:2752198/36‑1 (MQ=255)
                  cGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTaa                 <  1:7607606/36‑1 (MQ=255)
                    gCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAAc               <  1:3896069/36‑1 (MQ=255)
                     cGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACa              <  1:3860833/36‑1 (MQ=255)
                       tGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGc            <  1:2734891/36‑1 (MQ=255)
                        gCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCt           <  1:5331912/36‑1 (MQ=255)
                        gCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCt           <  1:7339025/36‑1 (MQ=255)
                          aGGGCGTACTGGCGGTTTTGTTCGTTAAACAGCTTg         <  1:5634314/36‑1 (MQ=255)
                           gggCGTACTGGCGGTTTTCTTCGTTAAACAGCTTgg        <  1:3080190/36‑1 (MQ=255)
                              cGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGAtt     <  1:4402970/36‑1 (MQ=255)
                              cGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGAtt     <  1:3501613/36‑1 (MQ=255)
                                tACTGGCGGTTTTCTTCGTTCAACAGCTTGGATtct   <  1:6049393/36‑1 (MQ=255)
                                 aCTGGCGGTTTTCTTCGTTAAACAGCTTGGATtctc  <  1:7418066/36‑1 (MQ=255)
                                   |                                 
GGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTC  >  REL606/3571626‑3571694

Base quality scores:  ATCG  < 3 ≤  ATCG  < 12 ≤  ATCG  < 24 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG