New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1207795 =	14 (0.170)	3 (0.040) +C	3/422	NT	20.1%	coding (295/630 nt)	stfP	e14 prophage; uncharacterized protein
	?	NC_000913	1209615 =	10 (0.120)					pseudogene (5/501 nt)	stfE	pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

ACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTA  >  NC_000913/1209452‑1209840                                                                                                                                                                    |                                                                                                                                                                                                                                  accacGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGGGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGCGCCTTCTCCCAAACCAACGGTTAGGAAAATGAAGAAAAAACAAGCAAAAGGCAACATTCCCCTTTTTAACAGGGGGGTTTAAAAGGCtt                                                                                                                                                        >  2:754605/1‑239 (MQ=255)                   cccagaggcagatgcttttgcggcattccacccccactgcatgattactgattgctcaataccaccaatcataaccggtctagctatgaaacaggtgcgtcctaatgtagccgttgctgcatttagttttgctccttctcccaacCCACCGTTTATGAAAATGACGAATTAACAAGCAAATGGCATCATTCCTGCTTTTACACGGGGGATTTAAAATGCTTATTGGCTATGTACGGGTa                                                                                                                                      >  2:860654‑M2/146‑239 (MQ=255)                          gTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaacaa                                                                                                                               <  1:68068/239‑1 (MQ=255)                                aaaGGCTTTGCCAGGTCCGGGGATTCTTCAGCCGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGAc                                                                                                                          <  1:754605/237‑1 (MQ=255)                                         gCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATACCAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGACTTAAAATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAg                                                                                                                >  1:495117/1‑239 (MQ=255)                                           cgtgtCCGGGTTTTTTTCAGCGTAAATAGCGGGACCATGGGATGTCAGCCACCATGGTGGCGGAGGGGCTGAAGGCCACGGAACCGGGAATCCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGACTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGAt                                                                                                              <  2:495117/235‑1 (MQ=255)                                                accccccctgaatgttttcttgttgctcaataccccccattcataccgggttagcttggaaacctgtgcgtcctaatgtagccgttggtgcattgagttttgctccgtctcccaacCCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCGACa                                                                                                         <  2:356858‑M2/123‑1 (MQ=255)                                                      actgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaacCCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCGGCTTTTACAAGGGGGATTTAACATGATTATTGGATATGTACACGTATCAAAAAATGAACAGAAACCAAATATAAAACGTaa                                                                                                   >  2:1073136‑M2/111‑239 (MQ=255)                                                      aTTCTTCAGCAGAAAAAGCGGCACCCTTGCATTTCAGCCAGGCTGGTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTa                                                                                                    <  1:784011/238‑1 (MQ=255)                                                                                                                                 ggACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAACGACCAGAACACAGATCTACAACTTAATGCGCTGCACTCTGCAGGATGCGAGCCGATTTTTGCAAACAAAAGAAGGGGCACAAACACCAAAAGGCCGGGCCt                        >  2:5118/1‑239 (MQ=255)                                                                                                                                             gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAAc                                                                                         <  1:285349/162‑1 (MQ=255)                                                                                                                                                      agCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACAtt   >  1:239131/1‑239 (MQ=255)                                                                                                                                                       gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTa  >  1:617509/1‑239 (MQ=255)                                                                                                                                                                    |                                                                                                                                                                                                                                  ACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTA  >  NC_000913/1209452‑1209840

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