New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_bme_pgi	= 1939098	NA (NA)	4 (0.100)	4/270	5.9	NA	intergenic (+7/+124)	pykA/lpxM	pyruvate kinase II/myristoyl‑acyl carrier protein (ACP)‑dependent acyltransferase
	?	NC_000913_3_bme_pgi	= 1939128	NA (NA)					noncoding (30/30 nt)	REP137 (repetitive extragenic palindromic) element; contains 1 REP sequences	REP137 (repetitive extragenic palindromic) element; contains 1 REP sequences
Rejected: Coverage evenness skew score above cutoff.

GTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCACCCAGGGCGACGTGATGAGTACCGTGGGTTCTACTAATACCACGCGTATTTTAACGGTAGAGTAAGTACGTTgccggatg‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913_3_bme_pgi/1938951‑1939098 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GCCGGATGTGGCGTTTTCGCCGCATCCGGCAACGTACTTACTCTACCGTTAAAATACGCGTGGTATTAGTAGAACCCACGGTACTCATCACGTCGCCCTGGGTGACAATCACCAGGTCACCAGACATCAAGTAACCTTTATCGCGCAGCA  <  NC_000913_3_bme_pgi/1939128‑1938979                                                                                                                                                                                                                                                                                                            GTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCACCCAGGGCGACGTGATGAGTACCGTGGGTTCTACTAATACCACGCGTATTTTAACGGTAGAGTAAGTACGTTG                                                                                                                                                       >  1:54979/1‑149                                                              ACCTGGTGATTGTCACCCAGGGCGACGTGATGAGTACCGTGGGTTCTACTAATACCACGCGTATTTTAACGGTAGAGTAAGTACGTTGCCGGATGTGGCGTTTTCGCCGCATCCGGCAACGTACTTACT                                                                                                              >  2:563728/1‑129                                                              ACCTGGTAATTGTCACCCAGGGCGACGTGATGAGTACCGTGGGTTCTACTAATACCACGCGTATTTTAACGGTAGAGTAAGTACGTTGCCGGATGTGGCGTTTTCGCCGCATCCGGCAACGTACTTACT                                                                                                              <  1:563728/129‑1                                                               CCTGGTGATTGTCACCCAGGGCGACGTGATGAGTACCGTGGGTTCTACTAATACCACGCGTATTTTAACGGTAGAGTAAGTACGTTGCC                                                                                                                                                     <  1:515366/89‑1                                                               CCTGGTGATTGTCACCCAGGGCGACGTGATGAGTACCGTGGGTTCTACTAATACCACGCGTATTTTAACGGTAGAGTAAGTACGTTGCC                                                                                                                                                     >  2:515366/1‑89                                                                       TTGTCACCCAGGGCGACGTGATGAGTACCGTGGGTTCTACTAATACCACGCGTATTTTAACGGTAGAGTAAGTACGTTGCCGGATGTGGCGTTTTCGCCGCATC                                                                                                                              >  1:112493/1‑104                                                                       TTGTCACCCAGGGCGACGTGATGAGTACCGTGGGTTCTACTAATACCACGCGTATTTTAACGGTAGAGTAAGTACGTTGCCGGATGTGGCGTTTTCGCCGCATC                                                                                                                              <  2:112493/104‑1                                                                                                                                                      CCGGATGTGGCGGTTTCGCCGCATCCGGCAACGTACTTACTCTACCGTTAAAATACGCGTGGTATTAGTAGAACCCACGGTACTCATCACGTCGCCCTGGGTGACAATCACCAGGTCACCAGACATCAAGTAACCTTTATCGCGCAGCA  <  2:2706/149‑1                                                                                                                                                                                                                                                                                                            GTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCACCCAGGGCGACGTGATGAGTACCGTGGGTTCTACTAATACCACGCGTATTTTAACGGTAGAGTAAGTACGTTgccggatg‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913_3_bme_pgi/1938951‑1939098 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GCCGGATGTGGCGTTTTCGCCGCATCCGGCAACGTACTTACTCTACCGTTAAAATACGCGTGGTATTAGTAGAACCCACGGTACTCATCACGTCGCCCTGGGTGACAATCACCAGGTCACCAGACATCAAGTAACCTTTATCGCGCAGCA  <  NC_000913_3_bme_pgi/1939128‑1938979

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