New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4584195 =	36 (0.970)	6 (0.170)	6/242	3.6	15.4%	coding (2567/3513 nt)	hsdR	endonuclease R Type I restriction enzyme
	?	NC_000913	4584232 =	32 (0.920)					coding (2530/3513 nt)	hsdR	endonuclease R Type I restriction enzyme
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CTGGAGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCACCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCCTTTCATCTGTTCGTACAGAATGCGGCTGCGTACTTTACGCAG  >  NC_000913/4584050‑4584331                                                                                                                                                  |                                                                                                                                         cTGGAGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCACCGGACGCATGgtgt                                                                                                                                       <  4:21467/149‑1 (MQ=255)    gAGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCACCGGACGCATGGTGTCGa                                                                                                                                    >  5:4964/1‑149 (MQ=255)             caccaGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCACCGGACGCATGGTGTCGACGCTCTCCa                                                                                                                           <  6:4964/149‑1 (MQ=255)                  gttgttCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCACCGGACGCATGGTGTCGACGCTCTCCAGc                                                                                                                         >  6:88883/1‑146 (MQ=255)                  gttgttCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCACCGGACGCATGGTGTCGACGCTCTCCAGc                                                                                                                         <  5:88883/146‑1 (MQ=255)                                tgCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCACCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGa                                                                                                         >  1:175515/1‑148 (MQ=255)                                   tCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCACCGGACGCATGgt                                                                                                                                         >  6:152881/1‑113 (MQ=255)                                   tCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCACCGGACGCATGgt                                                                                                                                         <  5:152881/113‑1 (MQ=255)                                             tGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGTATCCGTATTTTCTTTTACCAGCGTTTGCAGTTCCACCTTCTGGCGCACCACCACCGGACTCATGTTGTCGACGCTCTCCa                                                                                                                           <  8:114083/117‑1 (MQ=255)                                             tGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCACCGGACGCATGGTGTCGACGCTCTCCa                                                                                                                           >  7:114083/1‑117 (MQ=255)                                                    aTCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCACCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCATAAAGc                                                                                    >  8:17258/1‑149 (MQ=255)                                                        gCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCACCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGt                                                                                >  3:11477/1‑149 (MQ=255)                                                           tCGGTGATGTTATAGGGTTCTGAATCGGTAATTTCATTGACCAGCGTGTGCAGTTCCACCTTCGGGCGCACCACCACCGGACGCATGGTGTCGACGCGCTCCAGCGTGCTGTAGAAATCGACACAGTCAAAAATCTTAAAGCTGGgttt                                                                             <  4:14286/149‑1 (MQ=255)                                                              gTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCACCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAACTCTTAAAGCTGGTTTTAt                                                                           >  4:120842/1‑148 (MQ=255)                                                                  ttttATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCACCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTa                                                                            <  1:14105/143‑1 (MQ=255)                                                                  ttttATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCACCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTa                                                                            >  2:14105/1‑143 (MQ=255)                                                                         ggTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCACCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCggg                                                                <  2:175515/148‑1 (MQ=255)                                                                          gTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCACCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGac                                                                                                        >  4:90046/1‑107 (MQ=255)                                                                          gTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCACCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGac                                                                                                        <  3:90046/107‑1 (MQ=255)                                                                          gTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCACCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCa                                                              >  5:148295/1‑149 (MQ=255)                                                                           tttCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCACCGGACGCATGGTGTCGACGCtctc                                                                                                                             >  2:192812/1‑85 (MQ=255)                                                                           tttCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCACCGGACGCATGGTGTCGACGCtctc                                                                                                                             <  1:192812/85‑1 (MQ=255)                                                                           tttCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCACCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCa                                                              >  1:152764/1‑148 (MQ=255)                                                                              cTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCACCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATaa                                                           <  6:148295/148‑1 (MQ=255)                                                                              cTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCACCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATaa                                                           <  3:120842/148‑1 (MQ=255)                                                                                   tCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCACCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGt                                                     <  7:17258/149‑1 (MQ=255)                                                                                      gTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCACCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGgc                                                  >  7:239/1‑149 (MQ=255)                                                                                            tCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCACCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGc                                             <  4:11477/148‑1 (MQ=255)                                                                                              aTTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCACCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCgg                                                                 >  6:165747/1‑126 (MQ=255)                                                                                              aTTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCACCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCgg                                                                 <  5:165747/126‑1 (MQ=255)                                                                                                  aCCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCACCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTTTTCACCTCCGGGCATAAGCGCGTGGCGCGGCCTTTCa                                       <  5:3502/148‑1 (MQ=255)                                                                                                      gCGTTTGCAGTTCCACCTTCGGGCGCACCACCACCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCCTTTCATCTGt                                  <  2:152764/149‑1 (MQ=255)                                                                                                                      cTTCGGGCGCACCACCACCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCCTTTCATCTGTTCGTACAGAATGCGGc                  >  2:130584/1‑149 (MQ=255)                                                                                                                           ggCGCACCACCACCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCCTTTCATCTGTTCGTACAGAATGCGGCTGCg              >  5:187/1‑148 (MQ=255)                                                                                                                                        cGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCCTTTc                                        >  7:131097/1‑109 (MQ=255)                                                                                                                                        cGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCCTTTc                                        <  8:131097/109‑1 (MQ=255)                                                                                                                                                  gtgtCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCCTTTCATCTGTTCGTACAGAATGCGGCTGCGTACTTTACGCAg  <  6:116266/137‑1 (MQ=255)                                                                                                                                                  gtgtCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCCTTTCATCTGTTCGTACAGAATGCGGCTGCGTACTTTACGCAg  >  5:116266/1‑137 (MQ=255)                                                                                                                                                  |                                                                                                                                         CTGGAGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCACCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCCTTTCATCTGTTCGTACAGAATGCGGCTGCGTACTTTACGCAG  >  NC_000913/4584050‑4584331

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