New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_bme_pgi	= 1299498	0 (0.000)	19 (0.580)	16/272	1.2	100%	intergenic (+253/‑1684)	ychE/oppA	UPF0056 family inner membrane protein/oligopeptide ABC transporter periplasmic binding protein
	?	NC_000913_3_bme_pgi	1300698 =	0 (0.000)					intergenic (+1453/‑484)	ychE/oppA	UPF0056 family inner membrane protein/oligopeptide ABC transporter periplasmic binding protein

TAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAA  >  NC_000913_3_bme_pgi/1300558‑1300846                                                                                                                                             |                                                                                                                                                     ttacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATc                                                                                                                                              >  1:306450‑M2/141‑149 (MQ=255)             ataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatccagggcttgctgaagaataattgaaatcatattattaattccactgactttggtagaggaaagtgctaaATAATAATCAATTGTTAAatt                                                                                                                                  >  2:500962‑M2/129‑149 (MQ=255)              taaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAatta                                                                                                                                 >  1:334691‑M2/128‑149 (MQ=255)              taaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAatta                                                                                                                                 >  2:44883‑M2/128‑149 (MQ=255)                 ggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATtg                                                                                                                              >  1:580425‑M2/125‑149 (MQ=255)                              tgcttctccttagcatcatgtccttttcgctacatcataatcctgggcctgctgaaaaataattgaaatgatattattaattccactgcctttggtagaggaaagtaataaATAATAATCAATTGTTAAATTATTGTGCATTTCactac                                                                                                                 <  1:174904‑M2/38‑1 (MQ=255)                              tgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTaaa                                                                                                                                    <  2:4785‑M2/19‑1 (MQ=255)                              tgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTaaa                                                                                                                                    >  1:4785‑M2/112‑130 (MQ=255)                                       tcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattacttccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGt                                                                                                        <  1:550326‑M2/47‑1 (MQ=255)                                           catcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTATGCATTTCACTACTGGAACTGTAATc                                                                                                    >  1:502192‑M2/99‑149 (MQ=255)                                                        tcgccacatcataattctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaACAATAATCAATTGTTAAATTATTGTGTATTTCACTACTTGAACTGTAATCAGAAAAGATAGAc                                                                                       <  2:354329‑M2/64‑1 (MQ=255)                                                        tcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGAc                                                                                       <  2:349973‑M2/64‑1 (MQ=255)                                                                         tgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATctct                                                                      >  1:489742‑M2/69‑149 (MQ=255)                                                                                                    atattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGAtgtttgtt                                           <  1:500962‑M2/108‑1 (MQ=255)                                                                                                           taattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCAt                                    <  1:44883‑M2/115‑1 (MQ=255)                                                                                                                  actgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGCTTGACTTGAAAGATGTTTGTTAAGGCATGGATgcc                             >  2:12611‑M2/28‑148 (MQ=255)                                                                                                                  actgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCa                             >  2:376719‑M2/28‑149 (MQ=255)                                                                                                                            gtagaggaaagagctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGAt                   >  1:375851‑M2/18‑149 (MQ=255)                                                                                                                              agaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTc                 <  2:334691‑M2/134‑1 (MQ=255)                                                                                                                                             ataataatCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATaa  >  1:268861/1‑149 (MQ=255)                                                                                                                                             |                                                                                                                                                     TAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAA  >  NC_000913_3_bme_pgi/1300558‑1300846

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