New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 3434228 =39 (0.860)6 (0.130) 4/52 NT 15.0% noncoding (2/768 nt) IS1 repeat region
?REL606 = 4568881 29 (0.480)coding (374/2151 nt) yjiY predicted inner membrane protein

CATAAAGTCCTGTACCGCACCGGCCAGCACTACCCCCGCCAGCAGCCACAGCGTGCCGGGCAGGTAGCC  >  REL606/4568847‑4568915
                                  |                                  
cATAAAGTCCTGTACCGCACCGGCCAGCACTACCCc                                   >  1:1766593‑M2/1‑35 (MQ=255)
cATAAAGTCCTGTACCGCACCGGCCAGCACTACCCc                                   <  1:3319849‑M2/36‑2 (MQ=255)
 aTAAAGTCCTGTACCGCACCGGCCAGCACTACCCcc                                  >  1:7112631‑M2/1‑34 (MQ=255)
 aTAAAGTCCTGTACCGCACCGGCCAGCACTACCCcc                                  <  1:5645477‑M2/36‑3 (MQ=255)
 aTAAAGTCCTGTACCGCACCGGCCAGCACTACCCcc                                  >  1:4731923‑M2/1‑34 (MQ=255)
    aaGTCCTGTACCGCACCGGCCAGCACTACCCccgcc                               >  1:1620601‑M2/1‑31 (MQ=255)
    aaGTCCTGTACCGCACCGGCCAGCACTACCCccgcc                               >  1:3677276‑M2/1‑31 (MQ=255)
     aGTCCTGTACCGCACCGGCCAGCACTACCCccgcca                              >  1:4229261‑M2/1‑30 (MQ=255)
     aGTCCTGTACCGCACCGGCCAGCACTACCCccccca                              <  1:2332728‑M2/36‑7 (MQ=255)
      gTCCTGTACCGCACCGGCCAGCACTACCCccgccac                             <  1:6171601‑M2/36‑8 (MQ=255)
         cTGTACCGCACCGGCCAGCACTACCCCCGCcagcag                          <  1:1337275/36‑1 (MQ=255)
         cTGTACCGCACCGGCCAGCACTACCCCCGCcagcag                          <  1:5778576/36‑1 (MQ=255)
          tGTACCGCACCGGCCAGCACTACCCCCGCCagcagc                         >  1:1651605/1‑36 (MQ=255)
          tGTACCGCACCGGCCAGCACTACCCCCGCCagcagc                         >  1:4772316/1‑36 (MQ=255)
          tGTACCCCACCGGCCAGCACTACCCCCGCCagcagc                         <  1:7904022/36‑1 (MQ=255)
           gTACCGCACCGGCCAGCACTACCCCCGCCAGCAGcc                        >  1:3228419/1‑36 (MQ=255)
            tACCGCACCGGCCAGCACTACCCCCGCAAGCAGCca                       >  1:1267603/1‑36 (MQ=255)
             aCCGCACCGGCCAGCACTACCCCCGCCAGCAGCcac                      >  1:1521638/1‑36 (MQ=255)
             aCCGCACCGGCCAGCACTACCCCCGCCAGCAGCcac                      >  1:4008308/1‑36 (MQ=255)
               cccACCGGCCACCACTACCCCCGCCAGCAGCCACAg                    <  1:1347815/34‑1 (MQ=38)
                gCACCGGCCAGCACTACCCCCGCCAGCAGCCACAGc                   <  1:4354597/36‑1 (MQ=255)
                gCACCGGCCAGCACTACCCCCGCCAGCAGCCACAGc                   >  1:3785274/1‑36 (MQ=255)
                cccccGGCCACCACTACCCCCGCCAGCAGCCACAGc                   <  1:3365418/33‑1 (MQ=37)
                ccACCGCCCACCACAACCCCCGCCAGCAGCCACAGc                   <  1:6662422/35‑1 (MQ=21)
                   ccGGCCAGCACTACCCCCGCCAGCAGCCACAGCGTg                >  1:5963970/1‑36 (MQ=255)
                        cAGCACTACCCCCGCCAGCAGCCACAGCGTGCCggg           >  1:7397279/1‑36 (MQ=255)
                         accACTACCCCCGCCAGCAGCCACAGCGTGCCGGGc          <  1:3455960/34‑1 (MQ=255)
                         accACTACCCCCGCCAGCAGCCACAGCGTGCCGGGc          <  1:2249637/34‑1 (MQ=255)
                         aGCACTACCCCCGCCAGCAGCCACAGCGTGCCGGGc          >  1:7017990/1‑36 (MQ=255)
                         aGCACTACCCCCGCCAGCAGCCACAGCGTGCCGGGc          >  1:6152007/1‑36 (MQ=255)
                         aGCACTACCCCCGCCAGCAGCCACAGCGTGCCGGGc          <  1:5033217/36‑1 (MQ=255)
                          gCACTACCCCCGCCAGCAGCCACAGCGCGCCGGGCa         >  1:3885128/1‑36 (MQ=255)
                          gCACAACCCccgccagcacatctgccagcacgtagc         >  1:2100801‑M2/1‑9 (MQ=255)
                           cACTACCCCCGCCAGCAGCCACAGCGTGCCGGGCAg        >  1:7625222/1‑36 (MQ=255)
                           cACCACCCccgccagcacatctgccagcacgtagct        >  1:1842650‑M2/1‑8 (MQ=255)
                             ccccccccGCCAGCAGCCACAGCGTGCCGGGCAGGt      <  1:3828512/33‑1 (MQ=255)
                              tACCCCCGCCAGCAGCCACAGCGTGCCGGGCAGGTa     >  1:4543574/1‑36 (MQ=255)
                              cccccccGCCAGCAGCCACAGCGTGCCGGGCAGGTa     <  1:3312888/34‑1 (MQ=255)
                              catccccgccagcacatctgccagcacgtagctgtc     >  1:3208921‑M2/4‑5 (MQ=255)
                              aacCCccgccagcacatctgccagcacgtagctgtc     >  1:5668887‑M2/2‑5 (MQ=255)
                               aCCCCCGCCAGCAGCCACAGCGTGCCGGGCAGGTAg    >  1:4909731/1‑36 (MQ=255)
                                cccccgccagcacatctgccagcacgtagctgtcct   >  1:4239415‑M2/1‑3 (MQ=255)
                                cccccGCCAGCAGCCACAGCGTGCCGGGCAGGTAGc   <  1:6120245/36‑1 (MQ=255)
                                 ccccgccagcacatctgccagcacgtagctgtcctg  >  1:3474554‑M2/1‑2 (MQ=255)
                                 ccccGCCAGCAGCCACAGCGTGCCGGGCAGGTAGcc  >  1:1649122/1‑36 (MQ=255)
                                  |                                  
CATAAAGTCCTGTACCGCACCGGCCAGCACTACCCCCGCCAGCAGCCACAGCGTGCCGGGCAGGTAGCC  >  REL606/4568847‑4568915

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 18 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG