New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1207805	0 (0.000)	15 (0.620)	13/128	0.7	100%	coding (305/630 nt)	ycfK	e14 prophage; protein StfP
	?	NC_000913	= 1209602	0 (0.000)					pseudogene (54/537 nt)	stfE	e14 prophage; putative side tail fiber protein fragment

CGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAT  >  NC_000913/1209523‑1209679                                                                                |                                                                              cGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                              >  1:216861‑M2/1‑80 (MQ=255)             aTTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaac                                                                  >  2:216862‑M2/1‑68 (MQ=255)              tttCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaacc                                                                 >  2:216863‑M2/1‑67 (MQ=255)                    ccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggtt                                                           <  2:216864‑M2/81‑21 (MQ=255)                      agccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggtttt                                                         >  1:216865‑M2/1‑59 (MQ=255)                       gccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttc                                                        >  1:216866‑M2/1‑58 (MQ=255)                              ggCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagcc                                                 <  1:216645‑M2/81‑31 (MQ=255)                              ggCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagcc                                                 <  1:216646‑M2/81‑31 (MQ=255)                              ggCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagcc                                                 <  2:216647‑M2/81‑31 (MQ=255)                                  gAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgttt                                             >  2:216867‑M2/1‑47 (MQ=255)                                        cTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacag                                       <  1:216644‑M2/81‑41 (MQ=255)                                            aGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgc                                    <  1:216643‑M2/80‑44 (MQ=255)                                                 aCGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatcc                               <  2:216642‑M2/80‑49 (MQ=255)                                                     aaCCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgattt                          <  1:216641‑M2/81‑54 (MQ=255)                                                          ggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatat                     <  2:216640‑M2/81‑59 (MQ=255)                                                                    aGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacgg           >  1:216868‑M2/1‑13 (MQ=255)                                                                      gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggat         <  1:216639‑M2/81‑71 (MQ=255)                                                                             agagccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattctggcg  <  2:448292‑M2/81‑78 (MQ=255)                                                                                |                                                                              CGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAT  >  NC_000913/1209523‑1209679

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