New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 501044 =30 (0.560)3 (0.060) 3/62 NT 9.1% coding (238/261 nt) ybbD hypothetical protein
?REL606 4395045 = 30 (0.490)intergenic (‑120/‑62) yjfO/yjfP hypothetical protein/predicted hydrolase

TGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTTCCCTAACAAATGATGAT  >  REL606/501009‑501078
                                   |                                  
tGGATATATCAATATACACAAATGTAGAAGATGatt                                    >  1:4285608/1‑36 (MQ=255)
atgaagaaatccttagcacaagtttctgctggcggt                                    >  1:4164904‑M1/36‑36 (MQ=255)
atgaagaaatccttagcacaagtttctgctggcggt                                    <  1:7323356‑M1/1‑1 (MQ=255)
atgaagaaatccttagcacaagtctctgctggcggt                                    <  1:1838685‑M1/1‑1 (MQ=255)
 tgaagaaatccttagcacaagtttctgctggcggTa                                   >  1:264767‑M1/35‑36 (MQ=255)
 ggATATATCAATATAAACAAATGTAGAAGATGatta                                   <  1:3546522/36‑1 (MQ=255)
  gaagaaatccttagcacaagtttctgctggcggTAt                                  >  1:6536454‑M1/34‑36 (MQ=255)
  gaagaaatccttagcacaagtttctgctggcggTAt                                  <  1:6718331‑M1/3‑1 (MQ=255)
  gATATATCAATATACACAAATGTAGAAGATGattat                                  <  1:4733451/36‑1 (MQ=255)
   atatatCAATATACACAAATGTAGAAGATGattatt                                 <  1:2299111/36‑1 (MQ=255)
   aagaaatccttagcacaagtttctgctggcggTAtt                                 >  1:6780627‑M1/33‑36 (MQ=255)
     atatCAATATACACAAATGTAGAAGATGATTAtttt                               >  1:6453992/1‑36 (MQ=255)
     atatCAATATACACAAATGTAGAAGATGATTAtttt                               >  1:7629459/1‑36 (MQ=255)
      tatCAATATACACAAATGTAGAAGATGATTAttttt                              >  1:7139932/1‑36 (MQ=255)
      aaatccttagcacaagtttctgctggcggTAttttt                              <  1:366986‑M1/7‑1 (MQ=255)
        taaaTATACGCAAATGTAGAAGATGATTAttttttt                            <  1:1986200/34‑1 (MQ=38)
         cAATATACACAAATGTAGAAGATGATTATTTTTTTc                           >  1:7613139/1‑36 (MQ=255)
         cAATATACACAAATGTAGAAGATGATTATTTTTTTc                           >  1:120386/1‑36 (MQ=255)
           atatACACAAATGTAGAAGATGATTATTTTTTTCtt                         <  1:6855017/36‑1 (MQ=255)
            tataCACAAATGTAGAAGATGATTATTTTTTTCTTa                        <  1:1239957/36‑1 (MQ=255)
             ataCACAAATGTAGAAGATGATTATTTTTTTCTTAt                       <  1:739781/36‑1 (MQ=255)
             ataCACAAATGTAGAAGATGATTATTTTTTTCTTAt                       >  1:7870749/1‑36 (MQ=255)
             ataCACAAATGTAGAAGATGATTATTTTTTTCTTAt                       >  1:5075750/1‑36 (MQ=255)
              taCACAAATGTAGAAGATGATTATTTTTTTCTTAtt                      <  1:3249505/36‑1 (MQ=255)
              taCACAAATGTAGAAGATGATTATTTTTTTCTTAtt                      >  1:4040967/1‑36 (MQ=255)
               gcacaagtttctgctggcggTATTTTTTTTTTAttt                     <  1:5754430‑M1/16‑1 (MQ=255)
               acacAAATGTAGAAGATGATTATTTTTTTCTTAttt                     >  1:1984913/1‑36 (MQ=255)
                cacaAATGTAGAAGATGATTATTTTTTTCTTAtttt                    >  1:7916887/1‑36 (MQ=255)
                cacaAATGTAGAAGATGATTATTTTTTTCTTAtttt                    <  1:5887065/36‑1 (MQ=255)
                 acaagtttctgctggctgTATTTTTTTCTTAttttt                   <  1:901648‑M1/18‑1 (MQ=255)
                 acaAATGTAGAAGATGATTATTTTTTTCTTAttttt                   >  1:1899473/1‑36 (MQ=255)
                 acaAATGTAGAAGATGATTATTTTTTTCTTAttttt                   >  1:1582337/1‑36 (MQ=255)
                 acaAATGTAGAAGATGATTATTTTTTTCTTAttttt                   >  1:155367/1‑36 (MQ=255)
                   aaaTGTAGAAGATGATTATTTTTTTCTTATTTTTcc                 >  1:6889960/1‑36 (MQ=255)
                   aaaTGTAGAAGATGATTATTTTTTTCTTATTTTTcc                 >  1:1305723/1‑36 (MQ=255)
                         agaagaTGATTATTTTTTTCTTATTTTTCCCTAACa           <  1:5605805/36‑1 (MQ=255)
                            agaTGATTATTTTTTTCTTATTTTTCCCTAACAAat        >  1:3113685/1‑36 (MQ=255)
                             gatgatTATTTTTTTCTTATTTTTCCCTAACAAatg       <  1:4048285/36‑1 (MQ=255)
                              atgatTATTTTTTTCTTATTTTTCCCTAACAAatga      <  1:5288247/36‑1 (MQ=255)
                                 attattTTTTTCTTATTTTTCCCTAACAAacgatga   >  1:5896925/1‑36 (MQ=255)
                                  ttattTTTTTCTTATTTTTCCCTAACAAatgatgat  >  1:5529414/1‑36 (MQ=255)
                                   |                                  
TGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTTCCCTAACAAATGATGAT  >  REL606/501009‑501078

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 27 ≤  ATCG  < 34 ≤  ATCG  < 37 ≤  ATCG