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

TACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAA  >  REL606/1802303‑1802372
                                   |                                  
tACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGAt                                    >  1:6095162/1‑36 (MQ=255)
 aCAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATc                                   >  1:7483429/1‑36 (MQ=255)
  cAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATcc                                  >  1:2757314/1‑36 (MQ=255)
   aaCTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCt                                 >  1:3043507/1‑36 (MQ=255)
    aCTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTg                                >  1:261274/1‑36 (MQ=255)
      tGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTc                              >  1:7472275/1‑36 (MQ=255)
       ggttGAGAAGAGAAAAGTAAACCGACGATCTTGTcc                             >  1:2631270/3‑36 (MQ=21)
          tGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCAcc                          >  1:7956733/1‑36 (MQ=255)
            agaagagaaAAGAAAACCGCCGATCCTGTCCACCGc                        >  1:831865/1‑36 (MQ=255)
               agagaaAAGAAAACCGCCGATCCTGTCCACCGCAtt                     >  1:7578609/1‑36 (MQ=255)
               agagaaAAGAAAACCGCCGATCCTGTCCACCGCAtt                     >  1:3592009/1‑36 (MQ=255)
               agagaaAAGAAAACCGCCGATCCTGTCCACCGCAtt                     >  1:6382252/1‑36 (MQ=255)
                 aaaaaaaaaccccgccgaTCCTGTCCACCGCATTAc                   <  1:4053813‑M2/18‑1 (MQ=255)
                   aaaagaaaaCCGCCGATCCTGTCCACCGCATTACTg                 >  1:8165824/1‑36 (MQ=255)
                   aaaagaaaaCCGCCGATCCGGTACGCCGCGTTATTg                 >  1:4663697/1‑36 (MQ=12)
                   aaaaaatccccgccgaTCCTGTCCACCGCATTACTg                 <  1:3613825‑M2/20‑1 (MQ=255)
                    aaaaaaccccgccgaTCCTGTCCACCGCATTACTGc                <  1:6922917‑M2/21‑1 (MQ=255)
                    aaaaaaacccgccgaTCCTGTCCACCGCATTACTGc                <  1:5226555‑M2/21‑1 (MQ=255)
                      agaaaaCCGCCGATCCTGGCCACCGCATTACTGCaa              >  1:4326351/1‑36 (MQ=255)
                      aaacccccgccgaTCCTGTCCACCGCATTACTGCaa              <  1:5839144‑M2/23‑1 (MQ=255)
                       gaaaaCCGCCGATCCTGTCCACCGCATTACTGCaaa             >  1:2991258/1‑35 (MQ=255)
                       gaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAg             >  1:768196/1‑36 (MQ=255)
                        aaaaCCGCCGATCCTGTCCACCGCATTACTGCAAgg            >  1:8170503/1‑36 (MQ=255)
                         accccgccgaTCCTGTCCACCGCATTACTGCAAGGt           <  1:5380101‑M2/26‑1 (MQ=255)
                          ccccgccgaTCCTGTCCACCGCATTACTGCAAGGTa          <  1:3747616‑M2/27‑1 (MQ=255)
                          aaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTa          >  1:4416852/1‑36 (MQ=255)
                           cccgccgaTCCTGTCCACCGCATTACTGCAAGGTAg         <  1:8307634‑M2/28‑1 (MQ=255)
                            ccgccgaTCCTGTCCACCGCATTACTGCAAGGTAGt        <  1:1414916‑M2/29‑1 (MQ=255)
                              gccgaTCCTGTCCACCGCATTACTGCAAGGTAGTgg      <  1:1295553‑M2/31‑1 (MQ=255)
                              gccgaTCCTGTCCACCGCATTACTGCAAGGTAGTgg      >  1:3329560‑M2/6‑36 (MQ=255)
                              gccgaTCCTGTCCACCGCATTACTGCAAGGTAGTgg      >  1:5474912‑M2/6‑36 (MQ=255)
                               ccgaTCCTGTCCACCGCATTACTGCAAGGTAGTGGa     <  1:2891127‑M2/32‑1 (MQ=255)
                               ccgaTCCTGTCCACCGCATTACTGCAAGGTAGTGGa     <  1:6085655‑M2/32‑1 (MQ=255)
                               ccgaTCCTGTCCACCGCATTACTGCAAGGTAGTGGa     <  1:5497094‑M2/32‑1 (MQ=255)
                                cgaTCCTGTCCACCGCATTACTGCAAGGTATTGGAc    <  1:823881‑M2/33‑1 (MQ=255)
                                cgaTCCTGTCCACCGCATTACTGCAAGGTAGTGGAc    <  1:2234033‑M2/33‑1 (MQ=255)
                                 gaTCCTGTCCACCGCATTACTGCAAGGTAGTGGACa   >  1:2735555‑M2/3‑36 (MQ=255)
                                  aTCCTGTCCACCGCATTACTGCAAGGTAGTGGACaa  <  1:3749610‑M2/35‑1 (MQ=255)
                                  aTCCTGTCCACCGCATTACTGCAAGGTAGTGGACaa  >  1:8163698‑M2/2‑36 (MQ=255)
                                  aTCCTGTCCACCGCATTACTGCAAGGTAGTGGACaa  <  1:562560‑M2/35‑1 (MQ=255)
                                  aTCCTGTCCACCGCATTACTGCAAGGTAGTGGACaa  <  1:4965785‑M2/35‑1 (MQ=255)
                                   |                                  
TACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAA  >  REL606/1802303‑1802372

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