New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	= 1203261	10 (0.280)	8 (0.270)	7/138	1.4	49.8%	coding (305/630 nt)	BW25113_RS05990	hypothetical protein
	?	NZ_CP009273	= 1205058	8 (0.270)					coding (57/540 nt)	BW25113_RS06005	phage tail protein

GAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTT  >  NZ_CP009273/1204972‑1205146                                                                                       |                                                                                         gAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGcct                                                                                       <  1:551386‑M2/90‑4 (MQ=255)         ggCACCATGCCGTGTCAGACACTCTGTTGGCAGAAGGGGTGGATGGAAAGGAACCCGGGCACCAACAGGTAATGCAGAGCCTTCTCCCaa                                                                               <  1:523242/90‑1 (MQ=255)           cACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaac                                                                             <  2:691931‑M2/90‑14 (MQ=255)               aTTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAAc                                                                         >  2:542096/1‑90 (MQ=255)                       tCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcga                                                                 <  1:198162‑M2/90‑26 (MQ=255)                          gccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagag                                                              <  1:703578‑M2/90‑29 (MQ=255)                          gccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagag                                                              <  1:159305‑M2/90‑29 (MQ=255)                          gccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGaaaa                                                              >  2:236526/1‑90 (MQ=255)                                    tGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcac                                                    >  2:376941‑M2/1‑52 (MQ=255)                                            tGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaagggtttcgagagccgtttgcacagtgccat                                            <  1:774602‑M2/90‑47 (MQ=255)                                                   aGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcat                                     <  1:236526/90‑1 (MQ=255)                                                         cGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                               <  1:928621/90‑1 (MQ=255)                                                              ccGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaa                          >  2:198162‑M2/1‑26 (MQ=255)                                                                gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagag                                                              <  1:619025‑M2/52‑29 (MQ=255)                                                                gggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagag                                                              >  2:619025‑M2/1‑24 (MQ=255)                                                                            ggTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAtt            <  1:700643/90‑1 (MQ=255)                                                                                  gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACAt      <  2:285908/90‑1 (MQ=255)                                                                                  gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACAt      <  2:678989/90‑1 (MQ=255)                                                                                      agCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCtt  >  1:636668/1‑90 (MQ=255)                                                                                       |                                                                                         GAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTT  >  NZ_CP009273/1204972‑1205146

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