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)

CACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGG  >  REL606/1286671‑1286737
                                   |                               
cACTAAGCGAATTTTAGTTCACATAGACACTGCTTc                                 <  1:1626208/36‑1 (MQ=255)
agataatgcgaaaaaacagcctttccggtactctac                                 <  1:89882‑M1/1‑1 (MQ=255)
agataatgcgaaaaaacagcctttccggtactctac                                 >  1:4984889‑M1/36‑36 (MQ=255)
agataatgcgaaaaaacagcatttccggtactctac                                 <  1:1461148‑M1/1‑1 (MQ=255)
 gataatgcgaaaaaacagcctttccggtactctacg                                <  1:1368590‑M1/2‑1 (MQ=255)
    aatgcgaaaaaacagcctttccggtactctacggcg                             >  1:5724871‑M1/32‑36 (MQ=255)
    aatgcgaaaaaacagcctttccggtactctacggcg                             <  1:2845297‑M1/5‑1 (MQ=255)
    aaGCGAATTTTAGTTCACATAGACACTGCTTcggcg                             <  1:2479826/36‑1 (MQ=255)
      tgcgaaaaaacagcctttccggtactctaCGGCggg                           <  1:4737515‑M1/7‑1 (MQ=255)
      tgcgaaaaaacagcctttccggtactctaCGGCggg                           <  1:4603355‑M1/7‑1 (MQ=255)
       cGAATTTTAGTTCACATAGACACTGCTTCGGCgggg                          <  1:4748724/36‑1 (MQ=255)
       cGAATTTTAGTTCACATAGACACTGCTTCGGCgggg                          <  1:3594632/36‑1 (MQ=255)
        gAATTTTAGTTCACATAGACACTGCTTCGGCGGGGt                         <  1:3241774/36‑1 (MQ=255)
         gaaaaaacagcctttccggtactctaCggcgggttt                        <  1:132988‑M1/10‑4 (MQ=255)
         gaaaaaacagcctttccggtactctaCggcgggctt                        <  1:4749798‑M1/10‑4 (MQ=255)
          aaaaaacagcctttccggtactctaCGGCGGGGttt                       <  1:773203‑M1/11‑1 (MQ=255)
          aTTTTAGTTCACATAGACACTGCTTCGGCGGGGttt                       >  1:1265691/1‑36 (MQ=255)
          aTTTTAGTTCACATAGACACTGCTTCGGCGGGGttt                       >  1:4563735/1‑36 (MQ=255)
            tttAGTTCACATAGACACTGCTTCGGCGGGGttttt                     >  1:38932/1‑36 (MQ=255)
            aaaacagcctttccggtactctaCGGCGGGtttttg                     <  1:5653459‑M1/13‑2 (MQ=255)
             aaacagcctttccggtactctaCGGCGGGGtttttt                    <  1:1388866‑M1/14‑1 (MQ=255)
              aacagcctttccggtactctaCGGCGGGGttttttt                   <  1:4675942‑M1/15‑1 (MQ=255)
               acagcctttccggtactctaCGGCGGGttttttttt                  <  1:2890449‑M1/16‑2 (MQ=255)
               aGTTCACATAGACACTGCTTCGGCGGGGTTTTTTTa                  <  1:2731140/36‑1 (MQ=255)
                 ttCACATAGACACTGCTTCGGCGGGGTTTTTTTATg                <  1:3405869/36‑1 (MQ=255)
                      aTAGACACTGCTTCGGCGGGGTTTTTTTATGGGCAc           <  1:4521072/36‑1 (MQ=255)
                         gACACTGCTTCGGCGGGGTTTTTTTATGGGCACGGt        <  1:5298025/36‑1 (MQ=255)
                         gACACTGCTTCGGCGGGGTTTTTTTATGGGCACGGt        <  1:2603158/36‑1 (MQ=255)
                          acacTGCTTCGGCGGGGTTTTTTTATTGGCACGGTg       >  1:4885903/4‑36 (MQ=25)
                               gCTTCGGCGGGGGTTTTTTATGTGCACGGTGCgggg  >  1:2138923/1‑36 (MQ=37)
                               gCTTCGGCGGGGGTTTTTTATCGGCACGGTGTggtg  >  1:2969225/1‑34 (MQ=18)
                                   |                               
CACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGG  >  REL606/1286671‑1286737

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