New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1570598 =	130 (0.750)	7 (0.050) +CCGTTACCGTTAAAC	5/248	NT	4.2%	intergenic (‑109/+47)	gadC/gadB	glutamate:gamma‑aminobutyric acid antiporter/glutamate decarboxylase B, PLP‑dependent
	?	NC_000913	3666176 =	226 (1.300)					intergenic (‑366/+4)	gadX/gadA	acid resistance regulon transcriptional activator; autoactivator/glutamate decarboxylase A, PLP‑dependent
Rejected: Frequency below/above cutoff threshold.

TGATAACGTTTAACGGTAACGGTGTCCCGAAACGAACCCGTTTCGGGACA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/1570647‑1570598 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GTTATCAGGTGTGTTTAAAGCTGTTCTGCTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGT  >  NC_000913/3666176‑3666286                                                   |||||||||||||||                                                                                                                TGATAACGTTTAACGGTAACGGTGTCCCGAAACGAACCCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTA                            >  2:1046966/1‑150         TTTAACGGTAACGGTGTCCCGAAACGAACCCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCA                    >  2:1695388/1‑150         TTTAACGGTAACGGTGTCCCGAAACGAACCCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGGTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCA                    >  2:1695654/1‑150               GGTAACGGTGTCCCGAAACGAACCCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGACTTTGTAGTCTTCCAGCAACA              >  2:859256/1‑150                          CCCGAAACGAACCCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAG   >  2:1260739/1‑150                          CCCGAAACGAACCCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGAGCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCCAAG   >  1:549673/1‑150                           CCGAAACGAACCCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGT  >  2:339200/1‑150                                                   |||||||||||||||                                                                                                                TGATAACGTTTAACGGTAACGGTGTCCCGAAACGAACCCGTTTCGGGACA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/1570647‑1570598 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GTTATCAGGTGTGTTTAAAGCTGTTCTGCTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGT  >  NC_000913/3666176‑3666286

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