New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 5395374	32 (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

CGTAACTACTTCGCCAAATACTCCGGTGCAGCGCGCCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATCAAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAAGCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGCAACGGCGGAGCCTTGTCCACCTAAGCAGCCAACTCCATTACAATCCCC  >  NC_002947/5395246‑5395541                                                                                                                                                   |                                                                                                                                                      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:284369/150‑1 (MQ=255)                            gCAGCGCGCCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATCAAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAAGCTGGc                                                                                                                         <  2:623989/150‑1 (MQ=255)                               gcgcgcCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATCAAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTggcgg                                                                                                                                   <  1:452789/137‑1 (MQ=255)                               gcgcgcCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATCAAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTggcgg                                                                                                                                   >  2:452789/1‑137 (MQ=255)                                cgcgcCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATCAAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAAGCTGGCTaca                                                                                                                     >  1:732132/1‑150 (MQ=255)                                  cgcCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATCAAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAAGCTGGCTacaac                                                                                                                   <  2:844649/150‑1 (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                                                     >  2:371050/1‑149 (MQ=255)                                                                                                 aTCAAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAAGCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTtgggtg                                                     <  1:371050/149‑1 (MQ=255)                                                                                                         cctccAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAAGCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGCAAcggc                                            <  1:736365/150‑1 (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)                                                                                                                             aCCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAAGCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGCAACGGCGGAGCCTTGTCCACCTAagc                        >  1:570162/1‑150 (MQ=255)                                                                                                                                  cTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAAGCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTGGGGTGGGCAACGGCGGAGCCTTGTCCACCTAAGCAGCCa                   >  1:354810/1‑150 (MQ=255)                                                                                                                                        gCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAAGCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGCAACGGCGGAGCCTTGTCCACCTAAGCAGCCAACTCCa             <  2:38716/150‑1 (MQ=255)                                                                                                                                                  aaCTGGTTAAAGCCTTTGGCGGCAAAGCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGCAACGGCGGAGCCTTGTCCACCTAAGCAGCCAACTCCATTACAATccc   >  2:346554/1‑150 (MQ=255)                                                                                                                                                   aCTGGTTAAAGCCTTTGGCGGCCACGACGGCCACCCCAACTCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGCAACGGCGGAGCCTTGTCCACCTAAGCAGCCAACTCCATTACAATcccc  <  2:354810/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                      CGTAACTACTTCGCCAAATACTCCGGTGCAGCGCGCCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATCAAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAAGCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGCAACGGCGGAGCCTTGTCCACCTAAGCAGCCAACTCCATTACAATCCCC  >  NC_002947/5395246‑5395541

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.