New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1207790 =	2 (0.040)	25 (0.500)	23/568	1.2	91.1%	coding (290/630 nt)	stfP	e14 prophage; uncharacterized protein
	?	NC_000913	1209619 =	3 (0.060)					pseudogene (1/501 nt)	stfE	pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

GTTGCATCAGTCACAATCGTGGATCTTGACGGTAATCCATGAGCATGATCCTCCGTTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCTGGAAGCTGGATCGGCTGGGGCGTAGTATGCGGCATCTTGTCGTGCTGGTGGAGG  >  NC_000913/1209321‑1209919                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                             tgagctgcgcaccgtaataaaaacgaagtggttgtatatccccccgttgccactgcattgataggtgcgctgacattgtcataatttgcaacaacggcaacacaggcattattaaacgcaaccgggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgcttcttcttccaaaccaaCGt                                                                                                                                                                                                                                                                                                            >  1:417897‑M2/299‑301 (MQ=255)    gCATCAGTCACAATCGTGGATCTTGACGGTAATCCATGAGCATGATCCTCCGTTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTa                                                                                                                                                                                                                                                                                                         >  1:86562/1‑301 (MQ=255)         gcaccgtaataaaaacgaagtggttgtatatccccccgttgccactgcattgataggtgcgctgacattgtcataatttgcaacaacggcaacacaggcattattaaacgcaaccgggaatacaactgtatttccatcccccccagagccagatgctttggcggcattccacccccactgaatgattattgattgctcagtaccaccaatcataactggtacagattagaaaacggttcgtcctaatgtacccgttgcggtcttgagttcagctctttatctcaactcacCCTTTAGGaaa                                                                                                                                                                                                                                                                                                    >  2:468900‑M2/291‑301 (MQ=255)                                                             aaaggtgggctaacattgtcataagttgcaacaacggaaacacaggcattattaaacgcaaccgggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAttt                                                                                                                                                                                                                                                <  2:50059‑M2/63‑1 (MQ=255)                                                              cccccTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATCGCATCATTCCTGCTTTTACCAGGGGGATTTa                                                                                                                                                                                                                                               <  2:86561/300‑1 (MQ=255)                                                                    cgctgacattgtcataatttgcaacaacggcaacacaggcattattaaacgcaaccgggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGTGATTTAACAtga                                                                                                                                                                                                                                         >  1:482448‑M2/232‑300 (MQ=255)                                                                     gctgacattgtcataatttgcaacaacggcaacacaggcattattaaacgcaaccgggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaatacacccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGAAAATGTCATCATTACTGCTTTTAAAAGGGGGATTTAACATg                                                                                                                                                                                                                                          >  2:148004‑M2/231‑299 (MQ=255)                                                                              gtcataatttgcaacaacggcaacacaggcattattaaacgcaaccgggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgaggtttgctccttctcccaaactaaCGTTTCTGAAAATGAAGAAATAACAGACAAATTGCATCAGTCCTGCTTGTACCAGGGGGATTGACCAAGTTTattggttt                                                                                                                                                                                                                               >  2:427346‑M2/222‑298 (MQ=255)                                                                                    atgtgaaacaaaggcaacacaggcattattaaacgcaaccgggcataaaactgtagttccatcccccccagaggcagatgcttttgcggcattccacccacactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAcg                                                                                                                                                                                                                         <  1:348423‑M2/86‑1 (MQ=255)                                                                                                 gcaacacaggcattattaaacgcaaccgggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTTACATGCTTATTGGCTATGTACGCGTATCAACAAAt                                                                                                                                                                                                            >  1:31475‑M2/203‑301 (MQ=255)                                                                                                 gcaacacaggcattattaaacgcaaccgggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTTGCTATGTACGCGTATCAAcaaaa                                                                                                                                                                                                            >  1:284949‑M2/203‑300 (MQ=255)                                                                                                                  aaacgcaaccgggaatacaactgtatttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTCCGCGTATCAACAAATGACCAGAACACAGAt                                                                                                                                                                                             >  1:245655‑M2/186‑299 (MQ=255)                                                                                                                                     actgtatttccatcccccccagaggcagatgcttttgcggcactccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAAc                                                                                                                                                                        <  2:284945‑M2/135‑1 (MQ=255)                                                                                                                                            ttccatcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAg                                                                                                                                                                 <  1:202993‑M2/142‑1 (MQ=255)                                                                                                                                                 tcccccccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgacttttgctccttctcccaaaccaaCGTTTCTGAACATGAAGAAATAACAAGCACATGTCATCATTCCTCCTTTTACCAGGGGGATTTACCATGCTTATTGGCCTTGTCGGCGCACCAACAATTGCCCAGAACACAGACATACAACGTAATGCCCTGAACGGTGCAGGATg                                                                                                                                                             >  2:369791‑M2/155‑300 (MQ=255)                                                                                                                                                                 ctcccataccgtccttcccccccctctgaatgattactgattgcataatacctccaaacataaccggtatagctaataaaccggtgcgtcctaatttagccgttgctgcattgaattttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAg                                                                                                                                            <  2:255024‑M2/163‑1 (MQ=255)                                                                                                                                                                        ttccgccatccacccccacttaatgcttacttattttgctataccaccaatcaaaaccgttatagctatgaaacctgtgtgtcctaatgtcgccgttgctgcattgagttttgctccttctcccaaaccaaCTTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTTCTTTTACCAGGGGGATTTAACATGCTTATTTGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACATCGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGAt                                                                                                                                     <  2:25833‑M2/170‑1 (MQ=255)                                                                                                                                                                          cggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGCCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTTAACTGTGCAGGATGCGAgctgatttttgatgcgatcggaa                                                                                                                                   >  1:255027‑M2/130‑290 (MQ=255)                                                                                                                                                                            gaattccaccaccactgaatgattaccgattgctaaataccaccaatcataaccgggatagcgatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTTACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAg                                                                                                                                  <  2:31475‑M2/173‑1 (MQ=255)                                                                                                                                                                                                                        aatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACAAGAACACAGATCTACAACGTCATGCGCTGAACTGGGCAGGATGCGAGATGATTTTTGAAGAAAAGATAAGCGGCACAAAGTCCGAAAggccgggacagaaaaaaatcgacaggaa                                                                                     >  2:476563‑M2/84‑290 (MQ=255)                                                                                                                                                                                                                            cgccagccTGTTGGCGGGATGGCTGACGGCAACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTa                                                                                 <  1:210561/300‑1 (MQ=255)                                                                                                                                                                                                                               accggtttcatagctataccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCg                                                                              >  1:159621‑M2/77‑301 (MQ=255)                                                                                                                                                                                                                                  ggtaaggctatgaaaccgggtctgccgaaagtcgccgttgctgcattgagttgtgctccttctccaaaaccaaCGTTTATGAAAATGAAGAAAGAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGc                                                                            <  2:317303‑M2/227‑1 (MQ=255)                                                                                                                                                                                                                                         ctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGac                                                                     <  1:148006‑M2/234‑1 (MQ=255)                                                                                                                                                                                                                                               aaccggtgcgtcctaatgtagccgttgctggattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAGCTGCTCAGGACATTATCGGCAGGTGACACTCTgg                                                              <  1:468910‑M2/241‑1 (MQ=255)                                                                                                                                                                                                                                                ccccgttcttcctaatgtagccgtgtgtggattgagttttgctcctggtaccaaaccaaCGTTTATGAAAATGAAGAAATAACGAGTGAATGGCATCATTCCTGCTTTTACAGGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGGGAGCAGATTTTTGAAGAAAAGATAGGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGt                                                             <  2:417887‑M2/242‑1 (MQ=255)                                                                                                                                                                                                                                                                          atgcagtgagttttgctcctgatcccaaaccaaCGTTAATGAAAAGGAAGAAATAACAAGCAAATGGCAACATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCTGGAAGCTGGATCGGCTgggg                                    <  2:159619‑M2/267‑1 (MQ=255)                                                                                                                                                                                                                                                                              attgagttttgctccttcccccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGATTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACATATCTACAACGTAATGAGTTGACCTGTGCAGGATGCGAGTTGATATTTGACGACCAGATACGCGGCACAAGGCACGAAAGGCAGGGACGGAAACAACTGCCCAGAACATAATCGGCAGGAGACACCTCCGTTTTCTGGCACCTGTAtcggcggcggcggtg                               >  2:475553‑M2/30‑298 (MQ=255)                                                                                                                                                                                                                                                                                                         aaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGTTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGAAGAGAGCTGATTTCTGAAGAAAAGAAAAGCGGCACAAAGTCCGAAAGACCGGGACTGAAAAAACTGCTCAGGACACTATGGACCGGTGACACTCTGGATgtctgcaagctgtactggcgggggcgtagtagcggacatccgggagccctgtgggg    >  2:475943‑M2/3‑276 (MQ=255)                                                                                                                                                                                                                                                                                                           cGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTGACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGGTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCTGGAAGCTGGATCGGCTGGGGCGTAGTATGCGGCATCTTGTCGTGCTGGTggagg  <  1:369799/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                             GTTGCATCAGTCACAATCGTGGATCTTGACGGTAATCCATGAGCATGATCCTCCGTTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCTGGAAGCTGGATCGGCTGGGGCGTAGTATGCGGCATCTTGTCGTGCTGGTGGAGG  >  NC_000913/1209321‑1209919

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