New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 285973	21 (1.030)	3 (0.160)	3/274	2.0	13.5%	coding (579/1383 nt)	yagG	CP4‑6 prophage; putative sugar transporter
	?	NC_000913	= 286392	19 (1.020)					coding (998/1383 nt)	yagG	CP4‑6 prophage; putative sugar transporter
Rejected: Frequency below cutoff threshold.

GCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCGCGTTCAAGATGAT  >  NC_000913/285824‑286107                                                                                                                                                      |                                                                                                                                       gCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTG                                                                                                                                        >  2:121346‑M1/1‑150 (MQ=255)       cAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGctc                                                                                                                                  <  1:256567/150‑1 (MQ=255)         gTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTa                                                                                                                                >  2:161038/1‑150 (MQ=255)           ccTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGt                                                                                                                             <  1:271272/151‑1 (MQ=255)                  cttcttCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGctt                                                                                                                       >  1:223881/1‑150 (MQ=255)                    tcttcCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGGGGTGCTGCTCTACGTCTGcttctt                                                                                                                    >  1:5765/1‑151 (MQ=255)                                   cggGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAgc                                                                                                      <  1:1767/150‑1 (MQ=255)                                     ggCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAgcg                                                                                                     >  1:270652/1‑149 (MQ=255)                                         cgctcgctATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTacac                                                                                               <  1:301459/151‑1 (MQ=255)                                                 aTCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAgg                                                                                       <  1:97074/151‑1 (MQ=255)                                                      cGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTg                                                                                     <  1:101339/148‑1 (MQ=255)                                                        gCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGc                                                                                 >  2:27946/1‑150 (MQ=255)                                                                   ccGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGGGGGGgtattaaagacgaacaggaacagggggttgagggcaaaaatgagcgcgatgtgctccgccggggtga                                                                     >  2:107143‑M1/1‑84 (MQ=255)                                                                        ggTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCg                                                                >  1:158601/1‑151 (MQ=255)                                                                              catcatCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGAc                                                          <  1:197666/151‑1 (MQ=255)                                                                                    cGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGgtattaaagacgaacaggaacaggatgttgagggcaaaaatgagcgcgatgtgctccgccggggtgacgaaaatcagcagg                                                       >  1:64349‑M1/1‑67 (MQ=255)                                                                                      gCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAgc                                                   <  1:249696/150‑1 (MQ=255)                                                                                              gACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGgtattaaagacgaacaggaacaggatgttgagggcaaaaatgagcgcgatgtgctccgccggggtgacgaaaatcagcaggctgatcagcgagg                                          >  2:72792‑M1/1‑57 (MQ=255)                                                                                                               cTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGgcg                          >  1:2533/1‑150 (MQ=255)                                                                                                                 tCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAGAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGcat                       <  2:179291/151‑1 (MQ=255)                                                                                                                 tCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGcat                       <  1:92055/151‑1 (MQ=255)                                                                                                                     cgcCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCAtg                   <  1:287534/151‑1 (MQ=255)                                                                                                                          tgtgCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTgcgc              <  2:180761/151‑1 (MQ=255)                                                                                                                                     gggCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCGCGTTCAAgatg    <  2:195240/150‑1 (MQ=255)                                                                                                                                      ggCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCGCGTTCAAgatgat  >  2:170372/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                       GCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGGCTCGTCGGTGGCGAAAGACCTTAAGCTGCTGCTGGGCAACAGCCAGTGGCGCATCATGTGCGCGTTCAAGATGAT  >  NC_000913/285824‑286107

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