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

GGTTGCGCAGCAATATTCAGTGGTGAATTAATAGCCATAATTGCTGCCTCCCCTGCGTAATAAGGTTTG  >  REL606/1357518‑1357586
                                  |                                  
ggTTGCGCAGCAATATTCAGTGGTGAATTAATAGc                                    >  1:6203305/1‑35 (MQ=255)
ggTTGCGCAGCAATATTCAGTGGTGAATTAATAGc                                    >  1:1937271/1‑35 (MQ=255)
  ttGCGCAGCAATATTCAGTGGTGAATTACTAGACa                                  >  1:339956/1‑35 (MQ=37)
  aagaccagacagggcaaaacctaattttttagACa                                  <  1:8046496‑M2/3‑1 (MQ=255)
  aagaccagacagggcaaaacctaattttttagACa                                  <  1:15729689‑M2/3‑1 (MQ=255)
   tGCGCAGCAATATTCAGTGGTGAATTNATAGCCAt                                 >  1:54709/1‑35 (MQ=255)
   tGCGCAGCAATATTCAGTGGTGAATTAATAGCCAt                                 >  1:2017387/1‑35 (MQ=255)
   agaccagacagggcaaaacctaattttttagACAt                                 <  1:16166424‑M2/4‑1 (MQ=255)
    gcgcAGCAATATTCAGTGGTGAATTAATAGCCATa                                >  1:3569796/1‑35 (MQ=255)
    gaccagacagggcaaaacctaattttttagACATa                                >  1:11871270‑M2/31‑35 (MQ=255)
     cgcAGCAATATTCAGTGGTGAATTAATAGCCATaa                               >  1:801290/1‑35 (MQ=255)
     cgcAGCAATATTCAGTGGTGAATTAATAGCCATaa                               >  1:1699927/1‑35 (MQ=255)
     accagacagggcaaaacctaattttttagACATaa                               <  1:7903266‑M2/6‑1 (MQ=255)
      ccagacagggcaaaacctaattttttagACAtaaa                              <  1:9488975‑M2/7‑2 (MQ=255)
      ccagacagggcaaaacctaattttttagACAtaaa                              >  1:10173072‑M2/29‑34 (MQ=255)
      ccagacagggcaaaacctaattttttagACAtaaa                              >  1:14974846‑M2/29‑34 (MQ=255)
       cagcaATATTCAGTGGTGAATTAATAGCCATAAtt                             >  1:2403312/1‑35 (MQ=255)
         gcaATATTCAGTGGTGAATTAATAGCCATAATtgc                           >  1:5797004/1‑35 (MQ=255)
         gcaATATTCAGTGGTGAATTAATAGCCATAATtgc                           >  1:5618109/1‑35 (MQ=255)
          caATATTCAGTGGTGAATTAATAGCCATAATtgct                          >  1:1957909/1‑35 (MQ=255)
          caATATTCAGTGGTGAATTAATAGCCATAATtgct                          >  1:5893560/1‑35 (MQ=255)
          caATATTCAGTGGTGAATTAATAGCCATAATtgct                          >  1:3123343/1‑35 (MQ=255)
            atatTCAGTGGTGAATTAATAGCCATAATtgctgc                        >  1:2161501/1‑35 (MQ=255)
              atTCAGTGGTGAATTAATAGCCATAATTGCTGCCt                      >  1:3346148/1‑35 (MQ=255)
               ttCAGTGGTGAATTAATAGCCATAATTGCTGCCTc                     >  1:498729/1‑35 (MQ=255)
                  tgtggtaATTTATTAGCCATAATTGCTGCCTcccc                  <  1:1322934/34‑1 (MQ=255)
                  aGTGGTGAATTAATAGCCATAATTGCTGCCTcccc                  >  1:6196207/1‑35 (MQ=255)
                  aGTGGTGAATTAATAGCCATAATTGCTGCCTcccc                  >  1:5975310/1‑35 (MQ=255)
                  aGTGGTGAATTAATAGCCATAATTGCTGCCTcccc                  >  1:5269697/1‑35 (MQ=255)
                    tggtgAATTAATAGCCATAATTGCTGCCTCCCCTg                >  1:5921369/1‑35 (MQ=255)
                      gtgAATTAATAGCCATAATTGCTGCCTCCCCTGCg              >  1:5285246/1‑35 (MQ=255)
                      gtgAATTAATAGCCATAATTGCTGCCTCCCCTGCg              >  1:2129018/1‑35 (MQ=255)
                       taattttttagCCATAATTGCTGCCTCCCCTGCGt             <  1:2876031‑M2/24‑1 (MQ=255)
                       taatttattagCCATAATTGCTGCCTCCCCTGCGt             <  1:8281363‑M2/24‑1 (MQ=255)
                        gAATTAATAGCCATAATTGCTGCCTCCCCTGCGta            >  1:7376236/1‑35 (MQ=255)
                        aattttttagCCATAATTGCTGCCTCCCCTGCGta            <  1:4455462‑M2/25‑1 (MQ=255)
                         attttttagCCATAATTGCTGCCTCCCCTGCGtaa           <  1:11827370‑M2/26‑1 (MQ=255)
                         aaTTATTAGCCATAATTGCTGCCTCCCCTGCGtaa           <  1:4867057/35‑1 (MQ=255)
                          tttattagCCATAATTGCTGCCTCCCCTGCGtaat          <  1:12443962‑M2/27‑1 (MQ=255)
                           tttttagCCATAATTGCTGCCTCCCCTGCGtaata         <  1:10508763‑M2/28‑1 (MQ=255)
                              ttagCCATAATTTCTGCCTCCCCTGCGTAATAAgg      <  1:16078017‑M2/31‑1 (MQ=255)
                               tagCCATAATTGCTGCCTCCCCTGCGTAATAAGGt     <  1:16123999‑M2/32‑1 (MQ=255)
                               tagCCATAATTGCTGCCTCCCCTGCGTAATAAGGt     >  1:14026266‑M2/4‑35 (MQ=255)
                                agCCATAATTGCTGCCTCCCCTGCGTAATAAGGtt    <  1:14038044‑M2/33‑1 (MQ=255)
                                 gCCATAATTGCTGCCTCCCCTGCGTAATAAGGttt   <  1:16521953‑M2/34‑1 (MQ=255)
                                 gCCATAATTGCTGACTCCCCTGCGTAATAAGGttt   >  1:8761801‑M2/2‑35 (MQ=255)
                                  ccATAATTGCTGCCTCCCCTGCGTAATAAGGTTTg  <  1:2208360/35‑1 (MQ=255)
                                  ccATAATTGCTGCCTCCCCTGCGTAATAAGGTTTg  <  1:1391254/35‑1 (MQ=255)
                                  ccATAATTGCTGCCTCCCCTGCGTAATAAGGTTTg  <  1:7214099/35‑1 (MQ=255)
                                  ccATAATTGCTGCCTCCCCTGCGTAATAAGGTTTg  <  1:5136124/35‑1 (MQ=255)
                                  |                                  
GGTTGCGCAGCAATATTCAGTGGTGAATTAATAGCCATAATTGCTGCCTCCCCTGCGTAATAAGGTTTG  >  REL606/1357518‑1357586

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