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.

TCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCAGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCA  >  NC_000913/4585176‑4585465
                                                                                                                                             |                                                                                                                                                    
tCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCAGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCg                                                                                                                                                >  7:123982/1‑148 (MQ=255)
    aaaTGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCAGCGCCTGTTCGCCAAGAGAACGGCGGGCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACgg                                                                                                                                           <  3:35754/149‑1 (MQ=255)
            ttGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCAGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTtac                                                                                                                                    <  8:75398/148‑1 (MQ=255)
                aGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCAGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTtaccggta                                                                                                                               <  7:157351/149‑1 (MQ=255)
                     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)
                                                     cAGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGtttac                                                                                                                                     >  1:83131/1‑104 (MQ=255)
                                                     cAGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGtttac                                                                                                                                     <  2:83131/106‑3 (MQ=255)
                                                           cTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGAccggtaccggtaaaacccgtacggcaatcgcca                                                                                                                     >  1:29871‑M1/1‑83 (MQ=255)
                                                           cTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGAccggtaccggtaaaacccgtacggcaatcgcca                                                                                                                     <  2:29871‑M1/116‑34 (MQ=255)
                                                                 gCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGAccggtaccggtaaaacccgtacggcaatcgccatga                                                                                                                  <  3:40783‑M1/113‑37 (MQ=255)
                                                                 gCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGAccggtaccggtaaaacccgtacggcaatcgccatga                                                                                                                  >  4:40783‑M1/1‑77 (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                                                                                       <  7:18954/94‑1 (MQ=255)
                                                                                                               tGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGa                                                                                       >  8:18954/1‑94 (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)
                                                                                                                                             accggtaccggtaaaacccgtacggcaatcgccatgatgttccgcctgatccagtcccagcgttttaaacgcattctcttccttgtcgaccgccgttctcttggcgaacaggcgctgggcgcgtttgaagatacgcgtattaacggcga  <  6:48277‑M1/149‑149 (MQ=255)
                                                                                                                                             accggtaccggtaaaacccgtacggcaatcgccatgatgttccgcctgatccagtcccagcgttttaaacgcattctcttccttgtcgaccgccgttctcttggcgaacaggcgctgggcgcgtttgaagatacgcgtattaacggcga  >  8:118523‑M1/1‑1 (MQ=255)
                                                                                                                                             aTTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTc                                                     >  3:56855/1‑98 (MQ=255)
                                                                                                                                             aTTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTc                                                     <  4:56855/98‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
TCGAAAATGCTGTTGAAGGTGTCGCCGTTAATACGCGTATCTTCAAACGCGCCCAGCGCCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCA  >  NC_000913/4585176‑4585465

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.