New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4583137	51 (1.050)	4 (0.090)	4/270	6.7	7.9%	intergenic (‑89/+112)	hsdM/hsdR	DNA methyltransferase M/endonuclease R Type I restriction enzyme
	?	NC_000913	= 4583224	44 (0.940)					intergenic (‑176/+25)	hsdM/hsdR	DNA methyltransferase M/endonuclease R Type I restriction enzyme
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGC  >  NC_000913/4583082‑4583358                                                                                                                                               |                                                                                                                                       tCATCCAGCCGCAGAAGAAGAAATTGCCGAGGGTAATAGACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGGGCGTATAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAg                                                                                                                                        <  2:889646/143‑1 (MQ=255)     cGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCGAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGCGCCAGCCTg                                                                                                                                    >  2:147842/1‑143 (MQ=255)      gAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGa                                                                                                                                   >  1:807140/1‑143 (MQ=255)       aGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGaa                                                                                                                                  <  2:552336/143‑1 (MQ=255)         ccGCAGGAGAAGAAAGTGCGGAGGGTAATATGCGCAAGCTCAGTCAGGTTGCAGCCAGGCGGCAAGTGAGTGAAGCCCCGGTAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAgg                                                                                                                                <  2:117930/143‑1 (MQ=25)          cGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGa                                                                                                                               >  1:130364/1‑143 (MQ=255)          cGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGa                                                                                                                               >  1:320166/1‑143 (MQ=255)           gCACAAGAAGAAACTGCAGAGGGTAACATACACAAGATCATGCAGGTTGCATCAAGGCGGCAAGGGCGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGAt                                                                                                                              <  2:885924/143‑1 (MQ=37)              gaagaagaaATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGaa                                                                                                                           <  1:732352/143‑1 (MQ=255)                        ttGCCGAGGGTAATATACACAAAAGCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACgt                                                                                                                 <  2:832055/143‑1 (MQ=255)                           ccGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTc                                                                                                              >  1:743953/1‑143 (MQ=255)                               gggTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGc                                                                                                          <  2:611672/143‑1 (MQ=255)                                  tAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTATGTGTCCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAg                                                                                                       <  1:500240/143‑1 (MQ=255)                                   agtataCACAAAATCATTCAGGTTGCATCAAGGCGGCAAGGGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGc                                                                                                      <  2:608652/141‑1 (MQ=255)                                         cacaAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTcc                                                                                                <  1:140992/143‑1 (MQ=255)                                          acaAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTccg                                                                                               >  1:294488/1‑142 (MQ=255)                                            aaaaTCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAg                                                                                             >  1:662552/1‑143 (MQ=255)                                              aaTCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGat                                                                                           >  1:735242/1‑143 (MQ=255)                                                tCATGCGGGTTGCATCAGGGCGGGAAGTGAGTGAATCCTGGGGGGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGATGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGatat                                                                                         <  2:563220/143‑1 (MQ=38)                                                  aTTGAGGGTGCAGCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATaa                                                                                       <  2:329565/143‑1 (MQ=255)                                                      aGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGGGTAGACGTGTCAGGCCAGCTGGTCCCAGATATAATCGc                                                                                   >  1:40982/1‑143 (MQ=255)                                                        gTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTg                                                                                 >  1:614161/1‑143 (MQ=255)                                                            cATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAAtt                                                                             >  1:108176/1‑143 (MQ=255)                                                             ggcAACGCGGCAAGTGAGAGGATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAAttt                                                                            <  1:402467/141‑1 (MQ=255)                                                             aTCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAAttt                                                                            >  2:63658/1‑143 (MQ=255)                                                                 aGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGACGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGccc                                                                        >  2:769495/1‑143 (MQ=255)                                                                    cgggAAGTGAGTGAAGTGCGAGGATCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCcagc                                                                     <  2:576181/143‑1 (MQ=38)                                                                     ggcAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCGGATATAATCGCTGAATTTGCCcagca                                                                    <  2:803634/143‑1 (MQ=255)                                                                                     cgcGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATc                                                    <  1:674702/141‑1 (MQ=255)                                                                                     cccGGGCGCGTACAGAAGTACGTGCCCGGGGGGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATc                                                    <  1:730876/143‑1 (MQ=255)                                                                                     cccGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATc                                                    >  1:860500/1‑143 (MQ=255)                                                                                            gCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAg                                             <  1:547435/143‑1 (MQ=255)                                                                                               tACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGGCAGGCCAGCTCGTCGCAGATATAATCGCTGCATTTGCCGAGCAGGGGATCGCGCTTATCGCCACAGGtt                                          >  1:282535/1‑143 (MQ=255)                                                                                               tACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGACCGGGCTACAAGCGTCTGATATAATCGCTGCATTTGCCCCGCAGGGTATCGAGATTATCGTCAAAGGtt                                          >  1:202804/1‑143 (MQ=255)                                                                                                    aaGTACGTGACCGGGGTGGACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTg                                     <  1:875695/143‑1 (MQ=255)                                                                                                    aaGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTg                                     >  2:35110/1‑143 (MQ=255)                                                                                                      gTAGGTGACCGGGGTGAACGAGCGCAGCCAACGCACAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATAAAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTgca                                   <  2:400349/143‑1 (MQ=255)                                                                                                      gTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTgca                                   <  2:407027/143‑1 (MQ=255)                                                                                                      gTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTgca                                   <  1:624318/143‑1 (MQ=255)                                                                                                             aCCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGc                            <  2:67026/143‑1 (MQ=255)                                                                                                             aCCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGc                            >  1:246549/1‑143 (MQ=255)                                                                                                                  ggTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTcc                       <  2:662552/143‑1 (MQ=255)                                                                                                                          aGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGa               >  1:630585/1‑143 (MQ=255)                                                                                                                          aGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGa               >  1:521082/1‑143 (MQ=255)                                                                                                                             gcAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTCTCGTCACAGGTTCTTTGCAGCATCGCCTCCCcgccgcgccgc            >  1:808890/1‑142 (MQ=255)                                                                                                                              cAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATAGAATCGGTGAATTTGCCCAGCAGGGTATCGAGAGTCTCGTCAAAGGTTATTTGCAGCAGGGCCTTCGCGCCGAGACGGt           >  2:62978/1‑143 (MQ=255)                                                                                                                               aGCCAACGCAGAGGCAGCCTGCAGGATGAAGTGTATACATGTCAGGCCAGCTCGTCCCAGATATACTCGCTGAATTTGCGCAGCAGCGTATCGAGATTATCTTCAAAGGTGCTGTGCAGCATCGCCTTCCCTCCGCGACGGTg          >  2:190810/1‑143 (MQ=255)                                                                                                                                       cggagGCAGCCGGAAGGATGAGGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCACTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGc  <  2:212316/141‑1 (MQ=255)                                                                                                                                               |                                                                                                                                       TCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGC  >  NC_000913/4583082‑4583358

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