New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4620598 =	105 (1.050)	5 (0.050)	5/166	NT	5.3%	coding (996/1224 nt)	deoB	phosphopentomutase
	?	NC_000913	4620626 =	82 (0.880)					coding (1024/1224 nt)	deoB	phosphopentomutase
Rejected: Frequency below cutoff threshold.

AGTACCGGAATGTGTTCACGCGTGTGGTCAGTACCGGTCCAGGTCGGATCGCAACCGTGGTCAGCGGTGAGGATCAGGATGTCGTCATCGCGCAGCAGAGACATC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4620702‑4620598 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gacatcCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAA  >  NC_000913/4620626‑4620723                                                                                                                                                                                                             AGTACCGGAATGTGTTCACGCGTGTGGTCAGTACCGGTCCAGGTCGGATCGCAACCGTGGTCAGCGGTGAGGATCAGGATGTCGTCATCGCGCAGCAGAG                                                                                                         <  2:182644/100‑1   TACCGGAATGTGTTCACGCGTGTGGTCAGTACCGGTCCAGGTCGGATCGCAACCGTGGTCAGCGGTGAGGATCAGGATGTCGTCATCGCGCAGCAGAGAC                                                                                                       <  1:2475643/100‑1   TACCGGAATGTGTTCACGCGTGTGGTCAGTACCGGTCCAGGTCGGATCGCAACCGTGGTCAGCGGTGAGGATCAGGATGTCGTCATCGCGCAGCAGAGAC                                                                                                       <  2:2709382/100‑1    ACCGGAATGTGTTCACGCGTGTGGTCAGTACCGGTCCAGGTCGGATCGCAACCGTGGTCAGCGGTGAGGATCAGGATGTCGTCATCGCGCAGCAGAGACA                                                                                                      >  1:786577/1‑100    ACCGGAATGTGTTCACGCGTGTGGTCAGTACCGGTCCAGGTCGGATCGCAACCGTGGTCAGCGGTGAGGATCAGGATGTCGTCATCGCGCAGCAGAGACA                                                                                                      >  2:114983/1‑100     CCGGAATGTGTTCACGCGTGTGGTCAGTACCGGTCCAGGTCGGATCGCAACCGTGGTCAGCGGTGAGGATCAGGATGTCGTCATCGCGCAGCAGAGACAT                                                                                                     >  1:2030358/1‑100     CCGGAATGTGTTCACGCGTGTGGTCAGTACCGGTCCAGGTCGGATCGCAACCGTGGTCAGCGGTGAGGATCAGGATGTCGTCATCGCGCAGCAGAGACAT                                                                                                     >  2:377753/1‑100      CGGAATGTGTTCACGCGTGTGGTCAGTACCGGTCCAGGTCGGATCGCAACCGTGGTCAGCGGTGAGGATCAGGATGTCGTCGTCGCGCAGCAGAGACATC                                                                                                    <  2:1977071/100‑1      CGGAATGTGTTCACGCGTGTGGTCAGTACCGGTCCAGGTCGGATCGCAACCGTGGTCAGCGGTGAGGATCAGGATGTCGTCATCGCGCAGCAGAGACATC                                                                                                    <  2:242150/100‑1                                                                         ATCAGGATGTCGTCATCGCGCAGCAGAGACATCCTGATCCTCACCGCTGACCACGGTTGCGAACCGACATGGACCGGTACTGACCCCACGCGTGAACACA                                 >  2:2736713/1‑100                                                                            AGGATGTCGTCATCGCGCAGCAGAGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTC                              >  1:890350/1‑100                                                                             GGATGTCGTCATCGCGCAGCAGAGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCC                             >  1:1296123/1‑100                                                                               ATGTCGTCATCGCGCAGCAGAGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGG                           >  1:617586/1‑100                                                                                   CGTCATCGCGCAGCAGAGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACT                       >  2:466777/1‑100                                                                                                    GACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAG      <  1:1548680/100‑1                                                                                                    GACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAG      <  2:593772/100‑1                                                                                                      CATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTA    >  2:282612/1‑100                                                                                                      CATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTA    <  2:1759245/100‑1                                                                                                      CATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTA    <  1:666611/100‑1                                                                                                       ATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAA   <  2:1221967/100‑1                                                                                                       ATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTA    <  2:503324/99‑1                                                                                                        TCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGGAAGTAAA  >  2:671353/1‑100                                                                                                        TCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAA  >  2:2743764/1‑100                                                                                                                                                                                                             AGTACCGGAATGTGTTCACGCGTGTGGTCAGTACCGGTCCAGGTCGGATCGCAACCGTGGTCAGCGGTGAGGATCAGGATGTCGTCATCGCGCAGCAGAGACATC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4620702‑4620598 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gacatcCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAA  >  NC_000913/4620626‑4620723

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ 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.