New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 640118 =35 (0.690)6 (0.120) 3/62 NT 13.6% intergenic (+36/+19) cspE/ccrB cold shock protein E/camphor resistance protein CrcB
?REL606 = 2825478 41 (0.710)coding (741/1629 nt) ECB_02649 protein similar to L‑ribulokinase AraB

GATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTTGGTCTTTAGTGTGCGGTTGAGG  >  REL606/640083‑640153
                                   |                                   
cattagtcgttaatcctaactcttcagcagctgtgg                                     <  1:601547‑M1/1‑1 (MQ=255)
cattagtcgttaatcctaactcttcagcagctgtgg                                     <  1:3520744‑M1/1‑1 (MQ=255)
 attagtcgttaatcctaactcttcagcagctgtgGc                                    >  1:3484711‑M1/35‑36 (MQ=255)
 attagtcgttaatcctaactcttcagcagctgtgGc                                    >  1:3071309‑M1/35‑36 (MQ=255)
 aTACGTCAGCAAGAATTCAAAACCCGCTTAATCAGc                                    <  1:2942546/36‑1 (MQ=255)
  ttagtcgttaatcctaactcttcagcagctgtgGCg                                   <  1:4982480‑M1/3‑1 (MQ=255)
  ttagtcgttaatcctaactcttcagcagctgtgGCg                                   <  1:3763904‑M1/3‑1 (MQ=255)
  tACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCg                                   <  1:325118/36‑1 (MQ=255)
  tACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCg                                   <  1:7278325/36‑1 (MQ=255)
   tagtcgttaatcctaactcttcagcagctgtgGCgg                                  <  1:6662457‑M1/4‑1 (MQ=255)
   tagtcgttaatcctaactcttcagcagctgtgGCgg                                  <  1:6016900‑M1/4‑1 (MQ=255)
   tagtcgttaatcctaactcttcagcagctgtgGCgg                                  >  1:2624270‑M1/33‑36 (MQ=255)
     gTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGt                                <  1:2131089/36‑1 (MQ=255)
     gTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGt                                <  1:5499065/36‑1 (MQ=255)
      tCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGtt                               <  1:4793121/36‑1 (MQ=255)
      tCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGtt                               <  1:4157222/36‑1 (MQ=255)
       cAGCAAGAATTCAAAACCCGCTTAATCAGCGGGttt                              <  1:6206669/36‑1 (MQ=255)
       cAGCAAGAATTCAAAACCCGCTTAATCAGCGGGttt                              <  1:1904720/36‑1 (MQ=255)
        aGCAAGAATTCAAAACCCGCTTAATCAGCGGGtttt                             >  1:5905018/1‑36 (MQ=255)
        aGCAAGAATTCAAAACCCGCTTAATCAGCGGGtttt                             <  1:1609810/36‑1 (MQ=255)
         ttaatcctaactcttcagcagctgtgGCGGgttttg                            <  1:7043565‑M1/10‑2 (MQ=255)
         ttaatcctaactcttcagcagctgtgGCGGGttttt                            <  1:2164140‑M1/10‑1 (MQ=255)
         gCAAGAATTCAAAACCCGCTTAATCAGCGGGttttt                            <  1:2915122/36‑1 (MQ=255)
          cAAGAATTCAAAACCCGCTTAATCAGCGGGtttttt                           >  1:2411408/1‑36 (MQ=255)
           aaGAATTCAAAACCCGCTTAATCAGCGGGttttttt                          <  1:5641655/36‑1 (MQ=255)
            atcctaactcttcagcagctgtgGCGGGtttttttt                         <  1:4851749‑M1/13‑1 (MQ=255)
            aGAATTCAAAACCCGCTTAATCAGCGGGtttttttt                         <  1:5083807/36‑1 (MQ=255)
            aGAATTCAAAACCCGCTTAATCAGCGGGtttttttt                         >  1:4762158/1‑36 (MQ=255)
             tcctaactcttcagcagctgtgGCGGgttttttttc                        <  1:2611985‑M1/14‑2 (MQ=255)
             tcctaactcttcagcagctgtgGCGGGttttttttt                        <  1:116366‑M1/14‑1 (MQ=255)
             aaaTTCAAAACCCGCTTAATCAGCGGGttttttttt                        <  1:6679447/35‑1 (MQ=255)
              cctaactcttcagcagctgtgGCGGgtttttgtttc                       <  1:2756904‑M1/15‑2 (MQ=255)
              aaTTCAAAACCCGCTTAATCAGCGGGTTTTTTTTTg                       <  1:743930/36‑1 (MQ=255)
                  cAAAACCCGCTTAATCAGCGGGTTTTTTTTTGGTCt                   <  1:2426334/36‑1 (MQ=255)
                  cAAAACCCGCATAATCAGCGGGTTTTTTTTTGGTCt                   <  1:6445475/36‑1 (MQ=255)
                   aaaaCCCGCTTAATCAGCGGGTTTTTTTTTGGTCtt                  >  1:7139223/1‑36 (MQ=255)
                    aaccccGCTTAATCAGCGGGTTTTTTTTTGGTCttt                 <  1:1908274/33‑1 (MQ=255)
                    aaaCCCGCTTAATCAGCGGGTTTTTTTTTGGTCttt                 <  1:361764/36‑1 (MQ=255)
                    aaaCCCGCTTAATCAGCGGGTTTTTTTTTGGTCttt                 <  1:5659178/36‑1 (MQ=255)
                      aCCCGCTTAATCAGCGGGTTTTTTTTTGGTGTTTat               >  1:5503873/1‑35 (MQ=255)
                      aCCCGCTTAATCAGCGGGTTTTTTTTTGGTCTTTat               >  1:1539053/1‑35 (MQ=255)
                      aCCCGCTTAATCAGCGGGTTTTTTTTTGGTCTTTAg               >  1:6222422/1‑36 (MQ=255)
                       cccGCTTAATCAGCGGGTTTTTTTTTGGTCTTTAgt              <  1:1828890/36‑1 (MQ=255)
                       cccGCTTAATCAGCGGGTTTTTTTTTGGTCTTTAgt              >  1:4653727/1‑36 (MQ=255)
                        ccGCTTAATCAGCGGGTTTTTTTTTGGTCTTTAgtg             <  1:956291/36‑1 (MQ=255)
                         cGCTTAATCAGCGGGTTTTTTTTTGGTCTTTAgtgt            >  1:6656135/1‑36 (MQ=255)
                            ttAATCAGCGGGTTTTTTTTTGGTCTTTAGTGTGCg         <  1:3394098/36‑1 (MQ=255)
                             tAATCAGCGGGTTTTTTTTTTGTCTTTATTGTGCgg        >  1:1100072/1‑36 (MQ=37)
                             tAATCAGCGGGTTTTTTTTTGGTCTTTAGTGTGCgg        >  1:3996536/1‑36 (MQ=255)
                              aaTCAGCGGGTTTTTTTTTGGTCTTTAGTGTGCGGt       >  1:3067956/1‑36 (MQ=255)
                                   gCGGGTTTTTTTTTGGTCTTTATTGTGCGTTTGatg  >  1:3479514/1‑34 (MQ=37)
                                   GCGGGTTTTTTTTTGGTCTTTAGTGTTCGGTTGAgg  <  1:6163131‑M1/36‑1 (MQ=255)
                                   GCGGGTTTTTTTTTGGTCTTTAGTGTGCGGTTGAgg  >  1:1299932‑M1/1‑36 (MQ=255)
                                   |                                   
GATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTTGGTCTTTAGTGTGCGGTTGAGG  >  REL606/640083‑640153

Base quality scores:  ATCG  < 3 ≤  ATCG  < 8 ≤  ATCG  < 21 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG