New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 458531736 (0.990)4 (0.120) 4/226 2.8 9.5% coding (1445/3513 nt) hsdR endonuclease R Type I restriction enzyme
?NC_000913 = 4585344 42 (1.230)coding (1418/3513 nt) hsdR endonuclease R Type I restriction enzyme
Rejected: Frequency below/above cutoff threshold.

TCGCCGTTAATACGCGTATCTTCAAACGCGCCCAGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCC  >  NC_000913/4585197‑4585492
                                                                                                                                                   |                                                                                                                                                    
tCGCCGTTAATACGCGTATCTTCAAACGCGCCCAGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGt                                                                                                                                                     >  8:118523‑M2/1‑148 (MQ=255)
tCGCCGTTAATACGCGTATCTTCAAACGCGCCCAGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGt                                                                                                                                                     <  6:48277‑M2/149‑2 (MQ=255)
            cgcgTATCTTCAAACGCGCCCAGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGtcgccatcgcca                                                                                                                                          <  3:23572/148‑1 (MQ=255)
                      cAAACGCGCCCAGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGAt                                                                                                                                   <  2:103760/145‑1 (MQ=255)
                      cAAACGCGCCCAGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGAt                                                                                                                                   <  4:51053/145‑1 (MQ=255)
                         aCGCGCCCAGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCttgt                                                                                                                            >  1:54173/1‑149 (MQ=255)
                             gcCCAGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGcc                                                                                                                        >  4:3938/1‑149 (MQ=255)
                               ccAGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGcccc                                                                                                                      >  1:173044/1‑149 (MQ=255)
                                            gCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCTCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGcc                                                                                                          <  3:3938/148‑1 (MQ=255)
                                              cAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGa                                                                                                               <  8:11008/141‑1 (MQ=255)
                                              cAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGa                                                                                                               >  7:11008/1‑141 (MQ=255)
                                                            tCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGAGCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGAc                                                                                         <  2:54173/149‑1 (MQ=255)
                                                                   ggAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCt                                                                                  >  1:20532/1‑149 (MQ=255)
                                                                         gaATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCttgtt                                                                                                                           <  1:178274/102‑1 (MQ=255)
                                                                         gaATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCttgtt                                                                                                                           >  2:178274/1‑102 (MQ=255)
                                                                         gaATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCttg                                                                                                                             >  8:37957/1‑100 (MQ=255)
                                                                         gaATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCttg                                                                                                                             <  7:37957/100‑1 (MQ=255)
                                                                          aaTGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGtcgccatc                                                                                                                                              <  8:32787/82‑1 (MQ=255)
                                                                          aaTGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGtcgccatc                                                                                                                                              >  7:32787/1‑82 (MQ=255)
                                                                              cGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGataata                                                                       <  2:20532/149‑1 (MQ=255)
                                                                                    aaaCGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAg                                                                 <  8:123982/149‑1 (MQ=255)
                                                                                    aaaCGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAg                                                                 <  2:173044/149‑1 (MQ=255)
                                                                                          tGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGa                                                                                                                  >  8:18954/1‑94 (MQ=255)
                                                                                          tGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGa                                                                                                                  <  7:18954/94‑1 (MQ=255)
                                                                                                  gATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGcccc                                                                                                                      >  4:13698/1‑82 (MQ=255)
                                                                                                  gATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGcccc                                                                                                                      <  3:13698/82‑1 (MQ=255)
                                                                                                          ggAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAggg                                           <  5:48277/149‑1 (MQ=255)
                                                                                                                tcatGGCGATTGCCGTACGGGTTTTACCGGTACCGGtcgccatgatgttccgcctgatccagtcccagcgttttaaacgcattctcttccttgtcgaccgccgttctcttggc                                                                         >  4:40783‑M2/1‑36 (MQ=255)
                                                                                                                tcatGGCGATTGCCGTACGGGTTTTACCGGTACCGGtcgccatgatgttccgcctgatccagtcccagcgttttaaacgcattctcttccttgtcgaccgccgttctcttggc                                                                         <  3:40783‑M2/113‑78 (MQ=255)
                                                                                                                   tGGCGATTGCCGTACGGGTTTTACCGGTACCGGtcgccatgatgttccgcctgatccagtcccagcgttttaaacgcattctcttccttgtcgaccgccgttctcttggcgaacag                                                                   >  1:29871‑M2/1‑33 (MQ=255)
                                                                                                                   tGGCGATTGCCGTACGGGTTTTACCGGTACCGGtcgccatgatgttccgcctgatccagtcccagcgttttaaacgcattctcttccttgtcgaccgccgttctcttggcgaacag                                                                   <  2:29871‑M2/116‑84 (MQ=255)
                                                                                                                   tGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACg                                                                    <  6:29351/115‑1 (MQ=255)
                                                                                                                   tGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACg                                                                    >  5:29351/1‑115 (MQ=255)
                                                                                                                   tGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGAGTGCCTTTTCAACCGCGCGCACGGCATCGTCCTGATAATAACGCAGGCCCCGCTCGCACCTGCCAGGTTTATCGGc                                   >  8:119146/1‑148 (MQ=255)
                                                                                                                        aTTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTc                                                                                <  4:56855/98‑1 (MQ=255)
                                                                                                                        aTTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTc                                                                                >  3:56855/1‑98 (MQ=255)
                                                                                                                         ttGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCAc                            >  2:134083/1‑149 (MQ=255)
                                                                                                                         ttGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCAc                            >  2:136910/1‑149 (MQ=255)
                                                                                                                             cGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGAt                                                                                                                                   >  8:75706/1‑42 (MQ=255)
                                                                                                                             cGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGAt                                                                                                                                   <  7:75706/42‑1 (MQ=255)
                                                                                                                                   ggttggaccggtaccggtCGCCATCGCCAGCAGGATCTTGTGTTGCCCCTGGACGAGTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGtt                  <  7:119146/143‑1 (MQ=255)
                                                                                                                                    gTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCttgtt                                                                                                                           <  7:52244/43‑1 (MQ=255)
                                                                                                                                    gTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCttgtt                                                                                                                           >  8:52244/1‑43 (MQ=255)
                                                                                                                                         accggtaccggtCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGata                                                                          <  4:91954/87‑1 (MQ=255)
                                                                                                                                         accggtaccggtCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATGGCCTTTTCAACCGCGCGGACGGCATCTTCCTGata                                                                          >  3:91954/1‑87 (MQ=255)
                                                                                                                                                ccggtCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGa                                                                            >  8:115172/1‑78 (MQ=255)
                                                                                                                                                ccggtCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGa                                                                            <  7:115172/78‑1 (MQ=255)
                                                                                                                                                   gtCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGcc  >  5:37236/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TCGCCGTTAATACGCGTATCTTCAAACGCGCCCAGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCC  >  NC_000913/4585197‑4585492

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.