New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4583138	53 (1.050)	4 (0.080)	4/270	6.9	7.6%	intergenic (‑89/+112)	hsdM/hsdR	DNA methyltransferase M/endonuclease R Type I restriction enzyme
	?	NC_000913	= 4583225	47 (0.970)					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/4583083‑4583359                                                                                                                                               |                                                                                                                                       tCATCCAGCCGCAGAAGAAGAAATTGCCGAGGGTAATAGACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGGGCGTATAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAg                                                                                                                                        <  1:889646/143‑1 (MQ=255)     cGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCGAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGCGCCAGCCTg                                                                                                                                    >  1:147842/1‑143 (MQ=255)      gAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGa                                                                                                                                   >  2:807140/1‑143 (MQ=255)       aGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGaa                                                                                                                                  <  1:552336/143‑1 (MQ=255)         ccGCAGGAGAAGAAAGTGCGGAGGGTAATATGCGCAAGCTCAGTCAGGTTGCAGCCAGGCGGCAAGTGAGTGAAGCCCCGGTAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAgg                                                                                                                                <  1:117930/143‑1 (MQ=25)          cGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGa                                                                                                                               >  2:320166/1‑143 (MQ=255)          cGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGa                                                                                                                               >  2:130364/1‑143 (MQ=255)           gCACAAGAAGAAACTGCAGAGGGTAACATACACAAGATCATGCAGGTTGCATCAAGGCGGCAAGGGCGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGAt                                                                                                                              <  1:885924/143‑1 (MQ=37)              gaagaagaaATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGaa                                                                                                                           <  2:732352/143‑1 (MQ=255)                     aaaTTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCACCGCAGAGGCAGCCTGAAggctttcatggcga                                                                                                                    >  1:210142/1‑131 (MQ=38)                        ttGCCGAGGGTAATATACACAAAAGCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACgt                                                                                                                 <  1:832055/143‑1 (MQ=255)                           ccGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTc                                                                                                              >  2:743953/1‑143 (MQ=255)                               gggTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGc                                                                                                          <  1:611672/143‑1 (MQ=255)                                  tAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTATGTGTCCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAg                                                                                                       <  2:500240/143‑1 (MQ=255)                                   agtataCACAAAATCATTCAGGTTGCATCAAGGCGGCAAGGGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGc                                                                                                      <  1:608652/141‑1 (MQ=255)                                         cacaAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTcc                                                                                                <  2:140992/143‑1 (MQ=255)                                          acgagtgcgttCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTccc                                                                                               <  2:782984/134‑1 (MQ=255)                                          acaAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTccg                                                                                               >  2:294488/1‑142 (MQ=255)                                            aaaaTCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAg                                                                                             >  2:662552/1‑143 (MQ=255)                                              aaTCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGat                                                                                           >  2:735242/1‑143 (MQ=255)                                                tCATGCGGGTTGCATCAGGGCGGGAAGTGAGTGAATCCTGGGGGGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGATGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGatat                                                                                         <  1:563220/143‑1 (MQ=38)                                                  aTTGAGGGTGCAGCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATaa                                                                                       <  1:329565/143‑1 (MQ=255)                                                      aGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGGGTAGACGTGTCAGGCCAGCTGGTCCCAGATATAATCGc                                                                                   >  2:40982/1‑143 (MQ=255)                                                        gTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTg                                                                                 >  2:614161/1‑143 (MQ=255)                                                            cATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAAtt                                                                             >  2:108176/1‑143 (MQ=255)                                                             ggcAACGCGGCAAGTGAGAGGATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAAttt                                                                            <  2:402467/141‑1 (MQ=255)                                                             aTCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAAttt                                                                            >  1:63658/1‑143 (MQ=255)                                                                 aGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGACGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGccc                                                                        >  1:769495/1‑143 (MQ=255)                                                                    cgggAAGTGAGTGAAGTGCGAGGATCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCcagc                                                                     <  1:576181/143‑1 (MQ=38)                                                                     ggcAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCGGATATAATCGCTGAATTTGCCcagca                                                                    <  1:803634/143‑1 (MQ=255)                                                                                     cgcGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATc                                                    <  2:674702/141‑1 (MQ=255)                                                                                     cccGGGCGCGTACAGAAGTACGTGCCCGGGGGGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATc                                                    <  2:730876/143‑1 (MQ=255)                                                                                     cccGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATc                                                    >  2:860500/1‑143 (MQ=255)                                                                                         gtcacgtacttctgTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCa                                                >  2:380597/14‑143 (MQ=255)                                                                                            gCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAg                                             <  2:547435/143‑1 (MQ=255)                                                                                               tACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGGCAGGCCAGCTCGTCGCAGATATAATCGCTGCATTTGCCGAGCAGGGGATCGCGCTTATCGCCACAGGtt                                          >  2:282535/1‑143 (MQ=255)                                                                                               tACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGACCGGGCTACAAGCGTCTGATATAATCGCTGCATTTGCCCCGCAGGGTATCGAGATTATCGTCAAAGGtt                                          >  2:202804/1‑143 (MQ=255)                                                                                                    aaGTACGTGACCGGGGTGGACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTg                                     <  2:875695/143‑1 (MQ=255)                                                                                                    aaGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTg                                     >  1:35110/1‑143 (MQ=255)                                                                                                      gTAGGTGACCGGGGTGAACGAGCGCAGCCAACGCACAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATAAAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTgca                                   <  1:400349/143‑1 (MQ=255)                                                                                                      gTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTgca                                   <  2:624318/143‑1 (MQ=255)                                                                                                      gTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTgca                                   <  1:407027/143‑1 (MQ=255)                                                                                                             aCCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGc                            >  2:246549/1‑143 (MQ=255)                                                                                                             aCCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGc                            <  1:67026/143‑1 (MQ=255)                                                                                                                  ggTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTcc                       <  1:662552/143‑1 (MQ=255)                                                                                                                          aGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGa               >  2:630585/1‑143 (MQ=255)                                                                                                                          aGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGa               >  2:521082/1‑143 (MQ=255)                                                                                                                             gcAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTCTCGTCACAGGTTCTTTGCAGCATCGCCTCCCcgccgcgccgc            >  2:808890/1‑142 (MQ=255)                                                                                                                              cAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATAGAATCGGTGAATTTGCCCAGCAGGGTATCGAGAGTCTCGTCAAAGGTTATTTGCAGCAGGGCCTTCGCGCCGAGACGGt           >  1:62978/1‑143 (MQ=255)                                                                                                                               aGCCAACGCAGAGGCAGCCTGCAGGATGAAGTGTATACATGTCAGGCCAGCTCGTCCCAGATATACTCGCTGAATTTGCGCAGCAGCGTATCGAGATTATCTTCAAAGGTGCTGTGCAGCATCGCCTTCCCTCCGCGACGGTg          >  1:190810/1‑143 (MQ=255)                                                                                                                                       cggagGCAGCCGGAAGGATGAGGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCACTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGc  <  1:212316/141‑1 (MQ=255)                                                                                                                                               |                                                                                                                                       TCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGC  >  NC_000913/4583083‑4583359

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