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

ACCTTCAGTAAAAGTCTGGCCTCGATTGGCGGTTTTGCTGTCGGTTCGGAAGCGATGGAAGTGCTGCGCTA  >  REL606/3061086‑3061156
                                   |                                   
aacTTCAGTAAAAGTCTGGCCTCGATTGGCGGgttt                                     >  1:424417/3‑36 (MQ=38)
aCCTTCAGTAAAAGTCTGGCCTCGATTGGCGGtttt                                     <  1:14647/36‑1 (MQ=255)
aCCTTCAGTAAAAGTCTGGCCTCGATTGGCGGtttt                                     >  1:1746533/1‑36 (MQ=255)
 ccTTCAGTAAAAGTCTGGCCTCGATTGGCGGTTTTg                                    <  1:2931735/36‑1 (MQ=255)
  cTTCAGTAAAAGTCTGTCCTCGATTGGCGGgttttt                                   >  1:2510156/1‑34 (MQ=37)
   ttCAGTAAAAGTCTGGCCTCGATTGGCGGTTTTGCt                                  <  1:6633315/36‑1 (MQ=255)
    tCAGTAAAAGTCTGGCCTCGATTGGCGGTTTTGCTg                                 <  1:233823/36‑1 (MQ=255)
    tCAGTAAAAGTCTGGCCTCGATTGGCGGGTTTtttt                                 <  1:5236334‑M2/36‑5 (MQ=255)
     cAGTCAAAGTCTGGCCTCGATTGGCGGTTTTGCTGt                                <  1:3460338/36‑1 (MQ=255)
     cAGTAAAAGTCTGGCCTCGATTGGCGGTTTTttttt                                <  1:6287354‑M2/36‑6 (MQ=255)
     cAGTAAAAGTCTGGCCTCGATTGGCGGTTTTGCTGt                                <  1:7814542/36‑1 (MQ=255)
     cAGTAAAAGTCTGGCCTCGATTGGCGGTTTTGCTGt                                <  1:5947181/36‑1 (MQ=255)
       gTAAAAGTCTGGCCTCGATTGGCGGTTTTttttttt                              <  1:5472006‑M2/36‑8 (MQ=255)
         aaaaGTCTGGCCTCGATTGGCGGTTTTGCTGTCGGt                            <  1:6977965/36‑1 (MQ=255)
         aaaaGTCTGGCCTCGATTGGCGGTTTTGCTGTCGGt                            <  1:6032221/36‑1 (MQ=255)
         aaaaGTCTGGCCTCGATTGGCGGTTTTGCTGTCGGt                            <  1:496597/36‑1 (MQ=255)
         aaaaGTCTGGCCTCGATTGGCGGATTTTCTTTTGGt                            >  1:5963314/1‑36 (MQ=19)
          aaaGTCTGGCCTCGATTGGCGGTTTTGCTGTCGGtt                           <  1:2541101/36‑1 (MQ=255)
          aaaGTCTGGCCTCGATTGGCGGTTTTGCTGTCGGtt                           <  1:1750742/36‑1 (MQ=255)
            aGTCTGGCCTCGATTGGCGGTTTTGCTGTCGGTTCg                         <  1:3330274/36‑1 (MQ=255)
             gTCTGGCCTCGATTGGCGGGTTTtttttttgttttt                        <  1:3368538‑M2/36‑14 (MQ=255)
              tCTGGCCTCGATTGGCGGTTTTGCTGTCGGTTCGGa                       <  1:1810996/36‑1 (MQ=255)
               cTGGCCTCGATTGGCGGTTTTtttttctgtttttga                      <  1:232239‑M2/36‑16 (MQ=255)
                tGGCCTCGATTGGCGGTTTTGCTGTCGGTTCGGAAg                     <  1:5852641/36‑1 (MQ=255)
                 ggCCTCGATTGGCGGTTTTGCTGTCGGTTCGGAAGc                    <  1:8009545/36‑1 (MQ=255)
                   ccTCGATTGGCGGTTTTGCTGTCGGTTCGGAAGCGa                  <  1:4726655/36‑1 (MQ=255)
                     tCGATTGGCGGTTTTGCTGTCGGTTCGGAAGCGATg                <  1:1893162/36‑1 (MQ=255)
                     tCGATTGGCGGTTTTGCTGTCGGTTCGGAAGCGATg                <  1:6954189/36‑1 (MQ=255)
                     tCGATTGGCGGTTTTGCTGTCGGTTCGGAAGCGATg                <  1:6104092/36‑1 (MQ=255)
                      cGATTGGCGGTTTTGCTGTCGGTTCGGAAGCGATgg               <  1:789090/36‑1 (MQ=255)
                          tGGCGGTTTTGCTGTCGGTTCGGAAGCGATGGAAGt           <  1:1413287/36‑1 (MQ=255)
                               gTTTTGCTGTCGGTTCGGAAGCGATGGAAGtgctgc      <  1:3333250/36‑1 (MQ=255)
                               gTTTTGCTGTCGGTTCGGAAGCGATGGAAGtgctgc      >  1:1949046/1‑36 (MQ=255)
                                 tttGCTGTCGGTTCGGAAGCGATGGAAGTGCTgcgc    <  1:3765836/36‑1 (MQ=255)
                                  ttGCTGTCGGTTCGGAAGCGATGGAAGTGCTGCGCt   <  1:371143/36‑1 (MQ=255)
                                   tGCTGTCGGTTCGGAAGCGATGGAAGTGCTGCGCtc  <  1:3695475/36‑2 (MQ=255)
                                   tGCTGTCGGTTCGGAAGCGATGGAAGTGCTGCGCTa  <  1:789900/36‑1 (MQ=255)
                                   |                                   
ACCTTCAGTAAAAGTCTGGCCTCGATTGGCGGTTTTGCTGTCGGTTCGGAAGCGATGGAAGTGCTGCGCTA  >  REL606/3061086‑3061156

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