New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	251809 =	136 (1.130)	7 (0.060)	6/552	11.1	5.2%	coding (912/1056 nt)	dinB	DNA polymerase IV
	?	NC_000913	251873 =	123 (1.050)					coding (976/1056 nt)	dinB	DNA polymerase IV
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

AACCGCAACCGGTTGATCAATAAAGTATTTAGCTGGGTTTTTACGTAACTCAGTGATATTGACCGATTTTTCAGCGAGAATTCGATGCATACAGTGATACCCTCATAATAATGCACACCAGAATATACATAATAGTATACATCATAATCCCAGCACCAGTTGTCTTTCCATTTGCGGGTCAAGCAACGTCACATGCAGCCCCACCAGACGCACACCGCGCCCGCCGCGGCGTTCATCCCAGGTTTTACGCGCGGTGGCGATTAGATCAGCTTTATTCAGCCGCGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/252094‑251809 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gccgcggcGGGCGCGGTGTGCGTCTGGTGGGGCTGCATGTGACG  >  NC_000913/251873‑251908                                                                                                                                                                                                                                                                                                                                    AACCGCAACCGGTTGATCAATAAAGTATTTAGCTGGGTTTTTACGTAACTCAGTGATATTGACCGATTTTTCAGCGAGAATTCGATGCATACAGTGATACCCTCATAATAATGCACACCAGAATATACATAATAGTATACATCATAATCCCAGCACCAGTTGTCTTTCCATTTGCGGGTCAAGCAACGTCACATGCAGCCCCACCAGACGCACACCGCGCCCGCCGCGGCGTTCATCCCAGGTTTTACGCGCGGTGGCGATTAGATCAGCTTTATTCAGCCGCGGC                                      >  1:559874/1‑286   gcgaAACCGATTGATAAGTAAAGGAGTTAGCTGGGTTTTTACGTAACTCAGTGATATTGACCGATTTTTCATCGAGGATTCGATGCATACAGTGATACCCTCATAATAATGCACACCAGAATATACATAATAGTATACATCATAATCCCAGCACCAGTTGTCTTTCCATTTGCGGGTCAAGCAACGTCACATGCAGCCCCACCAGACGCACACCGCGCCCGCCGCGGCGTTCATCCCAGGTTTTACGCGCGGTGGCGATTAGATCAGCTTTATTCAGCCGCGGCGGGCG                                 <  2:78752/285‑1                                                                                                       TCATAATAATGCACACCAGAATATACATAATAGTATACATCATAATCCCAGCACCAGTTGTCTTTCCATTTGCGGGTCAAGCAACGTCACATGCAGCCCCACCAGACGCACACCGCGCCCGCCGCGGCGTTCATCCCAGGTTTTACGCGCGGTGGCGATTAGATCAGCTTTATTCAGCCGCGGCGGGCGCGGT                             >  1:37670/1‑193                                                                                                       TCATAATAATGCACACCAGAATATACATAATAGTATACATCATAATCCCAGCACCAGTTGTCTTTCCATTTGCGGGTCAAGCAACGTCACATGCAGCCCCACCAGACGCACACCGCGCCCGCCGCGGCGTTCATCCCAGGTTTTACGCGCGGTGGCGATTAGATCAGCTTTATTCAGCCGCGGCGGGCGCGGT                             <  2:37670/193‑1                                                                                                                                    ATAGTATACATCATAATCCCAGCACCAGTTGTCTTTCCATTTGCGGGTCAAGCAACGTCACATGCAGCCCCACCAGACGCACACCGCGCCCGCCGCGGCGTTCATCCCAGGTTTTACGCGCGGTGGCGATTAGATCAGCTTTATTCAGCCGCGGCGGGCG                                 <  1:301647/160‑1                                                                                                                                    ATAGTATACATCATAATCCCAGCACCAGTTGTCTTTCCATTTGCGGGTCAAGCAACGTCACATGCAGCCCCACCAGACGCACACCGCGCCCGCCGCGGCGTTCATCCCAGGTTTTACGCGCGGTGGCGATTAGATCAGCTTTATTCAGCCGCGGCGGGCG                                 >  2:301647/1‑160                                                                                                                                               CATAATCCCAGCACCAGTTGTCTTTCCATTTGCGGGTCAAGCAACGTCACATGCAGCCCCACCAGACGCACACCGCGCCCGCCGCGGCGTTCATCCCAGGTTTTACGCGCGGTGGCGATTAGATCAGCTTTATTCAGCCGCGGCGGGCGCGGTGTGCGTCTGGTGGGGCTGCATGTGACG  <  1:24203/180‑1                                                                                                                                               CATAATCCCAGCACCAGTTGTCTTTCCATTTGCGGGTCAAGCAACGTCACATGCAGCCCCACCAGACGCACACCGCGCCCGCCGCGGCGTTCATCCCAGGTTTTACGCGCGGTGGCGATTAGATCAGCTTTATTCAGCCGCGGCGGGCGCGGTGTGCGTCTGGTGGGGCTGCATGTGACG  >  2:24203/1‑180                                                                                                                                                                                                                                                                                                                                    AACCGCAACCGGTTGATCAATAAAGTATTTAGCTGGGTTTTTACGTAACTCAGTGATATTGACCGATTTTTCAGCGAGAATTCGATGCATACAGTGATACCCTCATAATAATGCACACCAGAATATACATAATAGTATACATCATAATCCCAGCACCAGTTGTCTTTCCATTTGCGGGTCAAGCAACGTCACATGCAGCCCCACCAGACGCACACCGCGCCCGCCGCGGCGTTCATCCCAGGTTTTACGCGCGGTGGCGATTAGATCAGCTTTATTCAGCCGCGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/252094‑251809 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gccgcggcGGGCGCGGTGTGCGTCTGGTGGGGCTGCATGTGACG  >  NC_000913/251873‑251908

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