New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	4132612 =	20 (1.200)	5 (0.300)	4/284	1.6	24.5%	coding (245/687 nt)	PP_3634	GNAT family acetyltransferase
	?	NC_002947	4132646 =	11 (0.660)					coding (279/687 nt)	PP_3634	GNAT family acetyltransferase

CAGTACGGCGTCGGCATGGAAGTCCACCAGTACCTGCAAGCGCTGGACGGCACCCGCGGTCTGGCTGTGGCGCTTACCCTGGCGCTGGATGCCGAAGCACCGGACCAGTTGCTGGGCTTTGCGCTGTCGCTGCCAGCCCAGGACGATGGGCAGGCCTGTGCATTGGCCTTTCTGGCCGTGTGCGCCAGCCATCGCCGCCAGGGCATCGCCCGTGCCCTGCTGGGCGACCTGCGCTCACGCCATGCCTGCGTGGAGCTGAACGCCTTTGC  >  NC_002947/4132503‑4132771                                                                                                                                                |                                                                                                                              cAGTACGGCGTCGGCATGGAAGTCCACCAGTACCTGCAAGCGCTGGACGGCACCCGCGGTCTGGCTGTGGCGCTTACCCTGGCGCTGGATGCCGAAGCACCGGACCAGTTGCTGGGCTTTGCGCTGTCGCTGCCAGCCCAGGACGATggg                                                                                                                         >  2:371023/1‑150 (MQ=255)                  ggAAGTCCACCAGTACCTGCAAGCGCTGGACGGCACCCGCGGTCTGGCTGTGGCGCTTACCCTGGCGCTGGATGCCGAAGCACCGGACCAGTTGCTGGGCTTTGCGCTGTCGCTGCCAGCCCAGGACGATGGGCAGGCCTGTGCATTGGc                                                                                                        >  2:255565/1‑150 (MQ=255)                       tCCACCAGTACCTGCAAGCGCTGGACGGCACCCGCGGTCTGGCTGTGGCGCTTACCCTGGCGCTGGATGCCGAAGCACCGGACCAGTTGCTGGGCTTTGCGCTGTCGCTGCCAGCCCAGGACGATGGGCAGGCCTGTGCATTGGCCTTTc                                                                                                   >  2:141639/1‑150 (MQ=255)                            gcccagcagggcacgggcgatgccctggcggcgatggctggcgcacacggccagaaaggccaatgcacaggcctgcccatcgTGCTGGGCTTTGCGCTGTCGCTGCCAGCCCAGGACGATggg                                                                                                                         >  2:152281‑M1/83‑123 (MQ=255)                            gcccagcagggcacgggcgatgccctggcggcgatggctggcgcacacggccagaaaggccaatgcacaggcctgcccatcgTGCTGGGCTTTGCGCTGTCGCTGCCAGCCCAGGACGATggg                                                                                                                         <  1:152281‑M1/41‑1 (MQ=255)                                  cTGCAAGCGCTGGACGGCACCCGCGGTCTGGCTGTGGCGCTTACCCTGGCGCTGGATGCCGAAGCACCGGACCAGTTGCTGGGCTTTGCGCTGTCGCTGCCAGCCCAGGACGATGGGCAGGCCTGTGCATTGGCCTTTCTGGCCGTGTgc                                                                                        >  2:197403/1‑150 (MQ=255)                                               aCGGCACCCGCGGTCTGGCTGTGGCGCTTACCCTGGCGCTGGATGCCGAAGAACCGGACCAGTTGCTGGGCTTTGCGCTGTCGCTGCCAGCCCAGGACGATGGGCAGGCCTGTGCATTGGCCTTTCTGGCCGTGTGCGCCAGCCATcgcc                                                                           <  2:185168/150‑1 (MQ=255)                                                   cACCCGCGGTCTGGCTGTGGCGCTTACCCTGGCGCTGGATGCCGAAGCACCGGACCAGTTGCTGGGCTTTGCGCTGTCGCTGCCAGCCCAGGACGATGGGCAGGCCTGTGCATTGGCCTTTCTGGCCGTGTGCGCCAGCCATCGCCGCCa                                                                       <  2:303115/150‑1 (MQ=255)                                                        gcgGTCTGGCTGTGGCGCTTACCCTGGCGCTGGATGCCGAAGCACCGGACCAGTTGCTGGGCTTTGCGCTGTCGCTGCCAGCCCAGGACGATGGGCAGGCCTGTGCATTGGCCTTTCTGGCCGTGTGCGCCAGCCATCGCCGCCAGGGCa                                                                  <  1:198116/150‑1 (MQ=255)                                                          ggTCTGGCTGTGGCGCTTACCCTGGCGCTGGATGCCGAAGCACCGGACCAGTTGCTGGGCTTTGCGCTGTCGCTGCCAGCCCAGGACGATGGGCAGGCCTGTGCATTGGCCTTTCTGGCCGTGTGCGCCAGCCATCGCCGCCAGGGCATc                                                                <  1:78851/150‑1 (MQ=255)                                                                                         aTGCCGAAGCACCGGACCAGTTGCTGGGCTTTGCGCTGTCGCTGCCAGCCCAGGACGATGGGCAGGCCTGTGCATTGGCCTTTCTGGCCGTGTGCGCCAGCCATCGCCGCCAGGGCATCGCCCGTGCCCTGCTGGGCGACCTGCGCTCAc                                 <  1:371023/150‑1 (MQ=255)                                                                                                       cccatcgtcctgggctggcagcgacagcgcaaagcccagcaCGATGGGCAGGCCTGTGCATTGGCCTTTCTGGCCGTGTGCGCCAGCCATCGCCGCCAGGGCATCGCCCGTGCCCTGCTGGGc                                              >  2:152281‑M2/42‑123 (MQ=255)                                                                                                       cccatcgtcctgggctggcagcgacagcgcaaagcccagcaCGATGGGCAGGCCTGTGCATTGGCCTTTCTGGCCGTGTGCGCCAGCCATCGCCGCCAGGGCATCGCCCGTGCCCTGCTGGGc                                              <  1:152281‑M2/82‑1 (MQ=255)                                                                                                                  gggctggcagcgacagcgcaaagcccaggaCGATGGGCAGGCCTGTGCATTGGCCTTTCTGGCCGTGTGCGCCAGCCATCGCCGCCAGGGCATCGCCCGTGCCCTGCTGGGCGACCTGCGCTCACGCCATGCCTGCGTGGAGCTGAACGc        >  1:26484‑M2/31‑150 (MQ=255)                                                                                                                   ggCTTTGCGCTGTCGCTGCCAGCCCAGGACGATGGGCAGGCCTGTGCATTGGCCTTTCTGGCCGTGTGCGCCAGCCATCGCCGCCAGGGCATCGCCCGTGCCCTGCTGGGCGACCTGCGCTCACGCCATGCCTGCGTGGAGCTGAACGcc       <  1:197403/150‑1 (MQ=255)                                                                                                                      attGCGCTGTCGCTGCCAGCCCAGGACGATGGGCAGGCCTGTGCATTGGCCTTTCTGGCCGTGTGCGCCAGCCATCGCCGCCAGGGCATCGCCCGTGCCCTGCTGGGCGACCTGCGCTCACGCCATGCCTGCGTGGAGCTGAACGCCttt    <  2:197903/149‑1 (MQ=255)                                                                                                                        gcagcgacagcgcaaagcccaggaCGATGGGCAGGCCTGTGCATTGGCCTTTCTGGCCGTGTGCGCCAGCCATCGCCGCCAGGGCATCGCCCGTGCCCTGCTGGGCGACCTGCGCTCACGCCATGCCTGCGTGGAGCTGAACGCCTTTGc  >  2:229785‑M2/25‑150 (MQ=255)                                                                                                                        gcagcgacagcgcaaagcccaggaCGATGGGCAGGCCTGTGCATTGGCCTTTCTGGCCGTGTGCGCCAGCCATCGCCGCCAGGGCATCGCCCGTGCCCTGCTGGGCGACCTGCGCTCACGCCATGCCTGCGTGGAGCTGAACGCCTTTGc  >  2:229800‑M2/25‑150 (MQ=255)                                                                                                                                                |                                                                                                                              CAGTACGGCGTCGGCATGGAAGTCCACCAGTACCTGCAAGCGCTGGACGGCACCCGCGGTCTGGCTGTGGCGCTTACCCTGGCGCTGGATGCCGAAGCACCGGACCAGTTGCTGGGCTTTGCGCTGTCGCTGCCAGCCCAGGACGATGGGCAGGCCTGTGCATTGGCCTTTCTGGCCGTGTGCGCCAGCCATCGCCGCCAGGGCATCGCCCGTGCCCTGCTGGGCGACCTGCGCTCACGCCATGCCTGCGTGGAGCTGAACGCCTTTGC  >  NC_002947/4132503‑4132771

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