New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	CP009273	= 1203261	8 (0.290)	10 (0.430)	8/194	NT	60.9%	coding (305/630 nt)	stfP	e14 prophage; uncharacterized protein
	?	CP009273	= 1205058	6 (0.260)					pseudogene (18/501 nt)	stfE	pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

ATCCATGAGCATGATCCTCCGTTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACAC  >  CP009273/1204811‑1205184                                                                                                                                                                                                                                                        |                                                                                                                               aTCCATGAGCTGTATCCTCCGTTGCATACCCCTGAATATTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAAGAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGTTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgaga                                                                                                     >  1:336010‑M2/1‑248 (MQ=255)                              cccTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGACCCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgat                                                                        <  1:191384‑M2/275‑57 (MQ=255)                                            aaTAGCGCGACCTTTATCAATCCCCCGCCNGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCATCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAATGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcc                                                                                                                        >  1:525207‑M2/1‑205 (MQ=255)                                                                                        aaggaaactcaCCAAGGAAAGCAGGCAATGTATTTGTCGGATAAGCCTGTGCCAGTTCCGGGTATTCTGCAGCAGAAAAAGCGGCACCATGGCATTTCAGACAGCCTGGTGGCGGAGTGGCTGAAGGCAACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattattgcggcttaacagcagcgcgcgaagtgcggtaagtaactgg              >  2:191380‑M2/5‑161 (MQ=255)                                                                                                                                      cTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTa                                                             >  2:99485/1‑182 (MQ=255)                                                                                                                                          gCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAg                                                         <  1:99485/182‑1 (MQ=255)                                                                                                                                                                                      tttCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGa      >  2:571636/1‑189 (MQ=255)                                                                                                                                                                                         cagccagcCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATga                                                                                                 >  1:298643/1‑95 (MQ=255)                                                                                                                                                                                          agccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAacac  <  1:571647/189‑1 (MQ=255)                                                                                                                                                                                             cagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaagaa                                                                                             <  2:298638/95‑1 (MQ=255)                                                                                                                                                                                                                        cACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggatt                                                         <  1:114458‑M2/104‑72 (MQ=255)                                                                                                                                                                                                                          cGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttc                                                                                                         <  1:368190‑M2/54‑24 (MQ=255)                                                                                                                                                                                                                            gAACCGTGACATCAACAGGTAATGCAGAGccttctcccaaaccaaggtgttcgagagccgttttctcagtgccatccgatttgatatcgcaaaacggtttcttg                                                     >  2:114457‑M2/1‑29 (MQ=255)                                                                                                                                                                                                                              aCCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgaga                                                                                                     >  2:368185‑M2/1‑27 (MQ=255)                                                                                                                                                                                                                                  ggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttc                                                                                                         <  1:366531‑M2/46‑24 (MQ=255)                                                                                                                                                                                                                                      acCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgaga                                                                                                     >  2:366526‑M2/1‑19 (MQ=255)                                                                                                                                                                                                                                                  tGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgc                                    <  2:16796‑M2/99‑93 (MQ=255)                                                                                                                                                                                                                                                      gagccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcga                                >  1:16796‑M2/1‑3 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                               ATCCATGAGCATGATCCTCCGTTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACAC  >  CP009273/1204811‑1205184

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