New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 = 539537432 (1.010)3 (0.100) 3/228 NT 8.8% coding (2268/2364 nt) hsdR type I restriction modification system endonuclease
?NC_002947 = 5395392 33 (1.130)coding (2286/2364 nt) hsdR type I restriction modification system endonuclease

AGTGAAAAAGCGTAACTACTTCGCCAAATACTCCGGTGCAGCGCGCCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATCAAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAAGCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGCAACGGCGGAGCCTTGTCCACCTAAGCA  >  NC_002947/5395236‑5395520
                                                                                                                                          |                                                                                                                                                  
aGTGAAAAAGCGTAACTACTTCGCCAAATACTCCGGTGCAGCGCGCCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACCTCGAGGACATCAAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCAcc                                                                                                                                         <  2:631927/150‑1 (MQ=255)
          cgcaaCTACTTCGCCAAATACTCCGGTGCAGCGCGCCATCTGCTGGAAGCACTGCTGGATAAATACGCCGACACCCGCATTGAGCACATCGAGGACATCAAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTg                                                                                                                               <  2:307822/147‑1 (MQ=255)
           gTAACTACTTCGCCAAATACTCCGGTGCAGCGCGCCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATCAAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAActgg                                                                                                                              <  2:683154‑M2/150‑5 (MQ=255)
                          aaTACTCCGGTGCAGCGCGCCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATCAAAATCCTCCAGCTCGATCCCTTTAACCAGTTCGATGGGTGCACCGAGctggttaaagggatcgagc                                                                                                               >  2:369609‑M2/1‑131 (MQ=255)
                                  ggTGCAGCGCGCCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATCAAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAAGCt                                                                                                       <  1:628143/150‑1 (MQ=255)
                                   gTGCAGCGCGCCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATCAAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAActggttaaa                                                                                                                         >  1:111148‑M2/1‑122 (MQ=255)
                                   gTGCAGCGCGCCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATCAAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAActggttaaa                                                                                                                         <  2:111148‑M2/131‑10 (MQ=255)
                                   gTGCAGCGCGCCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATCAAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGACAAAGCTg                                                                                                      >  1:752895/1‑150 (MQ=255)
                                    tGCAGCGCGCCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATCAAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGACAAAGCTgg                                                                                                     <  2:752895/150‑1 (MQ=255)
                                     gCAGCGCGCCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATCAAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAAGCTGGc                                                                                                    <  2:623989/150‑1 (MQ=255)
                                     gCAGCGCGCCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATCAAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAAGCTGGc                                                                                                    <  2:284369/150‑1 (MQ=255)
                                        gcgcgcCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATCAAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTggcgg                                                                                                              >  2:452789/1‑137 (MQ=255)
                                        gcgcgcCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATCAAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTggcgg                                                                                                              <  1:452789/137‑1 (MQ=255)
                                         cgcgcCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATCAAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAAGCTGGCTaca                                                                                                >  1:732132/1‑150 (MQ=255)
                                           cgcCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATCAAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAAGCTGGCTacaac                                                                                              <  2:844649/150‑1 (MQ=255)
                                              cAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATCAAAATCCTCCAGCTCGATCCCTTTAACCAGttcgatgggtgcaccgagctggttaaaggga                                                                                                                     >  2:712027‑M1/1‑93 (MQ=255)
                                              cAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATCAAAATCCTCCAGCTCGATCCCTTTAACCAGttcgatgggtgcaccgagctggttaaaggga                                                                                                                     <  1:712027‑M1/124‑32 (MQ=255)
                                              cAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATCAAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAAGCTGGCTACAACaaa                                                                                           >  2:118397/1‑150 (MQ=255)
                                                       gAAGTGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATCAAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAAGCTGGCTACAACAAAGCTATCCAt                                                                                  <  2:708722/150‑1 (MQ=255)
                                                       gAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATCAAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAAGCTGGCTACAACAAAGCTATCCAt                                                                                  <  1:779854/150‑1 (MQ=255)
                                                                  ggATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATCAAAATCCTCCGGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAAGCTGGCTACAACAAAGCTATCCATGAACTGGAAGa                                                                       >  2:337739/1‑150 (MQ=255)
                                                                                acacCGGCATTGAGCACATCGAGGACATCAAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAAGCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCa                                                         <  1:339904/150‑1 (MQ=255)
                                                                                acacCGGCATTGAGCACATCGAGGACATCAAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAAGCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCa                                                         <  2:744161/150‑1 (MQ=255)
                                                                                acacCGGCATTGAGCACATCGAGGACATCAAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAAGCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCa                                                         <  2:247610/150‑1 (MQ=255)
                                                                                    cGGCATTGAGCACATCGAGGACATCAAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAAGCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTa                                                     >  2:736365/1‑150 (MQ=255)
                                                                                    cGGCATTGAGCACATCGAGGACATCAAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAAGCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTa                                                     >  2:149747/1‑150 (MQ=255)
                                                                                            aGCACATCGAGGACATCAAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAAGCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGctc                                             >  1:771755/1‑150 (MQ=255)
                                                                                               acaTCGAGGACATCAAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAAGCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGtt                                                      >  2:292783/1‑138 (MQ=255)
                                                                                               acaTCGAGGACATCAAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAAGCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGtt                                                      <  1:292783/138‑1 (MQ=255)
                                                                                                     aGGACATCAAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAAGCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTtg                                    <  1:118397/150‑1 (MQ=255)
                                                                                                          aTCAAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAAGCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTtgggtg                                <  1:371050/149‑1 (MQ=255)
                                                                                                          aTCAAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAAGCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTtgggtg                                >  2:371050/1‑149 (MQ=255)
                                                                                                                  cctccAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAAGCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGCAAcggc                       <  1:736365/150‑1 (MQ=255)
                                                                                                                        gCTCGATCCCTTTAACCAGctcggtgcacccatcgaactggttaaagggatcgagctggaggattttgatgtcctcgatgtgctcaatgccggtgtcggcgtatttatccagcagcgcttccagcacctggcgcgctgcaccggagtatt                 >  2:369609‑M1/1‑19 (MQ=255)
                                                                                                                        gCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAAGCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGCAACGGCGGAGcc                 <  2:732132/150‑1 (MQ=255)
                                                                                                                          tCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAAGCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGCAACGGCGGAGCCtt               <  1:683154/150‑1 (MQ=255)
                                                                                                                            gATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAAGCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGCAACGGCGGAGCCTTGt             <  1:337739/150‑1 (MQ=255)
                                                                                                                              tCCCTTTAACCAGCTCGGTGCACCCATCGAActggttaaagggatcgagctggaggattttgatgtcctcgatgtgctcaatgccggtgtcggcgtatttatccagcagcgcttccagcacctg                                     <  1:712027‑M2/124‑94 (MQ=255)
                                                                                                                              tCCCTTTAACCAGCTCGGTGCACCCATCGAActggttaaagggatcgagctggaggattttgatgtcctcgatgtgctcaatgccggtgtcggcgtatttatccagcagcgcttccagcacctg                                     >  2:712027‑M2/1‑31 (MQ=255)
                                                                                                                              tCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAAGCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGCAACGGCGGAGCCTTGTcc           >  1:129770/1‑150 (MQ=255)
                                                                                                                                 cTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAAGCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGCAACGGCGGAGCCTTGTCCAcc        >  1:807744/1‑150 (MQ=255)
                                                                                                                                  tttAACCAGttcgatgggtgcaccgagctggttaaagggatcgagctggaggattttgatgtcctcgatgtgctcaatgccggtgtcggcgtatttatccagcagcgcttccagcacctggcgcgctgcac                          <  2:111148‑M1/131‑123 (MQ=255)
                                                                                                                                  tttAACCAGttcgatgggtgcaccgagctggttaaagggatcgagctggaggattttgatgtcctcgatgtgctcaatgccggtgtcggcgtatttatccagcagcgcttccagcacctggcgcgctgcac                          >  1:111148‑M1/1‑9 (MQ=255)
                                                                                                                                      aCCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAAGCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGCAACGGCGGAGCCTTGTCCACCTAagc   >  1:570162/1‑150 (MQ=255)
                                                                                                                                       ccAGttcgatgggtgcaccgagctggttaaagggatcgagctggaggattttgatgtcctcgatgtgctcaatgccggtgtcggcgtatttatccagcagcgcttccagcacctggcgcgctgcaccggagtatttggcgaagtagttac  <  2:683154‑M1/150‑147 (MQ=255)
                                                                                                                                          |                                                                                                                                                  
AGTGAAAAAGCGTAACTACTTCGCCAAATACTCCGGTGCAGCGCGCCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATCAAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAAGCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGCAACGGCGGAGCCTTGTCCACCTAAGCA  >  NC_002947/5395236‑5395520

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