New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4513364	139 (0.860)	6 (0.040)	5/248	NT	4.3%	coding (46/999 nt)	fecC	ferric citrate ABC transporter permease
	?	NC_000913	= 4513381	138 (0.900)					coding (29/999 nt)	fecC	ferric citrate ABC transporter permease
Rejected: Frequency below/above cutoff threshold.

ACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4513213‑4513364 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gtgctgcTGTGGGGGCTTCCCGTTGCAGCACTTATTATTATTTTCTGGCTGAGTCTGTTTTGCTACTCGGCCATTCCTGTTTCCGGAGCAGATGCAACCCGCGCCCTGCTGCCTGGACACACGCCAACGCTACCAGAAGCGCTGGTGCAAAACCT  <  NC_000913/4513381‑4513234                                                                                                                                                                                                                                                                                                              ACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGC                                                                                                                                                         <  3:41452/149‑1   CAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTG                                                                                                                                                       >  2:1389872/1‑149   CAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTG                                                                                                                                                       >  3:1423011/1‑149                               ACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCTGTGGGGGCTTCCCGTTGCAGCA                                                                                                                               <  3:780642/145‑1                               ACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCTGTGGGGGCTTCCCGTTGCAGCA                                                                                                                               >  4:780642/1‑145                                                               AGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCTGTGGGGGCTTCCCGTTGCAGCACTTAT                                                                                                                          >  1:1176035/1‑118                                                               AGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCTGTGGGGGCTTCCCGTTGCAGCACTTAT                                                                                                                          <  2:1176035/118‑1                                                                   AGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCTGTGGGGGCTTCCCGTTGCAGCA                                                                                                                               <  1:1265365/109‑1                                                                   AGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCTGTGGGGGCTTCCCGTTGCAGCA                                                                                                                               >  2:1265365/1‑109                                                                                                                                                        CTGTGGGGGCTTCCCGTTGCAGCACTTATTATTATTTTCTGGCTGAGTCTGTTTTGCTACTCGGCCATTCCTGTTTCCGGAGCAGATGCAACCCGCGCCCTGCTGCCTGGACACACGCCAACGCTACCAGAAGCGCTGGTGCAAAACCT  <  1:2125/149‑1                                                                                                                                                                                                                                                                                                              ACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4513213‑4513364 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gtgctgcTGTGGGGGCTTCCCGTTGCAGCACTTATTATTATTTTCTGGCTGAGTCTGTTTTGCTACTCGGCCATTCCTGTTTCCGGAGCAGATGCAACCCGCGCCCTGCTGCCTGGACACACGCCAACGCTACCAGAAGCGCTGGTGCAAAACCT  <  NC_000913/4513381‑4513234

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