New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 1286705 =26 (0.540)8 (0.170) 4/56 NT 27.1% intergenic (+29/‑154) narI/ECB_01206 nitrate reductase 1, gamma (cytochrome b(NR)) subunit/hypothetical protein
?REL606 3092577 = 17 (0.280)coding (756/927 nt) yqhG hypothetical protein

TCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGA  >  REL606/1286670‑1286740
                                   |                                   
tggtcagcgtaaaacgaacatgaatgcgggcggaat                                     <  1:7079218‑M1/1‑1 (MQ=255)
tggtcagcgtaaaacgaacatgaatgcgggcggaat                                     >  1:6285585‑M1/36‑36 (MQ=255)
tCACTAAGCGAATTTTAGTTCACATAGACCCTGCtt                                     <  1:1913736/36‑1 (MQ=255)
tCACTAAGCGAATTTTAGTTCACATAGACCCTGCtt                                     <  1:5693240/36‑1 (MQ=255)
 ggtcagcgtaaaacgaacatgaatgcgggcggaaTc                                    >  1:768829‑M1/35‑36 (MQ=255)
 ggtcagcgtaaaacgaacatgaatgcgggcggaaTc                                    <  1:5892520‑M1/2‑1 (MQ=255)
  gtcagcgtaaaacgaacatgaatgcgggcggaaTcg                                   >  1:2044552‑M1/34‑36 (MQ=255)
  aCTTAGCGAATTTTAGTTCACATAGACCCTGCTTcg                                   <  1:3458169/36‑1 (MQ=255)
  aCTAAGCGAATTTTAGTTCACATAGACCCTGCTTcg                                   <  1:4440347/36‑1 (MQ=255)
  aCTAAGCGAATTTTAGTTCACATAGACCCTGCTTcg                                   >  1:4616512/1‑36 (MQ=255)
   tcagcgtaaaacgaacatgaatgcgggcggaaTcgg                                  <  1:3447958‑M1/4‑1 (MQ=255)
   tcagcgtaaaacgaacatgaatgcgggcggaaTcgg                                  <  1:3125037‑M1/4‑1 (MQ=255)
   tcagcgtaaaacgaacatgaatgcgggcggaaTcgg                                  >  1:7946372‑M1/33‑36 (MQ=255)
   cTAACCGAATTTTAGTTCACATAGACCCTGCTTcgg                                  >  1:7087427/1‑36 (MQ=255)
    cagcgtaaaacgaacatgaatgcgggcggaaTcggc                                 >  1:2562114‑M1/32‑36 (MQ=255)
    cagcgtaaaacgaacatgaatgcgggcggaaTcggc                                 >  1:3492753‑M1/32‑36 (MQ=255)
     agcgtaaaacgaacatgaatgcgggcggaatcggct                                >  1:5504185‑M1/31‑35 (MQ=255)
     aaGCGAATTTTAGTTCACATAGACCCTGCTTcggcg                                >  1:2838910/1‑36 (MQ=255)
     aaGCGAATTTTAGTTCACATAGACCCTGCTTcggcg                                >  1:5627908/1‑36 (MQ=255)
      gcgtaaaacgaacatgaatgcgggcggaaTcggcgg                               <  1:6402519‑M1/7‑1 (MQ=255)
      aGCGAATTTTAGTTCACATAGACCCTGCTTcggcgg                               <  1:270427/36‑1 (MQ=255)
        gtaaaacgaacatgaatgcgggcggaaTCGGCgggg                             <  1:5302233‑M1/9‑1 (MQ=255)
        cGAATTTTAGTTCACATAGACCCTGCTTCGGCgggg                             <  1:4171540/36‑1 (MQ=255)
        cGAATTTTAGTTCACATAGACCCTGCTTCGGCgggg                             <  1:4347837/36‑1 (MQ=255)
         gAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGt                            <  1:3165770/36‑1 (MQ=255)
           aaacgaacatgaatgcgggcggaaTCGGCGGGTttt                          <  1:153809‑M1/12‑1 (MQ=255)
           aaacgaacatgaatgcgggcggaaTCGGCGGGGttt                          <  1:1020722‑M1/12‑1 (MQ=255)
           aTTTTAGTTCACATAGACCCTGCTTCGGCGGGGttt                          >  1:4461375/1‑36 (MQ=255)
            ttttAGTTCACATAGACCCTGCTTCGGCGGGGtttt                         <  1:4287109/36‑1 (MQ=255)
             tttAGTTCACATAGACCCTGCTTCGGCGGGGttttt                        <  1:5445801/36‑1 (MQ=255)
             tttAGTTCACATAGACCCTGCTTCGGCGGGGttttt                        >  1:7633006/1‑36 (MQ=255)
             tttAGTTCACATAGACCCTGCTTCGGCGGGGttttt                        <  1:7337948/36‑1 (MQ=255)
             acgaacatgaatgcgggcggaaTCGGCGGGTttttt                        <  1:6360370‑M1/14‑1 (MQ=255)
             acgaacatgaatgcgggcggaaTCGGCGGGGttttt                        <  1:1354095‑M1/14‑1 (MQ=255)
              ttAGTTCACATAGACCCTGCTTCGGCGGGGtttttt                       >  1:5302626/1‑36 (MQ=255)
               tAGTTCACATAGACCCTGCTTCGGCGGGGttttttt                      >  1:4108454/1‑36 (MQ=255)
               gaacatgaatgcgggcggaaTCGGCGGGgttttttc                      <  1:7754492‑M1/16‑2 (MQ=255)
               gaacatgaatgcgggcggaaTCGGCGGGGttttttt                      <  1:2993200‑M1/16‑1 (MQ=255)
                aGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTa                     <  1:4579619/36‑1 (MQ=255)
                aGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTa                     >  1:1942530/1‑36 (MQ=255)
                aGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTa                     >  1:2463309/1‑36 (MQ=255)
                 gTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTAt                    <  1:7027970/36‑1 (MQ=255)
                  ttCACATAGACCCTGCTTCGGCGGGGTTTTTTTATg                   <  1:4423557/36‑1 (MQ=255)
                  catgaatgcgggcggaaTCGGCGgggtttttttgtg                   <  1:3935195‑M1/19‑4 (MQ=255)
                   tCACATAGACCCTGCTTCGGCGGGGTTTTTTTAttg                  >  1:3216507/1‑34 (MQ=255)
                   tCACATAGACCCTGCTTCGGCGGGGTTTTTTTATgg                  <  1:4594645/36‑1 (MQ=255)
                    cacaTAGACCCTGCTTCGGCGGGGTTTTTTTATggg                 <  1:3217442/36‑1 (MQ=255)
                      caTAGCCCCTGCTTCGGCGGGGTTTTTTTATGGGCa               <  1:7319404/36‑1 (MQ=255)
                       aTAGACCCTGCTTCGGCGGGGTTTTTTTATTGGCAc              >  1:3820814/1‑36 (MQ=255)
                        tAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACg             <  1:5504026/36‑1 (MQ=255)
                         aGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACgg            >  1:4637085/1‑36 (MQ=255)
                            cccTGCTTCGGCGGGGTTTTTTTTTTTGCACGGTgt         >  1:5666796/1‑35 (MQ=17)
                                gCTTCGGCGGGGTTTTTTTATTTTCACGGTGCgggg     >  1:625560/1‑36 (MQ=18)
                                 tttCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGt    <  1:7805197/35‑1 (MQ=255)
                                  ttCGGCGGGGTTTTTTTATCGGCACGGTGCGGGGTg   >  1:5591124/1‑36 (MQ=255)
                                   tCGGCGGGGGTTTTTTATTGGACCGGTTCGGGGGGa  >  1:2804621/1‑36 (MQ=9)
                                   |                                   
TCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGA  >  REL606/1286670‑1286740

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 22 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG