New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	CP009273	= 361823	268 (1.000)	6 (0.020) +TTAATTAA	5/262	NT	2.2%	coding (831/834 nt)	mhpR	mhp operon transcriptional activator
	?	CP009273	4033991 =	299 (1.120)					noncoding (115/120 nt)	rrfA	5S ribosomal RNA of rrnA operon
Rejected: Frequency below/above cutoff threshold.

GACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCGCCATTCGATGGTGTCAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  CP009273/361699‑361823 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑AGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGAC  >  CP009273/4033991‑4034079                                                                                                                              ||||||||                                                                                          gaaACCGGCATAATCTGCGACATCGTATAACTTTACTGGTTTCACATTCACCACCATGAATTGACTCTCTTCCGTGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAATTAATTAAATGCATCAAATAAAACG                                                                          >  2:3320960/4‑150      cgtCATACTCTGCTAAATAGTATAACGTTACTGGTTTCAAATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCAGGCCATACCGCTAAAGGTTTTGAACCATTCGATGGTGTCAATTAATTAAAGGCATCAAATAAAACGAAAGG                                                                     <  1:3273954/147‑1       ggaATACTCTGCGACTTCTTATAACTTTTCTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCTATGGTGTAAATTAATTAAAGTCATCAAATAAAACGAAAGGC                                                                    <  2:3351409/147‑1                TGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAATTAATTAAAGGGATCAAATAAAACGAAGGGCTCAGTCGGA                                                           >  1:3658773/1‑150                                              ATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAATTAATTAAAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTGTCGTGTTATCTGTTGGTT                             >  1:4302584/1‑150                                                                         cgcGCGCCACCACGCCATCCCGCGCAAGGTTTGGCACCAGTCGATGGTGACAATTAATTAAAGCCATCAAATCACACGAAAGGCTCAGTCGCAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGAC  <  1:881381/147‑1                                                                                                                              ||||||||                                                                                          GACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCGCCATTCGATGGTGTCAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  CP009273/361699‑361823 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑AGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGAC  >  CP009273/4033991‑4034079

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