New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_pae_pgi	2792954 =	20 (1.230)	3 (0.200)	3/238	1.9	14.3%	coding (220/1281 nt)	gabT	4‑aminobutyrate aminotransferase, PLP‑dependent
	?	NC_000913_3_pae_pgi	2793007 =	17 (1.110)					coding (273/1281 nt)	gabT	4‑aminobutyrate aminotransferase, PLP‑dependent
Rejected: Frequency below/above cutoff threshold.

TAATCTCGCACAGCTCCAGATACGGCTCGTAAGCCAGCACCTGGAAGCAGGTGTGCGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913_3_pae_pgi/2793011‑2792954 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tgtgcgaGATTATGAATCAGAAGGTGCCGGGCGATTTCGCCAAGAAAACGCTGCTGGTTACGACCGGTTCCGAAGCGGTGGAAAACGCGGTAAAAATCGCCCGCGCCGCCACCAAACGTAGCGGCACCATCGCTTTTAGCGGCGCGTATCACGG  >  NC_000913_3_pae_pgi/2793007‑2793153                                                                                                                                                                                                               TAATCTCGCACAGCTCCAGATACGGCTCGTAAGCCAGCACCTGGAAGCAGGTGTGCGAGATTATGAATCAGAAGGTGCCGGGCGATTTCGCCAAGAAAACGCTGCTGGTTACGACCGGTTCCGAAGCGGTGGAAAACGCGGTAAAAATCG                                                         >  1:182556/1‑150              CTCCAGATACGGCTCGTAAGCCAGCACCTGGAAGCAGGTGTGCGAGATTATGAATCAGAAGGTGCCGGGCGATTTCGCCAAGAAAA                                                                                                            <  1:128148/86‑1              CTCCAGATACGGCTCGTAAGCCAGCACCTGGAAGCAGGTGTGCGAGATTATGAATCAGAAGGTGCCGGGCGATTTCGCCAAGAAAA                                                                                                            >  2:128148/1‑86                                                        CGAGATTATGAATCAGAAGGTGCCGGGCGATTTCGCCAAGAAAACGCTGCTGGTTACGACCGGTTCCGAAGCGGTGGAAAACGCGGTAAAAATCGCCCGCGCCGCCACCAAACGTAGCGGCACCATCGCTTTTAGCGGCGCGTATCACGG  >  4:100414/1‑150                                                         GAGATTATGAATCAGAAGGTGCCGGGCGATTTCGCCAAGAAAACGCTGCTGGTTACGACCGGTTCCGAAGCGGTG                                                                            <  1:13121/75‑1                                                         GAGATTATGAATCAGAAGGTGCCGGGCGATTTCGCCAAGAAAACGCTGCTGGTTACGACCGGTTCCGAAGCGGTG                                                                            >  2:13121/1‑75                                                                                                                                                                                                               TAATCTCGCACAGCTCCAGATACGGCTCGTAAGCCAGCACCTGGAAGCAGGTGTGCGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913_3_pae_pgi/2793011‑2792954 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tgtgcgaGATTATGAATCAGAAGGTGCCGGGCGATTTCGCCAAGAAAACGCTGCTGGTTACGACCGGTTCCGAAGCGGTGGAAAACGCGGTAAAAATCGCCCGCGCCGCCACCAAACGTAGCGGCACCATCGCTTTTAGCGGCGCGTATCACGG  >  NC_000913_3_pae_pgi/2793007‑2793153

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