New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	3303564 =	53 (1.030)	3 (0.070)	3/134	5.0	5.8%	coding (798/1368 nt)	dcuD	predicted transporter
	?	REL606	3303584 =	50 (1.080)					coding (818/1368 nt)	dcuD	predicted transporter
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATAACTTGCGCGAAAC  >  REL606/3303512‑3303659                                                                         |                                                                            tgATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCtggt                                                                          >  1:220047/1‑76 (MQ=255) tgATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCtggt                                                                          >  1:657581/1‑76 (MQ=255)      ccGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCtggtggtgg                                                                     <  2:34048/76‑1 (MQ=255)       cGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTggtggtggt                                                                    >  2:1005384/1‑76 (MQ=255)         ttAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGa                                                                  >  1:492438/1‑76 (MQ=255)          tAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGAt                                                                 <  1:1333171/76‑1 (MQ=255)          tAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCCGATGCAGTCAGAACTGCATCTGGTGGTGGGGAt                                                                 >  1:73349/1‑76 (MQ=255)           aaTCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATg                                                                <  2:647875/76‑1 (MQ=255)             tCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGtt                                                              <  1:1397849/76‑1 (MQ=255)             tCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGtt                                                              <  1:1645805/76‑1 (MQ=255)               cTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTAc                                                            <  1:1298203/76‑1 (MQ=255)                 gATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTg                                                          <  2:699071/76‑1 (MQ=255)                 gATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGGGGTGAGGTTGCTg                                                          >  2:240185/1‑76 (MQ=255)                   tGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAg                                                        <  1:1552222/76‑1 (MQ=255)                      tcggctcgCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGttt                                                     <  1:680208/76‑1 (MQ=255)                      tcggctcgCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGttt                                                     <  1:777774/76‑1 (MQ=255)                      tcggctcgCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGttt                                                     >  2:1548903/1‑76 (MQ=255)                       cggctcgCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTg                                                    >  2:86649/1‑76 (MQ=255)                       cggctcgCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTg                                                    >  2:60086/1‑76 (MQ=255)                       cggctcgCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGGGGTGATGTTACTGAGTTTg                                                    >  2:896324/1‑76 (MQ=255)                            cgCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACtgt                                               <  1:29563/76‑1 (MQ=255)                             gCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACtgtg                                              >  1:1145968/1‑76 (MQ=255)                             gCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACtgtg                                              <  2:1526113/76‑1 (MQ=255)                                gTTCCTCGCCCACGTCGGGCTGATGCAGTCCGAACTGCATCTGGTGGGGGTGATGTTACTGAGTTTGACTGTGACg                                           >  2:914835/1‑76 (MQ=255)                                gTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACg                                           >  2:1318705/1‑76 (MQ=255)                                gTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACg                                           >  2:772954/1‑76 (MQ=255)                                gTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACg                                           >  1:900152/1‑76 (MQ=255)                                gTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACg                                           >  2:529258/1‑76 (MQ=255)                                gTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACg                                           >  2:866571/1‑76 (MQ=255)                                gTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACg                                           >  2:249813/1‑76 (MQ=255)                                gTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACg                                           >  2:1625647/1‑76 (MQ=255)                                gTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACg                                           >  2:1149339/1‑76 (MQ=255)                                gTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGAGATGTTACTGAGCTTGACTGTGACg                                           >  2:1609375/1‑76 (MQ=255)                                gTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGGGGTGATGTTACTGAGTTTGACTGTGACg                                           >  2:1471338/1‑76 (MQ=255)                                      cGCCCACGTCGGGCTGATGCAGTCAGAACTGCCTCTGGCGGTGGCGATGTTACTGAGTTTGCCTGCGCCGAcgttt                                     >  2:185804/1‑76 (MQ=255)                                      cGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGAtgttgtgt                                     >  2:1657220/1‑74 (MQ=255)                                      cGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGAtgttt                                     >  2:791475/1‑76 (MQ=255)                                      cGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGAtgttt                                     >  2:976079/1‑76 (MQ=255)                                          cACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTttgttg                                 <  1:959906/76‑1 (MQ=255)                                               cGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGttc                            <  2:586141/76‑1 (MQ=255)                                                   cTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTcc                        <  1:1140740/76‑1 (MQ=255)                                                     gATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGc                      >  2:1358507/1‑76 (MQ=255)                                                      aTGCAGTCAGAACTTCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCa                     <  2:477438/76‑1 (MQ=255)                                                      aTGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCa                     <  2:1325122/76‑1 (MQ=255)                                                      aTGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCAGCa                     <  2:9419/76‑1 (MQ=255)                                                         cAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGc                  <  2:1246247/76‑1 (MQ=255)                                                         cAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGc                  <  2:99131/76‑1 (MQ=255)                                                            tCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATa               <  2:707977/76‑1 (MQ=255)                                                             cAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATaa              >  1:956379/1‑76 (MQ=255)                                                              aGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATAAc             <  2:133586/76‑1 (MQ=255)                                                                  cTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATAACTTgc         <  2:73349/76‑1 (MQ=255)                                                                   tGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATAACTTgcg        <  1:1585087/76‑1 (MQ=255)                                                                      aTCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATAACTTGCGCGa     >  1:548255/1‑76 (MQ=255)                                                                      aTCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATAACTTGCGCGa     <  2:436091/76‑1 (MQ=255)                                                                         tggtggtggTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATAACTTGCGCGAAAc  <  1:821834/76‑1 (MQ=255)                                                                         tggtggtggTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATAACTTGCGCGAAAc  <  2:223266/76‑1 (MQ=255)                                                                         |                                                                            TGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATAACTTGCGCGAAAC  >  REL606/3303512‑3303659

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