New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4396901 =	198 (0.990)	3 (0.020)	3/274	NT	1.5%	coding (837/1338 nt)	amiB	N‑acetylmuramoyl‑l‑alanine amidase II
	?	NC_000913	4396939 =	197 (1.030)					coding (875/1338 nt)	amiB	N‑acetylmuramoyl‑l‑alanine amidase II
Rejected: Frequency below/above cutoff threshold.

TACCCGCTGGGAATGACCGAACTGTAAATCCAGCACCGCCTGGCTTAAATAGGGGTCAGACTGACTGTTCGCCAGCACATCACCCGCCCCACCCAGCAGCTCCGACTGTTTCTCGTGCTGCTCCAGCCAGCTGGCCATTTCACTGTTGGCGCGACGGTTAGAGAGCACCCATACGGAAGCGCCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4397084‑4396901 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcgccaACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGTGATGTGCTGGCGAACAGTCAGTCTGACCCCTATTTAAGCCAGGCGGTGCTGGATTTACAGTTCGGTCATTCCCAGCGGGTAGG  >  NC_000913/4396939‑4397086                                                                                                                                                                                                                                                                                                                                              TACCCGCTGGGAATGACCGAACTGTAAATCCAGCACCGCCTGGCTTAAATAGGGGTCAGACTGACTGTTCGCCAGCACATCACCCGCCCCACCCAGCAGCTCCGACTGTTTCTCGTGCTGCTCCAGACAGCTGGCCATTTCACTGTTGGCGCGACGGTTAGAGAGCACCCATACGGAAGCGCCA                                                                                                                                                      >  4:315943/1‑184                               CAGCACCGCCTGGCTTAAATAGGGGTCAGACTGACTGTTCGCCAGCACATCACCCGCCCCACCCAGCAGCTCCGACTGTTTCTCGTGCTGCTCCAGCCAGCTGGCCATTTCACTGTTGGCGCGACGGTTAGAGAGCACCCATACGGAAG                                                                                                                                                           <  6:417546/149‑1                                 GCACCGCCTGGCTTAAATAGGGGTCAGACTGACTGTTCGCCAGCACATCACCCGCCCCACCCAGCAGCTCCGACTGTTTCTCGTGCTGCTCCAGCCAGCTGGCCATTTCACTGTTGGCGCGACGGTTAGAGAGCACCCATACGGAAGCG                                                                                                                                                         <  7:441915/149‑1                                 GCACCGCCTGGCTTAAATAGGGGTCAGACTGACTGTTCGCCAGCACATCACCCGCCCCACCCAGCAGCTCCGACTGTTTCTCGTGCTGCTCCAGCCAGCTGGCCATTTCACTGTTGGCGCGACGGTTAGAGAGCACCCATACGGAAGCG                                                                                                                                                         >  5:1071089/1‑149                                    CCGCCTGGCTTAAATAGGGGTCAGACTGACTGTTCGCCAGCACATCACCCGCCCCACCCAGCAGCTCCGACTGTTTCTCGTGCTGCTCCAGCCAGCTGGCCATTTCACTGTTGGCGCGACGGTTAGAGAGCACCCATACGGAAGCGCCA                                                                                                                                                      >  5:1323717/1‑149                                    CCGCCTGGCTTAAATAGGGGTCAGACTGACTGTTCGCCAGCACATCACCCGCCCCACCCAGCAGCTCCGACTGTTTCTCGTGCTGCTCCAGCCAGCTGGCCATTTCACTGTTGGCGCGACGGTTAGAGAGCACCCATACGGAAGCGCCA                                                                                                                                                      <  8:520051/149‑1                                                                                                                                            TCACTGTTGGCGCGACGGTTAGAGAGCACCCATACGGAAGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGTGATGTGCTGGCGAACAGTCAGTCTGACCCCTATTTAAGC                                              >  6:1295562/1‑149                                                                                                                                                       GCGACGGTTAGAGAGCACCCATACGGAAGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGA                                                                                                        >  5:230383/1‑80                                                                                                                                                       GCGACGGTTAGAGAGCACCCATACGGAAGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGA                                                                                                        <  6:230383/80‑1                                                                                                                                                                                   GCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGTGATGTGCTGGCGAACAGTCAGTCTGACCCCTATTTAAGCCAGGCGGTGCTGGATTTACAGTTCGGTCATTCCCAGCGG       >  6:292058/1‑149                                                                                                                                                                                     GCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGTGATGTGCTGGCGAACAGTCA                                                                 >  7:821088/1‑89                                                                                                                                                                                     GCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGTGATGTGCTGGCGAACAGTCA                                                                 <  8:821088/89‑1                                                                                                                                                                                     GCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGG                                                                                                         <  4:220301/49‑1                                                                                                                                                                                      CCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGTGATGTGCTGGCGAACAGTCAGTCTGACCCCTATTTAAGCCAGGCGGTGCTGGATTTACAGTTCGGTCATTCCCAGCGGGTA    <  8:214557/149‑1                                                                                                                                                                                        AACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGTGATGTGCTGGCGAACAGTCAGTCTGACCCCTATTTAAGCCAGGCGGTGCTGGATTTACAGTTCGGTCATTCCCAGCGGGTAGG  <  7:718721/149‑1                                                                                                                                                                                        AACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGTGATGTGCTGGCGAACAGTCAGTCTGACCCCTATTTAAGCCAGGCGGTGCTGGATTTACAGTTCGGTCATTCCCAGCGGGTAGG  >  7:786486/1‑149                                                                                                                                                                                        AACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGTGATGTGCTGGCGAACAGTCAGTCTGACCCCTATTTAAGCCAGGCGGTGCTGGATTTACAGTTCGGTCATTCCCAGCGGGTAGG  >  8:1045388/1‑149                                                                                                                                                                                        AACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGTGATGTGCTGGCGAACAGTCAGTCTGACCCCTATTTAAGCCAGGCGGTGCTGGATTTACAGTTCGGTCATTCCCAGCGGGTAGG  >  8:442023/1‑149                                                                                                                                                                                        AACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGTGATGTGCTGGCGAACAGTCAGTCTGACCCCTATTTAAGCCAGGCGGTGCTGGATTTACAGTTCGGTCATTCCCAGCGGGTAGG  <  8:514801/149‑1                                                                                                                                                                                                                                                                                                                                              TACCCGCTGGGAATGACCGAACTGTAAATCCAGCACCGCCTGGCTTAAATAGGGGTCAGACTGACTGTTCGCCAGCACATCACCCGCCCCACCCAGCAGCTCCGACTGTTTCTCGTGCTGCTCCAGCCAGCTGGCCATTTCACTGTTGGCGCGACGGTTAGAGAGCACCCATACGGAAGCGCCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4397084‑4396901 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcgccaACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGTGATGTGCTGGCGAACAGTCAGTCTGACCCCTATTTAAGCCAGGCGGTGCTGGATTTACAGTTCGGTCATTCCCAGCGGGTAGG  >  NC_000913/4396939‑4397086

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