New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 2448493 =23 (0.420)15 (0.270) 3/64 NT 39.5% intergenic (+24/+26) nupC/yfeA nucleoside (except guanosine) transporter/predicted diguanylate cyclase
?REL606 = 2774196 NA (NA)noncoding (1343/1343 nt) IS186 repeat region

CTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAG  >  REL606/2448458‑2448526
                                   |                                 
tcagatggttcaacccttaagttagcgcttatggga                                   <  1:4832192‑M1/1‑1 (MQ=255)
tcagatggttcaacccttaagttagcgcttatggga                                   >  1:5183058‑M1/36‑36 (MQ=255)
tcagatggttcaacccttaagttagcgcttatggga                                   >  1:6667605‑M1/36‑36 (MQ=255)
tcagatggttcaacccttaagttagcgcttatggga                                   <  1:6824994‑M1/1‑1 (MQ=255)
tcagatggttcaacccttaagttagcgcttatggga                                   >  1:1080774‑M1/36‑36 (MQ=255)
cTGGTGCTGTAAGACAATACATAAAAAAGCCGGGGa                                   <  1:1732939/36‑1 (MQ=255)
 cnnatggttcaacccttaagttngcgcttatgggAt                                  >  1:2071339‑M1/35‑36 (MQ=255)
 catatggttcaacccttaagttagcgcttatgggAt                                  >  1:2155472‑M1/35‑36 (MQ=255)
 cagatggttcaacccttaagttagcgcttatgggAt                                  >  1:4705316‑M1/35‑36 (MQ=255)
 cagatggttcaacccttaagttagcgcttatgggAt                                  <  1:2053423‑M1/2‑1 (MQ=255)
 cagatggttcaacccttaagttagcgcttatgggAt                                  <  1:1091181‑M1/2‑1 (MQ=255)
 cagatggttcaacccttaagttagcgcttatgggAt                                  <  1:3454589‑M1/2‑1 (MQ=255)
   gTGCTGTAAGACCATACATAAAAAAGCCGGGGATaa                                <  1:6175339/36‑1 (MQ=255)
    tGCTGTAAGACCATAAATAAAAAAGCCGGGGATAAt                               <  1:6865338/36‑1 (MQ=255)
      cTGTAAGACCATACATAAAAAAGCCGGGGATAATTc                             <  1:4384022/36‑1 (MQ=255)
        gTAAGACCATACATAAAAAAGCCGGGGATAATTacc                           >  1:2406523/1‑36 (MQ=255)
         tcaacccttaagttagcgcttatgggATAATTcccc                          <  1:1316355‑M1/10‑1 (MQ=255)
         tAAGACCATACATAAAAAAGCCGGGGATAATTcccc                          >  1:1792896/1‑36 (MQ=255)
         tAAGACCATACATAAAAAAGCCGGGGATAATTcccc                          >  1:4174910/1‑36 (MQ=255)
            gACCATACATAAAAAAGCCGGGGATAATTCCCCGGt                       <  1:170809/36‑1 (MQ=255)
              ccATACATAAAAAAGCCGGGGATAATTCCCCGGttt                     >  1:7062898/1‑36 (MQ=255)
               cttaagttagcgcttatgggATAATTCCCCGGtttt                    >  1:5452716‑M1/21‑36 (MQ=255)
               catacataAAAAAGCCGGGGATAATTCCCCGGtttt                    <  1:972148/36‑1 (MQ=255)
               catacataAAAAAGCCGGGGATAATTCCCCGGtttt                    <  1:2138667/36‑1 (MQ=255)
               catacataAAAAAGCCGGGGATAATTCCCCGGtttt                    <  1:4779782/36‑1 (MQ=255)
                ttaagttagcgcttatgggATAATTCCCCGGttttt                   >  1:891755‑M1/20‑36 (MQ=255)
                ttaagttagcgcttatgggATAATTCCCCGGttttt                   >  1:79624‑M1/20‑36 (MQ=255)
                 tacataAAAAAGCCGGGGATAATTCCCCGGtttttt                  >  1:6896833/1‑36 (MQ=255)
                  acataAAAAAGCCGGGGATAATTCCCCGGttttttt                 >  1:5959545/1‑36 (MQ=255)
                  acataAAAAAGCCGGGGATAATTACCCGGttttttt                 >  1:2696692/1‑36 (MQ=255)
                  acataAAAAAGCCCGGGATAATTCCCCGGttttttt                 >  1:2138016/1‑36 (MQ=255)
                   cataAAAAAGCCGGGGATAATTCCCCGGTTTTTTTa                <  1:553901/36‑1 (MQ=255)
                    ataAAAAAGCCGGGGATAATTCCCCGGTTTTTTTAc               <  1:7811671/36‑1 (MQ=255)
                         aaaGCCGGGGATAATTCCCCGGTTTTTTTACGCCTg          >  1:1834004/1‑36 (MQ=255)
                         aaaGCCGGGGATAATTCCCCGGTTTTTTTACGCCTg          >  1:5113257/1‑36 (MQ=255)
                           aGCCGGGGATAATTCCCCGGTTTTTTTACGCCTGtt        >  1:2534862/1‑36 (MQ=255)
                           aGCCGGGGATAATTCCCCGGTTTTTTTACGCCTGtt        >  1:4491948/1‑36 (MQ=255)
                             ccGGGGATAATTCCCCGGTTTTTTTACGCCTGTTaa      <  1:1164998/36‑1 (MQ=255)
                              cGGGGATAATTCCCCGGTTTTTTTACGCCTGTTAAt     >  1:7322326/1‑36 (MQ=255)
                                gggATAATTCCCCGGTTTTTTTACGCCTGTTAATCa   <  1:530462‑M1/33‑1 (MQ=255)
                                 ggATAATTCCCCGGTTTTTTTACGCCTGTTAATCAg  <  1:2318872‑M1/34‑1 (MQ=255)
                                   |                                 
CTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAG  >  REL606/2448458‑2448526

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 24 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG