New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1570641 =	107 (0.570)	6 (0.040) +ATAACGTTGTAAA	5/212	NT	9.5%	intergenic (‑152/+4)	gadC/gadB	glutamate:gamma‑aminobutyric acid antiporter/glutamate decarboxylase B, PLP‑dependent
	?	NC_000913	3666120 =	22 (0.120)					intergenic (‑310/+60)	gadX/gadA	acid resistance regulon transcriptional activator; autoactivator/glutamate decarboxylase A, PLP‑dependent

CAAAGCCTCCCTGAAATATCTCAGCGATCACCCGAAACTGCAGGGTATTGCCCAACAGAACAGCTTTAAACATACCTGATAAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/1570723‑1570641 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑AACCGAATGCCCAGCTGTTTTTTTAAAGGCTGG  >  NC_000913/3666120‑3666152                                                                                    |||||||||||||                                  CAAAGCCTCCCTGAAATATCTCAGCGATCACCCGAAACTGCAGGGTATTGCCCAGCAGAACAGCTTTAAACACACCTGATAACATAACGTTGTAAAAACCGAATGCCCAGCTGTTTTTTTAAAGGCTGG  >  3:705512/1‑129 CAAAGCCTCCCTGAAATATCTCAGCGATCACCCGAAACTGCAGGGTATTGCCCAGCAGAACAGCTTTAAACACACCTGATAACATAACGTTGTAAAAACCGAATGCCCAGCTGTTTTTTTAAAGGCTGG  <  4:705512/129‑1        TCCCTGAAATATCTCAGCGATCACCCGAAACTGCAGGGTATTGCCCAGCAGAACAGCTTTAAACACACCTGATAACATAACGTTGTAAAAACCGAATGCCCAGCTGTTTTTTTAAAGGCTGG  <  1:659723/122‑1        TCCCTGAAATATCTCAGCGATCACCCGAAACTGCAGGGTATTGCCCAGCAGAACAGCTTTAAACACACCTGATAACATAACGTTGTAAAAACCGAATGCCCAGCTGTTTTTTTAAAGGCTGG  >  2:659723/1‑122                   TCTCAGCGATCACCCGAAACTGCAGGGTATTGCCCAGCAGAACAGCTTTAAACACACCTGATAACATAACGTTGTAAAAACCGAATGCCCAGCTGTTTTTTTAAAGGCT    <  3:456038/109‑1                   TCTCAGCGATCACCCGAAACTGCAGGGTATTGCCCAGCAGAACAGCTTTAAACACACCTGATAACATAACGTTGTAAAAACCGAATGCCCAGCTGTTTTTTTAAAGGCT    >  4:456038/1‑109                                                                                    |||||||||||||                                  CAAAGCCTCCCTGAAATATCTCAGCGATCACCCGAAACTGCAGGGTATTGCCCAACAGAACAGCTTTAAACATACCTGATAAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/1570723‑1570641 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑AACCGAATGCCCAGCTGTTTTTTTAAAGGCTGG  >  NC_000913/3666120‑3666152

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