New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 60665121 (0.400)6 (0.120) 3/60 NT 20.3% coding (568/957 nt) fepB iron‑enterobactin transporter subunit
?REL606 3650765 = 26 (0.430)intergenic (+46/+141) cspA/hokA major cold shock protein/toxic polypeptide, small

GCGGTATAGACAATGGCAGTGACCGGCTGCGGCGGTAATTTGATTTGCTCTTTCGCCGCTGCCAGTTG  >  REL606/606617‑606684
                                  |                                 
gcggTATAGACAATGGCAGTGACCGGCTGggggggg                                  >  1:1161056/1‑35 (MQ=37)
 cggTATAGACAATGGCAGTGACCGGCTGGGGCGGgg                                 >  1:3560276‑M1/1‑34 (MQ=255)
 cggTATAGACAATGGCAGTGACCGGCTGCGGGGGgg                                 >  1:4133588‑M1/1‑34 (MQ=255)
 cggTATAGACAATGGCAGTGACCGGCTGCGGCGGgg                                 >  1:1729233‑M1/1‑34 (MQ=255)
    tataGACAATGGCAGTGACCGGCTGCGGCGGTAatt                              <  1:7094683/36‑1 (MQ=255)
    tataGACAATGGCAGTGACCGGCTGCGGAGGggatt                              >  1:6669102‑M1/1‑31 (MQ=255)
     ataGACAATGGCAGTGACCGGCTGCGGCGGTAattt                             <  1:1279153/36‑1 (MQ=255)
      taGACAATGGCAGTGACCGGCTGCGGCGGTAatttg                            <  1:7164174/36‑1 (MQ=255)
      taGACAATGGCAGTGACCGGCTGCGGCGGTAatttg                            <  1:1580023/36‑1 (MQ=255)
       aGACAATGGCAGTGACCGGCTGCGGCGGTAatttga                           <  1:3717127/36‑1 (MQ=255)
       aGACAATGGCAGTGACCGGCAGCGGCGGTAatttga                           <  1:1959162/36‑1 (MQ=255)
        gACAATGGCAGTGACCGGCTGCGGCGGTAatttgat                          <  1:1916978/36‑1 (MQ=255)
        gACAATGGCAGTGACCGGCTGCGGCGGTAatttgat                          <  1:1440395/36‑1 (MQ=255)
          cAATGGCAGTGACCGGCTGCGGCGGgggtttttttt                        >  1:708645‑M1/1‑25 (MQ=255)
           aaTGGCAGTGACCGGCTGCGGCGGTAatttgatttg                       <  1:2909956/36‑1 (MQ=255)
            aTGGCAGTGACCGGCTGCGGCGGTAATTTGATTTGc                      <  1:3048196/36‑1 (MQ=255)
             tGGCAGTGACCGGCTGCGGCGGggattttttttttt                     >  1:5698149‑M1/1‑22 (MQ=255)
              ggCAGTGACCGGCTGCGGCGGTAATTTGATTTGctc                    <  1:3361697/36‑1 (MQ=255)
              ggCAGTGACCGGCTGCGGCGGTAATTTGATTTGctc                    <  1:7116138/36‑1 (MQ=255)
              ggCAGTGACCGGCTGCGGCGGTAATTTGATTTGctc                    <  1:4265604/36‑1 (MQ=255)
                cAGTGACCGGCTGCGGCGGTAATTTGATTTGCTCtt                  <  1:381443/36‑1 (MQ=255)
                cAGTGACCGGCTGCGGCGGTAATTTGATTTGCTCtt                  <  1:4917050/36‑1 (MQ=255)
                        ggCTGCGGCGGTAATTTGATTTGCTCTTTCGCCGCt          <  1:5735049/36‑1 (MQ=255)
                        ggCTGCGGCGGTAATTTGATTTGCTCTTTCGCCGCt          <  1:5811042/36‑1 (MQ=255)
                          cTGCGGCGGTAATTTTTTTTTTTCTTTTTCCGCttc        >  1:6362606/1‑34 (MQ=0)
                             cggcggTAATTTGATTTGCTCTTTCGCCGCTGCCAg     <  1:1336395/36‑1 (MQ=255)
                              ggcggggatttttttatttgttttcaggaaataaat    >  1:7577879‑M1/1‑5 (MQ=255)
                              ggcggggatttttttatttgttttcaggaaataaat    <  1:899528‑M1/36‑32 (MQ=255)
                               gcggggatttttttatttgttttcaggaaataaata   <  1:2502404‑M1/36‑33 (MQ=255)
                               gcggTAATTTGATTTGCTCTTTCGCCGCTGCCAGtt   <  1:296711/36‑1 (MQ=255)
                                cggTAATTTGATTTGCTCTTTCGCCGCTGCCAGTTg  >  1:1180387/1‑36 (MQ=255)
                                cggTAATTTGATTTGCTCTTTCGCCGCTGCCAGTTg  >  1:2846509/1‑36 (MQ=255)
                                  |                                 
GCGGTATAGACAATGGCAGTGACCGGCTGCGGCGGTAATTTGATTTGCTCTTTCGCCGCTGCCAGTTG  >  REL606/606617‑606684

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