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

AAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAAT  >  REL606/3650730‑3650798
                                   |                                 
gcggtatagacaatggcagtgacccgctgcggcggg                                   >  1:4393111‑M2/36‑36 (MQ=255)
aaaGCACAGAATCTAATACCCCAGCCATTTGGCggg                                   <  1:3242866/36‑1 (MQ=25)
aaaGCACAGAATCGAAGATCCCTGCCATTTGGCggg                                   >  1:4368525/1‑36 (MQ=255)
 aaGCACAGAATCTAAGATCCCTGCCATTTGGCgggg                                  <  1:1099399/36‑1 (MQ=255)
 aaGCACAGAATCTAAGATCCCTGCCATTTGGCgggg                                  <  1:7050470/36‑1 (MQ=255)
  aGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGa                                 <  1:4905019/36‑1 (MQ=255)
  aGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGa                                 >  1:293529/1‑36 (MQ=255)
   gCACAGAATCTAAGATCCCTGCCATTTGGCGGGGAt                                <  1:6760452/36‑1 (MQ=255)
    tatagacaatggcagtgaccggctgcggcggGGGtt                               >  1:2817442‑M2/32‑36 (MQ=255)
     atagacaatggcagtgaccggctgcggcggGGAatt                              >  1:1317761‑M2/31‑36 (MQ=255)
     acaGAATCTAAGATCCCTGCCATTTGGCGGGGAttt                              <  1:2279028/36‑1 (MQ=255)
     acaGAATCTAAGATCCCTGCCATTTGGCGGGGAttt                              <  1:6913603/36‑1 (MQ=255)
      caGAATCTAAGATCCCTGCCATTTGGCGGGGAtttt                             <  1:6289239/36‑1 (MQ=255)
       agacaatggcagtgaccggctgcggcggGGAttttt                            >  1:1050899‑M2/29‑36 (MQ=255)
         acaatggcagtgaccggctgcggcggGGAttttttt                          >  1:6605894‑M2/27‑36 (MQ=255)
         acaatggcagtgaccggctgcggcggGGAttttttt                          >  1:4873211‑M2/27‑36 (MQ=255)
         aaTCTAAGATCCCTGCCATTTGGCGGGGAttttttt                          <  1:5891812/36‑1 (MQ=255)
           tCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTAt                        <  1:6444101/36‑1 (MQ=255)
              aaGATCCCTGCCATTTGGCGGGGATTTTTTTATTTg                     >  1:1738608/1‑36 (MQ=255)
                gATCCCTGCCATTTGGCGGGGATTTTTTTATTTTtt                   >  1:89467/1‑36 (MQ=255)
                gATCCCTGCCATTTGGCGGGGATTTTTTTATTTGtt                   <  1:1918430/36‑1 (MQ=255)
                gATCCCTGCCATTTGGCGGGGANTTTTTTATTTGtt                   <  1:4290395/36‑1 (MQ=255)
                 aTCCCTGCCATTTGGCGGGGATTTTTTTATTTGttt                  >  1:7147377/1‑36 (MQ=255)
                  tCCCTGCCATTTGGCGGGGATTTTTTTATTTGtttt                 >  1:3666155/1‑36 (MQ=255)
                    ccTGCCCTTTGGCGGGGATTTTTTTATTTGTTTTCa               <  1:4398231/36‑1 (MQ=255)
                    ccTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCa               <  1:5769403/36‑1 (MQ=255)
                       gCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGa            <  1:3922491/36‑1 (MQ=255)
                        ccATTTGGCGGGGATTTTTTTATTTGTTTTCAGGaa           >  1:7054236/1‑36 (MQ=255)
                         cATTTGGCGGGGATTTTTTTATTTGTTTTCATGaaa          >  1:3104219/1‑36 (MQ=255)
                          aTTTGGCGGGTATTTTTTTATTTGTTTTCAGGaaat         <  1:3405924/36‑1 (MQ=255)
                          aTTTGGCGGGGATTTTTTTATTTGTTTTCAGGaaat         >  1:1259022/1‑36 (MQ=255)
                           tttGGCGGGGATTTTTTTATTTTTTTTCAGGaaata        >  1:2738554/1‑36 (MQ=255)
                           tttGGCGGGGATTTTTTTATTTGTTTTCAGGaaata        <  1:6175728/36‑1 (MQ=255)
                            ttGGCGGGGATTTTTTTATTTGTTTTCAGGaaataa       >  1:452723/1‑36 (MQ=255)
                              ggcggGGATTTTTTTATTTGTTTTCAGGaaataaat     >  1:984275‑M2/6‑36 (MQ=255)
                                cggGGATTTTTTTATTTGTTTTCAGGAAataaataa   >  1:5236683‑M2/4‑36 (MQ=255)
                                 ggGGATTTTTTTATTTGTTTTCAGGAAATAaataat  <  1:6140404‑M2/34‑1 (MQ=255)
                                   |                                 
AAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAAT  >  REL606/3650730‑3650798

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