New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 5397806	52 (1.230)	3 (0.080)	3/240	NT	5.9%	coding (818/1731 nt)	hsdS	type I restriction modification system specificity protein
	?	NC_002947	= 5397858	47 (1.200)					coding (870/1731 nt)	hsdS	type I restriction modification system specificity protein

ACACCGCATCGTTGCCAAAGTCGAAGAACTAATGGCCCTATGCGACCGCTTGGAAGCCCAGCAGGCCGACGCCGAGTCCGCCCACGTCCAACTGGTGCAGGCGATGCTCGACAGCCTGACCCAAGCCATCGACGCCGCCGACTTCGCGACCAGTTGGCAGCGCCTAGCCGAGCACTTTCACACCCTGTTCACAAACGAGTTTGCCATCGACGCCCTCAAGAAAACCCTCTTGCAACTTGCCGTGATGGGAAAGCTTGTCCCGCAAGACGTTACTGACGAATCCGCCAGCGAA  >  NC_002947/5397711‑5398002                                                                                                                                                    |                                                                                                                                                 acacCGCTTCGTTTCCAATGTCGAAGAAGTAATGTGCCTATGCGACCGGTTGGAAGCGAAGCAGGCCGACGCGTAGTGCGCCGACGTACATCTGTTGCAGGCGATGCTCGACAGCCTGACGCAAGCCATCGAGGCCGCGGACTTcgcgt                                                                                                                                                 <  2:52718/149‑2 (MQ=255)            ttGCCAAAGTCGAAGAACTAATGGCCCTATGCGACCGCTTGGAAGCCCAGCAGGCCGACGCCGAGTCCGCCCACGTCCAACTGGTGCAGGCGATGCTCGACAGCCTGACCCAAGCCATCGACGCCGCCGACTTCGGGACCAGTTGGCAg                                                                                                                                      >  5:50481/1‑149 (MQ=255)               ccAAAGTCGAAGAACTAATGGCCCTATGCGACCGCTTGGAAGCCCAGCAGGCCGACGCCGAGTCCGCCCACGTCCAACTGGTGCAGGCGATGCTCGACAGCCTGACCCAAGCCATCGACGCCGCCGACTTCGCGACCAGTTGGCAgcg                                                                                                                                    <  4:190777/148‑1 (MQ=255)                 aaaGTCGAAGAACTAATGGCCCTATGCGACCGCTTGGAAGCCCAGCAGGCCGACGCCGAGTCAGACCACGTCCAAATGGTGCAGGCGAAGCTCGACAGCCTGACCCAAGCCATCGACGCCGCCGACTTCGCGACCAGTTGGCAGCGcc                                                                                                                                  >  5:74705/1‑148 (MQ=255)                   aGTCGAAGAACTAATGGCCCTATGCGACCGCTTGGAAGCCCAGCAGGCCGACGCCGAGTCCGCCCACGTCCAACTGGTGCAGGCGATGCTCGACAGCCTGACCCAAGCCATCGACGCCGCCGACTTCGCGACCAGTTGGCAGCGCCTAg                                                                                                                               <  2:211192/149‑1 (MQ=255)                       gaagaaCTAATGGCCCTATGCGACCGCTTGGAAGCCCAGCAGGCCGACGCCGAGTCCGCCCACGTCCAACTGGTGCAGGCGATGCTCGACAGCCTGACCCAAGCCATCGACGCCGCCGACTTCGCGACCAGTTGGCAGCGCCTAGcc                                                                                                                             >  8:13735/1‑147 (MQ=255)                       gaagaaCTAATGGCCCTATGCGACCGCTTGGAAGCCCAGCAGGCCGACGCCGAGTCCGCCCACGTCCAACTGGTGCAGGCGATGCTCGACAGCCTGACCCAAGCCATCGACGCCGCCGACTTCGCGACCAGTTGGCAGCGCCTAGCCGa                                                                                                                           <  7:177847/149‑1 (MQ=255)                        aagaaCTAATGGCCCTATGCGACCGCTTGGAAGCCCAGCAGGCCGACGCCGAGTCCGCCCACGTCCAACTGGTGCAGGCGATGCTCGACAGCCTGACCCAAGCCATCGACGCCGCCGACTTCGCGACCAGTTGGCAGCGCCTAGCCGAg                                                                                                                          >  5:115934/1‑149 (MQ=255)                        aagaaCTAATGGCCCTATGCGACCGCTTGGAAGCCCAGCAGGCCGACGCCGAGTCCGCCCACGTCCAACTGGTGCAGGCGATGCTCGACAGCCTGACCCAAGCCATCGACGCCGCCGACTTCGCGACCAGTTGGCAGCGCCTAGCCGAg                                                                                                                          <  7:13735/149‑1 (MQ=255)                         agaaCTAATGGCCCTATGCGACCGCTTGGAAGCCCAGCAGGCCGACGCCGAGTCCGCCCACGTCCAACTGGTGCAGGCGATGCTCGACAGCCTGACCCAAGCCATCGACGCCGCCGACTTCGCGACCAGTTGGCAGCGCCTAGCCGAGc                                                                                                                         <  3:48444/149‑1 (MQ=255)                          gaaCTAATGGCCCTATGCGACCGCTTGGAAGCCCAGCAGGCCGACGCCGAGTCCGCCCACGTCCAACTGGTGCAGGCGATGCTCGACAGCCTGACCCAAGCCATCGACGCCGCCGACTTCGGGACCAGTTGGCAGCGCCTAGCCGAGCa                                                                                                                        <  6:50481/149‑1 (MQ=255)                            aCTAATGGCCCTATGCGACCGCTTGGAAGCCCAGCAGGCCGACGCCGAGTCCGCCCACGTCCAACTGGTGCAGGCGATGCTCGACAGCCTGACCCAAGCCATCGACGCCGCCGACTTCGCGACCAGTTGGCAGCGCCTAGCCGAGCAc                                                                                                                       <  6:115934/148‑1 (MQ=255)                              tAATGGCCCTATGCGACCGCTTGGAAGCCCAGCAGGCCGACGCCGAGTCCGCCCACGTCCAACTGGTGCAGGCGATGCTCGACAGCCTGACCCAAGCCATCGACGCCGCCGACTTCGCGACCAGTTGGCAgcgc                                                                                                                                   <  6:203221/134‑1 (MQ=255)                              tAATGGCCCTATGCGACCGCTTGGAAGCCCAGCAGGCCGACGCCGAGTCCGCCCACGTCCAACTGGTGCAGGCGATGCTCGACAGCCTGACCCAAGCCATCGACGCCGCCGACTTCGCGACCAGTTGGCAgcgc                                                                                                                                   >  5:203221/1‑134 (MQ=255)                                aTGGCCCTATGCGACCGCTTGGAAGCCCAGCAGGCCGACGCCGAGTCCGCCCACGTCCAACTGGTGCAGGCGATGCTCGACAGCCTGACCCAAGCCATCGACGCCGCCGACTTCGCGACCAGTTGGCAGCGCCTAGCCGAGCACttt                                                                                                                    <  1:211151/147‑1 (MQ=255)                                       tATGCGACCGCTTGGAAGCCCAGCAGGCCGACGCCGAGTCCGCCCACGTCCAACTGGTGCAGGCGATGCTCGACCGACTGCCCCAAGCCATCGCCGCCGCCGACTGCGCGACCAGTTGGCAGCTCCTAGCCGCGCACTTTCACCCCCTg                                                                                                           >  3:156440/1‑149 (MQ=255)                                            gACCGCTTGGAAGCCCAGCAGGCCGAAGCCGAGTCCGCCCACGTCCAACTGGTGCAGGCGATGCTCGACCGCCTGACCCAAGCCATCGACGCCGCCGACTTCGAGACCAGTTGGCAGCGCCTAGCCGAGCACTTTCACACCCTGTTca                                                                                                       >  7:45318/1‑148 (MQ=255)                                               cGCTTGGAAGCCCAGCAGGCCGACGCCGAGTCCGCCCACGTCCAACTGGTGCAGGCGATGCTCGACAGCCTGACCCAAGCCATCGACGCCGCCGACTTCGCGACCAGTTGGCAGCGCCTAGCCGAGCACTTTCACACCCTGTTCACaaa                                                                                                   >  5:89188/1‑149 (MQ=255)                                                 cTTGGAAGCCCAGCAGGCCGACGCCGAGTCCGCCCACGTCCAACTGGTGCAGGCGATGCTCGACAGCCTGACCCAAGCCATCGACGCCGCCGACTTCGCGACCAGTTGGCAGCGCCTAGCCGAGCACTTTCACACCCTGTTCACAAAc                                                                                                  >  5:115101/1‑148 (MQ=255)                                                    ggAAGCCCAGCAGGCCGACGCCGAGTCCGCCCACGTCCAACTGGTGCAGGCGATGCTCGACAGCCTGACCCAAGCCATCGACGCCGCCGACTTCGCGACCAGTTGGCAGCGcc                                                                                                                                  >  8:75267/1‑113 (MQ=255)                                                    ggAAGCCCAGCAGGCCGACGCCGAGTCCGCCCACGTCCAACTGGTGCAGGCGATGCTCGACAGCCTGACCCAAGCCATCGACGCCGCCGACTTCGCGACCAGTTGGCAGCGcc                                                                                                                                  <  7:75267/113‑1 (MQ=255)                                                    ggAAGCCCAGCAGGCCGACGCCGAGTCCGCCCACGTCCAACTGGTGCAGGCGATGCTCGACAGCCTGACCCAAGCCATCGACGCCGCCGACTTCGCGACCAGTTGGCAGCGCCTAGCCGAGCACTTTCACACCCTGTTCACAAACGAGt                                                                                              >  1:35064/1‑149 (MQ=255)                                                       aGCCCTGCAGGCCGACGCCGAGTCCGCCCACGTCCAAGTGGTGCAGGCGATGCTCGACAGCCTGACCCAAGCCATCGACGCCGCCGACTTCGCGACCAGTTGGCAGCGCCTAGCCGAGCACTTTCACACCCTGTTCACAAACGAGTTTg                                                                                           <  5:218861/149‑1 (MQ=255)                                                             gcagGCCGACGCCGAGTCCGCCCACGTCCAACTGGTGCAGGCGATGCTCGACAGCCTGACCCAAGCCATCGACGCCGCCGACTTCGCGACCAGTTGGCAGCGCCTAGCCGAGCACTTTCACACCCTGt                                                                                                          <  1:265728/128‑1 (MQ=255)                                                             gcagGCCGACGCCGAGTCCGCCCACGTCCAACTGGTGCAGGCGATGCTCGACAGCCTGACCCAAGCCATCGACGCCGCCGACTTCGCGACCAGTTGGCAGCGCCTAGCCGAGCACTTTCACACCCTGt                                                                                                          >  2:265728/1‑128 (MQ=255)                                                                  ccGACGCCGAGTCCGCGCACGTCCAAATGGTGCAGGCGATGCTCGACAGCCTGACCCAAGCCATCGACGCCGCCGACTTCGCGACCAGTTGGCAGCGCCTAGCCGAGCACTTTCACACCCTGTTCACAAACGAGTTTGCCATCGACGcc                                                                                <  8:45318/149‑1 (MQ=255)                                                                     aCGCCGAGTCCGCCCACGTCCAACTGGTGCAGGCGATGCTCGACAGCCTGACCCAAGCCATCGACGCCGCCGACTTCGCGACCAGTTGGCAGCGCCTAGCCGAGCACTTTCACACCCTGTTCACAAACGAGTTTGCCATCGACGCCCTc                                                                             <  6:115101/149‑1 (MQ=255)                                                                               cGCCCACGTCCAACTGGTGCAGGCGATGCTCGACAGCCTGACCCAAGCCATCGACGCCGCCGACTTCGCGACCAGTTGGCAGCGCCTAGCCGAGCACTTTCACACCCTGTTCACAAACGAGTTTGCCATCGACGCCCTCAAGAAAAccc                                                                   >  6:273043/1‑149 (MQ=255)                                                                               cGCCCACGTCCAACTGGTGCAGGCGATGCTCGACAGCCTGACCCAAGCCATCGACGCCGCCGACTTCGCGACCAGTTGGCAGCGCCTAGCCGAGCACTTTCACACCCTGTTCACAAACGAGTTTGCCATCGACGCCCTCAAGAAAAccc                                                                   >  3:214335/1‑149 (MQ=255)                                                                                      gTCCAACTGGTGCAGGCGATGCTCGACAGCCTGACCCAAGCCATCGACGCCGCCGACTTCGCGACCAGTTGGCAGCGCCTAGCCGAGCACTTTCACACCCTGTTCACAAACGAGTTTGCCATCGACGCCCTCAAGAAAACCCTCTTGc                                                             >  7:106033/1‑148 (MQ=255)                                                                                        ccAACTGGTGCAGGCGATGCTCGACAGCCTGACCCAAGCCATCGACGCCGCCGACTTCGCGaccagttggacgtgggcggactcggcgtcggcctgctgggcttccaagcggtcgcatagggccattagttcttcgactttggcaacga                                                          <  1:241861‑M2/149‑89 (MQ=255)                                                                                         cAACTGGTGCAGGCGATGCTCGACAGCCTGACCCAAGCCATCGACGCCGCCGACTTCGCGACCAGTTGGCAGCGCCTAGCCGAGCACTTTCACACCCTGTTCACAAACGAGTTTGCCATCGACGCCCTCAAGAAAACCCTCTTGCAACt                                                         >  3:14932/1‑149 (MQ=255)                                                                                          aaCTGGTGCAGGCGATGCTCGACAGCCTGACCCAAGCCATCGACGCCGCCGACTTCGCGACCAGTTGGCAGCGCCTAGCCGAGCACTTTCACACCCTGTTCACAAACGAGTTTGCTATCGACGCCCTCAAGAAAACCCTCTTGCAACtt                                                        >  4:5672/1‑149 (MQ=255)                                                                                           aCTGGGGCATGCGATGCTCGACCGCCTGACCCACGCCCTCGACGCCGCCGACTTCGAGACCAGTGGGCAGCGCCTAGCCGAGCACTTTCACACCCTGTTCACAAACGAGTTTGCCATCGACGCCCTCAAGAAAACCCTCTTGCAACttt                                                       <  4:156440/149‑2 (MQ=255)                                                                                               gTGCAGGCGATGCTCGACAGCCTGACCCAAGCCATCGACGCCGCCGACTTCGCGaccagtttgacgtgggcggactcggcgtcggcctgctgggcttccaagcggtcgcgtagggccattagttcttcgactttggc                                                               >  2:166076‑M2/1‑54 (MQ=255)                                                                                               gTGCAGGCGATGCTCGACAGCCTGACCCAAGCCATCGACGCCGCCGACTTCGCGaccagttggacgtgggcggactcggcgtcggcctgctgggcttccaagcggtcgcgtagggccattagttcttcgactttggc                                                               <  1:166076‑M2/137‑84 (MQ=255)                                                                                               gTGCAGGCGATGCTCGACAGCCTGACCCAAGCCATCGACGCCGCCGACTTCGCGACCAGTTGGCAGCGCCTAGCCGAGCACTTTCACACCCTGTTCACAAACGAGTTTGCCATCGACGCCCTCAAGAAAACCCTCTTGCAACTTGcc                                                     <  5:273043/147‑1 (MQ=255)                                                                                                  cAGGCGATGCTCGACAGCCTGACCCAAGCCATCGACGCCGCCGACTTCGCGACCAGTTGGCAGCGCCTAGCCGAGCACTTTCACACCCTGTTCACAAACGAGTTTGCCATCGACGCCCTCAAGAAAACCCTCTTGCAACTTGCCGTGAt                                                <  2:35064/149‑1 (MQ=255)                                                                                                              gACAGCCTGACCCAAGCCATCGACGCCGCCGACTTCGCGACCAGTTGGCAGCGCCTAGCCGAGCACTTTCACACCCTGTTCACAAACGAGTTTGCCATCGACGCCCTCAAGAAAACCCTCTTGCAACTTGCCGTGATGGGAAAGCtt                                      <  6:74705/147‑1 (MQ=255)                                                                                                                 aGCCTGACCCAAGCCATCGACGCCGCCGACTTCGCGACCAGTTGGCAGCGCCTAGCCGAGCACTTTCACACCCTGTTCACAAACGAGTTTGCCATCGACGCCCTCAAGAAAACCCTCTTGCAACTTGCCGTGATGGGAAAGCTTGTccc                                 <  4:14932/149‑1 (MQ=255)                                                                                                                   ccTGACCCAAGCCATCGACGCCGCCGACTTCGCGACCAGTTGGCAGCGCCTAGCCGAGCACTTTCACACCCTGTTCACAAACGAGTTTGCCATCGACGCCCTCAAGAAAACCCTCTTGCAACTTGCCGTGATGGGAAAGCTTGTCCCg                                >  1:9570/1‑148 (MQ=255)                                                                                                                   ccTGACCCAAGCCATCGACGCCGCCGACTTCGCGACCAGTTGGCAGCGCCTAGCCGAGCACTTTCACACCCTGTTCACAAACGAGTTTGCCATCGACGCCCTCAAGAAAACCCTCTTGCAACTTGCCGTGATGGGAAAGCTTGTCCCg                                >  4:261830/1‑148 (MQ=255)                                                                                                                   ccTGACCCAAGCCATCGACGCCGCCGACTTCGCGACCAGTTGGCAGCGCCTAGCCGAGCACTTTCACACCCTGTTCACAAACGAGTTTGCCATCGACGCCCTCAAGAAAACCATCTTGCAACTTGCCGTGATGGGAAAGCTTGTCCCg                                >  3:95603/1‑148 (MQ=255)                                                                                                                     tGACCCAAGCCATCGACGCCGCCGACTTCGCGACCAGTTGGCAGCGCCTAGCCGAGCACTTTCACACCCTGTTCACAAACGAGTTTGCCATCGACGCCCTCAAGAAAACCCTCTTGCAACTTGCCGTGATGGGAAAGCTTGTCCCGc                               >  2:174244/1‑147 (MQ=255)                                                                                                                       aCCCAAGCCATCGACGCCGCCGACTTCGCGACCAGTTGGCAGCGCCTAGCCGAGCACTTTCACACCCTGTTCACAAACGAGTTTGCCATCGACGCCCTCAAGAAAACCCTCTTGCAACTTGCCGTGATGGGAAAGCTTGTCCCGCAAGa                           <  3:210819/149‑1 (MQ=255)                                                                                                                       aCCCAAGCCATCGACGCCGCCGACTTCGCGACCAGTTGGCAGCGCCTAGCCGAGCACTTTCACACCCTGTTCACAAACGAGTTTGCCATCGACGCCCTCAAGAAAACCCTCTTGCAACTTGCCGTGATGGGAAAGCTTGTCCCGCAAGa                           >  4:210819/1‑149 (MQ=255)                                                                                                                           aaGCCATCGACGCCGCCGACTTCGCGACCAGTTGGCAGCGCCTAGCCGAGCACTTTCACACCCTGTTCACAAACGAGTTTGCCATCGACGCCCTCAAGAAAACCCTCTTGc                                                             >  4:76830/1‑111 (MQ=255)                                                                                                                           aaGCCATCGACGCCGCCGACTTCGCGACCAGTTGGCAGCGCCTAGCCGAGCACTTTCACACCCTGTTCACAAACGAGTTTGCCATCGACGCCCTCAAGAAAACCCTCTTGc                                                             <  3:76830/111‑1 (MQ=255)                                                                                                                            aGCCATCGACGCCGCCGACTTCGCGACCAGTTGGCAGCGCCTAGCCGAGCACTTTCACACCCTGTTCACAAACGAGTTTGCCATCGACGCCCTCAAGAAAACCCTCTTGCAACTTGCCGTGATGGGAAAGCTTGTCCCGCAAGACGTTa                      <  4:214335/149‑1 (MQ=255)                                                                                                                                aTCGACGCCGCCGACTTCGCGACCAGTTGGCAGCGCCTAGCCGAGCACTTTCACACCCTGTTCACAAACGAGTTTGCTATCGACGCCCTCAAGAAAACCCTCTTGCAACTTGCCGTGATGGGAAAGCTTGTCCCGCAAGACGTTACTGa                  <  3:5672/149‑1 (MQ=255)                                                                                                                                        cgccgACTTCGCGACCAGTTGGCAGCGCCTAGCCGAGCACTTTCACACCCTGTTCACAAACGAGTTTGCCATCGACGCCCTCAAGAAAACCCTCTTGCAACTTGCCGTGATGGGAAAGCTTGTCCCGCAAGACGTTACTGACGAATCCg          >  2:6527/1‑149 (MQ=255)                                                                                                                                               ttCGCGACCAGTTGGCAGCGCCTAGCCGAGCACTTTCACACCCTGTTCACAAACGAGTTTGCCATCGACGCCCTCAAGAAAACCCTCTTGCAACTTGCCGTGATGGGAAAGCTTGTCCCGCAAGACGTTACTGACGAATcc           <  6:53657/141‑1 (MQ=255)                                                                                                                                               ttCGCGACCAGTTGGCAGCGCCTAGCCGAGCACTTTCACACCCTGTTCACAAACGAGTTTGCCATCGACGCCCTCAAGAAAACCCTCTTGCAACTTGCCGTGATGGGAAAGCTTGTCCCGCAAGACGTTACTGACGAATcc           >  5:53657/1‑141 (MQ=255)                                                                                                                                                 cgcgACCAGTTGGCAGCGCCTAGCCGAGCACTTTCACACCCTGTTCACAAACGAGTTTGCCATCGACGCCCTCAAGAAAACCCTCTTGCAACTTGCCGTGATGGGAAAGCTTGTCCCGCAAGACGTTACTGACGAATCCGCCAGCGaa  <  5:213326/148‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                 ACACCGCATCGTTGCCAAAGTCGAAGAACTAATGGCCCTATGCGACCGCTTGGAAGCCCAGCAGGCCGACGCCGAGTCCGCCCACGTCCAACTGGTGCAGGCGATGCTCGACAGCCTGACCCAAGCCATCGACGCCGCCGACTTCGCGACCAGTTGGCAGCGCCTAGCCGAGCACTTTCACACCCTGTTCACAAACGAGTTTGCCATCGACGCCCTCAAGAAAACCCTCTTGCAACTTGCCGTGATGGGAAAGCTTGTCCCGCAAGACGTTACTGACGAATCCGCCAGCGAA  >  NC_002947/5397711‑5398002

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