New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 1462251 =34 (0.730)6 (0.130) 6/54 NT 11.4% intergenic (‑13/‑328) mokB/trg regulatory peptide/methyl‑accepting chemotaxis protein III, ribose and galactose sensor receptor
?REL606 = 3413922 59 (0.970)intergenic (‑199/‑103) yhfA/crp hypothetical protein/DNA‑binding transcriptional dual regulator

ATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATAT  >  REL606/3413887‑3413957
                                   |                                   
aTGTATGCAAAGGACGTCACATTACCGTGCAGTACa                                     <  1:5125404/36‑1 (MQ=255)
aTGTATGCAAAGGACGTCACATTACCGTGCAGTACa                                     <  1:5033054/36‑1 (MQ=255)
aTGTATGCAAAGGACGTCACATTACCGTGCAGTACa                                     <  1:3149652/36‑1 (MQ=255)
aTGTATGCAAAGGACGTCACATTACCGTGCAGTACa                                     >  1:3492836/1‑36 (MQ=255)
 tGTATGCAAAGGACGTCACATTACCGTGCAGTACAg                                    <  1:2990216‑M2/36‑2 (MQ=255)
 tGTATGCAAAGGACGTCACATTACCGTGCAGTACAg                                    <  1:1261078‑M2/36‑2 (MQ=255)
 tGTATGCAAAGGACGTCACATTACCGTGCAGTACAg                                    <  1:406893‑M2/36‑2 (MQ=255)
  gTATGCAAAGAACGTCACATTACCGTGCAGTACAgt                                   >  1:6787242‑M2/1‑34 (MQ=255)
  gTATGAAAAGGACGTCACATTACCGTGCAGTACAgt                                   >  1:5145902‑M2/1‑34 (MQ=255)
   tATGCAAAGGACGTCACATTACCGTGCAGTACAgtt                                  <  1:2993974‑M2/36‑4 (MQ=255)
   tATGCAAAGGACGTCACATTACCGTGCAGTACAgtt                                  <  1:5697581‑M2/36‑4 (MQ=255)
   tATGCAAAGGACGTCACATTACCGTGCAGTACAgtt                                  >  1:7561292‑M2/1‑33 (MQ=255)
   tATGCAAAGGACGTCACATTACCGTGCAGTACAgtt                                  >  1:7560250‑M2/1‑33 (MQ=255)
   tATGCAAAGGACGTCACATTACCGTGCAGTACAgtt                                  <  1:2117504‑M2/36‑4 (MQ=255)
    aTGCAAAGGACGTCACATTACCGTGCAGTACAgttg                                 >  1:787965‑M2/1‑32 (MQ=255)
     tGCAAAGGACGTCACATTACCGTGCAGTACAgttga                                >  1:6179844‑M2/1‑31 (MQ=255)
     tGCAAAGGACGTCACATTACCGTGCAGTACAgttga                                <  1:5265776‑M2/36‑6 (MQ=255)
     tGCAAAGGACGTCACATTACCGTGCAGTACAgttga                                <  1:2406417‑M2/36‑6 (MQ=255)
     tGCAAAGGACGTCACATTACCGTGCAGTACAgttga                                <  1:3873236‑M2/36‑6 (MQ=255)
     tGCAAAGGACGTCACATTACCGTGCAGTACAgttga                                >  1:6904806‑M2/1‑31 (MQ=255)
     tGCAAAGGACGTCACATTACCGTGCAGTACAgttga                                <  1:865735‑M2/36‑6 (MQ=255)
      gCAAAGGACGTCACATTACCGTGCAGTACAgttgat                               <  1:1377552‑M2/36‑7 (MQ=255)
      gCAAAGGACGTCACATTACCGTGCAGTACAgttgat                               <  1:1388365‑M2/36‑7 (MQ=255)
       cAAAGGACGTCACATTACCGTGCAGTACAgttgata                              <  1:2118056‑M2/36‑8 (MQ=255)
        aaaGTACGTCACATTACCGTGCAGTACAgttgatag                             >  1:1875404‑M2/1‑28 (MQ=255)
        aaaGGACGTCACATTACCGTGCAGTACAgttgatag                             >  1:2428069‑M2/1‑28 (MQ=255)
         aaGGACGTCACATTACCGTGCAGTACAGTTGATAGc                            <  1:1268636/36‑1 (MQ=255)
         aaGGACGTCACATTACCGTGCAGTACAGTTGATAGc                            >  1:6151573/1‑36 (MQ=255)
         aaGGACGTCACATTACCGTGCAGTACAGTTGATAGc                            <  1:752640/36‑1 (MQ=255)
         aaGGACGGCACATTACCGTGCAGTACAGTTGATAGc                            >  1:969205/1‑36 (MQ=255)
          aGGACGTCACATTACCGTGCAGTACAGTTGATAGcc                           <  1:5505963/36‑1 (MQ=255)
          aGGACGTCACATTACCGTGCAGTACAGTTGATAGcc                           >  1:531955/1‑36 (MQ=255)
          aGGACGGCAAATTACCGTGCAGTACAGTTGATAGcc                           >  1:6298657/1‑36 (MQ=37)
           ggACGTCACATTACCGTGCAGTACAGTTGATAGccc                          <  1:4436109/36‑1 (MQ=255)
           ggACGTCACATTACCGTGCAGTACAGTTGATAGccc                          <  1:112085/36‑1 (MQ=255)
            gACGTCACATTACCGTGCAGTACAGTTGATAGcccc                         <  1:6880215/36‑1 (MQ=255)
            gACGTCACATTACCGTGCAGTACAGTTGATAGcccc                         <  1:380890/36‑1 (MQ=255)
            gACGTCACATTACCGTGCAGTACAGTTGATAGcccc                         <  1:5912925/36‑1 (MQ=255)
            gACGTCACATTACCGTGCAGTACAGTTGATAGcccc                         <  1:1449424/36‑1 (MQ=255)
            gACGTCACATTACCGTGCAGTACAGTTGATAGcccc                         >  1:2957341/1‑36 (MQ=255)
             aCGTCACATTACCGTGCAGTACAGTTGATAGCCCCt                        <  1:5401977/36‑1 (MQ=255)
               gTCACATTACCGTGCAGTACAGTTGATAGCCCCTTc                      >  1:6563978/1‑36 (MQ=255)
               gTCACATTACCGTGCAGTACAGTTGATAGCCCCTTc                      >  1:7431057/1‑36 (MQ=255)
               gTCACATTACCGTGCAGTACAGTTGATAGCCCCTTc                      >  1:3898951/1‑36 (MQ=255)
                tCACATTACCGTGCAGTACAGTTGATAGCCCCTTcc                     <  1:4925811/36‑1 (MQ=255)
                tCACATTACCGTGCAGTACAGTTGATAGCCCCTTcc                     <  1:398202/36‑1 (MQ=255)
                tCACATTACCGTGCAGTACAGTTGATAGCCCCTTcc                     <  1:6069335/36‑1 (MQ=255)
                tCACATTACCGTGCAGTACAGTTGATAGCCCCTTcc                     <  1:2719519/36‑1 (MQ=255)
                tCACATTACCGTGCAGTACAGTTGATAGCCCCTTcc                     <  1:2397837/36‑1 (MQ=255)
                tCACATTACCGTGCAGTACAGTTGATAGCCCCTTcc                     <  1:6968578/36‑1 (MQ=255)
                tCACATTACCGTGCAGTACAGTTGATAGCCCCTTcc                     <  1:2086957/36‑1 (MQ=255)
                 cacaTTACCGTGCAGTACAGTTGATAGCCCCTTccc                    <  1:4107094/36‑1 (MQ=255)
                 cacaTTACCGTGCAGTACAGTTGATAGCCCCTTccc                    <  1:843992/36‑1 (MQ=255)
                  acaTTACCGTGCAGTACAGTTGATAGCCCCTTCCCa                   <  1:1085241/36‑1 (MQ=255)
                  acaTTACCGTGCAGTACAGTTGATAGCCCCTTCCCa                   >  1:2145441/1‑36 (MQ=255)
                  acaTTACCGTGCAGTACAGTTGATAGCCCCTTCCCa                   >  1:3191351/1‑36 (MQ=255)
                   caTTACCGTGCAGTACAGTTGATAGCCCCTTCCCAg                  <  1:4563751/36‑1 (MQ=255)
                   caTTACCGTGCAGTACAGTTGATAGCCCCTTCCCAg                  <  1:6040772/36‑1 (MQ=255)
                   caTTACCGTGCAGTACAGTTGATAGCCCCTTCCCAg                  <  1:4468205/36‑1 (MQ=255)
                     ttACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGt                <  1:2612266/36‑1 (MQ=255)
                     ttACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGt                <  1:4655484/36‑1 (MQ=255)
                     ttACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGt                <  1:4484454/36‑1 (MQ=255)
                     ttACCGGGCAGTACAGTTGATAGCCCCTTCCCAGGt                <  1:5005981/36‑1 (MQ=255)
                      tACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTa               <  1:6587927/36‑1 (MQ=255)
                       aCCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAg              >  1:6626515/1‑36 (MQ=255)
                        ccGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGc             <  1:4684398/36‑1 (MQ=255)
                         cGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCg            >  1:3711512/1‑36 (MQ=255)
                         cGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCg            <  1:5702476/36‑1 (MQ=255)
                          gTGCAGTACAgttgatagaaatatcggtcggggcct           <  1:6812794‑M2/36‑27 (MQ=255)
                           tGCAGTACAgttgatagaaatatcggtcggggcctt          >  1:1833950‑M2/1‑9 (MQ=255)
                           tGCAGTACAGTTGATAGCCCCTTCCCTGGTAGCggg          <  1:1717395/36‑1 (MQ=255)
                           tGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCggg          >  1:5786197/1‑36 (MQ=255)
                           tGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCggg          <  1:206713/36‑1 (MQ=255)
                           tGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCggg          >  1:5526900/1‑36 (MQ=255)
                           tGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCggg          >  1:1470062/1‑36 (MQ=255)
                             cAGTACAgttgatagaaatatcggtcggggccttcg        >  1:2363815‑M2/1‑7 (MQ=255)
                             cAGTACAgttgatagaaatatcggtcggggccttcg        <  1:4894178‑M2/36‑30 (MQ=255)
                             cAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGaa        >  1:3161190/1‑36 (MQ=255)
                              aGTACAGTTTATAGCCCCTTCCCAGGTAGCGGGAAg       <  1:1102435/36‑1 (MQ=255)
                              aGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAg       <  1:4294772/36‑1 (MQ=255)
                               atACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGc      <  1:2956718/35‑1 (MQ=255)
                                tgcAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCa     <  1:1860407/34‑1 (MQ=255)
                                tACAgttgatagaaatatcggtcggggccttcgtct     >  1:5900701‑M2/1‑4 (MQ=255)
                                tACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCa     >  1:3700905/1‑36 (MQ=255)
                                tACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCa     <  1:5776907/36‑1 (MQ=255)
                                 aCAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCat    <  1:3651674/36‑1 (MQ=255)
                                 aCAGTTCATAGCCCCTTCCCAGGTAGCGGGAAGCat    <  1:3990766/36‑1 (MQ=255)
                                  cAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCata   <  1:2182599/36‑1 (MQ=255)
                                  cAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCata   >  1:801712/1‑36 (MQ=255)
                                  cAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCata   <  1:4028377/36‑1 (MQ=255)
                                  aaGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCata   <  1:5447777/35‑1 (MQ=255)
                                   agttgatagaaatatcggtcggggccttcgtctttc  <  1:1508845‑M2/36‑36 (MQ=255)
                                   aGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCatat  <  1:2538301/36‑1 (MQ=255)
                                   |                                   
ATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATAT  >  REL606/3413887‑3413957

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 26 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG