New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2288919 =	NA (NA)	32 (0.940)	31/484	0.2	100%	noncoding (1195/1195 nt)	IS5	repeat region
	?	NC_000913	= 4577768	0 (0.000)					coding (191/1047 nt)	mcrC	Type IV methyl‑directed restriction enzyme EcoKMcrBC subunit

ACTGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATAT‑CTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTG  >  NC_000913/2288698‑2289142                                                                                                                                                                                                                               |                                                                                                                                                                                                                                aCTGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATAT‑CTTTTTCCAGGACAACGCCATCAACTGTCTCACCATACGCAACAGTTGATGCCAATGCTATGGGGGAAGATaaggaagg                                                                                                                                                                                                                            <  2:223281/227‑1 (MQ=255) aCTGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATAT‑CTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGCATATaaggaagg                                                                                                                                                                                                                            >  1:223281/1‑227 (MQ=255)                             cAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATAT‑CTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATttctt                                                                                                                                                                         <  1:44486/250‑1 (MQ=255)                              aaTGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATAT‑CTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCa                                                                                                                                                                       <  1:13365/251‑1 (MQ=255)                                                    aTTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATAT‑CTTTTTCCAGGACAACACCATCCACTGTCTCACCATACTCAAAAGTTGATGCCAATGCTATGGGTGAATATAAGGAAAGTGCGAATAAACGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTTTTCATGTTTTAGCCGATTTTTTTTcc                                                                                                                                                 >  2:33914/1‑251 (MQ=39)                                                       tCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATAT‑CTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCCACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGt                                                                                                                                              >  1:2214/1‑251 (MQ=255)                                                                         ccAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATAT‑CTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGcc                                                                                                                            <  2:250828/251‑1 (MQ=255)                                                                                cTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGCCGACCCCATACACCAGCTGGATAT‑TTTTTTCCAGCAAAACCCCATAAAATCACTCACCCTACTCAAAAGTTGATTCCACTACTATGGGAGAATATAAGGAAGGTACGAATAAGCGGGGCAATTCTTATCGCCTGCCCTAGAACTTTCATTTCTTCATGTTTAAGACGGTTTTTTAACCCGTAAATGCCTTTAATCAGCCTCttt                                                                                                                       >  1:210955/1‑249 (MQ=255)                                                                                    tataTGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCCTAAACCAGCTGGATATACTTTTTC‑AGGACAACACCATCAACTGTCTCACCCTACTCCACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGGATAATCGGGGGAATTCTTTTCGGCTGACTCCGTCATTTCATTTCTTCGTGTTTGGGCCGATTTTTTCTCCCGTAAATGCCTTTTATCAGCCTCTTTATAccc                                                                                                                 >  1:2483/1‑250 (MQ=255)                                                                                                            gattTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATAT‑CTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtcatt                                                                                                                                                                                  <  1:302723/162‑1 (MQ=255)                                                                                                            gattTTCGTATTATTGGCGGGCCCCAAAAACAGCCGGATAT‑CTTTTTCCCGGACAACCCCCTCCACTGTCTCACCATACTCAACAGATGATGCCAATGATATGGGGGGATATAAGGAAAGGGCGAAAAAGCCGGGGAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTTAGCCGACTTTTTCTCCCGTAAATACCTTGTATCCGCCTTTTTTTGCCGTTTCTTTGCCATTTAAAGCTTTAt                                                                                         >  2:116772/1‑251 (MQ=16)                                                                                                               tttCGTATTATTGGCGGTCCCATAAACCAGCTGGATAT‑CTTGTTCCAGGACAACACCATCAACTGGCTCACCATACTCAACAGTTGATGCCAATGCTATGTGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACGCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATcc                                                                                       <  2:359831/250‑1 (MQ=255)                                                                                                                               gTCCCATAAACCAGCTGGATAT‑CTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTg                                                                          >  2:35976/1‑247 (MQ=255)                                                                                                                                 cccATAAACCAGCTGGATAT‑CTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGAtct                                                                    >  1:317894/1‑251 (MQ=255)                                                                                                                                             cTGGATAT‑CTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACg                                                        >  2:213579/1‑251 (MQ=255)                                                                                                                                                         tttCCAGGACAACACCATCAACTGTCTCACAATCCTCAACAGTTGATGCCAATGCTATGGGTGAATATCAGGAAGGTGCGAATAAGCGGGGAACTTCTTCTCAGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGt                                             >  1:338145/1‑251 (MQ=255)                                                                                                                                                                   aacaCCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGaa                                    >  2:47114/1‑250 (MQ=255)                                                                                                                                                                         aTCAACTGTCTCCCCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGtt                             <  1:212229/251‑1 (MQ=255)                                                                                                                                                                          tCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGATTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTg                            >  2:281272/1‑251 (MQ=255)                                                                                                                                                                                    cACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAG‑AAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCGCGCCCGGTGACGTATCAGGTGGTCCGCCCGAAACAGGTTGGCCAGCGTg                   <  1:143759/249‑1 (MQ=14)                                                                                                                                                                                              aaCAGTTGATGCCAATACTATGGGTGAATATAAGGCAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGc        >  1:353990/1‑251 (MQ=255)                                                                                                                                                                                                    tGATCCCAAAGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCAGTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGGAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTTGCCAGCGTGAATAACATCGCCAGTTg  >  2:280840/1‑251 (MQ=255)                                                                                                                                                                                                                               |                                                                                                                                                                                                                                ACTGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATAT‑CTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTG  >  NC_000913/2288698‑2289142

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