New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_vch_pgi	1207790 =	22 (1.200)	5 (0.310)	5/256	1.5	23.0%	coding (290/630 nt)	stfP	e14 prophage; uncharacterized protein
	?	NC_000913_3_vch_pgi	1209619 =	14 (0.860)					pseudogene (1/501 nt)	stfE	pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

GGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGG  >  NC_000913_3_vch_pgi/1207656‑1207937                                                                                                                                       |                                                                                                                                                    ggCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAg                                                                                                                                       >  1:164995/1‑149 (MQ=255)   cATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGa                                                                                                                                     <  1:176663/149‑1 (MQ=255)    aTCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGCTTGGCAGAAggct                                                                                                                                    >  2:227711/1‑147 (MQ=255)    aTCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTGGGCGATATCAAATCGGATGGCAGTGTCCACACGGCTCACGAAAACATCGCTTTAGGACAAgccc                                                                                                                                    >  2:227810/1‑146 (MQ=17)         attacgttgtagatctgtggtctggtcatttgttgatacgggtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCaaaa                                                                                                                               >  1:59414‑M1/127‑149 (MQ=255)              gttgtatatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCaa                                                                                                                          >  1:78365‑M1/122‑149 (MQ=255)                                 gtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCACTGCAGCAACGGCTACATTa                                                                                                       <  2:173666‑M1/47‑1 (MQ=255)                                              tacCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGtt                                                                                          <  2:164995/149‑1 (MQ=255)                                                   cACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTgg                                                                                                                                             <  2:84479/93‑1 (MQ=255)                                                   cACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTgg                                                                                                                                             >  1:84479/1‑93 (MQ=255)                                                        cgcgcgCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTggg                                                                                                                                            >  2:9349/1‑89 (MQ=255)                                                        cgcgcgCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTggg                                                                                                                                            <  1:9349/89‑1 (MQ=255)                                                                        gCCGCAAGACTCCGGTTGGCGTTAGCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtggt                                                                <  2:283960/149‑1 (MQ=255)                                                                        aaatccccctggtaatagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtggt                                                                >  2:75919‑M1/64‑149 (MQ=255)                                                                              aGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTAtt                                                          <  1:174089/149‑1 (MQ=255)                                                                                 agccGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAg                                                       <  1:264850/147‑1 (MQ=255)                                                                                   ccGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCa                                                     <  1:280952/149‑1 (MQ=255)                                                                                           aggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGGTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTaa                                             <  1:150464‑M1/105‑1 (MQ=255)                                                                                                 cAAAACGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTc                                       <  1:274096/149‑1 (MQ=255)                                                                                                  aaaTCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAAAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGAATTGAGCAATCAGTAATCATTCa                                      >  1:261218/1‑149 (MQ=255)                                                                                                   aaTCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCACTCAGTAATCATTCAg                                     >  1:160426/1‑149 (MQ=255)                                                                                                        gATGGCATTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTgggg                                <  1:288258/149‑1 (MQ=255)                                                                                                        gATGGCAACGGGCAAACGGCACCCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTgggg                                <  2:261218/149‑1 (MQ=255)                                                                                                                         gggtctcGAAAACCTGGGGTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGTTATTGAGCAATCAGTAAGCATTCTGTGGGGGTGGAATGCCGCAAAAg               >  2:32198/4‑149 (MQ=255)                                                                                                                         ggCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGGTGGGGTGGAGTGCCGCAAAAg               >  2:257509/1‑149 (MQ=255)                                                                                                                                      cTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAGTCTTTCAGTGGGGGTGGAGTGCTGCAGAAGCATCTGCATCTgg  >  1:23001/1‑149 (MQ=255)                                                                                                                                      cTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAAATTAGGACGCACCGGTTTCATAGCTGTACCGGTTATGATTGGTGGTATTGAGCAATCCGTAATCAGTCAGTGGGGGTGGAATGCCGTAAAAGCATCTGCCTCTgg  >  1:246855/1‑149 (MQ=255)                                                                                                                                       |                                                                                                                                                    GGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGG  >  NC_000913_3_vch_pgi/1207656‑1207937

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