New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 16728NA (NA)4 (0.100) 4/60 NT 15.4% noncoding (1343/1343 nt) IS186 repeat region
?REL606 2448493 = 22 (0.490)intergenic (+24/+26) nupC/yfeA nucleoside (except guanosine) transporter/predicted diguanylate cyclase

TGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCT  >  REL606/2448459‑2448528
                                  |                                   
tGGTGCTGTAAGACCATACATAAAAAAGCCGGGGAt                                    >  1:5518672/1‑36 (MQ=255)
 ggTGCTGTAAGACCATACATAAAAAAGCCGGGGATa                                   >  1:4514692/1‑36 (MQ=255)
    gCTGTAAGACCATACATAAAAAAGCCGGGGATAAtt                                <  1:4427017/36‑1 (MQ=255)
      gttcaaccgttaagttagcgcttatgggATAATTcc                              <  1:5585479‑M2/8‑1 (MQ=255)
       gTAAGACCATACATAAAAAAGCCGGGGATAATTcac                             <  1:3297566/36‑3 (MQ=255)
        tcaacccttaagttagcgcttatgggATAATTcccc                            >  1:1968580‑M2/27‑36 (MQ=255)
        tAAGACCATACATAAAAAAGCCGGGGATAATTcccc                            >  1:1284227/1‑36 (MQ=255)
         aaGACCATACATAAAAAAGCCGGGGATAATTCCCCg                           >  1:1962433/1‑36 (MQ=255)
         aaGACCATACATAAAAAAGCCGGGGATAATTCCCCg                           >  1:4947474/1‑36 (MQ=255)
         aaGACCATACATAAAAAAGCCGGGGATAATTACCCg                           >  1:5551546/1‑36 (MQ=255)
            aCCATACATAAAAAAGCCGGGGATAATTCCCCGGtt                        >  1:3742838/1‑36 (MQ=255)
             ccATACATAAAAAAGCCGTGGATAATTCCCCGGttt                       <  1:4468196/36‑1 (MQ=255)
              catacataAAAAAGCCGGGGATAATTCCCCGGtttt                      <  1:5646177/36‑1 (MQ=255)
              catacataAAAAAGCCGGGGATAATTCCCCGGtttt                      >  1:1974345/1‑36 (MQ=255)
              catacataAAAAAGCCGGGGATAATTCCCCGGtttt                      <  1:3636548/36‑1 (MQ=255)
                tacataAAAAAGCCGGGGATCATTCCCCGGtttttt                    >  1:2334400/1‑36 (MQ=255)
                tacataAAAAAGCCGGGGATAATTCCCCGGtttttt                    >  1:5528173/1‑36 (MQ=255)
                 acataAAAAAGCCGGGGATAATTCCCCGGttttttt                   >  1:3365141/1‑36 (MQ=255)
                 aagttagcgcttatgggATAATTCCCCGGttttttt                   <  1:5170960‑M2/19‑1 (MQ=255)
                  cataAAAAAGCCGGGGATAATTCCCCGGTTTTTTTa                  <  1:2099502/36‑1 (MQ=255)
                     aaaaaaGCCGGGGATAATTCCCCGGTTTTTTTACGc               >  1:2193830/1‑36 (MQ=255)
                      agcgcttatgggATAATTCCCCGGTTTTTTTACGcc              <  1:3033985‑M2/24‑1 (MQ=255)
                      aaaaaGCCGGGGATAATTCCCCGGTTTTTTTACgtc              >  1:3569360/1‑34 (MQ=255)
                       aaaaGCCGGGGATAATTCCCCGGTTTTTTTACGCCt             >  1:385586/1‑36 (MQ=255)
                       aaaaGCCGGGGATAATTCCCCGGTTTTTTTACGCCt             <  1:4423992/36‑1 (MQ=255)
                           gCCGGGGATAATTCCCCGGTTTTTTTACGCCTGTTa         <  1:3139919/36‑1 (MQ=255)
                            ccGGGGATAATTCCCCGGTTTTTTTACGCCTGTTaa        <  1:5557498/36‑1 (MQ=255)
                              ggggATAATTCCTCGGTTTTTTTACGCCTGTTAATc      <  1:13535/36‑1 (MQ=255)
                                 gATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGc   <  1:1393054‑M2/35‑1 (MQ=255)
                                  ATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCt  <  1:3987004‑M2/36‑1 (MQ=255)
                                  |                                   
TGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCT  >  REL606/2448459‑2448528

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 17 ≤  ATCG  < 29 ≤  ATCG  < 37 ≤  ATCG