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

TTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCACTCAGCAATGCG  >  REL606/4038494‑4038564
                                   |                                   
ttttttAGTCATAAAAAACCTGCCAGCGATGGCAgg                                     >  1:557553/1‑36 (MQ=255)
gcgttcagtggcgtttgtggatggaccagtggccgg                                     >  1:3077283‑M2/36‑36 (MQ=255)
 cgttcagtggcgtttgtggatggaccagtggcagGt                                    >  1:3948815‑M2/35‑36 (MQ=255)
  ttttAGTCATAAAAAACCTGCCAGCGATGGCAGGtt                                   <  1:7480114/36‑1 (MQ=255)
  ttttAGTCATAAAAAACCTGCCAGCGATGGCAGGtt                                   <  1:2614567/36‑1 (MQ=255)
  gttcagtggcgtttgtggatggaccagtggcagGtt                                   >  1:5881481‑M2/34‑36 (MQ=255)
   tttAGTCATAAACAACCTGCCAGCGATGGCAGGttt                                  >  1:2584283/1‑36 (MQ=255)
    ttAGTCATAAAAAACCTGCCAGCGATGGCAGGtttt                                 <  1:4526088/36‑1 (MQ=255)
    ttAGTCATAAAAAACCTGCCAGCGATGGCAGGtttt                                 <  1:5776216/36‑1 (MQ=255)
    tcagtggcgtttgtggatggaccagtggcagGtttt                                 >  1:4260886‑M2/32‑36 (MQ=255)
     tAGTCATAAAAAACCTGCCAGCGATGGCAGGttttt                                >  1:7608426/1‑36 (MQ=255)
     tAGTCATAAAAAACCTGCCAGCGATGGCAGGttttt                                <  1:6032311/36‑1 (MQ=255)
     tAGTCATAAAAAACCTGCCAGCGATGGCAGGttttt                                <  1:8232043/36‑1 (MQ=255)
     tAGTCATAAAAAACCTGCCAGCGATGGCAGGttttt                                >  1:5306086/1‑36 (MQ=255)
      agtggcgtttgtggatggaccagtggcagTtttttt                               >  1:6224466‑M2/30‑36 (MQ=255)
       gtggcgtttgtggatggaccagtggcagGttttttt                              >  1:4870508‑M2/29‑36 (MQ=255)
       gtggcgtttgtggatggaccagtggcagGttttatt                              >  1:5171608‑M2/29‑36 (MQ=255)
       gTCATAAAAAACCTGCCAGCGATGGCAGGttttttt                              <  1:4936859/36‑1 (MQ=255)
        tggcgtttgtggatggaccagtggcagTtttttttt                             >  1:2233725‑M2/28‑36 (MQ=255)
        tggcgtttgtggatggaccagtggcagGttttgttt                             >  1:5026594‑M2/28‑36 (MQ=255)
        tCATAAAAAACCTGCCAGCGATGGCAGGtttttttt                             <  1:914976/36‑1 (MQ=255)
         ggcgtttgtggatggaccagtggcagGtttattttt                            >  1:3960250‑M2/27‑36 (MQ=255)
           tAAAAAACCTGCAAGCGATGGCAGGTTTTTTTTTGc                          >  1:2914041/1‑36 (MQ=255)
            aaaaaaCCTGCCAGCGATGGCAGGTTTTTTTTTGcc                         >  1:1499798/1‑36 (MQ=255)
             aaaaaCCTGCCAGCGATGGCAGGTTTTTTTTTGCCt                        >  1:2969467/1‑36 (MQ=255)
             aaaaaCCTGCCAGCGATGGCAGGTTTTTTTTTGCCt                        <  1:1330479/36‑1 (MQ=255)
             aaaaaCCTGCCAGCGATGGCAGGTTTTTTTTTGCCt                        <  1:18555/36‑1 (MQ=255)
              aaacCCTGGCATCGATGTGAGG.TTTTTTTTGCCTGt                      <  1:7063645/36‑1 (MQ=2)
              aaaaCCTGCCAGCGATGGCAGGTTTTTTTTTTCCTg                       >  1:8155156/1‑36 (MQ=255)
               acaCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGt                      <  1:3543022/34‑1 (MQ=255)
               aaaCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGt                      >  1:6592568/1‑36 (MQ=255)
                aaCCTGCCAGCGATGGCTGGTTTTTTTTTGCCTGTc                     >  1:1250786/1‑36 (MQ=255)
                aaCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTc                     >  1:6649004/1‑36 (MQ=255)
                aaCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTc                     >  1:1147561/1‑36 (MQ=255)
                aaCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTc                     <  1:2365172/36‑1 (MQ=255)
                aaCCTGCCAGCGATGGCAGGTTTTTTTTTGACTGTc                     >  1:1917100/1‑36 (MQ=255)
                 aCCTGCCAGCGATGGCAGGTTTTTTTTTTCCTTTCg                    >  1:3292731/1‑36 (MQ=21)
                  ccTGCCAGCGATGGCCGGTTTTTTTTTGCCTGTCGt                   <  1:4160592/36‑1 (MQ=255)
                  ccTGCCAGCGATGGATGGTTTTTTTTTGCCTGTCGt                   <  1:6735294/36‑1 (MQ=25)
                   cTGCCAGCGATGGCAGGTTTTTTTTTTCCTGTCGTg                  >  1:7715486/1‑36 (MQ=255)
                     gCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAg                >  1:5624489/1‑36 (MQ=255)
                      cgaGCGATGGCAGTTTTTTTTTTGCCTGTCGTGAGc               <  1:869127/34‑1 (MQ=37)
                      ccAGCGATGGCAGGTTTTTTTTTTCCTGTCGTGAGc               >  1:2830609/1‑36 (MQ=255)
                      ccAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGc               <  1:6930215/36‑1 (MQ=255)
                      ccAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGc               <  1:2293567/36‑1 (MQ=255)
                      ccAGCGATGGCAGGATTTTTTTTGCCTGTCGTGAGc               >  1:4275852/1‑36 (MQ=255)
                       cAGCGATGGCAGGTTTTTTTTTGCCTTTCGTGAGCa              >  1:3354406/1‑36 (MQ=255)
                       cAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCa              >  1:385365/1‑36 (MQ=255)
                       cAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCa              >  1:3864834/1‑36 (MQ=255)
                       cAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCa              <  1:4210437/36‑1 (MQ=255)
                       cAGCGATGGCAGGGTTTTTTTGGCATGTCGTGAGCa              >  1:3253390/1‑36 (MQ=25)
                        aGCGATGGCAGGTTTTTTTTTGACTGTCTTGAGCAc             >  1:1325689/1‑36 (MQ=37)
                          ctaTGGCAGGTTTTTTTTTGCCTGTCGTGAGCACTc           <  1:4566658/34‑1 (MQ=255)
                           gATGGCAGGTTTTTTTTTTCCTGTCGTGAGCACTCa          >  1:266989/1‑36 (MQ=255)
                           gATGGCAGGTTTTTTTTTGCCTGTCGTGAGCACTCa          >  1:6571378/1‑36 (MQ=255)
                           gATGGCAGGTTTTTTTTTGCCTGTCGTGAGCACTCa          <  1:1378292/36‑1 (MQ=255)
                           gATGGCAGGTTTTTTTTTGCCTGTCGTGAGCACTCa          >  1:8429546/1‑36 (MQ=255)
                           gATGGCAGGTTTTTTTTTGCCTGTCGTGAGAAATCa          >  1:2389525/1‑36 (MQ=37)
                            aTGGCAGGTTTTTTTTTGCCTGTCGTGAGCACTCAg         <  1:2761124/36‑1 (MQ=255)
                             tggcagGTTTTTTTTTGCCTGTCGTGAGCACTCAGc        <  1:2391998‑M2/30‑1 (MQ=255)
                             tggcagGTTTTTTTTTGCCTGTCGTGAGCACTCAGc        <  1:5888286‑M2/30‑1 (MQ=255)
                              ggcagGTTTTTTTTTGCCTGTCGTGAGCACTCAGCa       >  1:8019640‑M2/6‑36 (MQ=255)
                              ggcagGTTTTTTTTTGCCTGTCGTGAGCACTCAGCa       <  1:5665139‑M2/31‑1 (MQ=255)
                              ggcagGTTTTTTTTTGCCTGTCGTGAGCACTCAGCa       <  1:4448957‑M2/31‑1 (MQ=255)
                               gcagGTTTTTTTTTGCCTGTCGTGAGCACTCAGCaa      >  1:5811516‑M2/5‑36 (MQ=255)
                                cagGTTTTTTTTTGCCTGTCGTGAGCACTCAGCAAt     >  1:967812‑M2/4‑36 (MQ=255)
                                 agGTTTTTTTTTGCCTGTCGTGAGCACTCAGCAATg    >  1:1749806‑M2/3‑36 (MQ=255)
                                  gGTTTTTTTTTGCCTGTCGTGAGCACTCAGCAATGc   >  1:3557365‑M2/2‑36 (MQ=255)
                                   gTTTTTTTTTGCCTGTCGTGAGCACTCAGCAATGCg  <  1:7458555/36‑1 (MQ=255)
                                   gTTTTTTTTTGCCTGTCGTGAGCACTCAGCAATGCg  <  1:5076599/36‑1 (MQ=255)
                                   |                                   
TTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCACTCAGCAATGCG  >  REL606/4038494‑4038564

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