New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	CP009273	1203246 =	3 (0.160)	8 (0.470)	8/244	NT	73.9%	coding (290/630 nt)	stfP	e14 prophage; uncharacterized protein
	?	CP009273	1205075 =	3 (0.180)					pseudogene (1/501 nt)	stfE	pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

CGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCT  >  CP009273/1204935‑1205205                                                                                                                                             |                                                                                                                                   cGGAAAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGATTGGCTGAAGTTCACGGTTCCGTGATATCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGa                                                                                                                            <  1:104257/149‑1 (MQ=255)                  tccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                                                                                           <  2:33173‑M2/26‑1 (MQ=255)                                    ttactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                                                                                         <  1:197012‑M2/44‑1 (MQ=255)                                       ctgattgctaaatacccccattcgtaagcagcatggctatgaatccggtgcgtcctaatgtagccgttgctgcattgagttttgcttcttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCt                                                                                      <  1:119643‑M2/47‑1 (MQ=255)                                                           ctcataaccgctatagcgatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctccccaaccaaCGTTTATGAAAATGAAGTAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACa                                                                  <  1:39218‑M2/67‑1 (MQ=255)                                                                            cGGAGTGGCTGAAGGCCACGGAACCGGGTCACCAACAGGTAATGAAAAACATTATCCCAAACCAACGATTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACAAGCTAATTGGCTATGTa                                                 >  1:111766/1‑149 (MQ=255)                                                                                      gtgcgtcctaatgtagccgttgctgcattaagttttggtccttctcccaaaccaaCGTTTATGTAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaa                                       >  1:212296‑M2/56‑149 (MQ=255)                                                                                                               cattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCCGGGGGATTTAACATGCTTATTGGCTATGAACGCGTATCAACAAATGACCAGAACACAGATCTACa              >  2:162676‑M2/31‑149 (MQ=255)                                                                                                                    agttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTa         <  1:246225‑M2/124‑1 (MQ=255)                                                                                                                      ttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAAt       >  2:43963‑M2/24‑149 (MQ=255)                                                                                                                      tGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAAt       >  1:292583/1‑149 (MQ=255)                                                                                                                           atccttctcccaaaccagCGTTTATGAAGATGAAGAAATAACAAGCCAATGTCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTTTCATCAAATGACCAGAACACAGATCTACAATGTGATGCGCt  <  2:308608‑M2/131‑1 (MQ=255)                                                                                                                                             |                                                                                                                                   CGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCT  >  CP009273/1204935‑1205205

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