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

CAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAA  >  CP009273/1203100‑1203378                                                                                                                                                   |                                                                                                                                     catcctgcacagttcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTgg                                                                                                                                   >  1:202806‑M1/147‑150 (MQ=255)            aGGCATCCGGTGCCAGCGTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAG                                                                                                                        <  1:30924‑M1/150‑1 (MQ=255)            aGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGGTTGGGAGAAG                                                                                                                        >  1:81721‑M1/1‑150 (MQ=255)               cagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAg                                                                                                                     <  1:54577‑M1/18‑1 (MQ=255)                                                   aGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGa                                                                                     <  2:124219/146‑1 (MQ=255)                                                       actgacCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagc                                                                             <  1:244450‑M1/150‑43 (MQ=255)                                                                    cgcgcgCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaaca                                                                <  1:41066‑M1/150‑56 (MQ=255)                                                                      cgcgcTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGACAACCTTGGTTTGGGAGAAGGctcggcactacctgttggtgtcccggttccgtggccttcagccactccgccaa                                                                  >  2:152312‑M1/1‑93 (MQ=255)                                                                             gcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATg                                                       <  1:112852‑M1/80‑1 (MQ=255)                                                                                tgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatatcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGAtt                                                    <  1:169732‑M1/83‑1 (MQ=255)                                                                                 gttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGAGCCGGTTTCATAGCTATACCGGTCATGGTtg                                                   <  1:199659‑M1/84‑1 (MQ=255)                                                                                                                aTCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTGGGGAGAAGGctctgcattacccgtcggcgtcccggttccgtggccttcagccactccgccaacacgctggctgaaatgcaatggtgccgctttttctgctgaag                        >  2:151512‑M1/1‑51 (MQ=255)                                                                                                                aTCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTGTGGGAGAAGGctctgcactacatgttggagtctcggttccgtggccttcagccactccgccaacaggctgggccaaatgcgatggtgccggtttttctgctgagg                        >  2:151542‑M1/1‑51 (MQ=255)                                                                                                                 tCGGGTGGCACTCGGTGAACGGCTCGCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACCCAAGGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTCTGATTTGTGTTATTGAGCAATCAGTAATCATTCa                       <  2:73040/146‑1 (MQ=255)                                                                                                                   ggATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatcccc                 >  1:119698‑M1/1‑48 (MQ=255)                                                                                                                   ggATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaataccc                 >  1:225054‑M1/1‑48 (MQ=255)                                                                                                                        cttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTggggg                <  2:6050‑M1/119‑1 (MQ=255)                                                                                                                          tgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTg              >  2:48312‑M1/26‑146 (MQ=255)                                                                                                                                  cttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAATGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTGAGTGGGGGTGGAATGccgccaa  <  1:5408‑M1/133‑4 (MQ=255)                                                                                                                                                   |                                                                                                                                     CAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAA  >  CP009273/1203100‑1203378

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