New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 12231147 (0.940)9 (0.180) 6/58 NT 17.6% coding (227/855 nt) ampE predicted inner membrane protein
?REL606 1529589 = 37 (0.620)coding (19/1698 nt) sfcA malate dehydrogenase, (decarboxylating, NAD‑requiring) (malic enzyme)

AGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTAACTCACTCTTTT  >  REL606/1529554‑1529622
                                   |                                 
aGGGCCAGCGTAAGGGATATAAAGCGAACGCTGttt                                   >  1:3394534/1‑36 (MQ=255)
 tgttcaacgttcccacgctactggtgtggctggttt                                  >  1:6981541‑M2/35‑36 (MQ=255)
 tgttcaacgttcccacgctactggtgtggctggttt                                  >  1:5495281‑M2/35‑36 (MQ=255)
 gggCCAGCGTAAGGGATATAAAGCGAACGCTGtttt                                  >  1:93459/1‑36 (MQ=255)
 gggCCAGCGTAAGGGATATAAAGCGAACGCTGtttt                                  <  1:1677098/36‑1 (MQ=255)
  ggCCAGCGTAAGGGATATAAAGCGAACGCTGttttt                                 <  1:733312/36‑1 (MQ=255)
  ggCCAGCGTAAGGGATATAAAGCGAACGCTGttttt                                 <  1:4365517/36‑1 (MQ=255)
   gCCAGCGTAAGGGATATAAAGCGAACGCTGtttttt                                >  1:6800885/1‑36 (MQ=255)
   gCCAGCGTAAGGGATATAAAGCGAACGCTGtttttt                                <  1:1517570/36‑1 (MQ=255)
   gCCAGCGTAAGGGATATAAAGCGAACGCTGtttttt                                >  1:227632/1‑36 (MQ=255)
    ccAGCGTAAGGGATATAAAGCGAACGCTGttttttt                               <  1:600811/36‑1 (MQ=255)
     caacgttcccacgctactggtgtggctggttttttg                              >  1:1970347‑M2/31‑36 (MQ=255)
     cAGCGTAAGGGATATAAAGCGAACGCTGTTtttttg                              >  1:1029766/1‑36 (MQ=255)
      tgCGTAAGGGATATAAAGCGAACGCTGTTtttttgt                             <  1:6447068/35‑1 (MQ=255)
      aacgttcccacgctactggtgtggctgtttttttgt                             >  1:3543888‑M2/30‑36 (MQ=255)
       gCGTAAGGGATATAAACCGAACGCTGTTtttttgtt                            >  1:4734908/1‑36 (MQ=255)
       acgttcccacgctactggtgtggctgtttttttttt                            >  1:2765457‑M2/29‑36 (MQ=255)
        cGTAAGGGATATAAAGCGAACGCTGTTtttttgttt                           >  1:2051824/1‑36 (MQ=255)
         gttcccacgctactggtgtggctgtttttttttttt                          >  1:1698595‑M2/27‑36 (MQ=255)
         gttcccacgctactggtgtggctgttttttggtttt                          >  1:4177602‑M2/27‑36 (MQ=255)
         gTAAGTGATATAAAGCGAACGCTGTTtttttgtttt                          <  1:2924268/36‑1 (MQ=255)
         gTAAGGGATATAAAGCGAACGCTGTTtttttgtttt                          <  1:1506048/36‑1 (MQ=255)
          ttcccacgctactggtgtggctgttttttttttttt                         >  1:186452‑M2/26‑36 (MQ=255)
          tAAGGGATATAAAGCGAACGCTGTTtttttgttttt                         >  1:6326198/1‑36 (MQ=255)
          tAAGGGATATAAAGCGAACGCTGTTtttttgttttt                         <  1:2677534/36‑1 (MQ=255)
          tAAGGGATATAAAGCGAACGCTGTTtttttgttttt                         <  1:5539329/36‑1 (MQ=255)
          tAAGGGATATAAAGCGAACGCTGTTtttttgttttt                         >  1:4491348/1‑36 (MQ=255)
           tcccacgctactggtgtggctgtttttttgtttttt                        >  1:2886918‑M2/25‑35 (MQ=255)
            aGGGATATAAAGCGAACGCTGTTtttttgtttttgt                       >  1:2003404/1‑35 (MQ=255)
            aGGGATATAAAGCGAACGCTGTTtttttgtttttgt                       >  1:3239245/1‑35 (MQ=255)
             gggATATCAAGCGAACGCTGTTTTTTTGTTTTTGGt                      <  1:1829899/36‑1 (MQ=255)
             gggATATAAAGCGAACGCTGTTTTTTTGTTTTTGGt                      >  1:3805935/1‑36 (MQ=255)
             gggATATAAAGCGAACGCTGTTTTTTTGTTTTTGGt                      >  1:1369602/1‑36 (MQ=255)
              ggATATAAAGCGAACGCTGTTTTTTTGTTTTTGGtt                     >  1:3721784/1‑36 (MQ=255)
               gATATAAAGCGATCGCTGTTTTTTTGTTTTTGGTTc                    <  1:6713513/36‑1 (MQ=255)
               gATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTc                    >  1:7035316/1‑36 (MQ=255)
                atatAAAGCGAACGCTGTTTTTTTGTTTTTGGTTcc                   <  1:6966586/36‑1 (MQ=255)
                atatAAAGCGAACGCTGTTTTTTTGTTTTTGGTTcc                   <  1:1738427/36‑1 (MQ=255)
                 tataAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCa                  <  1:1870690/36‑1 (MQ=255)
                  ataAAGCGAACGCTGTTTTTTTGTTTTTTGTTCCAt                 >  1:6938786/1‑36 (MQ=255)
                  ataAAGCGAACGCTGTTTTTTTGTTTTTGGTTCGAt                 >  1:6933483/1‑36 (MQ=255)
                   taAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATg                <  1:2968497/36‑1 (MQ=255)
                    aaaGCGAACTCTGTTTTTGGGTTTTTGGTTCCATGt               <  1:2221199/36‑1 (MQ=37)
                    aaaGCGAACGCTGTTTTTTTGTTTTTGTTTCCATGt               >  1:7189410/1‑36 (MQ=255)
                    aaaGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGt               <  1:7654532/36‑1 (MQ=255)
                      aGCGAACGCTGTTTTTTTGTTTTTTTTTCCATTTaa             >  1:2383220/1‑36 (MQ=25)
                      aGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTaa             >  1:7329827/1‑36 (MQ=255)
                      aGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTaa             >  1:5615478/1‑36 (MQ=255)
                        cGAACGCTGTTTTTTTGTTTTTGGTTCCATGTAact           >  1:872797/1‑36 (MQ=255)
                         gAACGCTGTTTTTTTGTTTTTGTTTCCATGTAactc          >  1:1370308/1‑36 (MQ=255)
                          aaCGCTGTTTTTTTGTTTTTGGTTCCATGTAactca         <  1:1793887/36‑1 (MQ=255)
                           aCGCTGTTTTTTTGTTTTTGGTTCCATGTAactcac        <  1:5392262/36‑1 (MQ=255)
                           aCGCTGTTTTTTTGTTTTTGGTTCCATGTAactcac        <  1:3729184/36‑1 (MQ=255)
                            cGCTGTTTTTTTGTTTTTGGTTCCATGTAactcact       >  1:7227896/1‑36 (MQ=255)
                            cGCTGTTTTTTTGTTTTTGGTTCCATGTAactcact       >  1:1649937/1‑36 (MQ=255)
                            cGCTGTTTTTTTGTTTTTGGTTCCATGTAactcact       <  1:3756256/36‑1 (MQ=255)
                             gctgttTTTTTGTTTTTTGTTCCATGTAactcactc      >  1:1600772‑M2/7‑36 (MQ=255)
                                gttTTTTTGTTTTTGGTTCCATGTAACTCACTCttt   >  1:6541722‑M2/4‑36 (MQ=255)
                                gttTTTTTGTTTTTGGTTCCATGTAACTCACTCttt   >  1:7344240‑M2/4‑36 (MQ=255)
                                gttTTTTTGTTTTTGGTTCCATGTAACTAACTCttt   >  1:1137152‑M2/4‑36 (MQ=255)
                                 ttTTTTTGTTTTTGGTTCCATGTAACTCACTCtttt  <  1:305289‑M2/34‑1 (MQ=255)
                                 tctttttgtttttgGTTGCATGTAACTCACTCtttt  <  1:4794738/34‑1 (MQ=38)
                                   |                                 
AGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTAACTCACTCTTTT  >  REL606/1529554‑1529622

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