New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	168704 =	0 (0.000)	7 (0.420)	5/276	1.1	100%	coding (1221/2244 nt)	fhuA	ferrichrome outer membrane transporter/phage receptor
	?	NC_000913	182989 =	0 (0.000)					coding (527/1158 nt)	cdaR	DNA‑binding transcriptional activator CdaR

TGCTGAAAACAAAACCTCGCAAAACAGCGTTTATGGTTACGGCGTCTGCTCCGATCCGGCGAATGCTTACAGCAAACAGTGTGCGGCATTAGCGCCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCTGTACAATCCGGTGAATACCGATTTCGACTTCAATGCCAAAGATCCGGCAAACTCCGGCCCTTACCGCATTCTGAATAAACAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCTGGTCACCCTAGGCGG  >  NC_000913/168488‑168934                                                                                                                                                                                                                         |                                                                                                                                                                                                                                       cgcgaaaacgcagctggccgtactctttcatgcgggtaatcagttgttcaactcgattacgatgatcttctgcataccagcgcccaacagagttccacgccggtttcaacaccaccatttcggttagcgagacaatcgccaccagattattacgctcgggcgtagtcagcgcgttttgcagttgttgtaactccgccattgcgctgtccacgccaaGCTGACCGGTGTCGACTTTATGCGTATGCGTAATg                                                                                                                                                                                                      <  1:75733‑M1/35‑1 (MQ=255)                         ctgtcatgcgggtaatcagttgttcaactcgcttacgatgatcttctgcatcccagcgcccaaaagagttcaacgccggtttcaacaccaccatttcggttagcgagacaatcgccaccagatgagtacgctcgggcgtagtcagcgcgttgtgcagttgttgtaacgccgccttggcgctgtccacgccaaGCTGACCGGTGTc                                                                                                                                                                                                                            >  1:387910‑M1/193‑205 (MQ=255)                                                                                                       gtttcaacacccccatttcggttagcgcgacaatcgccaccagattattacgctcgggcgtagtcagcgcgttttgcagttgttgtaactccgccattgcgctgtccacgccaaGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCTGTACAATCCGGTGAATACCGATTTCGACTTCAATGc                                                                                                                 <  1:136822‑M1/120‑1 (MQ=255)                                                                                                                           gttagcgagacaatcgccaccagattattacgctcgggcgtagtcagcgcggtttgcatgtgttgtaactccgccagtgcgctgtccacgccaaGCTGACCGGTGTCGACTGTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCAc                                                                                                                                                               >  1:365010‑M1/95‑168 (MQ=255)                                                                                                                                                                gcgtagtcagcgcgttttgcagttgttgtaactccgccattgcgctgtccacgccaaGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCTGTACAATCCGGTGAATACCGATGTCGACTTCAATGCCAAAGATCCGGCAAACTCCGGCCCTTACCGCATTCTGAATAAACAGAAACAAACGGGCGTTTATGTTCAGGATc                                       >  1:445265‑M1/58‑251 (MQ=255)                                                                                                                                                                    agtcagcgcgttttgcagttgttgtaactccgccattgcgctgtccacgccaaGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCTGTACAATCCGGTGAATACCGATTTCGACTTCAATGCCAAAGATCCGGCAAACTCCGGCCCTTACCGCATTCTGAATAAACAGAAACAAACGGGCGTTTATGTTCAGGATCagn                                    >  1:40830‑M1/54‑249 (MQ=255)                                                                                                                                                                                                     nattgcgctgtccacgccaaGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCTGTACAATCCGGTGAATACCGATTTCGACTTCAATGCCAAAGATCCGGCAAACTCCGGCCCTTACCGCATTCTGAATAAACAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCTGGTCACCCTAGGCgg  <  1:359113‑M1/231‑1 (MQ=255)                                                                                                                                                                                                                         |                                                                                                                                                                                                                                       TGCTGAAAACAAAACCTCGCAAAACAGCGTTTATGGTTACGGCGTCTGCTCCGATCCGGCGAATGCTTACAGCAAACAGTGTGCGGCATTAGCGCCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCTGTACAATCCGGTGAATACCGATTTCGACTTCAATGCCAAAGATCCGGCAAACTCCGGCCCTTACCGCATTCTGAATAAACAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCTGGTCACCCTAGGCGG  >  NC_000913/168488‑168934

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