New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 626190	38 (1.080)	7 (0.210)	7/260	3.1	15.5%	coding (338/816 nt)	mtfA	glucose‑regulated peptidase
	?	NC_002947	= 626200	40 (1.190)					coding (328/816 nt)	mtfA	glucose‑regulated peptidase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GCATGTCGAGCTTGTGCGCCAGCTCATGGATGACCAGGTTGTAAGCCTCCCAGCCACCACTGGCCTGCACGCCGTTCCAGGCCAGGATGATCGGCCCCTGTTGCCAGGCTTCGCCGCTGTGTTCGCCGTCCCACACATGCTCCACGCCAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/626041‑626190 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgccacCGCGATGCCAGTGGCGTGGAGCATGTGT  <  NC_002947/626200‑626173                                                                                                                                                                                    GCATGTCGAGCTTGTGCGCCAGCTCATGGATGACCAGGTTGTAAGCCTCCCAGCCACCACTGGCCTGCACGCCGTTCCAGGCCAGGATGATCGGCCCCTGTTGCCAGGCTTCGCCGCTGTGTTCGCCGTCCCACACATGCTCCACGCCAC                              >  1:312454/1‑150 GCATGTCGAGCTTGTGCGCCAGCTCATGGATGACCAGGTTGTAAGCCTCCCAGCCACCACTGGCCTGCACGCCGTTCCAGGCCAGGATGATCGGCCCCTGTTGCCAGGCTTCGCCGCTGTGTTCGCCGTCCCACACATGCTCCACGCCAC                              >  1:394535/1‑150                             GATGACCAGGTTGTAAGCCTCCCAGCCACCACTGGCCTGCACGCCGTTCCAGGCCAGGATGATCGGCCCCTGTTGCCAGGCTTCGCCGCTGTGTTCGCCGTCCCACACATGCTCCACGCCACCGCGATGCCAGTGGCGTGGAGCATGTGT  <  1:28652/150‑1                                   CAGGTTGTAAGCCTCCCAGCCACCACTGGCCTGCACGCCGTTCCAGGCCAGGATGATCGGCCCCTGTTGCCAGGCTTCGCCGCTGTGTTCGCCGTCCCACACATGCTCCACGCCACCGCGATGCCAGTGGC               <  1:200399/131‑1                                   CAGGTTGTAAGCCTCCCAGCCACCACTGGCCTGCACGCCGTTCCAGGCCAGGATGATCGGCCCCTGTTGCCAGGCTTCGCCGCTGTGTTCGCCGTCCCACACATGCTCCACGCCACCGCGATGCCAGTGGC               >  2:200399/1‑131                                                                                           TCGGCCCCTGTTGCCAGGCTTCGCCGCTGTGTTCGCCGTCCCACACATGCTCCACGCCACCGCGATGCCAGTGGCGTGGA          >  1:58901/1‑80                                                                                           TCGGCCCCTGTTGCCAGGCTTCGCCGCTGTGTTCGCCGTCCCACACATGCTCCACGCCACCGCGATGCCAGTGGCGTGGA          <  2:58901/80‑1                                                                                                         CAGGCTTCGCCGCTGTGTTCGCCGTCCCACACATGCTCCACGCCACCGCGATGCCAGTG                 >  1:212304/1‑59                                                                                                         CAGGCTTCGCCGCTGTGTTCGCCGTCCCACACATGCTCCACGCCACCGCGATGCCAGTG                 <  2:212304/59‑1                                                                                                                                                                                    GCATGTCGAGCTTGTGCGCCAGCTCATGGATGACCAGGTTGTAAGCCTCCCAGCCACCACTGGCCTGCACGCCGTTCCAGGCCAGGATGATCGGCCCCTGTTGCCAGGCTTCGCCGCTGTGTTCGCCGTCCCACACATGCTCCACGCCAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/626041‑626190 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgccacCGCGATGCCAGTGGCGTGGAGCATGTGT  <  NC_002947/626200‑626173

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