New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3783607 =	120 (1.340)	6 (0.070)	6/250	9.4	5.1%	coding (55/1020 nt)	gpsA	glycerol‑3‑phosphate dehydrogenase (NAD+)
	?	NC_000913	3783630 =	112 (1.330)					coding (32/1020 nt)	gpsA	glycerol‑3‑phosphate dehydrogenase (NAD+)
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CCGCAACTGAACCTTGCGCCGGTTAACTTCGATGCGCTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCGTACGGCACCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3783759‑3783607 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cggcaccgATCACAGTCATTGAAGCATTACGTTGGTTCATCAGGCATCCTGATGTTCTTCAGTACCTTCGCCAGCCTGCTGCTGCAAATAGTTCATGAACAGCGCATCGAAGTTAACCGGCGCAAGGTTCAGTTGCGGGAATGTACCGCG  >  NC_000913/3783630‑3783771                                                                                                                                                                                                                                                                                                         CCGCAACTGAACCTTGCGCCGGTTAACTTCGATGCGCTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCGTACGGC                                                                                                                                                    >  5:338326/1‑149  CGCAACTGAACCTTGCGCCGGTTAACTTCGATGCGCTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCGTACGGCA                                                                                                                                                   >  4:156545/1‑149      ACTGAACCTTGCGCCGGTTAACTTCGATGCGCTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCGTACGGCACCG                                                                                                                                                >  7:157520/1‑148                                                                                                                          GACTGTGATCGGTGCCGGCTCGTACGGCACCGATCACAGTCATTGAAGCATTACGTTGGTTCATCAGGCATCCTGATGTTCTTCAGTACCTTCGCCAGCCTGCTGCTGCAAA                                                                <  1:179668/112‑1                                                                                                                          GACTGTGATCGGTGCCGGCTCGTACGGCACCGATCACAGTCATTGAAGCATTACGTTGGTTCATCAGGCATCCTGATGTTCTTCAGTACCTTCGCCAGCCTGCTGCTGCAAA                                                                >  2:179668/1‑112                                                                                                                           ACTGTGATCGGTGCCGGCTCGTACGGCACCGATCACAGTCATTGAAGCATTACGTTGGTTCATCAGGCATCCTGATGTTCTTCAGTACCTTCGCCAGCCTGCTGCTGC                                                                   <  7:166346/108‑1                                                                                                                           ACTGTGATCGGTGCCGGCTCGTACGGCACCGATCACAGTCATTGAAGCATTACGTTGGTTCATCAGGCATCCTGATGTTCTTCAGTACCTTCGCCAGCCTGCTGCTGC                                                                   >  8:166346/1‑108                                                                                                                                  TCGGTGCCGGCTCGTACGGCACCGATCACAGTCATTGAAGCATTACGTTGGTTCATCAGGCATCCTGATGTTCTTCAGTACCTTCGCCAGCCTGCTGCTGCAAATAGTTCATGAACAGCGC                                               >  5:370683/1‑121                                                                                                                                  TCGGTGCCGGCTCGTACGGCACCGATCACAGTCATTGAAGCATTACGTTGGTTCATCAGGCATCCTGATGTTCTTCAGTACCTTCGCCAGCCTGCTGCTGCAAATAGTTCATGAACAGCGC                                               <  6:370683/121‑1                                                                                                                                                   GGCACCGATCACAGTCATTGAAGCATTACGTTGGTTCATCAGGCATCCTGATGTTCTTCAGTACCTTCGCCAGCCTGCTGCTGCAAATAGTTCATGAACAGCGCATCGAAGTTAACCGGCGCAAGGTTCAGTTGCGGGAATGTACCGCG  >  6:220912/1‑149                                                                                                                                                                                                                                                                                                         CCGCAACTGAACCTTGCGCCGGTTAACTTCGATGCGCTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCGTACGGCACCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3783759‑3783607 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cggcaccgATCACAGTCATTGAAGCATTACGTTGGTTCATCAGGCATCCTGATGTTCTTCAGTACCTTCGCCAGCCTGCTGCTGCAAATAGTTCATGAACAGCGCATCGAAGTTAACCGGCGCAAGGTTCAGTTGCGGGAATGTACCGCG  >  NC_000913/3783630‑3783771

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