New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1164698 =	32 (1.040)	4 (0.150)	4/114	2.2	11.5%	coding (604/1026 nt)	nagZ	beta N‑acetyl‑glucosaminidase
	?	NC_000913	1165027 =	34 (1.290)					coding (933/1026 nt)	nagZ	beta N‑acetyl‑glucosaminidase
Rejected: Frequency below/above cutoff threshold.

ACCGCTCGCCGGACGCGGATCAACATCACTGTAGATCACATGCGCAGGCATAATGGCGTCGAGTTTATTTTCGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/1164772‑1164698 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑attttcgcgACAGGAACTGATGGACTCGGCTCGCTGGAAAGCGATCAGCACCCGTCTGAATCAGTTACATG  >  NC_000913/1165027‑1165088                                                                                                                                           ACCGCTCGCCGGACGCGGATCAACATCACTGTAGATCACATGCGCAGGCATAATGGCGTCGAGTTTAT                                                                       >  1:122401/1‑68 ACCGCTCGCCGGACGCGGATCAACATCACTGTAGATCACATGCGCAGGCATAATGGCGTCGAGATTAT                                                                       >  2:118346/1‑68                       ACATCACTGTAGATCACATGCGCAGGCATAATGGCGTCGAGTTCATTTTCGCGACAGGAACTGATGGA                                                 <  1:257813/68‑1                                 AGATCACATGCGCAGGCATAATGGCGTCGAGTTTATTTTCGCGACAGGAACTGATGGACTCGGCTCGC                                       <  2:159794/68‑1                                            GCAGGCATAATGGCGTCGAGTTTATTTTCGCGACAGGAACTGATGGACTCGGCTCGCTGGAAAGCGAT                            >  1:844284/1‑68                                                      TGGCGTCGAGTTCATTTTCGCGACAGGAACTGATGGACTCGGCTCGCTGGAAAGCGATCAGCACCCGT                  <  2:618067/68‑1                                                                     TTTCGCGACAGGAACTGATGGACTCGGCTCGCTGGAAAGCGATCAGCACCCGTCTGAATCAGTTACAT   >  1:1137299/1‑68                                                                      TTCGCGACAGGAACTGATGGACTCGGCTCGCTGGAAAGCGATCAGCGCCCGTCTGAATCAGTTACATG  >  1:219331/1‑68                                                                                                                                           ACCGCTCGCCGGACGCGGATCAACATCACTGTAGATCACATGCGCAGGCATAATGGCGTCGAGTTTATTTTCGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/1164772‑1164698 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑attttcgcgACAGGAACTGATGGACTCGGCTCGCTGGAAAGCGATCAGCACCCGTCTGAATCAGTTACATG  >  NC_000913/1165027‑1165088

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