New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_vch_pgi	= 1207805	18 (0.930)	5 (0.290)	5/258	1.6	27.0%	coding (305/630 nt)	stfP	e14 prophage; uncharacterized protein
	?	NC_000913_3_vch_pgi	= 1209602	11 (0.640)					pseudogene (18/501 nt)	stfE	pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

CACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCA‑ACAAATGACCAGAACACAGATCTACAACGTAAT  >  NC_000913_3_vch_pgi/1209454‑1209744                                                                                                                                                     |                                                                                                                                                cacGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGCAACCGGGACACCAACAGGTAATGCagag                                                                                                                                                 >  2:28633/1‑149 (MQ=255)             gCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCGCCATTGCATTTCAGCCAGCCTGTTGGCGTAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaa                                                                                                                                     <  2:319175‑M2/149‑13 (MQ=255)                      ttGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGttt                                                                                                                            <  1:51631/149‑1 (MQ=255)                            gctaaGCCTTGTCCAGTTCCGGTTACTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgaga                                                                                                                      >  2:178859‑M2/3‑122 (MQ=255)                                 gCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaa                                                                                                                 >  2:114966/1‑149 (MQ=255)                                 gCCTTTGCCAGTTCCGGGTATTCGTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaa                                                                                                                 <  2:193701/149‑1 (MQ=255)                                                gggTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCaaa                                                                                                  >  1:288826/1‑149 (MQ=255)                                                                aaaaaGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCt                                                                                  <  1:328044/149‑1 (MQ=255)                                                                             gttgCATGTCGGCCAGCATGGTGGCGGAGTGGCTGGAGGCCACGGAACCGGGACACCAAAAGGTAATGCAGAGccttctcccaaaccaggtgtgtcgagggcagtttgctcagtgccatccgatttgatatcgccaaacggattcttgc                                                                     >  2:145950‑M2/2‑73 (MQ=255)                                                                               tGCTTTTCATCCAGCTTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGTACACCAACAGGTAATGCAGAGCCTTCTCTCAAACCAACGTTTATGAAAATGAATAAATAACAAGCAAATGGCATCATTCCTGCTTTTACTAGGGGGAtt                                                                   <  2:128861/149‑1 (MQ=255)                                                                                         ccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCtt                                                         <  2:288687/149‑1 (MQ=255)                                                                                                  tGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTAt                                                <  1:112466/149‑1 (MQ=255)                                                                                                              tGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacgaaagcagccggcctaacagcggcgcgcgaagtccggtaagta                                    <  1:8472‑M2/149‑110 (MQ=255)                                                                                                                        ggAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCA‑ACAAATGAcc                         <  1:114966/149‑1 (MQ=255)                                                                                                                               ggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgcccc                   <  1:282246‑M2/149‑127 (MQ=255)                                                                                                                                      aaCAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCA‑ACAAGTGACCAGAACACAGACCTa           >  2:90744/1‑149 (MQ=255)                                                                                                                                         aGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggc         <  1:185892‑M2/149‑137 (MQ=255)                                                                                                                                               tGCAGAGCCTTCTCCCAAACCAACTTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGAGGATTTAACATGATAATTGGCAATCTA‑ACGCTTCAGACAAATGACCAGACAACAGATCTACAACATAAt  >  2:70823/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                CACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCA‑ACAAATGACCAGAACACAGATCTACAACGTAAT  >  NC_000913_3_vch_pgi/1209454‑1209744

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