New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1570641 =	77 (1.050)	5 (0.080) +ATAACGTTGTAAA	5/268	7.2	5.7%	intergenic (‑152/+4)	gadC/gadB	glutamate:gamma‑aminobutyric acid antiporter/glutamate decarboxylase B, PLP‑dependent
	?	NC_000913	3666120 =	106 (1.450)					intergenic (‑310/+60)	gadX/gadA	acid resistance regulon transcriptional activator; autoactivator/glutamate decarboxylase A, PLP‑dependent
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GAACTGTTGCTGGAAGACTACAAAGCCTCCCTGAAATATCTCAGCGATCACCCGAAACTGCAGGGTATTGCCCAACAGAACAGCTTTAAACATACCTGATAAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/1570743‑1570641 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑AACCGAATGCCCAGCTGTTTTTTTAAAGGCTGGGCATTCGGTTTTTAC  >  NC_000913/3666120‑3666167                                                                                                        |||||||||||||                                                 GAACTGTTGCTGGAAGACTACAAAGCCTCCCTGAAATATCTCAGCGATCACCCGAAACTGCAGGGTATTGCCCAGCAGAACAGCTTTAAACACACCTGATAACATAACGTTGTAAAAACCGAATGCCCAGCTGTTTTTTTAAAGGCGGG                 <  2:330020/149‑1   ACTGTTGCTGGAAGAGTACTAAGCCTCCCTGAAATATCTCAGCGATCACCCGAAACTGCAGGGTATTGCCCAGCAGAACAGCTTTAAACACACCTGATAACATAACGTTGTAAAAACCGAATGCCCAGCTGTTTTTTTAAAGGCTGGGC               <  2:628478/149‑1      GTTGCTGGAAGACTACAAAGCCTCCCTGAAATATCTCAGCGATCACCCGAAACTGCAGGGTATTGCCCAGCAGATCAGCTTTAAACACACCTGATAACATAACGTTGTAAAAACCGAATGCCCAGCTGTTTTTTTAAAGGCTGGGCATT            <  2:785027/149‑1        TGCTGGAAGACTACAAAGCCTCCCTGAAATATCTCAGCGATCACCCGAAACTGCAGGGTATTGCCCAGCAGAACAGCTTTAAACACACCTGATAACATAACGTTGTAAAAACCGAATGCCCAGCTGTTTTTTTAAAGGCTGGGCATTCG          <  2:1075890/149‑1                GACTACAAAGCCTCCCTGAAATATCTCCGCGATCACCCGAAACTGCAGGGTATTGCCCAGCAGAACAGCTTTAAACACACCTGATAACATAACGTTGTAAAAACCGAATGCCCAGCTGTTTTTTTAAAGGCTGGGCATTCGGTTTTTAC  <  1:184090/149‑1                                                                                                        |||||||||||||                                                 GAACTGTTGCTGGAAGACTACAAAGCCTCCCTGAAATATCTCAGCGATCACCCGAAACTGCAGGGTATTGCCCAACAGAACAGCTTTAAACATACCTGATAAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/1570743‑1570641 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑AACCGAATGCCCAGCTGTTTTTTTAAAGGCTGGGCATTCGGTTTTTAC  >  NC_000913/3666120‑3666167

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