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

AAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCG  >  REL606/3650730‑3650800
                                   |                                   
aaaGCACAGAATCTAAGATCCCTGCCATTTGGCtgg                                     >  1:3918117/1‑36 (MQ=255)
aaaGCACAGAATCTAAGATCCCTGCCATTTGGCggg                                     >  1:5723753/1‑36 (MQ=255)
aaaGCACAGAATCTAAGATCCCTGCCATTTGGCggg                                     <  1:432929/36‑1 (MQ=255)
 cggtatagacaatggcagtgaccggctggggcgggg                                    >  1:3560276‑M2/35‑36 (MQ=255)
 cggtatagacaatggcagtgaccggctgcggggggg                                    >  1:4133588‑M2/35‑36 (MQ=255)
 cggtatagacaatggcagtgaccggctgcggcgggg                                    >  1:1729233‑M2/35‑36 (MQ=255)
 aaGCACAGAATCTAAGATCCCTGCCATTTGGCgggg                                    <  1:6273008/36‑1 (MQ=255)
   gCACAGAATCTAAGATCCCTGCCATTTGGCGGGGAt                                  >  1:4508351/1‑36 (MQ=255)
    tatagacaatggcagtgaccggctgcggaggGGAtt                                 >  1:6669102‑M2/32‑36 (MQ=255)
    cacaGAATCTAAGATCCCTGCCATTTGGCGGGGAtt                                 >  1:1796165/1‑36 (MQ=255)
    cacaGAATCTAAGATCCCTGCCATTTGGCGGGGAtt                                 >  1:5138109/1‑36 (MQ=255)
     acaGAATCTAAGATCCCTGCCATTTGGCGGGGAttt                                >  1:916012/1‑36 (MQ=255)
     acaGAATCTAAGATCCCTGCCATTTGGCGGGGAttt                                >  1:77029/1‑36 (MQ=255)
     acaGAATCTAAGATCCCTGCCATTTGGCGGGGAttt                                >  1:3420533/1‑36 (MQ=255)
       aGAATCTAAGATCCCTGCCATTTGGCGGGGAttttt                              >  1:3152072/1‑36 (MQ=255)
       aGAATCTAAGATCCCTGCCATTTGGCGGGGAttctt                              <  1:1362782/36‑1 (MQ=255)
        gAATCTAAGATCCCTGCCATTTGGCGGGGAtttttt                             <  1:3830508/36‑1 (MQ=255)
         aaTCTAAGATCCCTGCCATTTGGCGGGGAttttttt                            >  1:7640588/1‑36 (MQ=255)
         aaTCTAAGATCCCTGCCATTTGGCGGGGAttttttt                            <  1:3331549/36‑1 (MQ=255)
          caatggcagtgaccggctgcggcggGGgtttttttt                           >  1:708645‑M2/26‑35 (MQ=255)
           tCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTAt                          >  1:5115491/1‑36 (MQ=255)
            cTAAGATCCCTGCCATTTGGCGGGGATTTTTTTAtt                         >  1:4047372/1‑36 (MQ=255)
             tggcagtgaccggctgcggcggGGATTTTTTTTttt                        >  1:5698149‑M2/23‑36 (MQ=255)
               aGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGt                      >  1:5226994/1‑36 (MQ=255)
                gATCCCTGCCATTTGGCGGGGATTTTTTTATTTGtt                     <  1:5646209/36‑1 (MQ=255)
                 aTCCCTGCCATTTGGCGGGGATTTTTTTATTTGttt                    >  1:365900/1‑36 (MQ=255)
                   cccTGCCATTTGGCGGGGATTTTTTTATTTGTTTTc                  <  1:5744522/36‑1 (MQ=255)
                   cccTGCCATTTGGCGGGGATTTTTTTATTTGTTTTc                  <  1:5178154/36‑1 (MQ=255)
                    ccTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCa                 <  1:7904638/36‑1 (MQ=255)
                     cTGCCATTTAGCGGGGATTTTTTTATTTGTTTTCAg                >  1:6031499/1‑36 (MQ=255)
                       gCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGa              <  1:7525793/36‑1 (MQ=255)
                        ccATTTGGCGGGGATTTTTTTATTTGTTTTCAGGaa             >  1:7686094/1‑36 (MQ=255)
                        ccATTTGGCGGGGATTTTTTTATTTGTTTTCAGGaa             >  1:4725041/1‑36 (MQ=255)
                          aTTTGGCGGGGATTTTTTTATTTGTTTTCAGGaaat           >  1:2143767/1‑36 (MQ=255)
                           tttNGCGGGGATTTNNTTATTNNTTTTCAGGaaata          <  1:3379985/36‑1 (MQ=37)
                           tttGGCGGGGATTTTTTTATTTGTTTTCAGGaaata          <  1:7822659/36‑1 (MQ=255)
                           tttGGCGGGGATTTTTTTATTTGTTTTCAGGaaata          <  1:7214100/36‑1 (MQ=255)
                             tGGCGGGGATTTTTTTATTTGTTTTCAGGaaataaa        >  1:3622381/1‑36 (MQ=255)
                             tGGCGGGGATTTTTTTATTTGATTTCAGGaaataaa        >  1:7004537/1‑36 (MQ=255)
                              ggcggGGATTTTTTTATTTGTTTTCAGGaaataaat       >  1:7577879‑M2/6‑36 (MQ=255)
                              ggcggGGATTTTTTTATTTGTTTTCAGGaaataaat       <  1:899528‑M2/31‑1 (MQ=255)
                               gcggGGATTTTTTTATTTGTTTTCAGGAaataaata      <  1:2502404‑M2/32‑1 (MQ=255)
                                   ggATTTTTTTATTTGTTTTCAGGAAATAAATAatcg  <  1:1658247/36‑1 (MQ=255)
                                   ggATTTTTTTATTTGTTTTCAGGAAATAAATAatcg  <  1:4710454/36‑1 (MQ=255)
                                   |                                   
AAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCG  >  REL606/3650730‑3650800

Base quality scores:  ATCG  < 3 ≤  ATCG  < 14 ≤  ATCG  < 26 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG