New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_vch_pgi	1207790 =	17 (0.770)	5 (0.260)	5/258	2.0	21.2%	coding (290/630 nt)	stfP	e14 prophage; uncharacterized protein
	?	NC_000913_3_vch_pgi	1209619 =	22 (1.130)					pseudogene (1/501 nt)	stfE	pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

CAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGG  >  NC_000913_3_vch_pgi/1207644‑1207938                                                                                                                                                   |                                                                                                                                                     cgtcctggagagttcggcgccttacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataaccatgttaactcccactagtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTg                                                                                                                                                    >  2:160899‑M1/147‑149 (MQ=255) cAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCCCGAATACGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTg                                                                                                                                                    <  1:172178/149‑1 (MQ=255) cAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTGAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTg                                                                                                                                                    <  1:44908/149‑1 (MQ=255)  atcctgcacagttcatcggattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTgg                                                                                                                                                   >  2:119569‑M1/146‑149 (MQ=255)      tgcacagttcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTg                                                                                                                                               >  1:342877‑M1/142‑149 (MQ=255)        gtggAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTggg                                                                                                                                             <  2:247324/149‑1 (MQ=255)                                      gggCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAATGGCTCTAGCAAACCTTGGTTGGGGAGAAGGAGCAGAGCTCAATGCAGCAACGGc                                                                                                               >  1:45414/1‑149 (MQ=21)                                              gcatttgttgatatgcgtatatatgcaataagcaagtaaattccccctcgtaagagcaggattgatgccatgtgcttgttatttcttcattttcataaccgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGACTACATTAg                                                                                                       >  2:56657‑M1/102‑149 (MQ=255)                                                           acCTCACTTCGCGCGCTGCTGTTGAGCCGCTAGAATCCGTTTAGAAATATGAAAGCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGttt                                                                                          <  1:51956/149‑1 (MQ=25)                                                           acCGCACTTCACGCGCTAATGTTAAGCCGCAAGAATCCGTTTGCAGAAGTCCGTTCGAGTGGCACTGTGCAAACGACTCTCGAAAACCTTGGTTTGGGTGATGGGACAGAACTCAAGGCAGCAACGGCTACATTGGGACGCACGGGgtt                                                                                          >  2:258261/1‑149 (MQ=16)                                                              gtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGACTACATTCGGACGCACCGGTTTCAt                                                                                       <  2:13956‑M1/64‑1 (MQ=255)                                                              gCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGGGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCAt                                                                                       >  1:63989/1‑149 (MQ=255)                                                                  ttCGTGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATAGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCt                                                                                   <  2:210301/149‑1 (MQ=255)                                                                      ccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATAc                                                                               <  2:218488‑M1/72‑1 (MQ=255)                                                                             gcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTCTACCGGTTAt                                                                        <  2:223904‑M1/79‑1 (MQ=255)                                                                              tgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCCCCGGTTTCATAGCTATACCGGTTATg                                                                       >  1:42134/1‑149 (MQ=255)                                                                                   aGCCGCAAGAACCCGTTTGGCGCTATCAAATCGGATGGCACTGTGCAACCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCATTGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtgg                                                                  <  1:131177/149‑1 (MQ=255)                                                                                            aaTCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTAttaa                                                         <  2:262304/149‑3 (MQ=255)                                                                                              tCCGTGTGGCGATATCAAATCGGATGGCACGGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAAATCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGc                                                       <  2:123285/149‑1 (MQ=255)                                                                                                   ttGGCGATATCAAATCGGAGGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGCGAGTCGGTAAGACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCa                                                  <  1:381880/149‑1 (MQ=38)                                                                                                       cGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTaa                                              <  2:369469/149‑1 (MQ=255)                                                                                                                    gATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTgggg                                 >  1:356783/1‑149 (MQ=255)                                                                                                                    gATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTgggg                                 <  2:275895/149‑1 (MQ=255)                                                                                                                               tttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCg                      <  1:190144‑M1/129‑1 (MQ=255)                                                                                                                                                  cTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCTGAGACGGTTGAATGCCGCAAAAGCATCTGCCTCTgg   <  2:180361/149‑1 (MQ=255)                                                                                                                                                   TTGGCTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTCGGACGCACCGGTTTCATAGCTAGACGGGGTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAACACATCTGCCTCTggg  <  2:189071‑M1/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                     CAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGG  >  NC_000913_3_vch_pgi/1207644‑1207938

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