New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 1095862 =53 (0.520)3 (0.050) 3/28 NT 10.4% coding (274/1509 nt) putP proline:sodium symporter
?REL606 1585457 = 21 (0.350)coding (495/666 nt) marC predicted transporter

TGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGAAAAAGATCAGCGGCGGCGCAACCATCAGTACCC  >  REL606/1585423‑1585491
                                  |                                  
tGCGTAGGCTTCCCCACAAAATAACCGCCACCAgg                                    >  1:16761658/1‑35 (MQ=255)
 gCGTAGGCTTCCCCACAAAATAACCGCCACCAGGa                                   >  1:10448215/1‑35 (MQ=255)
  cGTAGGCTTCCCCACAAAATAACCGCCACCAGGaa                                  >  1:8798399/1‑35 (MQ=255)
    tAGGCTTCCCCACAAAATAACCGCCACCAGGaaaa                                >  1:7188178/1‑35 (MQ=255)
    tAGGCTTCCCCACAAAATAACCGCCACCAGGaaaa                                >  1:5162291/1‑35 (MQ=255)
     atgCTTCCCCACAAAATAACCGCCACCAGGaaaaa                               >  1:14488756/3‑35 (MQ=255)
     aGGCTTCCCCACAAAATAACCGCCACCAGGaaaaa                               >  1:9369376/1‑35 (MQ=255)
      ggCTTCCCCACAAAATAACCGCCACCAGGAAAAAg                              >  1:16576621/1‑35 (MQ=255)
       gCTTCCCCACAAAATAACCGCCACCAGGAAAAAGa                             >  1:9339638/1‑35 (MQ=255)
        cTTCCCCACAAAATAACCGCCACCAGGAAAAAGAt                            >  1:17652683/1‑35 (MQ=255)
         ttCCCCACAAAATAACCGCCACCAGGAAAAAGATc                           >  1:16495114/1‑35 (MQ=255)
         ttCCCCACAAAATAACCGCCACCAGGAAAAAGATc                           >  1:287180/1‑35 (MQ=255)
         ttCCCCACAAAATAACCGCCACCAGGAAAAAGATc                           >  1:6634573/1‑35 (MQ=255)
          tCCCCACAAAATAACCGCCACCAGGAAAAAGATCa                          >  1:16538583/1‑35 (MQ=255)
             ccACAAAATAACCGCCACCAGGAAAAAGATCAgcg                       >  1:15964482/1‑35 (MQ=255)
             ccACAAAATAACCGCCACCAGGAAAAAGATCAgcg                       >  1:2626620/1‑35 (MQ=255)
              cacaAAATAACCGCCACCAGGAAAAAGATCAgcgg                      >  1:5074925/1‑35 (MQ=255)
               acaAAATAACCGCCACCAGGAAAAAGATCAgcggc                     >  1:11678017/1‑35 (MQ=255)
                caacccccccgccaccagGAAAAAGATCAgcggcg                    <  1:4683427‑M2/17‑1 (MQ=255)
                caAAATAACCGCCACCAGGAAAAAGATCAgcggcg                    >  1:19242872/1‑35 (MQ=255)
                caAAATAACCGCCACCAGGAAAAAGATCAgcggcg                    >  1:16479076/1‑35 (MQ=255)
                caAAATAACCGCCACCAGGAAAAAGATCAgcggcg                    >  1:4905742/1‑35 (MQ=255)
                caAAATAACCGCCACCAGGAAAAAGATCAgcggcg                    >  1:10581483/1‑35 (MQ=255)
                 aacccccccgccaccagGAAAAAGATCAgcggcgg                   <  1:7920503‑M2/18‑1 (MQ=255)
                 aaaaTAACCGCCACCAGGAAAAAGATCAgcggcgg                   >  1:20104199/1‑35 (MQ=255)
                 aaaaTAACCGCCACCAGGAAAAAGATCAgcggcgg                   >  1:20440634/1‑35 (MQ=255)
                  acccccccgccaccagGAAAAAGATCAgcggcggc                  <  1:8492430‑M2/19‑1 (MQ=255)
                  aaaTAACCGCCACCAGGAAAAAGATCAgcggcggg                  >  1:9665850/1‑34 (MQ=255)
                  aaaTAACCGCCACCAGGAAAAAGATCAgcggcggc                  >  1:2568527/1‑35 (MQ=255)
                   cccccccgccaccagGAAAAAGATCAgcggcggcg                 <  1:8565600‑M2/20‑1 (MQ=255)
                   accccccgccaccagGAAAAAGATCAgcggcggcg                 <  1:2207295‑M2/20‑1 (MQ=255)
                   aaTAACCGCCACCAGGAAAAAGATCAgcggcgggg                 >  1:7224101/1‑33 (MQ=255)
                   aaTAACCGCCACCAGGAAAAAGATCAgcggcgggg                 >  1:6423221/1‑33 (MQ=255)
                    ccccccgccaccagGAAAAAGATCAGCGGCGgcgc                <  1:303035‑M2/21‑1 (MQ=255)
                    ccccccgccaccagGAAAAAGATCAGCGGCGgcgc                <  1:4393222‑M2/21‑1 (MQ=255)
                    aTAACCGCCACCAGGAAAAAGATCAgcggcgggac                >  1:275115/1‑32 (MQ=255)
                      ccccgccaccagGAAAAAGATCAGCGGCGGCGCaa              <  1:14390364‑M2/23‑1 (MQ=255)
                      aaCCGCCACCAGGAAAAAGATCAGCGGCGGGGCaa              >  1:14544668/1‑35 (MQ=255)
                       cccgccaccagGAAAAAGATCAGCGGCGGCGCAAc             <  1:10089573‑M2/24‑1 (MQ=255)
                        ccGCCACCAGGAAAAAGATCAGCGGCGGGGCGAcc            >  1:16188164/1‑35 (MQ=37)
                         cgccaccagGAAAAAGATCAGCGGCGGCGCAGCCa           <  1:20301505‑M2/26‑1 (MQ=255)
                         cgccaccagGAAAAAGATCAGCGGCGGCGCAACCa           <  1:14405645‑M2/26‑1 (MQ=255)
                         cgccaccagGAAAAAGATCAGCGGCGGCGAAACCa           <  1:13096617‑M2/26‑1 (MQ=255)
                          gccaccagGAAAAGGATCAGCGGCGGCGCAACCAt          <  1:1235556‑M2/27‑1 (MQ=255)
                          gccaccagGAAAAAGATCAGCGGCGGCGCAACCAt          <  1:14666411‑M2/27‑1 (MQ=255)
                          gccaccagGAAAAAGATCAGCGGCGGCGCAACCAt          <  1:14051809‑M2/27‑1 (MQ=255)
                           ccaccaGGAAAAAGATCAGCGGCGGCGGAAACacc         >  1:8312310/1‑33 (MQ=25)
                              ccagGAAAAAGATCAGCGGCGGCGCAGCCATCAGt      <  1:16332014‑M2/31‑1 (MQ=255)
                              ccaGGAAAAAGATCAGCGGCGGGGGAACCAACAGt      >  1:10961236/1‑35 (MQ=25)
                               cagGAAAAAGATCAGCGGCGGCGCAACCATCAGTa     <  1:12950621‑M2/32‑1 (MQ=255)
                                 gGAAAAAGATCAGCGGCGGCGCAACCATCAGTAcc   <  1:14201917‑M2/34‑1 (MQ=255)
                                  gAAAAAGATCAGCGGCGGCGCAACCATCAGTAccc  <  1:12463806/35‑1 (MQ=255)
                                  |                                  
TGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGAAAAAGATCAGCGGCGGCGCAACCATCAGTACCC  >  REL606/1585423‑1585491

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.