New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 15973030 (0.590)15 (0.300) 10/62 2.7 36.1% coding (5/585 nt) yadN predicted fimbrial‑like adhesin protein
?REL606 1357552 = 23 (0.410)coding (4/1386 nt) ycjJ putrescine importer

CATAAGACCAGACAGGGCAAAACCTAATTTTTTAGACATAAAAATCCTTGAATAAAATTCCATTAGTCA  >  REL606/159696‑159764
                                  |                                  
cATAAGACCAGACAGGGCAAAACCTAATTTTTTAG                                    >  1:3359731‑M1/1‑35 (MQ=255)
cATAAGACCAGACAGGGCAAAACCTAATTTTTTAG                                    <  1:15868541‑M1/35‑1 (MQ=255)
cATAAGACCAGACAGGGCAAAACCTAATTTTTTAG                                    <  1:8782756‑M1/35‑1 (MQ=255)
cATAAGACCAGACAGGGCAAAACCTAATTTTTTAG                                    >  1:9838940‑M1/1‑35 (MQ=255)
cATAAGACCAGACAGGGCAAAACCTAATTTTTTAG                                    <  1:9751879‑M1/35‑1 (MQ=255)
 aTAAGACCAGACAGGGCAAAACCTAATTTTTTAGa                                   <  1:2177613/35‑1 (MQ=255)
 aTAAGACCAGACAGGGCAAAACCTAATTTTTTAGa                                   >  1:4826030/1‑35 (MQ=255)
 aTAAGACCAGACAGGGCAAAACCTAATTTTTTAGa                                   >  1:6321572/1‑35 (MQ=255)
 aTAAGACCAGACAGGGCAAAACCTAATTTTTTAGa                                   <  1:7237469/35‑1 (MQ=255)
 aTAAGACCAGACAGGGCAAAACCTAATTTTTTAGa                                   <  1:15189692‑M1/35‑2 (MQ=255)
  tAAGACCAGACAGGGCAAAACCTAATTTTTTAGac                                  <  1:11841233‑M1/35‑3 (MQ=255)
  tAAGACCAGACAGGGCAAAACCTAATTTTTTAGAc                                  <  1:4666109/35‑1 (MQ=255)
   aaGACCAGACAGGGCAAAACCTAATTTTTTAGaca                                 <  1:15729689‑M1/35‑4 (MQ=255)
   aaGACCAGACAGGGCAAAACCTAATTTTTTAGaca                                 <  1:8046496‑M1/35‑4 (MQ=255)
   aaGACCAGACAGGGCAAAACCTAATTTTTTAGACa                                 <  1:1185431/35‑1 (MQ=255)
    aGACCAGACAGGGCAAAACCTAATTTTTTAGacat                                <  1:16166424‑M1/35‑5 (MQ=255)
     gACCAGACAGGGCAAAACCTAATTTTTTAGacata                               >  1:11871270‑M1/1‑30 (MQ=255)
     gACCAGACAGGGCAAAACCTAATTTTTTAGACATa                               >  1:7613808/1‑35 (MQ=255)
     gACCAGACAGGGCAAAACCTAATTTTTTAGACATa                               >  1:7785034/1‑35 (MQ=255)
     gACCAGACAGGGCAAAACCTAATTTTTTAGACATa                               <  1:3373087/35‑1 (MQ=255)
      aCCAGACAGGGCAAAACCTAATTTTTTAGacataa                              <  1:7903266‑M1/35‑7 (MQ=255)
       ccAGACAGGGCAAAACCTAATTTTTTAGacataaa                             <  1:9488975‑M1/35‑8 (MQ=255)
       ccAGACAGGGCAAAACCTAATTTTTTAGacataaa                             >  1:10173072‑M1/1‑28 (MQ=255)
       ccAGACAGGGCAAAACCTAATTTTTTAGacataaa                             >  1:14974846‑M1/1‑28 (MQ=255)
        cAGACAGGGCAAAACCTAATTTTTTAGACATaaaa                            >  1:6747960/1‑35 (MQ=255)
        cAGACAGGGCAAAACCTAATTTTTTAGACATaaaa                            >  1:6259416/1‑35 (MQ=255)
         aGACAGGGCAAAACCTAATTTTTTAGACATaaaaa                           >  1:6471186/1‑35 (MQ=255)
         aGACAGGGCAAAACCTAATTTTTTAGACATaaaaa                           <  1:3117665/35‑1 (MQ=255)
          ggcAGGGCAAAACCTAATTTTTTAGACATAAAAAt                          <  1:4628684/33‑1 (MQ=255)
          gACAGGGCAAAACCTAATTTTTTAGACATAAAAAt                          >  1:3287453/1‑35 (MQ=255)
            cAGGGCAAAACCTAATTTTTTAGACATAAAAATcc                        <  1:1358833/35‑1 (MQ=255)
               ggCAAAACCTNATTTTTTAGACATAAAAATCCTTg                     <  1:2455531/35‑1 (MQ=255)
               ggCAAAACCTAATTTTTTAGACATAAAAATCCTTg                     >  1:6393975/1‑35 (MQ=255)
               ggCAAAACCTAATTTTTTAGACATAAAAATCCTTg                     >  1:7415182/1‑35 (MQ=255)
               ggCAAAACCTAATTTTTTAGACATAAAAATCCTTg                     >  1:4905055/1‑35 (MQ=255)
                    aaCCTAATTTTTTAGACATAAAAATCCTTGAATaa                >  1:5140108/1‑35 (MQ=255)
                     aCCTAATTTTTTAGACATAAAAATCCTTGAATaaa               >  1:6783975/1‑35 (MQ=255)
                      ccTAATTTTTTAGACATAAAAATCCTTGAATaaaa              <  1:5933636/35‑1 (MQ=255)
                      ccTAATTTTTTAGACATAAAAATCCTTGAATaaaa              <  1:2852217/35‑1 (MQ=255)
                       cTAATTTTTTAGACATAAAAATCCTTGAATAAAAt             >  1:7460842/1‑35 (MQ=255)
                       cTAATTTTTTAGACATAAAAATCCTTGAATAAAAt             >  1:3668625/1‑35 (MQ=255)
                        tAATTTTTTAGccataattgctgcctcccctgcgt            <  1:2876031‑M1/35‑25 (MQ=255)
                        tAATTTATTAGccataattgctgcctcccctgcgt            <  1:8281363‑M1/35‑25 (MQ=255)
                         aaTTTTTTAGccataattgctgcctcccctgcgta           <  1:4455462‑M1/35‑26 (MQ=255)
                         aaTTTTTTAGACATAAAAATCCTTGAATAAAATTc           >  1:4880515/1‑35 (MQ=255)
                         aaTTTTTTAGACATAAAAATCCTTGAATAAAATTc           <  1:2977415/35‑1 (MQ=255)
                          aTTTTTTAGccataattgctgcctcccctgcgtaa          <  1:11827370‑M1/35‑27 (MQ=255)
                           ttttttAGACATAAAAATCCTTGAATAAAATTCCa         <  1:2519274/35‑1 (MQ=255)
                           tttattAGccataattgctgcctcccctgcgtaat         <  1:12443962‑M1/35‑28 (MQ=255)
                            tttttAGccataattgctgcctcccctgcgtaata        <  1:10508763‑M1/35‑29 (MQ=255)
                            tttttAGACATAAAAATCCTTGAATAAAATTCCAt        <  1:7779236/35‑1 (MQ=255)
                              tttAGACATAAAAATCCTTGAATAAAATTCCATTa      >  1:4196529/1‑35 (MQ=255)
                              tttAGACATAAAAATCCTTGAATAAAATTCCATTa      <  1:1716550/35‑1 (MQ=255)
                               ttAGccataatttctgcctcccctgcgtaataagg     <  1:16078017‑M1/35‑32 (MQ=255)
                               ttAGACATAAAAATCCTTGAATAAAATTCCATTAg     <  1:5810778/35‑1 (MQ=255)
                                tggACATAAAAATCCTTGAATAAAATTCCATTAGt    >  1:153801/3‑35 (MQ=255)
                                tAGccataattgctgcctcccctgcgtaataaggt    <  1:16123999‑M1/35‑33 (MQ=255)
                                tAGccataattgctgcctcccctgcgtaataaggt    >  1:14026266‑M1/1‑3 (MQ=255)
                                tAGACATAAAAATCCTTGAATAAAATTCCATTAGt    <  1:6241503/35‑1 (MQ=255)
                                 aGccataattgctgcctcccctgcgtaataaggtt   <  1:14038044‑M1/35‑34 (MQ=255)
                                 aGACATAAAAATCCTTGAATAAAATTCCATTAGTc   >  1:5090313/1‑35 (MQ=255)
                                  gccataattgctgcctcccctgcgtaataaggttt  <  1:16521953‑M1/35‑35 (MQ=255)
                                  gccataattgctgactcccctgcgtaataaggttt  >  1:8761801‑M1/1‑1 (MQ=255)
                                  |                                  
CATAAGACCAGACAGGGCAAAACCTAATTTTTTAGACATAAAAATCCTTGAATAAAATTCCATTAGTCA  >  REL606/159696‑159764

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 24 ≤  ATCG  < 31 ≤  ATCG  < 33 ≤  ATCG