New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 1539621 =78 (0.760)8 (0.100) 4/38 NT 13.6% coding (51/2400 nt) dos cAMP phosphodiesterase, heme‑regulated
?REL606 = 3363882 40 (0.490)coding (771/1290 nt) rrmB 16S rRNA m5C967 methyltransferase, S‑adenosyl‑L‑methionine‑dependent

ACGGTGAACACATTTTGGATCTTTGTGCCGCCCCCGGCGGTAAAACAACGCATATCCTTGAGGTGGCAC  >  REL606/3363848‑3363916
                                  |                                  
aCGGTGAACACATTTTGGATCTTTGTGCCGccccc                                    >  1:9409499/1‑35 (MQ=255)
aCGGTGAACACATTTTGGATCTTTGTGCCGccccc                                    >  1:8870683/1‑35 (MQ=255)
aCGGTGAACACATTTTGGATCTTTGTGCCGccccc                                    <  1:18806760/35‑1 (MQ=255)
aCGGTGAACACATTTTGGATCTTTGTACCGccccc                                    >  1:12812523/1‑35 (MQ=255)
 cGGTGAACACATTTTGGATCTTTGTGCCGCCCCCg                                   <  1:17123629‑M2/35‑2 (MQ=255)
 cGGTGAACACATTTTGGATCTTTGTGCCGCCCCCg                                   >  1:8033970‑M2/1‑34 (MQ=255)
 cGGTGAACACATTTTGGATCTTTGTGCCGCCCCCg                                   >  1:10368893‑M2/1‑34 (MQ=255)
 cGGTGAACACATTTTGGATCTTTGTGCCGCCCCCg                                   <  1:18748705‑M2/35‑2 (MQ=255)
  ggTGAACACATTTTGGATCTTTGTGCCGCCCCCgg                                  >  1:4610040‑M2/1‑33 (MQ=255)
  ggTGAACACATTTTGGATCTTTGTGCCGCCCCCgg                                  <  1:19526442‑M2/35‑3 (MQ=255)
  ggTGAACACATTTTGGATCTTTGTGCCGCCCCCgg                                  <  1:13307573‑M2/35‑3 (MQ=255)
   gTGAACACATTTTGGATCTTTGTGCCGCCCCCggc                                 >  1:9953363‑M2/1‑32 (MQ=255)
   gTGAACACATTTTGGATCTTTGTGCCGCCCCCggc                                 <  1:10184431‑M2/35‑4 (MQ=255)
   gTGAACACATTTTGGATCTTTGTGCAGCCCCCggc                                 >  1:5518761‑M2/1‑32 (MQ=255)
    tGAACACATTTTGGATCTTTGTGCCGCCCCCggcg                                <  1:4148125‑M2/35‑5 (MQ=255)
    tGAACACATTTTGGATCTTTGTGCCGCCCCCggcg                                <  1:11383578‑M2/35‑5 (MQ=255)
     gAACACATTTTGGATCTTTGTGCCGCCCCCggcgg                               <  1:19729640‑M2/35‑6 (MQ=255)
     gAACACATTTTGGATCTTTGTGCCGCCCCCggcgg                               <  1:18136557‑M2/35‑6 (MQ=255)
      aaCACATTTTGGATCTTTGTGCCGCCCCCggcggg                              <  1:11898487‑M2/35‑7 (MQ=255)
      aaCACATTTTGGATCTTTGTGCCGCCCCCggcggg                              <  1:8532131‑M2/35‑7 (MQ=255)
      aaCACATTTTGGATCTTTGTGCCGCCCCCGGCGGt                              <  1:19641699/35‑1 (MQ=255)
      aaCACATTTTGGATCTTTGTGCCGCCCCCGGCGGt                              <  1:18500960/35‑1 (MQ=255)
      aaCACATTTTGGATCTTTGTGCCGCCCCCGGCGGt                              <  1:781882/35‑1 (MQ=255)
      aaCACATTTTGGATCTTTGTGCCGCCCCCGGCGGt                              <  1:9930961/35‑1 (MQ=255)
       acacATTTTGGATCTTTGTGCCGCCCCCGGCGGTa                             <  1:6561409/35‑1 (MQ=255)
        cacaTTTTGGATCTTTGTGCCGCCCCCGGCGGTaa                            <  1:11877991/35‑1 (MQ=255)
        cacaTTTTGGATCTTTGTGCCGCCCCCGGCGGTaa                            <  1:3793039/35‑1 (MQ=255)
         acaTTTTGGATCTTTGTGCCGCCCCCGGCGGTaaa                           <  1:2725948/35‑1 (MQ=255)
         acaTTTTGGATCTTTGTGCCGCCCCCGGCGGTaaa                           <  1:4668366/35‑1 (MQ=255)
          caTTTTGGATCTTTGTGCCGCCCACGGCGGTaaaa                          <  1:11836777/35‑1 (MQ=255)
           tttttGGATCTTTGTGCCGCCCCCGGCGGTAAaac                         <  1:18801850/34‑1 (MQ=255)
           aTTTTGGATCTTTGTGCCGCCCCCggcggggaaaa                         <  1:2635236‑M2/35‑12 (MQ=255)
            ttttGGATCTTTGTGCCGCCCCCggcggggaaaaa                        <  1:16223364‑M2/35‑13 (MQ=255)
            ttttGGATCTTTGTGCCGCCCCCggcggggaaaaa                        <  1:18958961‑M2/35‑13 (MQ=255)
            ttttGGATCTTTGTGCCGCCCCCGGCGGTAAaaca                        <  1:3561353/35‑1 (MQ=255)
            ttttGGATCTTTGTGCCGCCCCCGGCGGTAAaaca                        <  1:15255954/35‑1 (MQ=255)
             tttGGATCTTTGTGCCGCCTCCGGCGGGGAaacaa                       >  1:11347078/1‑35 (MQ=25)
             tttGGATCTTTGTGCCGCCCCCggcggggaaacaa                       <  1:15364491‑M2/35‑14 (MQ=255)
             tttGGATCTTTGTGCCGCCCCCggcggggaaaaaa                       <  1:20202255‑M2/35‑14 (MQ=255)
             tttGGATCTTTGTGCCGCCCCCggcggggaaaaaa                       <  1:527614‑M2/35‑14 (MQ=255)
              ttGGATCTTTGTGCCGCCCCCggcggggaaaaaaa                      <  1:8707507‑M2/35‑15 (MQ=255)
              ttGGATCTTTGTGCCGCCCCCGGCGGTAAaacaac                      <  1:5369592/35‑1 (MQ=255)
              ttGGATCTTTGTGCCGCCCCCGGCGGTAAaacaac                      <  1:9558200/35‑1 (MQ=255)
               tGGATCTTTGTGCCGCCCCCGGCGGTAAAACAACg                     <  1:11662412/35‑1 (MQ=255)
               tGGATCTTTGTGCCGCCCCCGGCGGTAAAACAACg                     <  1:17319016/35‑1 (MQ=255)
               tGGATCTTTGTGCCGCCCCCGGCGGTAAAACAACg                     <  1:11273451/35‑1 (MQ=255)
               tGGATCTTTGTGCCGCCCCCGGCGGTAAAACAACg                     <  1:15309569/35‑1 (MQ=255)
                ggATCTTTGTGCCGCCCCCggcggggaaaaaaagc                    <  1:13586159‑M2/35‑17 (MQ=255)
                ggATCTTTGTGCCGCCCCCggcggggaaaaaaaaa                    <  1:16714079‑M2/35‑17 (MQ=255)
                ggATCTTTGTGCCGCCCCCGGCGGTAAAACAACGc                    <  1:8088421/35‑1 (MQ=255)
                ggATCTTTGTGCCGCCCCCGGCGGTAAAACAACGc                    <  1:2136446/35‑1 (MQ=255)
                ggATCTTTGTGCCGCCCCCGGCGGGGAAACAACGc                    >  1:13635753/1‑35 (MQ=37)
                 gATCTTTGTGCCGCCCCCGGCGGTAAAACAACGCa                   <  1:3125965/35‑1 (MQ=255)
                 gATCTTTGTGCCGCCCCCGGCGGTAAAACAACGCa                   <  1:3574867/35‑1 (MQ=255)
                 gATCTTTGTGCCGCCCCCGGCGGGGAAAAAACGCa                   >  1:14470780/1‑35 (MQ=25)
                  aTCTTTGTGCCGCCCCCGGCGGTAAAACAACGCat                  <  1:9377722/35‑1 (MQ=255)
                  aTCTTTGTGCCGCCCCCGGCGGTAAAACAACGCat                  <  1:6390102/35‑1 (MQ=255)
                  aTCTTTGTGCCGCCCCCGGCGGGGAAAAAACACac                  >  1:19106490/1‑34 (MQ=16)
                   tCTTTGTGCCGCCCCCGGCGGTAAAACAACGCata                 <  1:10242873/35‑1 (MQ=255)
                   tCTTTGTGCCGCCCCCGGCGGTAAAACAACGCata                 <  1:19420872/35‑1 (MQ=255)
                   tCTTTGTGCCGCCCCCGGCGGTAAAACAACGCata                 <  1:17877747/35‑1 (MQ=255)
                   tCTTTGTGCCGCCCCCGGCGGGGAAACAAAACata                 >  1:18919106/1‑35 (MQ=16)
                    cTTTGTGCCGCCCCCGGCGGTAAAACAACGCatat                <  1:15115115/35‑1 (MQ=255)
                     tttGTGCCGCCCCCGGCGGTAAAACAACGCATATc               <  1:3482667/35‑1 (MQ=255)
                     tttGTGCCGCCCCCGGCGGTAAAACAACGCATATc               <  1:3746224/35‑1 (MQ=255)
                      ttGTGCCGCCCCCGGCGGTAAAACAACGCATATcc              <  1:18271403/35‑1 (MQ=255)
                      ttGTGCCGCCCCCGGCGGTAAAACAACGCATATcc              <  1:3247452/35‑1 (MQ=255)
                      ttGTGCCGCCCCCGGCGGTAAAACAACGCATATcc              <  1:12782378/35‑1 (MQ=255)
                       tgtgCCGCCCCCGGCGGTAAAACAACGCATATCCt             <  1:2034476/35‑1 (MQ=255)
                        gtgCCGCCCCCGGCGGTAAAACAACGCATATACtt            <  1:4278263/35‑1 (MQ=255)
                         tgCCGCCCCCGGCGGTAAAACAACGCATATCCTTg           <  1:5637387/35‑1 (MQ=255)
                         tgCCGCCCCCGGCGGTAAAACAACGCATATCCTTg           <  1:5668711/35‑1 (MQ=255)
                            cgccCCCGGCGGTAAAACAACGCATATCCTTGAgg        <  1:8162086/35‑1 (MQ=255)
                            cgccCCCGGCGGTAAAACAACGCATATCCTTGAgg        <  1:9983483/35‑1 (MQ=255)
                            cgccCCCGGCGGTAAAACAACGCATATCCTTGAgg        <  1:6667130/35‑1 (MQ=255)
                            cgccCCCGGCGGTAAAACAACGCATATCCTTGAgg        <  1:12270260/35‑1 (MQ=255)
                             gccCCCGGCGGTAAAACAACGCATATCCTTGAGGt       <  1:8045666/35‑1 (MQ=255)
                             gccCCCGGCGGTAAAACAACGCATATCCTTGAGGt       <  1:8746610/35‑1 (MQ=255)
                             gccCCCGGCGGTAAAACAACGCATATCCTTGAGGt       <  1:17993016/35‑1 (MQ=255)
                              cccccGGCGGTAAAACAACGCATATACTTGAGGTg      <  1:19201367/35‑1 (MQ=255)
                               ccccGGCGGTAAAACAACGCATATCCTTGAGGTgg     <  1:687390/35‑1 (MQ=255)
                               ccccGGCGGTAAAACAACGCATATCCTTGAGGTgg     <  1:29843/35‑1 (MQ=255)
                               ccccGGCGGTAAAACAACGCATATCCTTGAGGTgg     <  1:15903404/35‑1 (MQ=255)
                               ccccGGCGGTAAAACAACGCATATCCTTGAGGTgg     <  1:10179034/35‑1 (MQ=255)
                                cccGGCGGTAAAACAACGCATATCCTTGAGGTGGc    <  1:14806275/35‑1 (MQ=255)
                                cccGGCGGTAAAACAACGCATATCCTTGAGGTGGc    <  1:9660256/35‑1 (MQ=255)
                                 ccGGCGGTAAAACAACGCATATCCTTGAGGTGGCa   <  1:5696735/35‑1 (MQ=255)
                                 ccGGCGGTAAAACAACGCATATCCTTGAGGTGGCa   <  1:14063539/35‑1 (MQ=255)
                                  cggcggTAAAACAACGCATATCCTTGAGGTGGCAc  <  1:7402254/35‑1 (MQ=255)
                                  cggcggTAAAACAACGCATATCCTTGAGGTGGCAc  <  1:15130690/35‑1 (MQ=255)
                                  |                                  
ACGGTGAACACATTTTGGATCTTTGTGCCGCCCCCGGCGGTAAAACAACGCATATCCTTGAGGTGGCAC  >  REL606/3363848‑3363916

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ 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.