New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 230935238 (0.730)5 (0.100) 3/60 NT 15.2% coding (795/1203 nt) yfaY competence damage‑inducible protein A
?REL606 3944634 = 18 (0.300)intergenic (+23/‑124) proM/aslB tRNA‑Pro/predicted regulator of arylsulfatase activity

GTGAATTGCTCGCTCAACGTCAGGCTGAACTGGCGGTTTTGCAATTCGCGACTGATCTGCGCGGGCAGTCC  >  REL606/2309317‑2309387
                                   |                                   
gTGAATTGCTCGCTCAACGTCAGGCTGAACTGGCGG                                     >  1:2014409‑M1/1‑36 (MQ=255)
gTGAATTGCTCGCTAAACGTCAGGCTGAACTGGCGG                                     >  1:6001966‑M1/1‑36 (MQ=255)
 tGAATTGCTCGCTCAACGTCAGGCTGAACTGGCGGt                                    <  1:2690875/36‑1 (MQ=255)
  gAATTGCTCGCTCAACGTCAGGCTGAACTGGCGGgt                                   >  1:5178745‑M1/1‑34 (MQ=255)
   aaTTGCTCGCTCAACGTCAGGCTGAACTGGCGGttt                                  <  1:3177114/36‑1 (MQ=255)
    aTTGCTCGCTCAACGTCAGGCTGAACTGGCGGtttt                                 <  1:5234583/36‑1 (MQ=255)
     ttGCTCGCTCAACGTCAGGCTGAACTGGCGGTTTTg                                <  1:1416449/36‑1 (MQ=255)
      tGCTCGCTCAACGTCAGGCTGAACTGGCGGTTTTGc                               <  1:2968698/36‑1 (MQ=255)
       gctcgctcAACGTCAGGCTGAACTGGCGGtgttttt                              >  1:423406‑M1/1‑29 (MQ=255)
       gctcgctcAACGTCAGGCTGAACTGGCGGTTTTGCa                              <  1:4169013/36‑1 (MQ=255)
       gctcgctcAACGTCAGGCTGAACTGGCGGTTTTGCa                              <  1:787281/36‑1 (MQ=255)
        ctcgctcAACGTCAGGCTGAACTGGCGGggtttttt                             >  1:2032517‑M1/1‑28 (MQ=255)
        ctcgctcAACGTCAGGCTGAACTGGCGGgctttttt                             >  1:6613612‑M1/1‑28 (MQ=255)
        ctcgctcAACGTCAGGCTGAACTGGCGGTTTTGCaa                             <  1:6104788/36‑1 (MQ=255)
         tcgctcAACGTCAGGCTGAACTGGCGGTTTTGCAAt                            <  1:5805399/36‑1 (MQ=255)
         tcgctcAACGTCAGGCTGAAATGGCGGTTTTGCAAt                            <  1:6388588/36‑1 (MQ=255)
          cgctcAACGTCAGGCTGAACTGGCGGTTTTGCAAtt                           <  1:5805752/36‑1 (MQ=255)
           gctcAACGTCAGGCTGAGCTGGCGGTTTTGCAATTc                          <  1:687179/36‑1 (MQ=255)
           gctcAACGTCAGGCTGAACTGTCTGTTTTGCAATTc                          <  1:5513690/36‑1 (MQ=37)
           gctcAACGTCAGGCTGAACTGGCGGTTTTGCAATTc                          <  1:2927492/36‑1 (MQ=255)
           gctcAACGTCAGGCTGAACTGGCGGTTTTGCAATTc                          <  1:403548/36‑1 (MQ=255)
             tcAACGTCAGGCTGAACTGGCGGggttttttttttt                        >  1:3053623‑M1/1‑23 (MQ=255)
             tcAACGTCAGGCTGAACTGGCGGTTTTGCAATTcgc                        <  1:24846/36‑1 (MQ=255)
              cAACGTCAGGCTGAACTGGCGGTTTTGCAATTcgcg                       <  1:6771200/36‑1 (MQ=255)
              cAACGTCAGGCTGAACTGGCGGTTTTGCAATTcgcg                       <  1:2128751/36‑1 (MQ=255)
               aaCGTCAGGCTGAACTGGCGGTTTTGCAATTCGCGa                      <  1:262437/36‑1 (MQ=255)
               aaCGTCAGGCTGAACTGGCGGTTTTGCAATTCGCGa                      <  1:2036552/36‑1 (MQ=255)
               aaCGTCAGGCTGAACTGGCGGTTTTGCAATTCGCGa                      <  1:4378600/36‑1 (MQ=255)
                 cGTCAGGCTGAACTGGCGGTTTTGCAATTCGCGACt                    <  1:3837872/36‑1 (MQ=255)
                 cGTCAGGCTGAACTGGCGGGTTTGTAATTCTTGACt                    >  1:6452657/1‑36 (MQ=12)
                   tCAGGCTGAACTGGCGGTTTTGCAATTCGCGACTGa                  <  1:197545/36‑1 (MQ=255)
                    cAGGCTGAACTGGCGGTTTTGCAATTCGCGACTGAt                 <  1:4574155/36‑1 (MQ=255)
                    cAGGCTGAACTGGCGGTTTTGCAATTCGCGACTGAt                 <  1:3820609/36‑1 (MQ=255)
                    cAGGCTGAACTGGCGGTTTTGCAATTCGCGACTGAt                 <  1:2546308/36‑1 (MQ=255)
                     aGGCTGAACTGGCGGTTTTGCAATTCGCGACTGATc                <  1:4360760/36‑1 (MQ=255)
                     aGGCTGAACTGGCGGTTTTGCAATTCGCGACTGATc                <  1:5093226/36‑1 (MQ=255)
                       gCTGAACTGGCGGTTTTGCAATTCGCGACTGATCTg              <  1:5781388/36‑1 (MQ=255)
                        cTGAACTGGCGGTTTTGCAATTCGCGACTGATCTgc             <  1:1607040/36‑1 (MQ=255)
                        cTGAACTGGCGGTTTTGCAATTCGCGACTGATCTgc             <  1:1313452/36‑1 (MQ=255)
                         tGAACTGGCGGTTTTGCAATTCGCGACTGATCTgcg            <  1:1956752/36‑1 (MQ=255)
                          gAACTGGCGGTTTTGCAATTCGCGACTGATCTgcgc           <  1:4265961/36‑1 (MQ=255)
                           aaCTGGCGGTTTTGCAATTCGCGACTGATCTTcgcg          >  1:4929345/1‑36 (MQ=255)
                            aCTGGCGGTTTTGCAATTCGCGACTGATCTGCGCgg         <  1:4311122/36‑1 (MQ=255)
                            aCTGGCGGTTTTGCAATTCGCGACTGATCTGCGCgg         <  1:3976522/36‑1 (MQ=255)
                              tGGCGGTTTTGCAATTCGCGACTGATCTGCGCGGGc       <  1:2014657/36‑1 (MQ=255)
                                gCGGTTTTGCAATTCGCGACTGATCTGCGCGGGCAg     >  1:6328332/1‑36 (MQ=255)
                                 cGGggttttttgttttctgtgcatttcgtcactttc    >  1:5712339‑M1/1‑3 (MQ=255)
                                 cGGggttttttgttttctgtgcatttcgtcactttc    >  1:4387668‑M1/1‑3 (MQ=255)
                                  ggggttttttgttttctgtgcatttcgtccctttcc   <  1:7500632‑M1/36‑35 (MQ=255)
                                  ggggttttttgttttctgtgcatttcgtcactttcc   >  1:6394775‑M1/1‑2 (MQ=255)
                                  ggTTTTGCAATTCGCGACTGATCTGCGCGGGCAGTc   <  1:2256733/36‑1 (MQ=255)
                                   gggttttttgttttctgtgcatttcgtcactttccc  >  1:595873‑M1/1‑1 (MQ=255)
                                   gggttttttgttttctgtgcatttcgtcactttccc  <  1:2957668‑M1/36‑36 (MQ=255)
                                   gggttttttgttttctgtgcatttcgtcactttccc  >  1:905893‑M1/1‑1 (MQ=255)
                                   gTTTTGCAATTCGCGACTGATCTGCGCGGGCAGTcc  >  1:543800/1‑36 (MQ=255)
                                   gTTTTGCAATTCGCGACTGATCTGCGCGGGCAGGcc  <  1:4134981/36‑1 (MQ=255)
                                   |                                   
GTGAATTGCTCGCTCAACGTCAGGCTGAACTGGCGGTTTTGCAATTCGCGACTGATCTGCGCGGGCAGTCC  >  REL606/2309317‑2309387

Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 22 ≤  ATCG  < 31 ≤  ATCG  < 37 ≤  ATCG