New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 1980072 =32 (0.630)10 (0.200) 7/58 NT 22.0% intergenic (+48/‑295) yodA/yodB conserved metal‑binding protein/predicted cytochrome
?REL606 4038529 = 39 (0.630)intergenic (+41/‑176) bipA/yihL GTP‑binding protein/predicted DNA‑binding transcriptional regulator

TAACGCCTCGATGCCGCTTTAGTAAGTCATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGC  >  REL606/1980037‑1980106
                                   |                                  
tcaCGCCTCGATGCCGCTTTAGTAAGTCATCATAAc                                    <  1:5545678/34‑1 (MQ=255)
tAACGCCTCGATGCCGCTTTAGTAAGTCATCATAAc                                    <  1:1935549/36‑1 (MQ=255)
gcattgctgagtgctcacgacaggcaaaaaaaaacc                                    >  1:3557365‑M1/36‑36 (MQ=255)
 cattgctgagtgctcacgacaggcaaaaaaaaacCt                                   >  1:1749806‑M1/35‑36 (MQ=255)
 aaCGCCTCGATGCCGCTTTAGTAAGTCATCATAACt                                   >  1:6461170/1‑36 (MQ=255)
 aaCGCCTCGATGCCGCTTTAGTAAGTCATCATAACt                                   <  1:433290/36‑1 (MQ=255)
  attgctgagtgctcacgacaggcaaaaaaaaacCTg                                  >  1:967812‑M1/34‑36 (MQ=255)
   ttgctgagtgctcacgacaggcaaaaaaaaacCTGc                                 >  1:5811516‑M1/33‑36 (MQ=255)
   cGCCTCGATGCCGCTTTAGTAAGTCATCATAACTGc                                 >  1:721921/1‑36 (MQ=255)
    tgctgagtgctcacgacaggcaaaaaaaaacCTGcc                                >  1:8019640‑M1/32‑36 (MQ=255)
    tgctgagtgctcacgacaggcaaaaaaaaacCTGcc                                <  1:5665139‑M1/5‑1 (MQ=255)
    tgctgagtgctcacgacaggcaaaaaaaaacCTGcc                                <  1:4448957‑M1/5‑1 (MQ=255)
    gCCTCGATGCCGCTTTAGTAAGTCATCATAACTGcc                                >  1:7857793/1‑36 (MQ=255)
    gCCTCGATGCCGCTTTAGTAAGTCATCATAACTGca                                >  1:4966502/1‑35 (MQ=255)
     gctgagtgctcacgacaggcaaaaaaaaacCTGCCa                               <  1:2391998‑M1/6‑1 (MQ=255)
     gctgagtgctcacgacaggcaaaaaaaaacCTGCCa                               <  1:5888286‑M1/6‑1 (MQ=255)
     ccTCGATGCCGCTTTAGTAAGTCATCATAACTGCCa                               <  1:4459637/36‑1 (MQ=255)
     ccTCGATGCCGCTTTAGTAAGTCATCATAACTGCCa                               >  1:2154756/1‑36 (MQ=255)
     ccTCGATGCCGCTTTAGTAAGTAATCATAACTGcaa                               >  1:4210108/1‑34 (MQ=38)
         gATGCCGCTTTAGTAAGACATCATAACTGTCACTgg                           >  1:3028096/1‑36 (MQ=38)
           tGCCGCTTTAGTAAGTCATCATAACTGCCACTGGtc                         >  1:5168837/1‑36 (MQ=255)
               tttttAGAAAGTCATAAAACCTGCCACTGGtccatc                     <  1:6396179/34‑1 (MQ=16)
                 tttAGTAAGTCATCATAACTGCCACTGGtccatcca                   >  1:1711032/1‑36 (MQ=255)
                 tttAGTAAGTCATCATAACTGCCACTGGtccatcca                   >  1:1903835/1‑36 (MQ=255)
                 tttAGTAAGTCATCATAACTGCCACTGGtccatcca                   >  1:132675/1‑36 (MQ=255)
                 tttAATAAGTAATAATAACTGCCACTGGtccatcca                   <  1:4901420/36‑1 (MQ=25)
                  ttAGTAAGTCATCATAACTGCCACTGGTCCATCcac                  >  1:457587/1‑36 (MQ=255)
                  ttAGTAAGTCATCATAACTGCCACTGGTCCATCcac                  >  1:2909487/1‑36 (MQ=255)
                   tAGTAAGTCATCATAACTGCCACTGGTCCATCcaca                 >  1:1227423/1‑36 (MQ=255)
                   tAGTAAGTCATCATAACTGCCACTGGTCCATCcaca                 >  1:794788/1‑36 (MQ=255)
                      tAAGTCATCATAACTGCCACTGGTCCATCCACAAAc              >  1:4974207/1‑36 (MQ=255)
                      tAAGTCATCATAACTGCCACTGGTCCATCCACAAAc              >  1:1713375/1‑36 (MQ=255)
                       aaGTCATCATAACTGCCACTGGTCCATCCACAAACg             >  1:259585/1‑36 (MQ=255)
                        aGTCATCATAACTGCCACTGGGCCATACACAAACGc            >  1:5790658/1‑36 (MQ=37)
                         gTCATCATAACTGCCACTGGTCCATCCACAAACGcc           >  1:4850650/1‑36 (MQ=255)
                         aaaaataaacCTGCCACTGGTCCATCCACAAACGcc           >  1:3960250‑M1/11‑36 (MQ=255)
                          tcatcaTAACTGCCACTGGTCCATCCACAAACGCCa          >  1:1715204/1‑36 (MQ=255)
                          tcatcaTAACTGCCACTGGTCCATCCACAAACGCCa          >  1:311120/1‑36 (MQ=255)
                          tcatcaTAACTGCCACTGGTCCATCCACAAACGCCa          >  1:1041762/1‑36 (MQ=255)
                          tcatcaTAACTGCCACTGGTCCAGCCTCAAACGcta          >  1:4452680/1‑34 (MQ=37)
                          aaacaaaacCTGCCACTGGTCCATCCACAAACGCCa          >  1:5026594‑M1/10‑36 (MQ=255)
                          aaaaaaaaaCTGCCACTGGTCCATCCACAAACGCCa          >  1:2233725‑M1/10‑36 (MQ=255)
                           catcatAACTGCCACTGGTCCATCCACAAACGCCAc         >  1:3814145/1‑36 (MQ=255)
                           catcatAACTGCCACTGGTCCATCCACAAACGCCAc         >  1:7808064/1‑36 (MQ=255)
                           aataaaacCTGCCACTGGTCCATCCACAAACGCCAc         >  1:5171608‑M1/9‑36 (MQ=255)
                           aaaaaaacCTGCCACTGGTCCATCCACAAACGCCAc         >  1:4870508‑M1/9‑36 (MQ=255)
                            atcatAACTGCCACTGGTCCATCCACAAACGCCACt        >  1:4662693/1‑36 (MQ=255)
                            aaaaaaaCTGCCACTGGTCCATCCACAAACGCCACt        >  1:6224466‑M1/8‑36 (MQ=255)
                             tcatAACTGCCACTGGTCCATCCACAAACGCCACTg       >  1:7632170/1‑36 (MQ=255)
                             tcatAACTGCCACTGGTCCATCCACAAACGCCACTg       >  1:3741385/1‑36 (MQ=255)
                              catAACTGCCACTGGTCCATCCACAAACGCCACTGa      >  1:5079658/1‑36 (MQ=255)
                              aaaacCTGCCACTGGTCCATCCACAAACGCCACTGa      >  1:4260886‑M1/6‑36 (MQ=255)
                               aaaaCTGCCACTGGTCCATCCACAAACGCCACTGaa     <  1:7539781/34‑1 (MQ=255)
                                tAACTGCCACTGGTCCATCCACAAACGCCACTGAAc    >  1:2818490/1‑36 (MQ=255)
                                tAACTGCCACTGGTCCATCCACAAACGCCACTGAAc    >  1:1133150/1‑36 (MQ=255)
                                tAACTGCCACTGGTCCATCCACAAACGCCACTGAAc    >  1:2277438/1‑36 (MQ=255)
                                aacCTGCCACTGGTCCATCCACAAACGCCACTGAAc    >  1:5881481‑M1/4‑36 (MQ=255)
                                 acCTGCCACTGGTCCATCCACAAACGCCACTGAACg   >  1:3948815‑M1/3‑36 (MQ=255)
                                 aaCTGCCACTGGTCCATCCACAAACGCCACTGAACg   >  1:420889/1‑36 (MQ=255)
                                 aaCTGCCACTGGTCCATCCACAAACGCCACTGAACg   >  1:3974302/1‑36 (MQ=255)
                                  cCGGCCACTGGTCCATCCACAAACGCCACTGAACGc  >  1:3077283‑M1/2‑36 (MQ=255)
                                   |                                  
TAACGCCTCGATGCCGCTTTAGTAAGTCATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGC  >  REL606/1980037‑1980106

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