New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1207805	5 (0.210)	8 (0.350)	7/460	NT	77.4%	coding (305/630 nt)	stfP	e14 prophage, uncharacterized protein
	?	NC_000913	= 1209602	0 (0.000)					pseudogene (18/501 nt)	stfE	pseudogene, e14 prophage, side tail fiber protein fragment family,Phage or Prophage Related

GGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACA  >  NC_000913/1207533‑1208055                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                           ggttAAGTTAGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGCCCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAg                                                                                                                                                                                                                                                          <  1:26950/272‑1 (MQ=255)        tccagacaaccagaggttctcctgccgataatgtcctgagcagttttttcagtcccggcctttcggactttgtgccgcttatcttttcttcaaaaatcagctcgcatcctgcacagttcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaagccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcagtttcataaacgTTGGTTTGGGAGAAGGAGCAAAacg                                                                                                                                                                                                                                                   >  2:68045‑M1/251‑274 (MQ=255)         tCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTACCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTAAAGCCGCAAGAATCCGTGTGGCGAGATAAAAACGGATGGCAGTGTGAAAACGGCTATCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACt                                                                                                                                                                                                                                                   >  2:219393/1‑274 (MQ=255)             cgcgGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacct                                                                                                                                                                                                                                               >  2:37056‑M1/1‑261 (MQ=255)             cgcgGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCCCTAGAAAACCTTGGTTTGGGCGAAGGctcggcattccc                                                                                                                                                                                                                                                >  1:43371‑M1/1‑261 (MQ=255)                 ggTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGGGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGTCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGc                                                                                                                                                                                                                                                           <  2:43371‑M1/258‑2 (MQ=255)                      gtgtcccctgccgataagtgcctgtgcagttttttcagtcccggtctttcggactttgggccgcttagcttgtcttcacaaatcagctcgcattctgcacagttcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACgg                                                                                                                                                                                                                                     >  2:169645‑M1/237‑275 (MQ=255)                      gtattagagaccgattatgtcctgagcagttttttcagtcccggcctttcggactttgtgccgcttatcttgtcttcaaaaatcagctcgcatcctgcacagtgcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACgg                                                                                                                                                                                                                                     >  2:198927‑M1/237‑275 (MQ=255)                               gccgataatgtcctgagcagttttttcagtcccggcctttcggactttgtgccgcttatcttgtcttcaaaaatcagctcgcatcctgcacagttcagcgcattacgttttagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttattttttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATTCAGATAAGGCTAcattcg                                                                                                                                                                                                                            <  2:125421‑M1/48‑3 (MQ=255)                               gccgataatgtcctgagcagttttttcagtcccggcctttcggactttgtgccgcttatcttgtcttcaaaaatcagctcgcatcctgcacagttcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAg                                                                                                                                                                                                                            <  2:207604‑M1/48‑1 (MQ=255)                                           cccccAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagcc                                                                                                                                                                                                                 <  1:37056‑M1/274‑44 (MQ=255)                                            tgggcagttttttcagtcccggcctttcggactttgtgccgcttatcttgtcttcaaaaatcagctcgcatcctgcacagttcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGttt                                                                                                                                                                                                               >  1:198049‑M1/215‑275 (MQ=255)                                              agcagttttttcagtcccggcctttcggactttgtgccgcttatcttgtcttcaaaaatcagctcgcatcctgcacagttcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAatgcagga                                                                                                                                                                                                                                         <  1:198926‑M1/35‑3 (MQ=255)                                               cAGACCGGAACGCCTGCCACCGATCTGGATGATGACTACTTGGGCATGTTCCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCAt                                                                                                                                                                                                            <  2:4143/275‑1 (MQ=255)                                                      gCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATCTCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAGACTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaaaa                                                                                                                                                                                                     >  1:136315‑M1/1‑220 (MQ=255)                                                                              tgtgccgcttatcttgtcttcaaaaatcagctcgcatcctgcacagttcagcgcattacgttgtagatctgtgttatggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAg                                                                                                                                                                             >  1:207603‑M1/181‑275 (MQ=255)                                                                                                                  gCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaag                                                                                                                                         >  2:19630‑M1/1‑160 (MQ=255)                                                                                                                                                atctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAAt                                                                                                           <  2:38030‑M1/161‑1 (MQ=255)                                                                                                                                                          gcACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGCAATACAGTTGTAt                                                                                                 >  2:40468/1‑275 (MQ=255)                                                                                                                                                                       cttaccgagctcagggcgctgcggttaagccgcaAGAATCCGGTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgcctgatttacgtggtgagtttattcgtggct                                                                                     <  1:122360‑M1/259‑168 (MQ=255)                                                                                                                                                                                        gctgctGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgcctgatttacgtggtgagtttattcgtggctgggatgacgggcggggga                                                                   <  1:19630‑M1/275‑186 (MQ=255)                                                                                                                                                                                                                             atgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCAc                                                                                                                                                                                                                     <  2:2176‑M1/55‑1 (MQ=255)                                                                                                                                                                                                                                 catttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGt                                                                                                                                                                                                                 >  1:2176‑M1/34‑92 (MQ=255)                                                                                                                                                                                                                                                  aCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACggggg          >  2:192844/1‑274 (MQ=255)                                                                                                                                                                                                                                                     attttcaaaaacgGTGGTTGGGAGGAAGGAGCAAAACTCAATGCGACAACGGCTACATTAGGACGGACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGata      >  1:125420‑M1/14‑275 (MQ=255)                                                                                                                                                                                                                                                       tctcGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgcctgatttacgtggtgagtttattcgtggctgggatgacgggcgggggattgatacaggccgctctattttaagtattcaggggcatgcaaggcggactgatgttgcgggat    >  2:139578‑M1/1‑27 (MQ=255)                                                                                                                                                                                                                                                         tcGAAAACCTTGGTTGGGGAGAAGGGGAAAAACTCAAGGCAGCAACGGCTACATTAGGACGAACCGGTTTGATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACa  <  1:37893/275‑1 (MQ=255)                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                           GGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACA  >  NC_000913/1207533‑1208055

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