New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	5599510 =	42 (0.860)	4 (0.090)	4/250	NT	8.4%	coding (3972/5418 nt)	PP_4924	subtilase family serine protease
	?	NC_002947	5599554 =	47 (1.010)					coding (3928/5418 nt)	PP_4924	subtilase family serine protease

CAAGCTTCCACTCCCCAGCCGAGCGCTCTGCGCGATGATGGGTCGACATGAACCCATATTTGAACGCCCCTGAACGCGAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACAAAGTCGCAACCTGCCCCGCCGCAACCTGTTGAGCCACTG  >  NC_002947/5599365‑5599658                                                                                                                                                  |                                                                                                                                                     cAAGCTTCCACTCCCCAGCCGAGCGCTCTGCGCGATGATGGGTCGACATGAACCCATATTTGAACGCCCCTGAACGCGAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTgccg                                                                                                                                                   >  8:109324/1‑149 (MQ=255) accaaggccaatgaacggctgctgacagccagcagtgaaccagtggctcaacaggttgcggcggggcaggttgcgactttgtctttgaccctgccggctgatttgctggtggagcatgtggaggttcatgtccactcgatggtcggccg                                                                                                                                                   <  3:8554‑M1/4‑1 (MQ=255)        ccACTCCCCAGCCGAGCGCTCTGCGCGATGATGGGTCGACATGAACCCATATTTGAACGCCCCTGAACGCGAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGc                                                                                                                                             <  4:188249/148‑1 (MQ=255)        ccACTCCCCAGCCGAGCGCTCTGCGCGATGATGGGTCGACATGAACCCATATTTGAACGCCCCTGAACGCGAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGc                                                                                                                                             >  7:153766/1‑148 (MQ=255)          aCTCCCCAGCCGAGCGCTCTGCGCGATGATGGGTCGACATGAACCCATATTTGAACGCCCCTGAACGCGAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGa                                                                                                                                           >  2:136450/1‑148 (MQ=255)           cTCCCCAGCCGAGCGCTCTGCGCGATGATGGGTCGACATGAACCCATATTTGAACGCCCCTGAACGCGAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGaca                                                                                                                                         >  2:104472/1‑149 (MQ=255)                  gCCGAGCGCTCTGCGCGATGATGGGTCGACATGAACCCATATTTGAACGCCCCTGAACGCGAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGt                                                                                                                                  <  2:268153/149‑1 (MQ=255)                    cGAGCGCTCTGCGCGATGATGGGTCGACATGAACCCATATTTGAACGCCCCTGAACGCGAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGa                                                                                                                                <  2:11882/149‑1 (MQ=255)                           tctGCGCGATGATGGGTCGACATGAACCCATATTTGAACGCCCCTGAACGCGAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCAt                                                                                                                         >  2:230673/1‑149 (MQ=255)                           tctGCGCGATGATGGGTCGACATGAACCCATATTTGAACGCCCCTGAACGCGAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCACGCGCCACCAGCGTGAGCGTCAt                                                                                                                         >  7:260078/1‑149 (MQ=255)                                    tgatgGGTCGACATCAACCCATATTTGAACGCCCCTGAACGCGAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGATCCCTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAg                                                                                                                <  1:230673/149‑1 (MQ=255)                                         gggcGACATGAACCCATATTTGAACGCCCCTGAACGCGAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAggcggc                                                                                                           <  5:207827/146‑1 (MQ=255)                                          gTCGACATGAACCCATATTTGAACGCCCCTGAACGCGAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAggcggc                                                                                                           >  6:207827/1‑148 (MQ=255)                                                  gAACCCATATTTGAACGCCCCTGAACGCGAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATc                                                                                                   >  5:80653/1‑148 (MQ=255)                                                   acaggttgcggcggggcaggttgcgactttgtctttgaccctgccggctgatttgctggtggagcatgtggaggttgatgtccactcgatggtcgGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAgg                                                                                                               >  1:131006‑M1/96‑135 (MQ=255)                                                   acaggttgcggcggggcaggttgcgactttgtctttgaccctgccggctgatttgctggtggagcatgtggaggttgatgtccactcgatggtcgGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAgg                                                                                                               <  2:131006‑M1/40‑1 (MQ=255)                                                     cccATATTTGAACGCCCCTGAACGCGAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGt                                                                                               >  1:286121/1‑149 (MQ=255)                                                           tttGAACGCCCCTGAACGCGAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACAt                                                                                         >  4:261530/1‑149 (MQ=255)                                                              gAACGCCCCTGAACGCGAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCa                                                                                       <  2:286121/148‑1 (MQ=255)                                                               aaCGCCCCTGAACGCGAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAggcggc                                                                                                           >  2:305587/1‑127 (MQ=255)                                                               aaCGCCCCTGAACGCGAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAggcggc                                                                                                           <  1:305587/127‑1 (MQ=255)                                                                 cGCCCCTGAACGCGAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCt                                                                                   <  3:261530/149‑1 (MQ=255)                                                                   gccgtGAAGGCGAATCGCCCAGATGGTCATCACCAGAGCCCGGCGCCTTTCCTGTGCGTTCGAGCAGTACGCTGCGCGTGGGGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGAGATCTACCTcc                                                                                 <  8:260078/145‑1 (MQ=255)                                                                    cccTGAACGCGAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGa                                                                                            <  6:28156/137‑1 (MQ=255)                                                                    cccTGAACGCGAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGa                                                                                            >  5:28156/1‑137 (MQ=255)                                                                         aaCGCGAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATg                                                                            <  6:80653/148‑1 (MQ=255)                                                                             cgAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTcca                                                                       <  7:109324/149‑1 (MQ=255)                                                                              gAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTccac                                                                      >  4:216587/1‑149 (MQ=255)                                                                              gAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCAGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTccac                                                                      <  5:261409/149‑1 (MQ=255)                                                                                aTCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTccacca                                                                    >  4:292762/1‑149 (MQ=255)                                                                                        gATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCa                                                            <  8:153766/149‑1 (MQ=255)                                                                                               cggctgatttgctggtggagcatgtggaggttgatgtccactcgatggtcgGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAggc                                                                                                              >  4:39275‑M1/52‑92 (MQ=255)                                                                                               cggctgatttgctggtggagcatgtggaggttgatgtccactcgatggtcgGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAggc                                                                                                              <  3:39275‑M1/41‑1 (MQ=255)                                                                                                tcaCCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGCCACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACAAGCAAAACAGCCGGCAg                                                    >  2:266817/1‑149 (MQ=255)                                                                                                  aCCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCTAATCAGCCGGCAgg                                                   <  1:266817/148‑1 (MQ=255)                                                                                                       agCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCaaa                                             <  1:136450/149‑1 (MQ=255)                                                                                                       agCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCaaa                                             >  7:35104/1‑149 (MQ=255)                                                                                                       agCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCaaa                                             <  7:287915/149‑1 (MQ=255)                                                                                                         cccGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTcc                                                                                 >  8:66218/1‑111 (MQ=255)                                                                                                         cccGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTcc                                                                                 <  7:66218/111‑1 (MQ=255)                                                                                                            ggCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACa                                         >  8:168818/1‑148 (MQ=255)                                                                                                                                gAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACAAAGTCGCAACCTGCCccgccg                    >  2:46568/1‑149 (MQ=255)                                                                                                                                 aGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACAAAGTCGCAACCTGCCCcgccgc                   >  4:225236/1‑149 (MQ=255)                                                                                                                                    aGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACAAAGTCGCAACCTGCCCCGCCGCAAc                >  3:301990/1‑149 (MQ=255)                                                                                                                                        cGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACAAAGTCGCAACCTGCCCCGCCGCAACCTGt            <  1:46568/149‑1 (MQ=255)                                                                                                                                         gCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACAAAGTCGCAACCTGCCCCGCCGCAACCTGtt           >  2:237899/1‑149 (MQ=255)                                                                                                                                          cTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACAAAGTCGCACCCTGCCCCGCCGCAACCTGTTg          >  7:30447/1‑149 (MQ=255)                                                                                                                                           tGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACAAAGTCGCAACCTGCCCCGCCGCAACCTGtt           >  7:36368/1‑147 (MQ=255)                                                                                                                                             cgcgTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACAAAGTCGCAACCTGCCCCGCCGCAACCTGTTGa         <  7:168818/147‑1 (MQ=255)                                                                                                                                                gTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACAAAGTCGCAACCTGCCCCGCCGCAACCTGTTGCGCCAc    >  1:162043/1‑149 (MQ=255)                                                                                                                                                gTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACAAAGTCGCAACCTGCCCCGCCGCAACCTGTTGAGCCAc    <  1:104472/149‑1 (MQ=255)                                                                                                                                                  GCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACAAAGTCGCAACCTGCCCCGCCGCAACCTGTTGAGCCACTg  <  1:141244‑M1/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                     CAAGCTTCCACTCCCCAGCCGAGCGCTCTGCGCGATGATGGGTCGACATGAACCCATATTTGAACGCCCCTGAACGCGAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACAAAGTCGCAACCTGCCCCGCCGCAACCTGTTGAGCCACTG  >  NC_002947/5599365‑5599658

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