New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1672789 =	10 (0.110)	6 (0.070)	5/180	NT	41.4%	intergenic (+8/+31)	ydgD/mdtI	putative peptidase/multidrug efflux system transporter
	?	NC_000913	1672815 =	8 (0.100)					intergenic (+34/+5)	ydgD/mdtI	putative peptidase/multidrug efflux system transporter

CCGGTTTTCGCGACAAGCTGGATCAACTGTCGCAAAAATAATGTTCAAACGCTGCCCGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGATCCAACCTGCGGCTAACGTGGCGGCAATACCAAACCCGCCCCACAATGCATAA  >  NC_000913/1672741‑1672958                                                                           |                                                                                                                                                ccGGTTTTCGCGACAAGCTGGATCAACTGTCGCAAAAATAATGTTCAAACGCTGCCCGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAgg                                                                                               <  1:355480/125‑1 (MQ=255) ccGGTTTTCGCGACAAGCTGGATCAACTGTCGCAAAAATAATGTTCAAACGCTGCCCGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAgg                                                                                               >  2:355480/1‑125 (MQ=255)                                   aaaaTAATGTTCAAACGCTGCCCGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGc                                                                                     >  3:286644/1‑101 (MQ=255)                                   aaaaTAATGTTCAAACGCTGCCCGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGc                                                                                     <  4:286644/101‑1 (MQ=255)                                                    ctgcccgcgctgtcgggcagcgtCTTCATCAGGCAAGTTTc                                                                                                                                >  4:703517‑M2/24‑41 (MQ=255)                                                    ctgcccgcgctgtcgggcagcgtCTTCATCAGGCAAGTTTc                                                                                                                                <  3:703517‑M2/18‑1 (MQ=255)                                                    ctgcccgcgctgtcgggcagcgtCTTCATCAGGCAAGTTTCACCATGa                                                                                                                         <  4:592309‑M2/25‑1 (MQ=255)                                                    ctgcccgcgctgtcgggcagcgtCTTCATCAGGCAAGTTTCACCATGa                                                                                                                         >  3:592309‑M2/24‑48 (MQ=255)                                                         cgcgctgtcgggcagcgtCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGATCCAACCTGCGGc                                          <  3:439390‑M2/104‑1 (MQ=255)                                                         cgcgctgtcgggcagcgtCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGATCCAACCTGCGGc                                          >  4:439390‑M2/19‑122 (MQ=255)                                                           aCAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAg                                                                                      <  4:949145/76‑1 (MQ=255)                                                           aCAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAg                                                                                      >  3:949145/1‑76 (MQ=255)                                                           aCAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGa                                                       >  2:490973/1‑107 (MQ=255)                                                           aCAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGa                                                       <  1:490973/107‑1 (MQ=255)                                                              gcgcGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGa                                                       <  4:1115818/104‑1 (MQ=255)                                                              gcgcGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGa                                                       >  3:1115818/1‑104 (MQ=255)                                                                       gcgtCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGATCCAACCTGCGGCTAACGTGGCGGCAATACCAAACCCGCCCCACAATGCATaa  >  3:998838‑M2/5‑148 (MQ=255)                                                                           |                                                                                                                                                CCGGTTTTCGCGACAAGCTGGATCAACTGTCGCAAAAATAATGTTCAAACGCTGCCCGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGATCCAACCTGCGGCTAACGTGGCGGCAATACCAAACCCGCCCCACAATGCATAA  >  NC_000913/1672741‑1672958

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