New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1207790 =	8 (0.310)	13 (0.670)	10/100	0.5	76.4%	coding (290/630 nt)	stfP	e14 prophage; uncharacterized protein
	?	NC_000913	1209619 =	2 (0.100)					pseudogene (1/501 nt)	stfE	pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

TGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTT  >  NC_000913/1209559‑1209681                                                             |                                                               aaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATg                                                         >  1:201851‑M2/61‑68 (MQ=255)              taatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                            <  1:601911‑M2/21‑1 (MQ=255)              aaatgtagccgttgatgcattgagttttgctccttctcccaaacctaCGTTTATGAAAATGAAGTAAt                                            <  2:738098‑M2/21‑1 (MQ=255)                 tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAAc                                         >  1:773043‑M2/45‑68 (MQ=255)                 tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAAc                                         >  2:54150‑M2/45‑68 (MQ=255)                  gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACa                                        >  1:660697‑M2/44‑68 (MQ=255)                     gccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGc                                     <  2:408316‑M2/28‑1 (MQ=255)                                attgagttttgctcattctcccaaaccaaCGATTATGAAAATGAAGAAATAACAAGCAAATGGcatca                          >  1:932599‑M2/30‑68 (MQ=255)                                aGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatca                          <  1:823085/68‑1 (MQ=255)                                 ggTAATGCAGAGACTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatcat                         >  2:436771/1‑68 (MQ=255)                                  tgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                        >  1:533247‑M2/28‑68 (MQ=255)                                     gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                     >  2:846001‑M2/25‑68 (MQ=255)                                       tttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                   <  2:417835‑M2/46‑1 (MQ=255)                                       tttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                   <  2:501315‑M2/46‑1 (MQ=255)                                        ttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCt                  <  1:814122‑M2/47‑1 (MQ=255)                                                    ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg      <  2:201851‑M2/59‑1 (MQ=255)                                                     caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa     >  2:113523‑M2/9‑68 (MQ=255)                                                     caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa     <  1:915325‑M2/60‑1 (MQ=255)                                                      aaacctaCGTTTATGATAATGAAGATATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGTgag    <  1:860111‑M2/61‑2 (MQ=255)                                                        accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAttt  >  1:244769‑M2/6‑68 (MQ=255)                                                             |                                                               TGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTT  >  NC_000913/1209559‑1209681

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