New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 1109621 =29 (0.500)7 (0.120) 5/70 NT 17.9% coding (461/1359 nt) ycdT predicted diguanylate cyclase
?REL606 = 3061117 35 (0.610)coding (752/1173 nt) ECB_02854 conserved hypothetical protein

TGGCACCTTCAGTAAAAGTCTGGCCTCGATTGGCGGTTTTGCTGTCGGTTCGGAAGCGATGGAAGTGCTGC  >  REL606/3061082‑3061152
                                   |                                   
tGGCACCTTCAGTAAAAGTCTGGCCTCGATTGGCgg                                     <  1:465337/36‑1 (MQ=255)
tGGCACCTTCAGTAAAAGTCTGGCCTCGATTGGCgg                                     >  1:4428775/1‑36 (MQ=255)
tGGCACCTTCAGTAAAAGTCTGGCCTCGATTGGCgg                                     >  1:186761/1‑36 (MQ=255)
 ggCACCTTCAGTAAAAGTCTGGCCTCGATTGGCGGt                                    <  1:4516462/36‑1 (MQ=255)
 ggCACCTTCAGTAAAAGTCTGGCCTCGATTGGCGGg                                    <  1:3643962‑M2/36‑2 (MQ=255)
 ggCACCTTCAGTAAAAGTCTGGCCTCGATTGGAGGg                                    <  1:5045887‑M2/36‑2 (MQ=255)
  gatccTTCAGTAAAAGTCTGGCCTCGATTGGCGGtt                                   <  1:7404764/33‑1 (MQ=255)
  gCACCTTCAGTAAAAGTCTGGCCTCGATTGGCGGgt                                   <  1:1275845‑M2/36‑3 (MQ=255)
   cACCTTCAGTAAAAGTCTGGCCTCGATTGGCGGgtt                                  <  1:4906535‑M2/36‑4 (MQ=255)
    aCCTTCGGTAAAAGTCTGGCCTCGATTGGCGGtttt                                 <  1:6303362/36‑1 (MQ=255)
    aCCTTCAGTAAAAGTCTGGCCTGGATTGGCGGgttt                                 <  1:1683437‑M2/36‑5 (MQ=255)
    aCCTTCAGTAAAAGTCTGGCCTCGATTGGCGGtttt                                 <  1:1772474/36‑1 (MQ=255)
    aCCTTCAGTAAAAGTCTGGCCTCGATTGGCGGtttt                                 <  1:4770483/36‑1 (MQ=255)
     ccTTCAGTAAAAGTCTGGCCTCGATTGGCGGctttt                                <  1:4756890‑M2/36‑6 (MQ=255)
     ccTTCAGTAAAAGTCTGGCCTCGATTGGCGGTTTTg                                <  1:5640516/36‑1 (MQ=255)
     ccTTCAGTAAAAGTCTGGCCTCGATTGGCGGTTTTg                                <  1:49822/36‑1 (MQ=255)
       ttCAGTAAAAGTCTGGCCTCGATTGGCGGTTTTGCt                              <  1:6140264/36‑1 (MQ=255)
       ttCAGTAAAAGTCTGGCCTCGATTGGCGGTTTTGCt                              <  1:6219589/36‑1 (MQ=255)
        tCAGTAAAAGTCTGGCCTCGATTGGCGGTTTTGCTg                             <  1:402962/36‑1 (MQ=255)
         cAGTAAAAGTCTGGCCTCGATTGGCGGTTTTGCTGt                            <  1:203643/36‑1 (MQ=255)
         cAGTAAAAGTCTGGCCTCGATTGGCGGTTTTGCTGt                            <  1:4526652/36‑1 (MQ=255)
          aGTAAAAGTCTGGCCTCGATTGGCGGTTTTGCTGTc                           <  1:1627542/36‑1 (MQ=255)
            tAAAAGTCTGGCCTCGATTGGCGGTTTTGCTGTCgg                         <  1:750544/36‑1 (MQ=255)
             aaaaGTCTGGCCTCGATTGGCGGgtttgttttttgt                        <  1:5650929‑M2/36‑14 (MQ=255)
             aaaaGTCTGGCCTCGATTGGCGGTTTTGCTGTCGGt                        <  1:978038/36‑1 (MQ=255)
               aaGTCTGGCCTCGATTGGCGGTTTTGCTGTCGGTTc                      <  1:1215830/36‑1 (MQ=255)
               aaGTCTGGCCTCGATTGGCGGTTTTGCTGTCGGTTc                      <  1:5593140/36‑1 (MQ=255)
                aGTCTGGCCTCGATTGGCGGTTTTGCTGTCGGTTCg                     <  1:1235772/36‑1 (MQ=255)
                 gTCTGGCCTCGATTGGCGGTTTTGCTGTCGGTTCgg                    <  1:202025/36‑1 (MQ=255)
                 gTCTGGCCTCGATTGGCGGTTTTGCTGTCGGTTCgg                    <  1:2895790/36‑1 (MQ=255)
                  tcgggCCTCGATTGGCGGTTTTGCTGTCGGTTCGGa                   <  1:3098989/33‑1 (MQ=255)
                  tCTGGCCTCGATTGGCGGTTTTGCTGTCGGTTCGGa                   <  1:2242233/36‑1 (MQ=255)
                     ggCCTCGATTGGCGGTTTTGCTGTCGGTTCGGAAGc                <  1:4970653/36‑1 (MQ=255)
                      gCCTCGATTGGCGGTTTTGCTGTCGGTTCGGAAGCg               <  1:4580906/36‑1 (MQ=255)
                        cTCGATTGGCGGTTTTGCTGTCGGTTCGGAAGCGAt             <  1:2275409/36‑1 (MQ=255)
                         tCGATTGGCGGTTTTGCTGTCGGTTCGGAAGCGATg            <  1:3944921/36‑1 (MQ=255)
                          cGATTGGCGGTTTTGCTGTCGGTTCGGAAGCGATgg           <  1:7205960/36‑1 (MQ=255)
                          cGATTGGCGGTTTTGCTGTCGGTTCGGAAGCGATgg           <  1:1381974/36‑1 (MQ=255)
                           gATTGGCGGTTTTGCTGTCGGTTCGGAAGCGATGGa          <  1:4646237/36‑1 (MQ=255)
                            aTTGGCGGTTTTGCTGTCGGTTCGGAAGCGATGGaa         <  1:6879837/36‑1 (MQ=255)
                             ttGGCGGTTTTGCTGTCGGTTCGGAAGCGATGGAAg        <  1:2621473/36‑1 (MQ=255)
                              tGGCGGTTTTGCTGTCGGTTCGGAAGCGATGGAAGt       <  1:2136587/36‑1 (MQ=255)
                                 cGGTTTTGCTGTCGGTTCGGAAGCGATGGAAGtgct    >  1:5244558/1‑36 (MQ=255)
                                  ggTTTTGCTGTCGGTTCGGAAGCGATGGAAGtgctg   <  1:745176/36‑1 (MQ=255)
                                   gTTTTGCTGTCGGTTCGGAAGCGATGGAAGtgctgc  <  1:1852749/36‑1 (MQ=255)
                                   |                                   
TGGCACCTTCAGTAAAAGTCTGGCCTCGATTGGCGGTTTTGCTGTCGGTTCGGAAGCGATGGAAGTGCTGC  >  REL606/3061082‑3061152

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 22 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG