New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 1286706 =15 (0.400)9 (0.240) 4/58 NT 32.7% intergenic (+30/‑153) narI/ECB_01206 nitrate reductase 1, gamma (cytochrome b(NR)) subunit/hypothetical protein
?REL606 = 2780684 22 (0.490)intergenic (‑286/‑30) cysJ/ygcM sulfite reductase, alpha subunit, flavoprotein/6‑pyruvoyl tetrahydrobiopterin synthase (PTPS)

TAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGAT  >  REL606/2780649‑2780718
                                   |                                  
tAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTa                                    <  1:5076309/36‑1 (MQ=255)
 aGATAATGCGAAAAAACAGCCTTTCCGGTACTCTAc                                   <  1:89882‑M2/36‑2 (MQ=255)
 aGATAATGCGAAAAAACAGCCTTTCCGGTACTCTAc                                   >  1:4984889‑M2/1‑35 (MQ=255)
 aGATAATGCGAAAAAACAGCATTTCCGGTACTCTAc                                   <  1:1461148‑M2/36‑2 (MQ=255)
  gATAATGCGAAAAAACAGCCTTTCCGGTACTCTAcg                                  <  1:1368590‑M2/36‑3 (MQ=255)
   aTAATGCTAAAACACAGCCTTTCCGGTACTCTAcgg                                 <  1:2670430/36‑1 (MQ=37)
     aaTGCGAAAAAACAGCCTTTCCGGTACTCTAcggcg                               >  1:5724871‑M2/1‑31 (MQ=255)
     aaTGCGAAAAAACAGCCTTTCCGGTACTCTAcggcg                               <  1:2845297‑M2/36‑6 (MQ=255)
      aTGCGAAAAAACAGCCTTTCCGGAACTCTAaggcgg                              >  1:1336736/1‑36 (MQ=37)
       tGCGAAAAAACAGCCTTTCCGGTACTCTAcggcggg                             <  1:4737515‑M2/36‑8 (MQ=255)
       tGCGAAAAAACAGCCTTTCCGGTACTCTAcggcggg                             <  1:4603355‑M2/36‑8 (MQ=255)
       tGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGt                             <  1:4954137/36‑1 (MQ=255)
          gAAAAAACAGCCTTTCCTTTACTCTACGGCGGtttt                          <  1:4293882/36‑1 (MQ=37)
          gAAAAAACAGCCTTTCCGGTACTCTAcggcgggttt                          <  1:132988‑M2/36‑11 (MQ=255)
          gAAAAAACAGCCTTTCCGGTACTCTAcggcgggctt                          <  1:4749798‑M2/36‑11 (MQ=255)
           aaaaaaCAGCCTTTCCGGTACTCTAcggcggggttt                         <  1:773203‑M2/36‑12 (MQ=255)
            aaaaaCATCCTTTCCGGTACTCTACGGCGGTTTTAt                        <  1:1336838/36‑1 (MQ=255)
            aaaaaCAGCCTTTCCGGTACTCTACGGCGGattttt                        >  1:5855087/1‑34 (MQ=38)
            aaaaaCAGCCTTTCCGGTACTCTACGGCGGTTTTAt                        <  1:2711561/36‑1 (MQ=255)
             aaaaCAGCCTTTCCGGTACTCTAcggcgggtttttg                       <  1:5653459‑M2/36‑14 (MQ=255)
              aaaCAGCCTTTCCGGTACTCTAcggcggggtttttt                      <  1:1388866‑M2/36‑15 (MQ=255)
              aaaCAGCCTTTCCGGTACTCTACGGCGGTTTTATCg                      <  1:369760/36‑1 (MQ=255)
               aaCAGCCTTTCCGGTACTCTAcggcggggttttttt                     <  1:4675942‑M2/36‑16 (MQ=255)
               aaCAGCCTTTCCGGTACTCTACGGCGGTTTTATCGt                     <  1:171923/36‑1 (MQ=255)
                aCAGCCTTTCCGGTACTCTAcggcgggttttttttt                    <  1:2890449‑M2/36‑17 (MQ=255)
                 cAGTCTTTACTGTACTCTACGGCGGTTTTATCGTcc                   <  1:514731/36‑1 (MQ=21)
                 cAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTcc                   <  1:5275875/36‑1 (MQ=255)
                  aGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCt                  <  1:3045274/36‑1 (MQ=255)
                    ccTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGt                <  1:3290836/36‑1 (MQ=255)
                       ttCCGGTACTCTACGGCGGTTTTATCGTCCTGTaga             <  1:1960074/36‑1 (MQ=255)
                        tCCGGTACTCTACGGCGGTTTTATCGTCCTGTagag            <  1:3068904/36‑1 (MQ=255)
                           ggTACTCTACGGCGGTTTTATCGTCCTGTAGAGaaa         <  1:3120873/36‑1 (MQ=255)
                            gTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAAt        <  1:3251342/36‑1 (MQ=255)
                             tcctctACGGCGGTTTTATCTTCCTGTAGAGAAAtt       <  1:5323044/34‑1 (MQ=37)
                             tACTCTACGGCGGTTTTATCGTCCTGTAGAGAAAtt       <  1:1101455/36‑1 (MQ=255)
                                  tACGGCGGTTTTATCGTCCTGTAGAGAAATTatgat  <  1:3562574/36‑1 (MQ=255)
                                  tACGGCGGTTTTATCGTCCTGTAGAGAAATTatgat  <  1:2146379/36‑1 (MQ=255)
                                  tACGGCGGTTTTATCGTCCTGTAGAGAAATTatgat  <  1:876285/36‑1 (MQ=255)
                                   |                                  
TAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGAT  >  REL606/2780649‑2780718

Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 20 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG