New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 1286705 =35 (0.710)8 (0.160) 6/56 NT 25.0% intergenic (+29/‑154) narI/ECB_01206 nitrate reductase 1, gamma (cytochrome b(NR)) subunit/hypothetical protein
?REL606 1802338 = 13 (0.210)intergenic (‑136/+67) ydjR/spy hypothetical protein/envelope stress induced periplasmic protein

TCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGG  >  REL606/1286670‑1286735
                                   |                              
ttgtccactaccttgcagtaatgcggtggacaggat                                >  1:8163698‑M1/36‑36 (MQ=255)
ttgtccactaccttgcagtaatgcggtggacaggat                                <  1:3749610‑M1/1‑1 (MQ=255)
ttgtccactaccttgcagtaatgcggtggacaggat                                <  1:4965785‑M1/1‑1 (MQ=255)
ttgtccactaccttgcagtaatgcggtggacaggat                                <  1:562560‑M1/1‑1 (MQ=255)
 tgtccactaccttgcagtaatgcggtggacaggaTc                               >  1:2735555‑M1/35‑36 (MQ=255)
 cACTAAGCGAATTTTAGTTCACATAGACACTGCTTc                               <  1:3136583/36‑1 (MQ=255)
 cACTAAGCGAATTTTAGTTCACATAGACACTGCTTc                               <  1:6743314/36‑1 (MQ=255)
 cACTAAGCGAATTTTAGTTCACATAGACACTGCTTc                               <  1:1693917/36‑1 (MQ=255)
 cACTAAGCGAATTTTAGTTCACATAGACACTGCTTc                               <  1:5892558/36‑1 (MQ=255)
  gtccactaccttgcagtaatgcggtggacaggaTcg                              <  1:2234033‑M1/3‑1 (MQ=255)
  gtccaataccttgcagtaatgcggtggacaggaTcg                              <  1:823881‑M1/3‑1 (MQ=255)
  aCTAAGTGAATTTTAGTTCACATAGACACTGCTTcg                              <  1:6869881/36‑1 (MQ=38)
   tccactaccttgcagtaatgcggtggacaggaTcgg                             <  1:6085655‑M1/4‑1 (MQ=255)
   tccactaccttgcagtaatgcggtggacaggaTcgg                             <  1:2891127‑M1/4‑1 (MQ=255)
   tccactaccttgcagtaatgcggtggacaggaTcgg                             <  1:5497094‑M1/4‑1 (MQ=255)
   cTAAGCGAATTTTAGTTCACATAGACACTGCTTcgg                             <  1:3689718/36‑1 (MQ=255)
    tAAGCGAATTTTAGTTCACATAGACACTGCTTcggc                            >  1:4259560/1‑36 (MQ=255)
    tAAGCGAATTTTAGTTCACATAGACACTGCTTcggc                            <  1:7922416/36‑1 (MQ=255)
    ccactaccttgcagtaatgcggtggacaggaTcggc                            >  1:5474912‑M1/32‑36 (MQ=255)
    ccactaccttgcagtaatgcggtggacaggaTcggc                            <  1:1295553‑M1/5‑1 (MQ=255)
    ccactaccttgcagtaatgcggtggacaggaTcggc                            >  1:3329560‑M1/32‑36 (MQ=255)
     aaGCGAATTTTAGTTCACATAGACACTGCTTcggcg                           >  1:6208847/1‑36 (MQ=255)
     aaGCGAATTTTAGTTCACATAGACACTGCTTcggcg                           <  1:1817074/36‑1 (MQ=255)
      actaccttgcagtaatgcggtggacaggaTcggcgg                          <  1:1414916‑M1/7‑1 (MQ=255)
      aGCGAATTTTAGTTCACATAGACCCTGCTTcggcgg                          >  1:4259022/1‑36 (MQ=255)
       gCGAATTTTAGTTCACATAGACACTGCTTCGGCggg                         <  1:6530717/36‑1 (MQ=255)
       gCGAATTTTAGTTCACATAGACACTGCTTCGGCggg                         >  1:1468870/1‑36 (MQ=255)
       ctaccttgcagtaatgcggtggacaggaTCGGCggg                         <  1:8307634‑M1/8‑1 (MQ=255)
        taccttgcagtaatgcggtggacaggaTCGGCgggg                        <  1:3747616‑M1/9‑1 (MQ=255)
        cGAATTTTAGTTCACATAGACCCTGCTTCGGCgggg                        <  1:4898156/36‑1 (MQ=255)
         gAATTTTAGTTCACATAGACACTGCTTCGGCGGGGt                       <  1:2517595/36‑1 (MQ=255)
         accttgcagtaatgcggtggacaggaTCGGCGGGGt                       <  1:5380101‑M1/10‑1 (MQ=255)
          aaTTTTAGTTCACATAGACACTGCTTCGGCGGGGtt                      <  1:7665839/36‑1 (MQ=255)
          aaTTTTAGTTCACATAGACACTGCTTCGGCGGGGtt                      <  1:4755305/36‑1 (MQ=255)
           aTTTTAGTTCACATAGACACTGCTTCGGCGGGGttt                     <  1:2952910/36‑1 (MQ=255)
           aTTTTAGTTCACATAGACACTGCTTCGGCGGGGttt                     >  1:1082475/1‑36 (MQ=255)
            ttgcagtaatgcggtggacaggaTCGGCGGGGgttt                    <  1:5839144‑M1/13‑1 (MQ=255)
             tttAGTTCACATAGACCCTGCTTCGGCGGGGttttt                   >  1:3178623/1‑36 (MQ=255)
             tttAGTTCACATAGACACTGCTTCGGCGGGGttttt                   <  1:7409515/36‑1 (MQ=255)
              ttAGTTCACATAGACCCTGCTTCGGCGGGGtttttt                  >  1:3950724/1‑36 (MQ=255)
              ttAGTTCACATAGACACTGCTTCGGCGGGGtttttt                  >  1:7448567/1‑36 (MQ=255)
              gcagtaatgcggtggacaggaTCGGCGGGTtttttt                  <  1:5226555‑M1/15‑1 (MQ=255)
              gcagtaatgcggtggacaggaTCGGCGGGGtttttt                  <  1:6922917‑M1/15‑1 (MQ=255)
               tAGTTCACATAGACACTGCTTCGGCGGGGttttttt                 <  1:8157667/36‑1 (MQ=255)
               tAGTTCACATAGACACTGCTTCGGCGGGGttttttt                 <  1:1068933/36‑1 (MQ=255)
               cagtaatgcggtggacaggaTCGGCGGGGatttttt                 <  1:3613825‑M1/16‑1 (MQ=255)
                 gtaatgcggtggacaggaTCGGCGGggttttttttt               <  1:4053813‑M1/18‑3 (MQ=255)
                  ttCACATAGACACTGCTTCGGCGGGGTTTTTTTAtt              >  1:4642612/1‑35 (MQ=255)
                  ttCACATAGACACTGCTTCGGCGGGGTTTTTTTAtt              >  1:1378409/1‑35 (MQ=255)
                  ttCACATAGACACTGCTGCGGCGGGGTTTTTTTATg              <  1:4299734/36‑1 (MQ=37)
                   tCACATAGACACTGCTTCGGCGGGGTTTTTTTATgg             <  1:2489385/36‑1 (MQ=255)
                    cacaTAGACACTGCTTCGGCGGGGTTTTTTTATggg            <  1:960377/36‑1 (MQ=255)
                     acaTAGACACTGCTTCGGCGGGGTTTTTTTATTGGc           >  1:1421669/1‑36 (MQ=25)
                     acaTAGACACTGCTTCGGCGGGGTTTTTTTATGGGc           <  1:4966621/36‑1 (MQ=255)
                      caTAGACACTGCTTCGGCGGGGTTTTTTTATTGGCa          >  1:7308201/1‑36 (MQ=25)
                      caTAGACACTGCTTCGGCGGGGTTTTTTTATGGGCa          <  1:7870263/36‑1 (MQ=255)
                      caTAGACACTGCTTCGGCGGGGTTTTTTTATGGGCa          >  1:481629/1‑36 (MQ=255)
                       aTAGACACTGCTTCGGCGGGGTTTTTTTATTTGCAc         >  1:8169907/1‑36 (MQ=18)
                       aTAGACACTGCTTCGGCGGGGTTTTTTTATTGGAAc         >  1:2613685/1‑36 (MQ=18)
                        tAGACACTGCTTCGGCGGGGTTTTTTTTTGGGCagg        >  1:2482278/1‑34 (MQ=21)
                        tAGACACTGCTTCGGCGGGGTTTTTTTATTGGCACg        >  1:7961363/1‑36 (MQ=25)
                          gACACTGCTTCGGCGGGGTTTTTTTATGGGCACGGt      <  1:6212519/36‑1 (MQ=255)
                          gACACTGCTTCGGCGGGGTTTTTTTATGGGCACGGt      <  1:5604073/36‑1 (MQ=255)
                          gACACTGCTTCGGCGGGGTTTTTTTATGGGCACGGt      <  1:4831082/36‑1 (MQ=255)
                          gACACTGCTTCGGCGGGGTTTTTTTATGGGCACGGt      <  1:350876/36‑1 (MQ=255)
                            cacTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGc    <  1:4058646/34‑1 (MQ=255)
                             gcTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCg   <  1:7642699/35‑1 (MQ=255)
                              cTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCgg  >  1:7394097/1‑36 (MQ=255)
                                   |                              
TCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGG  >  REL606/1286670‑1286735

Base quality scores:  ATCG  < 3 ≤  ATCG  < 13 ≤  ATCG  < 24 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG