New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4591667	71 (0.510)	4 (0.030)	4/214	NT	5.5%	coding (11/1656 nt)	tsr	methyl‑accepting chemotaxis protein I, serine sensor receptor
	?	NC_000913	= 4591703	71 (0.540)					coding (47/1656 nt)	tsr	methyl‑accepting chemotaxis protein I, serine sensor receptor

TTTCCAGGCCGAAAATCTTGCATCGGTCCACAGGAAAGAGAAACCATGTTAAAACG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4591612‑4591667 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aaaacgGCCAAAACCAGCAGTAAGCTGGTCACAATTTTGATACGTTTTAACATGGTTTCTCTTTCCTGTGGACCGATGCAAGATTTTCGGCCTGGAAAGGAAAAACTTTATGAATATCCCGCGGAGATTACGCCGAGTGAATTTTATTCA  <  NC_000913/4591703‑4591560                                                                                                                                                                                                          TTTCCAGGCCGAAAATCTTGCATCGGTCCACAGGAAAGAGAAACCATGTTAAAACGGCCAAAACCAGCAGTAAGCTGGTCACAATTTTGATACGT                                                                                                           <  1:290623/95‑1 TTTCCAGGCCGAAAATCTTGCATCGGTCCACAGGAAAGAGAAACCATGTTAAAACGGCCAAAACCAGCAGTAAGCTGGTCACAATTTTGATACGT                                                                                                           >  2:290623/1‑95                      ATCGGTCCACAGGAAAGAGAAACCATGTTAAAACGGCCAAAACCAGCAGTAAGCTGGTCACAATTTTGATACGTT                                                                                                          >  3:677526/1‑75                      ATCGGTCCACAGGAAAGAGAAACCATGTTAAAACGGCCAAAACCAGCAGTAAGCTGGTCACAATTTTGATACGTT                                                                                                          <  4:677526/75‑1                                                    AAACGGCCAAAACCAGCAGTAAGCTGGTCACAATTTTGATACGTTTTAACATGGTTTCTCTTTCCTGTGGACCGATGCAAGATTTTCGGCCTGGAAAGGAAAAACTTTATGAATATCCCGCGGAGATTACGCCGAGTGAATTTTATTCA  <  1:123028/149‑1                                                                                                                                                                                                          TTTCCAGGCCGAAAATCTTGCATCGGTCCACAGGAAAGAGAAACCATGTTAAAACG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4591612‑4591667 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aaaacgGCCAAAACCAGCAGTAAGCTGGTCACAATTTTGATACGTTTTAACATGGTTTCTCTTTCCTGTGGACCGATGCAAGATTTTCGGCCTGGAAAGGAAAAACTTTATGAATATCCCGCGGAGATTACGCCGAGTGAATTTTATTCA  <  NC_000913/4591703‑4591560

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