New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3783137 =	24 (0.800)	4 (0.140)	4/208	2.9	18.0%	coding (525/1020 nt)	gpsA	glycerol‑3‑phosphate dehydrogenase (NAD+)
	?	NC_000913	3783176 =	14 (0.500)					coding (486/1020 nt)	gpsA	glycerol‑3‑phosphate dehydrogenase (NAD+)
Rejected: Frequency below/above cutoff threshold.

GACCTTTGCCGATGATCTCCAGCAGCTGCTGCACTGCGGCAAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3783179‑3783137 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cggcaaaGGTCTGATCGGTCGAGGCCAGCGAAATAGCTGTCGGTAAACCTGCCGCCAGTTCTTTCGCAAACGTTGGGCCAGAGATAACCGCCAGCGGAATTTGA  >  NC_000913/3783176‑3783272                                                                                                                                              GACCTTTGCCGATGATCTCCAGCAGCTGCTGCACTGCGGCAAAGGTCTGATCGGTCGAGGCCAGCGAAATAGCTGTCG                                                                >  5:4689/1‑78 GACCTTTGCCGATGATCTCCAGCAGCTGCTGCACTGCGGCAAAGGTCTGATCGGTCGAGGCCAGCGAAATAGCTGTCG                                                                <  6:4689/78‑1     TTTGCCGATGATCTCCAGCAGCTGCTGCACTGCGGCAAAGGTCTGATCGGTCGAGGCCAGCGAAATAGCTGTCGGTAAACCTGCCGCCAGTTCTTTCGCAAACGTTGGGCCAGAGATAACCGCCAGCGGAATTTGA  <  3:81910/136‑1     TTTGCCGATGATCTCCAGCAGCTGCTGCACTGCGGCAAAGGTCTGATCGGTCGAGGCCAGCGAAATAGCTGTCGGTAAACCTGCCGCCAGTTCTTTCGCAAACGTTGGGCCAGAGATAACCGCCAGCGGAATTTGA  >  4:81910/1‑136                                        CAAAGGTCTGATCGGTCGAGGCCAGCGAAATAGCTGTCGGTAAACCTGCCGCCAGTTCTTTCGCAAACGTTGGG                             >  4:64056/1‑74                                                                                                                                              GACCTTTGCCGATGATCTCCAGCAGCTGCTGCACTGCGGCAAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3783179‑3783137 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cggcaaaGGTCTGATCGGTCGAGGCCAGCGAAATAGCTGTCGGTAAACCTGCCGCCAGTTCTTTCGCAAACGTTGGGCCAGAGATAACCGCCAGCGGAATTTGA  >  NC_000913/3783176‑3783272

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.