New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1570641 =	109 (1.180)	3 (0.040) +ATAACGTTGTAAA	3/234	NT	2.7%	intergenic (‑152/+4)	gadC/gadB	glutamate:gamma‑aminobutyric acid antiporter/glutamate decarboxylase B, PLP‑dependent
	?	NC_000913	3666120 =	134 (1.450)					intergenic (‑310/+60)	gadX/gadA	acid resistance regulon transcriptional activator; autoactivator/glutamate decarboxylase A, PLP‑dependent
Rejected: Frequency below/above cutoff threshold.

TGGAAGACTACAAAGCCTCCCTGAAATATCTCAGCGATCACCCGAAACTGCAGGGTATTGCCCAACAGAACAGCTTTAAACATACCTGATAAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/1570733‑1570641 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑AACCGAATGCCCAGCTGTTTTTTTAAAGGCTGGGCATTCGGTTTTTACAACGTTATGTTATCAGGTGTGTTTAAAGCTGTTCTGCTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTT  >  NC_000913/3666120‑3666245                                                                                              |||||||||||||                                                                                                                               TGGAAGACTACAAAGCCTCCCTGAAATATCTCAGCGATCACCCGAAACTGCAGGGTATTGCCCAGCAGAACAGCTTTAAACACACCTGATAACATAACGTTGTAAAAACCGAATGCCCAGCTGTTTTTTTAAAGGCTGGGCATTCGG                                                                                       <  7:21814/147‑1 TGGAAGACTACAAAGCCTCCCTGAAATATCTCAGCGATCACCCGAAACTGCAGGGTATTGCCCAGCAGAACAGCTTTAAACACACCTGATAACATAACGTTGTAAAAACCGAATGCCCAGCTGTTTTTTTAAAGGCTGGGCATTCGG                                                                                       >  8:21814/1‑147                                                                                    ACCTGATAACATAACGTTGTAAAAACCGAATGCCCAGCCTTTAAAAAAACAGCTGGGCATTCGGTTTTTACAACGTTATGTTATCAGGTGTGTTTAAAGCTGTTCTGCTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTT  >  8:49631/1‑149                                                                                              |||||||||||||                                                                                                                               TGGAAGACTACAAAGCCTCCCTGAAATATCTCAGCGATCACCCGAAACTGCAGGGTATTGCCCAACAGAACAGCTTTAAACATACCTGATAAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/1570733‑1570641 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑AACCGAATGCCCAGCTGTTTTTTTAAAGGCTGGGCATTCGGTTTTTACAACGTTATGTTATCAGGTGTGTTTAAAGCTGTTCTGCTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTT  >  NC_000913/3666120‑3666245

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