New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3614604	17 (0.850)	5 (0.280)	3/118	1.8	23.7%	coding (938/1575 nt)	nikA	nickel/heme ABC transporter periplasmic binding protein
	?	NC_000913	= 3615035	17 (0.950)					coding (1369/1575 nt)	nikA	nickel/heme ABC transporter periplasmic binding protein

CTCTTAATTACGCGGTAAACAAAAAATCGCTGATTGATAACGCGTTGTATGGCACCCAGCAGGTCGCCGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3614535‑3614604 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tcgccgaTCTCTTTATCAATCAGCGGTTTGTCGGCTAATCCTTGCTGTGCCTGGAAGTCAGCGTGTGACG  <  NC_000913/3615035‑3614973                                                                                                                                       CTCTTAATTACGCTGTAAACAAAAAATCGCTGATTGATAACGCGTTGTATGGCACCCAGCAGGTCGCC                                                                   <  2:237516/68‑1 CTCTTAATTACGCGGTAAACAAAAAATCGCTGATTGATAACGCGTTGTATGGCACCCAGCAGGTCGCC                                                                   >  1:351206/1‑68 CTCTTAATTACGCGGTAAACAAAAAATCGCTGATTGATAACGCGTTGTATGGCACCCAGCAGGCCGCC                                                                   <  1:51821/68‑1               GTAAACAAAAAATCTCTGATTGATAACGCGTTGTATGGCACCCAGCAGGTCGCCGATCTCTTTATCAA                                                     <  2:286623/68‑1               GTAAACAAAAAATCGCTGATTGATAACGCGTTGTATGGCACCCAGCAGGTCGCCGATCTCTTTATCAA                                                     <  1:457339/68‑1                        AAATCGCTGATTGATAACGCGTTGTATGGCACCCAGCAGGTCGCCGATCTCTTTATCAA                                                     <  1:358902/59‑1                        AAATCGCTGATTGATAACGCGTTGTATGGCACCCAGCAGGTCGCCGATCTCTTTATCAA                                                     >  2:358903/1‑59                                       AACGCGTTGTATGGCACCCAGCAGGTCGCCGATCTCTTTATCAATCAGCGGTTTGTCGGCTAATCCTT                             >  1:378879/1‑68                                                                  GCCGATCTCTTTATCAATCAGCGGTTTGTCGGCTAATCCTTGCTGTGCCTGGAAGTCAGCGTGTGACG  >  2:526104/1‑68                                                                                                                                       CTCTTAATTACGCGGTAAACAAAAAATCGCTGATTGATAACGCGTTGTATGGCACCCAGCAGGTCGCCGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3614535‑3614604 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tcgccgaTCTCTTTATCAATCAGCGGTTTGTCGGCTAATCCTTGCTGTGCCTGGAAGTCAGCGTGTGACG  <  NC_000913/3615035‑3614973

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