New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1207805	6 (0.190)	15 (0.590)	14/128	0.8	75.4%	coding (305/630 nt)	ycfK	e14 prophage; protein StfP
	?	NC_000913	= 1209602	5 (0.200)					pseudogene (54/537 nt)	stfE	e14 prophage; putative side tail fiber protein fragment

GCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAC  >  NC_000913/1209522‑1209670                                                                                 |                                                                     gCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCagag                                                                      >  2:237413/1‑81 (MQ=255)  cGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                     <  1:237414‑M2/81‑2 (MQ=255)    gCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGcct                                                                   >  1:237415‑M2/1‑78 (MQ=255)      aCCATTGCATTTCAGCCAGCCTGTTGGCGGGGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttc                                                                 >  1:237416‑M2/1‑76 (MQ=255)         aTTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcc                                                              <  2:237417‑M2/81‑9 (MQ=255)            gCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctccca                                                            <  1:237418‑M2/80‑11 (MQ=255)              aTTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaac                                                         <  1:237461‑M2/81‑14 (MQ=255)               tttCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaacc                                                        >  1:237419‑M2/1‑67 (MQ=255)                   agccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaagg                                                    <  2:237420‑M2/81‑19 (MQ=255)                      cagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttt                                                 >  1:237421‑M2/1‑60 (MQ=255)                       agccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggtttt                                                <  2:237422‑M2/81‑23 (MQ=255)                        gccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttc                                               >  2:237423‑M2/1‑58 (MQ=255)                         ccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATg                                              <  1:237424/81‑1 (MQ=255)                            gTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAGTGCAGAGccttctcccaaaccaaggttttcga                                             <  2:237189‑M2/79‑26 (MQ=255)                             ttGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgaga                                           >  2:237425‑M2/1‑53 (MQ=255)                                gCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagcc                                        >  1:237426‑M2/1‑50 (MQ=255)                                  ggAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtt                                     >  1:237427‑M2/1‑48 (MQ=255)                                   gAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgttt                                    >  2:237429‑M2/1‑47 (MQ=255)                                   gAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaaga                                    <  2:237428/81‑1 (MQ=255)                                       ggCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgca                                 <  1:237187‑M2/80‑38 (MQ=255)                                       ggCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgca                                 <  1:237188‑M2/80‑38 (MQ=255)                                        gCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcac                                >  1:237430‑M2/1‑42 (MQ=255)                                              ggCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgcca                         <  1:237186‑M2/81‑46 (MQ=255)                                                   cGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGATTATGAAAATGAAGAAATAACAAGCAAATgg                    >  2:237463/1‑81 (MQ=255)                                                     gAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatt                  >  2:237431‑M2/1‑29 (MQ=255)                                                                   aCAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaac    <  2:237183‑M2/81‑67 (MQ=255)                                                                     aGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAc  <  2:237432/81‑1 (MQ=255)                                                                     aGGTAATGCAGAGCCTTCTCCCAAACCAACGATTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAc  <  2:237433/81‑1 (MQ=255)                                                                                 |                                                                     GCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAC  >  NC_000913/1209522‑1209670

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