New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3033670 =	104 (0.950)	6 (0.070) +CCCACAGAAATAAATTTCACTTTA	4/150	12.3	3.3%	coding (1505/1518 nt)	lysS	lysine tRNA synthetase, constitutive
	?	NC_000913	4353155 =	358 (3.280)					intergenic (+74/+45)	ghoT/lysU	toxin of GhoTS toxin‑antitoxin pair; membrane‑lytic protein; stimulator of persister cell formation/lysine tRNA synthetase, inducible
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below cutoff threshold.

CCAACAGCCATACCATCCGCGACGTTATTCTGTTCCCGGCGATGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3033715‑3033670 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGA  >  NC_000913/4353155‑4353230                                               ||||||||||||||||||||||||                                                                             ctAACAGCCATACTATTCGCGACGTTATTCTCTTCCCGGCGATGCGCCCACAGAAATAAATTTCACTTTAATGAACGAAGCAGTCGCCTGACTGCTTCGT                                                <  1:2305451/98‑1 ctAACAGCCATACTATTCGCGACGTTATTCTCTTCCCGGCGATGCGCCCACAGAAATAAATTTCACTTTAATGAACGAAGCAGTCGCCTGACTGCTTCGT                                                <  2:684530/98‑1  tAACAGCCATACTATTCGCGACGTTATTCTCTTCCCGGCGATGCGCCCACAGAAATAAATTTCACTTTAATGAACGAAGCAGTCGCCTGACTGCTTCGTT                                               <  2:1014200/99‑1   AACAGCCATACTATTCGCGACGTTATTCTCTTCCCGGCGATGCGCCCACAGAAATAAATTTCACTTTAATGAACGAAGCAGTCGCCTGACTGCTTCGTTC                                              <  1:2004811/100‑1       GCCATACTATTCGCGACGTTATTCTCTTCCCGGCGATGCGCCCACAGAAATAAATTTCACTTTAATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTA                                          <  2:925253/100‑1            actattCGCGACGTTATTCTCTTCCCGGCGATGCGCCCACAGAAATAAATTTCACTTTAATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTG                                     <  2:345593/94‑1                                                CCACAGAAATAAATTTCACTTTAATGAACGAAGCAGTCACCTGACTGCTTCGTTCATTAAAGAGAAATATATTTCTGAGGGCGCATCGCCGGGAAGAGA  >  2:1860489/1‑99                                               ||||||||||||||||||||||||                                                                             CCAACAGCCATACCATCCGCGACGTTATTCTGTTCCCGGCGATGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3033715‑3033670 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGA  >  NC_000913/4353155‑4353230

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