New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4578534	2 (0.060)	33 (1.050)	33/478	0.1	95.7%	coding (321/897 nt) coding (804/1380 nt)	mcrB mcrB	5‑methylcytosine‑specific restriction enzyme subunit McrB 5‑methylcytosine‑specific restriction enzyme subunit McrB
	?	NC_000913	4583087 =	1 (0.030)					intergenic (‑39/+162)	hsdM/hsdR	type I restriction enzyme EcoKI methylase subunit/type I restriction enzyme EcoKI endonuclease subunit

TGCTCCTGGCCGATGCGGGATTTCAGGTCATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCAT  >  NC_000913/4582849‑4583329                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                   gcgtagggcatagtcaacaacggccagagagcgatcggcagtattcattaaaccgatgatataaacattctccgggacatagaatcgttcttcatcgttttcggagtaggttaggggaacagaccagttttcacctcgtttatcatgttccattaacatcatcacttcgccaaatactttactgagattggcacgattgatttcatctataataaaaatatactttttctctggctgcGAGCCGCagaag                                                                                                                                                                                                                                         <  2:109745‑M2/12‑1 (MQ=255)         catagtcaacaacggccagagagcgatcggcagtattcattaaaccgatgatataaacattctccgggacatagaatcgttcttcatcgttttcggagtaggttaggggaacagaccagttttcacctcgtttatcatgttccattaacatcatcacttcgccaaatactttactgagattggcacgattgatttcatctataataaaaatatactttttctctggctgcGAGCCGAAGAAGAAGAAATTg                                                                                                                                                                                                                                >  1:61810‑M2/231‑251 (MQ=255)               caacaacggccagagagcgatcggcagtattcattaaaccgatgatataaacattctccgggacatagaatcgttcttcatcgttttcggagtaggttaggggaacagaccagttttcacctcgtttatcatgttccattaacatcatcacttcgccaaatactttactgagattggcacgattgatttcatctataataaaaatatactttttctctggctgcGAGCCGCAGAAGAAGAAATTGCCGAgg                                                                                                                                                                                                                          <  1:8785‑M2/27‑1 (MQ=255)                                                ttaaaccgatgatataaacattctccgggacatagaatcgttcttcatcgttttcggagtaggttaggggaacagaccagttttcacctcgtttatcatgttccattaacatcatcacttcgccaaatactttactgagattggcacgattgatttcatctataataaaaatatactttttctctggctgcGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCa                                                                                                                                                                                          >  2:295827‑M2/192‑250 (MQ=255)                                                    accgatgatataaacattctccgggacatagaatcgttcttcatcgttttcggagtaggttaggggaacagaccagttttcacctcgtttatcatgttccattaacatcatcacttcgccaaatactttactgagattggcacgattgatttcatctataataaaaatatactttttctctggctgcGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAggcg                                                                                                                                                                                     >  2:38470‑M2/188‑251 (MQ=255)                                                        atgatataaacattctccgggacatagaatcgttcttcatcgttttcggagtaggttaggggaacagaccagttttcacctcgtttatcatgttccattaacatcatcacttcgccaaatactttactgagattggcacgattgatttcatctataataaaaatatactttttctctggctgcGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCa                                                                                                                                                                                  >  2:245451‑M2/184‑250 (MQ=255)                                                                  cattctccgggacatagaatcgttcttcatcgttttcggagtaggttaggggaacagaccagttttcacctcgtttatcatgttccattaacatcatcacttcgccaaatactttactgagattggcacgattgatttcatctataataaaaatatactttttctctggctgcGAGCCGCAGCAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAagtgagtg                                                                                                                                                                          >  2:142462‑M2/174‑248 (MQ=255)                                                                  cattctccgggacatagaatcgttcttcatcgttttcggagtaggttaggggaacagaccagttttcacctcgtttatcatgttccattaacatcatcacttcgccaaatactttactgagattggcacgattgatttcatctataataaaaatatactttttctctggctgcGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAagtgagtg                                                                                                                                                                          <  1:142462‑M2/75‑1 (MQ=255)                                                                        ggTCTCTTTACACAGTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCTCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCgg                                                                                                                                                                 <  2:165111/251‑1 (MQ=255)                                                                          gggacatagaatcgttcttcatcgttttcggagtaggttaggggaacagaccagttttcacctcgtttatcatgttccattaacatcatcacttcgccaaatactttactgagattggcacgattgatttcatctataataaaaatatactttttctctggctgcGAGCCGCAGAAGAAGAAATTGCCGCGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTTAATCCCCggg                                                                                                                                                                >  2:315781‑M2/166‑250 (MQ=255)                                                                           ggacatagaatcgttcttcatcgttttcggagtaggttaggggaacagaccagttttcacctcgtttatcatgttccattaacatcatcacttcgccaaatactttactgagattggcacgattgatttcatctataataaaaatatactttttctctggctgcGAGCCGCAGAAGACGAAATTGCCGAGGGTAATATACACAACATCATTCCGGTTGCATCAAGGCGGCAAGGGAGTGAATCACCggg                                                                                                                                                                >  2:166430‑M2/165‑249 (MQ=255)                                                                               atagaatcgttcttcatcgttttcggagtaggttaggggaacagaccagttttcacctcgtttatcatgttccattaacatcatcacttcgccaaatactttactgagattggcacgattgatttcatctataataaaaatatactttttctctggctgcGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGt                                                                                                                                                           >  1:103195‑M2/161‑250 (MQ=255)                                                                                 agaatcgttcttcatcgttttcggagtaggttaggggaacagaccagttttcacctcgtttatcatgttccattgacatcatcacttcgccaaatactttactgcgatcggcacgattgatttcatctataataaaaatatactttttctctggctgcGAGCCGCAGACGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCCGGTTGCATCAAGGCTGACAGTGAATGAATCACCGGGCGCGTacg                                                                                                                                                        >  2:238122‑M2/159‑250 (MQ=255)                                                                                                cgttgtcggagtaggttaggggaacagaccagttttcacctcgtttatcatgttccattaacatcatcacttcgccaaatactttactgagattggcacaattgatgtcatctaaaataaaaatattctttttcactggctgcGTGCCGCAGAAGAAGAAAGTGCCGAGGGTAAGATACGCAAAATCATTCAGGTTGTATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCgg                                                                                                                                         <  2:290811‑M2/108‑1 (MQ=255)                                                                                                            aggttaggggaacagaccagttttcacctcgtttatcatgttccattaacatcatcacttcgccaaatactttactgagattggcacgattgatttcatctataataaaaatatactttttctctggctgcGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAgc                                                                                                                              <  1:166430‑M2/119‑1 (MQ=255)                                                                                                                    ggaaaagacctgttttcaccgcgtttataatgttccattaacatcatcacttcgccaaataatttactgagagtggcacgattgatttcatctataataaaaatatactttttctcgggctgcGAGCCGCAGAAGAAGAAATTGCCGATGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAAc                                                                                                                     <  1:64647‑M2/128‑1 (MQ=255)                                                                                                                            ccagttttcacctcgtttatcatgttccattaacatcatcacttcgccaaatactttactgagattggcacgattgatttcatctataataaaaatatactttttctctggctgcGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGc                                                                                                             <  2:103195‑M2/136‑1 (MQ=255)                                                                                                                               gttttcacctcgtttatcatgttccattaacatcatcacttcgccaaatactttactgagattggcacgattgatttcatctataataaaaatatactttttctctggctgcGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGc                                                                                                          <  1:66969‑M2/139‑1 (MQ=255)                                                                                                                                             tatcatgttccattaacatcatcacttcgccaaatactttactgagattggcacgattgatttcatctataataaaaatatactttttctctggctgcGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAg                                                                                             <  1:315781‑M2/152‑1 (MQ=255)                                                                                                                                              atcatgttccattaacatcatcacttcgccaaatactttactgagattggcatgattgatttcatctataataaaaatatactttttctctggctgcGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAgt                                                                                            <  1:90743‑M2/153‑1 (MQ=255)                                                                                                                                               tcatgttccattaacatcatcacttcgccaaatactttactgagattggcacgattgatttcatctataataaaaatatactttttctctggctgcGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCACCGCAGAGGCAGCCTGAAGGATGAAgtg                                                                                           >  1:80010‑M2/97‑250 (MQ=255)                                                                                                                                                    ttccattaacatcatcacttcgccaaatactttactgagattggcacgattgatttcatctataataaaaatatactttttctctggctgcGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACg                                                                                     >  1:9412‑M2/92‑251 (MQ=255)                                                                                                                                                                                actttactgagattggcacgattgatttcatctataataaaaatatcctttttctctggctgcGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATaa                                                          <  2:170345‑M2/187‑1 (MQ=255)                                                                                                                                                                                                  cgattgatttcatctataataaaaatatactttttctctggctgcGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCcagc                                        >  1:93146‑M2/46‑250 (MQ=255)                                                                                                                                                                                                             atctataataaaaatatactttttctctggctgcGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGa                                                                                                                                             <  2:211813‑M2/104‑1 (MQ=255)                                                                                                                                                                                                             atctataataaaaatatactttttctctggctgcGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGa                                                                                                                                             >  1:211813‑M2/35‑138 (MQ=255)                                                                                                                                                                                                             aactataataaaaatatactttttctctggctgcGAGCCTCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCgaga                            <  2:246982‑M2/217‑1 (MQ=255)                                                                                                                                                                                                               ctataataaaaatatactttttctctggctgcGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGAtt                          >  2:282691‑M2/33‑251 (MQ=255)                                                                                                                                                                                                                        aaatatactttttctctggctgcGAGCCGCAGACGAAGAAATTGCCGCGGTTAATATACACAAAATCATTCTGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTTACCGGGGTGAACGAGCGCAGCCAACGCCGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCCCGACCCAGATATAATCGCATAATTTTCCCAGCAGAGTATCGATATTTTCGTCaa                  >  2:223667‑M2/24‑250 (MQ=255)                                                                                                                                                                                                                                 ttttctctggctgcGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTg        >  1:255002‑M2/15‑251 (MQ=255)                                                                                                                                                                                                                                    tctcgggatgcGAGCGGCAGAAGAAGAAATTGCCGATGGTAATATACACAAAATCATGCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGGATACGTGTCAGGCCAGCTCGGCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGAGTATCGTCAAAGGTTCTTTgca      <  2:15945‑M2/239‑1 (MQ=255)                                                                                                                                                                                                                                       ctggctgcGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGCATTTGCCCAGCAGGGTCTCGAGATTATCGTCACAGGTTCTTTGCAGCAt  >  2:100445‑M2/9‑251 (MQ=255)                                                                                                                                                                                                                                           ctgcGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAAt                                                 >  2:208879‑M2/5‑200 (MQ=255)                                                                                                                                                                                                                                           ctgcGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAAt                                                 <  1:208879‑M2/196‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                   TGCTCCTGGCCGATGCGGGATTTCAGGTCATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCAT  >  NC_000913/4582849‑4583329

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