New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 31070538 (0.930)3 (0.070) 3/50 NT 8.6% coding (1279/1548 nt) yahF predicted acyl‑CoA synthetase with NAD(P)‑binding domain and succinyl‑CoA synthetase domain
?REL606 580473 = 26 (0.450)coding (857/3144 nt) cusA copper/silver efflux system, membrane component

CATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCC  >  REL606/580438‑580507
                                   |                                  
cATTGCCGAACTAAACGGCGAAGGCGAAGTggcggg                                    >  1:5633218/1‑36 (MQ=255)
cATTGCCGAACTAAACGGCGAAGGCGAAGTggcggg                                    >  1:1820793/1‑36 (MQ=255)
 aTTGCCGAACTAAACGGCGAAGGCGAAGTggcgggc                                   <  1:7490088/36‑1 (MQ=255)
    gCCGAACTAAACGGCGAAGGCGAAGTGGcgggcggg                                <  1:4269379/36‑1 (MQ=255)
    atcgaagaagcgcgcgatccagaagtggcggggggg                                >  1:5017828‑M2/32‑36 (MQ=255)
     ccGAACTAAACGGCGAAGGCGAAGTGGCGGGGgggg                               >  1:1249855/1‑36 (MQ=255)
      cGAACTAAACGGCGAAGGCGAACTGGCGGGCGGggt                              <  1:4612137/36‑1 (MQ=255)
       gaagaagcgcgcgatccagaagtggcggGGGGggtg                             >  1:1638501‑M2/29‑36 (MQ=255)
       gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggtg                             <  1:3782362/36‑1 (MQ=255)
       gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggtg                             <  1:3780599/36‑1 (MQ=255)
       gAACTAAACGGCGAAGGCGAAGTAGCGGGGgggggg                             >  1:3079899/1‑34 (MQ=25)
        aagaagcgcgcgatccagaagtggcggGGGGggtgg                            >  1:5662320‑M2/28‑36 (MQ=255)
        aaCTAAACGGCGAAGGCGAAGTGGCGGGGGGggggg                            >  1:2856051/1‑36 (MQ=37)
         aCTAAACGGCGAANGCGAAGTGGCGGGGGGgggggg                           >  1:2649912/1‑35 (MQ=37)
         aCTAAACGGCGAAGGCGAAGTGGCGGGGGGgggggg                           >  1:4123179/1‑35 (MQ=37)
         aCTAAACGGCGAAGGCGAAGTGGCGGGCGGggtggt                           <  1:101325/36‑1 (MQ=255)
           tAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGa                         <  1:2973448/36‑1 (MQ=255)
           tAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGa                         <  1:475501/36‑1 (MQ=255)
           tAAACGGCGAAGGCGAAGGGGCGGGCGGggcggaca                         >  1:1370995/1‑33 (MQ=25)
             gtcGGGGAAGGTGGAGTTGCGGGATGGGTGGTGATc                       <  1:234803/34‑1 (MQ=9)
             aaCGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATc                       <  1:660252/36‑1 (MQ=255)
                 gcgaaggcgaagTGGCGGGCGGGGTGGTGATCCTGc                   <  1:5557760/36‑1 (MQ=255)
                  cgaaggcgaagTGGCGGGCGGGGTGGTGATCCTGCg                  <  1:5316880/36‑1 (MQ=255)
                   gaaggcgaagTGGCGGGCGGGGTGGTGATCCTGCGa                 <  1:7511655/36‑1 (MQ=255)
                    aaggcgaagTGGCGGGCGGGGGGGGGGTCCTTCggt                >  1:3890804/1‑34 (MQ=16)
                     aggcgaggTGGCGGGAGGGGTGGTTATCCTGCGATc               <  1:5830896/36‑1 (MQ=21)
                      ggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATcc              <  1:1800197/36‑1 (MQ=255)
                       gcgaagTGGCGGGGGGGGTGGTGATACTGGGAtttg             >  1:1754117/1‑33 (MQ=19)
                        cgaagTGGCGGGCGGGGTGGTGGTGCTGCGATCCgg            >  1:3699445/1‑36 (MQ=37)
                         gaagtggcggGCGGGGTGGTGATCCTGCGATCCGGc           <  1:3116353‑M2/26‑1 (MQ=255)
                         gaagTGGCGGGCGGGGGGGTGATCCTTCGGTCCGGc           >  1:5931456/1‑36 (MQ=25)
                            gtggcggGCGGGGTGGTGATCCTGCGATCCGGCaaa        >  1:6410201‑M2/8‑36 (MQ=255)
                            gtggcggGCGGGGGGGTGATCCTGCGATCCGGCaaa        >  1:5364130‑M2/8‑36 (MQ=255)
                             gggcggGCGGGGTGGTGATCCTGCTATCCGGCaaaa       <  1:2549052‑M2/30‑1 (MQ=255)
                              ggcggGCGGGGTGGTGATCCTGCGATCCGGCaaaaa      <  1:3915771‑M2/31‑1 (MQ=255)
                               gcggGCGGGGTGGTGATCCTGCGCTCCGGCAAAAAc     >  1:3225286‑M2/5‑36 (MQ=255)
                               gcggGCGGGGTGGTGATCCTGCGATCCGGCAAAAAc     <  1:5189818‑M2/32‑1 (MQ=255)
                                cggGCGGGGTGGTGATCCTGCGATCCGGCAAAAACg    >  1:1176380‑M2/4‑36 (MQ=255)
                                 ggGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGc   >  1:1430831‑M2/3‑36 (MQ=255)
                                 ggGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGc   >  1:4986751‑M2/3‑36 (MQ=255)
                                 ggGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGc   <  1:3413056‑M2/34‑1 (MQ=255)
                                  ggcgggGTGGCGACCCTGCNATCCGGCAAAAACGcc  >  1:2630903/1‑36 (MQ=25)
                                  gGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGcc  >  1:5600625‑M2/2‑36 (MQ=255)
                                   |                                  
CATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCC  >  REL606/580438‑580507

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 17 ≤  ATCG  < 27 ≤  ATCG  < 37 ≤  ATCG