New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 1529587 =69 (1.120)5 (0.080) 4/70 NT 9.5% coding (21/1698 nt) sfcA malate dehydrogenase, (decarboxylating, NAD‑requiring) (malic enzyme)
?REL606 = 3061121 26 (0.420)coding (756/1173 nt) ECB_02854 conserved hypothetical protein

ACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTAACTCACTCTTTT  >  REL606/1529552‑1529622
                                   |                                   
aCAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGt                                     <  1:4004519/36‑1 (MQ=255)
aCAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGt                                     <  1:2350525/36‑1 (MQ=255)
  aGGGCCAGCGTAAGGGATATAAAGCGAACGCTGttt                                   <  1:1161064/36‑1 (MQ=255)
  aGGGCCAGCATAAGGGATATAAAGCGAACGCTGttt                                   <  1:843368/36‑1 (MQ=255)
   tcagtaaaagtctggcctcgattggcgggttttttt                                  <  1:5236334‑M1/4‑1 (MQ=255)
    gggcAGCGTAAGTCATATGAAGCGAACGCTGttttt                                 <  1:188839/33‑1 (MQ=19)
    ggCCAGCGTAAGGGATATAAAGCGAACGCTGttttt                                 >  1:3953254/1‑36 (MQ=255)
    ggCCAGCGTAAGGGATATAAAGCGAACGCTGttttt                                 >  1:3944440/1‑36 (MQ=255)
    ggCCAGCGTAAGGGATATAAAGCGAACGCTGttttt                                 <  1:6386615/36‑1 (MQ=255)
    cagtaaaagtctggcctcgattggcggttttttttt                                 <  1:6287354‑M1/5‑1 (MQ=255)
     gCCAGCGTAAGGGATATAAAGCGAACGCTGtttttt                                >  1:4170424/1‑36 (MQ=255)
      gtaaaagtctggcctcgattggcggttttttttttt                               <  1:5472006‑M1/7‑1 (MQ=255)
      ccAGCGTAAGGGATATAAAGCGAACGCTGttttttt                               >  1:4720193/1‑36 (MQ=255)
      ccAGCGTAAGGGATATAAAGCGAACGCTGttttttt                               >  1:5053596/1‑36 (MQ=255)
      ccAGCGTAAGGGATATAAAGCGAACGCTGttttttt                               <  1:264869/36‑1 (MQ=255)
      ccAGCGTAAGGGATATAAAGCGAACGCTGttttttt                               <  1:6392379/36‑1 (MQ=255)
       cAGCGTAAGGGATATAAAGCGAACGCTGTTtttttg                              <  1:7318323/36‑1 (MQ=255)
        aGCGTAAGGGATATAAAGCGAACGCTGTTtttttgt                             >  1:3637092/1‑36 (MQ=255)
        aGCGTAAGGGATATAAAGCGAACGCTGTTtttttgt                             >  1:7433171/1‑36 (MQ=255)
         gCGTAAGGGATATAAAGCGAACGCTGTTtttttgtt                            <  1:2266504/36‑1 (MQ=255)
           gTAAGGGATATAAAGCGAACGCTGTTtttttgtttt                          >  1:2620586/1‑36 (MQ=255)
           gTAAGGGATATAAAGCGAACGCTGTTtttttgtttt                          <  1:4149534/36‑1 (MQ=255)
           gTAAGGGATATAAAGCGAACGCTGTTtttttgtttt                          >  1:3391749/1‑36 (MQ=255)
            tAAGGGATATAAAGCGAACGCTGTTtttttgttttt                         <  1:3050140/36‑1 (MQ=255)
            tAAGGGATATAAAGCGAACGCTGTTtttttgttttt                         >  1:6516100/1‑36 (MQ=255)
            tAAGGGATATAAAGCGAACGCTGTTttttgtttttt                         >  1:315062/1‑36 (MQ=37)
            gtctggcctcgattggcgggtttTTtttttgttttt                         <  1:3368538‑M1/13‑1 (MQ=255)
             aaGGGATATAAAGCGAACGCTGTTtttttgtttttt                        >  1:1148433/1‑35 (MQ=255)
             aaGGGATATAAAGCGAACGCTGTTtttttgtttttt                        >  1:501244/1‑35 (MQ=255)
             aaGGGATATAAAGCGAACGCTGTTtttttgtttttg                        <  1:5925394/36‑1 (MQ=255)
              ctggcctcgattggcggttttTttttctgtttttga                       <  1:232239‑M1/15‑2 (MQ=255)
              aGGGATATAAAGCGAACGCTGTTtttttgtttttgt                       >  1:177120/1‑35 (MQ=255)
               gggATATAAAGCGAACGCTGTTTTTTTGTTTTTGGt                      >  1:3959549/1‑36 (MQ=255)
                 gATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTc                    >  1:5542934/1‑36 (MQ=255)
                 gATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTc                    <  1:5804783/36‑1 (MQ=255)
                   tataAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCa                  <  1:3690261/36‑1 (MQ=255)
                   tataAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCa                  <  1:7039133/36‑1 (MQ=255)
                    ataAAGCGAACGCTGTTTTTTTGTTTTTTGTTCCAt                 >  1:5013166/1‑36 (MQ=255)
                    ataAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCAt                 >  1:414934/1‑36 (MQ=255)
                    ataAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCAt                 <  1:5958615/36‑1 (MQ=255)
                     taAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATg                <  1:3531525/36‑1 (MQ=255)
                     taAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATg                <  1:2414448/36‑1 (MQ=255)
                     taAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATg                <  1:8144538/36‑1 (MQ=255)
                      aaaGCGAACGCTGTTTTTTTGTTTTTTTTTCCATGt               >  1:3571646/1‑36 (MQ=37)
                      aaaGCGAACGCTGTTTTTTTGTTTTTTGTTCCAttt               >  1:7905081/1‑34 (MQ=38)
                      aaaGCGAACGCTGTTTTTTTGTTTTTGTTTCCATGt               >  1:2015385/1‑36 (MQ=255)
                      aaaGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGt               <  1:667071/36‑1 (MQ=255)
                      aaaGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGt               >  1:4665484/1‑36 (MQ=255)
                      aaaGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGt               >  1:4067448/1‑36 (MQ=255)
                      aaaGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGt               >  1:3122281/1‑36 (MQ=255)
                      aaaGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGt               >  1:8176994/1‑36 (MQ=255)
                       aaGCGAACGCTGTTTTTTTGTTTTTTGTTCCATGTa              >  1:1669074/1‑36 (MQ=255)
                       aaGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTa              >  1:1505023/1‑36 (MQ=255)
                       aaGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTa              <  1:6628199/36‑1 (MQ=255)
                        aGCGACCGCTGTTTTTTTGTTTTTGGTTCCATGTaa             <  1:1549933/36‑1 (MQ=255)
                        aGCGAACGCTGTTTTTTTGTTTTTTGTTCCATGTaa             >  1:5714305/1‑36 (MQ=255)
                         gCGAACGCTGTTTTTTTGTTTTTTGTTCCATGTAac            >  1:440555/1‑36 (MQ=255)
                         gCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTAac            >  1:5719250/1‑36 (MQ=255)
                         gCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTAac            <  1:688976/36‑1 (MQ=255)
                         gCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTAac            <  1:601009/36‑1 (MQ=255)
                         gCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTAac            <  1:2546527/36‑1 (MQ=255)
                          cGAACGCTGTTTTTTTGTTTTTTGTTCCATGTAact           >  1:1731530/1‑36 (MQ=255)
                          cGAACGCTGTTTTTTTGTTTTTGGTTCCATGTAact           >  1:1231328/1‑36 (MQ=255)
                          cGAACGCTGTTTTTTTGTTTTTGGTTCCATGTAact           <  1:3984831/36‑1 (MQ=255)
                          cGAACGCTGTTTTTTTGTTTTTGGTTCCATGTAact           <  1:4234759/36‑1 (MQ=255)
                          cGAACGCTGTTTTTTTGTTTTTGGTTCCATGTAact           <  1:7813209/36‑1 (MQ=255)
                           gAACGCTGTTTTTTTGTTTTTGGTTCCATGTAactc          >  1:2278926/1‑36 (MQ=255)
                           gAACGCTGTTTTTTTGTTTTTGGTTCCATGTAactc          >  1:2901180/1‑36 (MQ=255)
                            aaCGCTGTTTTTTTGTTTTTGGTTCCATGTAactca         <  1:1125280/36‑1 (MQ=255)
                              cGCTGTTTTTTTGTTTTTGGTTCCATGTAactcact       <  1:6998432/36‑1 (MQ=255)
                              cGCTGTTTTTTTGTTTTTGGTTCCATGTAactcact       <  1:4688521/36‑1 (MQ=255)
                               gCTGTTTTTTTGTTTTTGGTTCCATGTAactcactc      <  1:7430268/36‑1 (MQ=255)
                               gCTGTTTTTTTGTTTTTGGTTCCATGTAactcactc      <  1:5903916/36‑1 (MQ=255)
                               gCTGTTTTTTTGTTTTTGGTTCCATGTAactcactc      <  1:5291389/36‑1 (MQ=255)
                               gCTGTTTTTTTGTTTTTGGTTCCATGTAactcactc      <  1:4690614/36‑1 (MQ=255)
                                cTGTTTTTTTGTTTTTGGTTCCATGTAACTCActct     >  1:7303168/1‑36 (MQ=255)
                                   tttttttGTTTTTGGTTCCATGTAACTCACTCtttt  >  1:2915638/1‑36 (MQ=255)
                                   tttttttGTTTTTGGTTCCATGTAACTCACTCtttt  <  1:914103/36‑1 (MQ=255)
                                   |                                   
ACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTAACTCACTCTTTT  >  REL606/1529552‑1529622

Base quality scores:  ATCG  < 3 ≤  ATCG  < 13 ≤  ATCG  < 25 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG