New junction evidence
  seq id position reads (cov) reads (cov) score skew annotation gene product
* ? REL606 72407 =24 (0.510)5 (0.110) 4/60 5.3 coding (446/1701 nt) araB ribulokinase
?REL606 1286712 = 49 (0.890)intergenic (+36/‑147) narI/ECB_01206 nitrate reductase 1, gamma (cytochrome b(NR)) subunit/hypothetical protein
Rejected: Position hash score below cutoff.

GCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGT  >  REL606/1286677‑1286745
                                   |                                 
gCGTATTTTAGTTCACATAGACACTGCTTCGGCggg                                   <  1:4165360/36‑1 (MQ=255)
gCGTATTTTAGTTCACATAGACACTGCTTCGGCggg                                   <  1:6120965/36‑1 (MQ=255)
gCGAATTTTAGTTCACATAGACACTGCTTCGGCggg                                   <  1:2576795/36‑1 (MQ=255)
gCGAATTTTAGTTCACATAGACACTGCTTCGGCggg                                   >  1:2703177/1‑36 (MQ=255)
gCGAATTTTAGTTCACATAGACACTGCTTCGGCggg                                   <  1:6295019/36‑1 (MQ=255)
 cGAATTTTAGTTCACATAGACACTGCTTCGGCgggg                                  <  1:7318502/36‑1 (MQ=255)
  gtaTTTTAGTTCACATAGACACTGCTTCGGCGGGGt                                 <  1:454/34‑1 (MQ=37)
  gAATTTTAGTTCACATAGACACTGCTTCGGCGGGGt                                 <  1:2646949/36‑1 (MQ=255)
  gAATTTTAGTTCACATAGACACTGCTTCGGCGGGGt                                 >  1:7313874/1‑36 (MQ=255)
  gAATTTTAGTTCACATAGACACTGCTTCGGCGGGGt                                 <  1:4536138/36‑1 (MQ=255)
   aaTTTTAGTTCGCATAGACACTGCTTCGGCGGGGtt                                <  1:4166011/36‑1 (MQ=25)
   aaTTTTAGTTCACATAGACACTGCTTCGGCGGGGtt                                <  1:1383104/36‑1 (MQ=255)
   aaTTTTAGTTCACATAGACACTGCTTCGCCGGGGtt                                >  1:6672761/1‑36 (MQ=25)
    aTTTTAGTTCACATAGACACTGCTTCGGCGGGGttt                               >  1:690117/1‑36 (MQ=255)
     ttttAGTTCACATAGACGCTGCTTCGGCGGGGtttt                              <  1:6382735/36‑1 (MQ=255)
     ttttAGTTCACATAGACACTGCTTCGGCGGGGtttt                              <  1:974897/36‑1 (MQ=255)
     ttttAGTTCACATAGACACTGCTTCGGCGGGGtttt                              <  1:187094/36‑1 (MQ=255)
      tttAGTTCACATAGACGCTGCTTCGGCGGGGttttt                             <  1:3639647/36‑1 (MQ=255)
       ttAGTTCACATAGACACTGCTTCGGCGGGGtttttt                            <  1:7060926/36‑1 (MQ=255)
         aGTTCACATAGACACTGCTTCGGCGGGGTTTTTTTa                          <  1:4948525/36‑1 (MQ=255)
          gTTCGCATAGACACTGCTTCGGCGGGGTTTTTTTAt                         <  1:6294349/36‑1 (MQ=21)
           ttCACATAAACACTGCTTCGGCTGGGTTTTTTTATg                        >  1:3145471/1‑36 (MQ=18)
            tCACATAGACACTGCTTCGGCGGGGTTTTTTTATgg                       <  1:647771/36‑1 (MQ=255)
             cacaTAGACACTGCTTCGGCGGGGTTTTTTTATggg                      <  1:6474932/36‑1 (MQ=255)
               caTAGACACTGCTTCGGCGGGGTTTTTTTATGGGCa                    <  1:5845514/36‑1 (MQ=255)
               caTAGACACTGCTTCGGCGGGGTTTTTTTATGGGCa                    >  1:3739890/1‑36 (MQ=255)
               caTAGACACTGCTTCGGCGGGGTTTTTTTATGGGCa                    <  1:6326730/36‑1 (MQ=255)
                  aGACACTGCTTCGGCGGGGTTTTTTTATTGGCTCgg                 >  1:5897368/1‑36 (MQ=18)
                  aGACACTGCTTCGGCGGGGTTTTTTTATGGGCTCgg                 >  1:4887121/1‑36 (MQ=21)
                  aGACACTGCTTCGGCGGGGTTTTTTTATGGGCGCgg                 >  1:7133379/1‑36 (MQ=21)
                  aGACACTGCTTCGGCGGGGTTTTTTTATGGGCACgg                 >  1:5220990/1‑36 (MQ=255)
                   gACACTGCTTCGGCGGGGTTTTTTTATGGGCACGGt                <  1:2190891/36‑1 (MQ=255)
                   gACACTGCTTCGGCGGGGTTTTTTTATGGGCACGGt                <  1:620856/36‑1 (MQ=255)
                   gACACTGCTTCGGCGGGGTTTTTTTATGGGCACGGt                <  1:5405476/36‑1 (MQ=255)
                   gACACTGCTTCGGCGGGGTTTTTTTATGGGCACGGt                <  1:4821635/36‑1 (MQ=255)
                        tGCTTCGGCGGGGTTTTTTTATGGGTGCGGTGCggg           >  1:5861636/1‑36 (MQ=21)
                        tGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCggg           <  1:5973600/36‑1 (MQ=255)
                        gtcTTCAGCTGGGCTTTTTTATGGGCACGGTGCggt           <  1:4923176/34‑2 (MQ=18)
                         gCTTCGGCGGGGTTTTTTTATGGGCACGGTGCgggg          <  1:5762832/36‑1 (MQ=255)
                          cTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGt         <  1:3733176/36‑1 (MQ=255)
                          cTTCGGCGGGGTATTTTTATGGGCACGGTCCgagga         <  1:4690389/36‑2 (MQ=21)
                            tCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGa       <  1:6006293/36‑1 (MQ=255)
                            tCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGa       <  1:5918264/36‑1 (MQ=255)
                            tCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGa       <  1:640196/36‑1 (MQ=255)
                             cggcggGGTTTTTTTATGGGCACGGTGCGGGGTggg      >  1:5887734/1‑34 (MQ=255)
                             cggcggGGTTTTTTTATGGGCACGGTGCGGGGTGAg      <  1:4864122/36‑1 (MQ=255)
                               ggggggGTTTGTTATGGGCACGGTGCGGGGTGAGtt    <  1:488720/34‑1 (MQ=25)
                                cggGGTTTTTTTATGGGCACGGTGCGGGGTGCGTTg   >  1:4697124/1‑36 (MQ=255)
                                 ggggTTTTTTTATGGGCACGGTGCGGGGTGAGTTGt  <  1:753203/36‑1 (MQ=255)
                                   |                                 
GCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGT  >  REL606/1286677‑1286745

Base quality scores:  ATCG  < 3 ≤  ATCG  < 8 ≤  ATCG  < 18 ≤  ATCG  < 29 ≤  ATCG  < 32 ≤  ATCG