New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 590471NA (NA)26 (0.610) 20/66 NT 68.4% noncoding (1977/1977 nt) IS150 repeat region
?REL606 = 3630716 12 (0.270)intergenic (‑171/‑305) ldrD/yhjV toxic polypeptide, small/predicted transporter

CAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGAACCTTTCCGGCTTGCCGCTGATTCTACGATTGCT  >  REL606/3630681‑3630751
                                   |                                   
cAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGG                                     <  1:3149755‑M2/36‑1 (MQ=255)
 aGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGt                                    <  1:792884‑M2/36‑2 (MQ=255)
  gggATAAAGCCCCCGCAACATTGCGCCTCACCGGtg                                   <  1:1693572‑M2/36‑3 (MQ=255)
   ggATAAAGCCCCCGCAACATTGCGCCTCACCGGtgt                                  >  1:4791161‑M2/1‑33 (MQ=255)
      tAAAGCCCCCGCAACATTGCGCCTCACCGGtgtact                               >  1:2784274‑M2/1‑30 (MQ=255)
       aaaGCCCCCGCAACATTGCGCCTCACCGGtgtactg                              >  1:3844889‑M2/1‑29 (MQ=255)
        aaGCCCCCGCAACATTGCGCCTCACCGGtgtactga                             <  1:3168903‑M2/36‑9 (MQ=255)
        aaGCCCCCGCAACATTGCGCCTCACCGGtgtactga                             <  1:2133692‑M2/36‑9 (MQ=255)
         aGCCCCCGCAACATTGCGCCTCACCGGtgtactgac                            >  1:3425274‑M2/1‑27 (MQ=255)
          gCCCCCGCACCATTGCGCCTCACCGGtgtactgacc                           >  1:386842‑M2/1‑26 (MQ=255)
           cccccGCAATATTGCGCCTCACCGGtgtactgaccc                          >  1:2903545‑M2/1‑25 (MQ=255)
           cccccGCAACATTGCGCCTCACCGGtgtactgaccc                          >  1:5340674‑M2/1‑25 (MQ=255)
            ccccGCAACATTGCGCCTCACCGGtgtactgacccc                         >  1:1627062‑M2/1‑24 (MQ=255)
             cccGCAACATTGCGCCTCACCGGtgtactgacccca                        >  1:806559‑M2/1‑23 (MQ=255)
              ccGCAACATTGCGCCTCACCGGtgtactgaccccaa                       <  1:5088475‑M2/36‑15 (MQ=255)
               cGCAACATTGCGCCTCACCGGtgtactgaccccaaa                      <  1:1531713‑M2/36‑16 (MQ=255)
                gCAACATTGCGCCTCACCGGtgtactgaccccaaaa                     <  1:440536‑M2/36‑17 (MQ=255)
                gCAACATTGCGCCTCACCGGtgtactgaccccaaaa                     <  1:332214‑M2/36‑17 (MQ=255)
                        gcgcCTCACCGGtgtactgaccccaaaaagttggac             >  1:400759‑M2/1‑12 (MQ=255)
                        gcgcCTCACCGGtgtactgaccccaaaaagttggac             <  1:420952‑M2/36‑25 (MQ=255)
                        gcgcCTCACCGGtgtactgaccccaaaaagttggac             >  1:2502611‑M2/1‑12 (MQ=255)
                           ccTCACCGGtgtactgaccccaaaaagttggacagt          >  1:3752164‑M2/1‑9 (MQ=255)
                           ccTCACCGGtgtactgaccccaaaaagttggacagt          >  1:744781‑M2/1‑9 (MQ=255)
                            cTCACCGGtgtactgaccccaaaaagttggacagtt         >  1:4424439‑M2/1‑8 (MQ=255)
                             tCACCGGtgtactgaccccaaaaagttggacagtta        <  1:5728154‑M2/36‑30 (MQ=255)
                               aCCGGtgtactgaccccaaaaaggtggacagttaaa      <  1:5919170‑M2/36‑32 (MQ=255)
                               aCCGGggtactgaccccaaaaagttggacagttaaa      <  1:843158‑M2/36‑32 (MQ=255)
                               aCCGGAACCTTTCCGGCTTGCCGCTGATTCTATGAt      >  1:3172029/1‑36 (MQ=12)
                               aCCGGAACCTTTCCGGCTTGCCGCTGATTCTACGAt      <  1:3876465/36‑1 (MQ=32)
                               aCCGGAACCTTTCCGGCTTGCCGCTGATTCTACGAt      >  1:2990890/1‑36 (MQ=32)
                                ccGGAGCCTTTCCGGCTTGCCGCTGATTCTACGAtt     <  1:1023445/36‑1 (MQ=255)
                                ccGGAACCTTTCCGGCTTGCCGCTGATTCTACGAtt     <  1:5522995/36‑1 (MQ=35)
                                  ggAACCTTTCCGGCTTGCCGATGATTCTACGATTGc   >  1:5435192/1‑36 (MQ=255)
                                   gAACCTTTCCGGTTTGCCGCTGATTCTACGATTgtt  >  1:1146507/1‑34 (MQ=11)
                                   gAACCTTTCCGGCTTGCCGCTGATTCTACGATTGCt  <  1:4768824/36‑1 (MQ=255)
                                   gAACCTTTCCGGCTTGCCGCTGATTCTACGATTGCt  <  1:3451596/36‑1 (MQ=255)
                                   gAACCTTTCCGGCTTGCCGCTGATTCTACGATTGCt  >  1:3154589/1‑36 (MQ=255)
                                   gAACCTTTCCGGCTTGCCGCTGATTCTACGATTGCt  >  1:2580102/1‑36 (MQ=255)
                                   gAACCTTTCCGGCTTGCCGCTGATTCTACGATTGCt  <  1:847161/36‑1 (MQ=255)
                                   |                                   
CAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGAACCTTTCCGGCTTGCCGCTGATTCTACGATTGCT  >  REL606/3630681‑3630751

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 20 ≤  ATCG  < 29 ≤  ATCG  < 37 ≤  ATCG