New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 550165	116 (0.970)	3 (0.030)	3/512	NT	2.2%	coding (539/816 nt)	ylbF	putative anaerobic allantoin catabolic oxamate carbamoyltransferase; DUF2877 family protein
	?	NC_000913	= 551959	147 (1.240)					coding (636/1068 nt)	purK	N5‑carboxyaminoimidazole ribonucleotide synthase
Rejected: Frequency below/above cutoff threshold.

GGGCGTTTCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/549876‑550165 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cggcGATTATGCAGGAGCTGGGCTATGTGGGCGTGATGGCGATGGAGTGTTTTGTCACCCCGCAAGGTCTGTTGATCAACGAACTGGCACCGCGTGTGCATAACAGCGGTCACTGGACACAAAACGGTGCCAGCATCAGCCAGTTTGAGCTGCATCTGCGGGCGATTACCGATCTGCCGTTACCGCAACCAGTGGTGAATAATCCGTCGGTGATGATCAATCTGATTGGTAGCGATGTGAATTATGACTGGCTGAAACTGCCGCTGGTGCATCTGCACTGGTACGACAAA  <  NC_000913/551959‑551674                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  GGGCGTTTCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGCGTCATGATGACACGCTGAGCGGCATGCTGCTGGCGG                                                                                                                                                                                                                                                                                                 >  1:121462/1‑289                                                                          TTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCGATTATGCAGGAGCTGGGCTATGTGGGCGTGATGGCGATGGAGTGTTTTGTCACCCCGCAAGGTCTGTTGA                                                                                                                                                                                                                         >  1:203713/1‑288                                                                                                                              GGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCGATTATGCAGGAGCTGGGCTATGTGGGCGTGATGGCGATGGAGTGTTTTGTCACCCCGCAAGGTCTGTTGATCAACGAACTGGCACCGCGTGTGCATAACAGCGGTCACTGGACACAAAACGGT                                                                                                                                                                    <  2:203713/289‑1                                                                                                                                                                                                                                                                                        tgGCTGGCGGCCATTATGCAGGGGCTGGGCTATGTGGTGGGGTTGGCGATGGAATGTTTTGTCACCCCCCAAGGTCTGTTGAGCAACGAACTGGCACCGCGGGTGCATAACAGCGGTCACTGGACACAAAACGGTGCCAGCATCAGCCAGTTTGAGCTGCATCTGCAGGCGATTACCGATCTGCCGTTACCGCAACCAGTGGTGAATAATCCGTCGGTGATTATCAATCTGATTGGTAGCGATGTGAATTATGACTGGCTGAAACTGCCGCTGGTGCATCTGCACTGG           <  2:239319/286‑1                                                                                                                                                                                                                                                                                        aTGCTGGCGGGGTTTATGAAGGTGCTGGGCTATGTGGGCGTAATGGCGATGTAGTGTTTTGTCACCCCGCAAGGTCTGTTGATCAACGAACTGGCACCGCGTGTGCATAACAGCGGTCACTGTACACAAAACGGTGCCAGCATCAGCCAGTTTGAGCTGCATCTGCGGTCGATTACCGATCTGCCGTTACCGCAACCAGTTGTGAATAATCCGTCGGTGATGGTCAATCTGATTGGTAGCGATGTGAATTATGACTGGCTGAAACTGCCGCTGGTGCATCTGCACTGG           <  1:452708/287‑1                                                                                                                                                                                                                                                                                                GGCGATTATGCAGGAGCTGGGCTATGTGGGCGTGATGGCGATGGAGTGTTTTGTCACCCCGCAAGGTCTGTTGATCAACGAACTGGCACCGCGTGTGCATAACAGCGGTCACTGGACACAAAACGGTGCCAGCATCAGCCAGTTTGAGCTGCATCTGCGGGCGATTACCGATCTGCCGTTACCGCAACCAGTGGTGAATAATCCGTCGGTGATGATCAATCTGATTGGTAGCGATGTGAATTATGACTGGCTGAAACTGCCGCTGGTGCATCTGCACTGGTACGACAAA  >  1:426670/1‑289                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  GGGCGTTTCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/549876‑550165 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cggcGATTATGCAGGAGCTGGGCTATGTGGGCGTGATGGCGATGGAGTGTTTTGTCACCCCGCAAGGTCTGTTGATCAACGAACTGGCACCGCGTGTGCATAACAGCGGTCACTGGACACAAAACGGTGCCAGCATCAGCCAGTTTGAGCTGCATCTGCGGGCGATTACCGATCTGCCGTTACCGCAACCAGTGGTGAATAATCCGTCGGTGATGATCAATCTGATTGGTAGCGATGTGAATTATGACTGGCTGAAACTGCCGCTGGTGCATCTGCACTGGTACGACAAA  <  NC_000913/551959‑551674

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