New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 5601631	36 (0.690)	4 (0.080)	4/246	NT	10.2%	coding (1851/5418 nt)	PP_4924	subtilase family serine protease
	?	NC_002947	= 5601675	37 (0.770)					coding (1807/5418 nt)	PP_4924	subtilase family serine protease

TCAGCAAGCTGTCGCCCGAGCCGCCCTTGATGAAAACTCCGCGTTGACCGCCGCTGGCGATGATGTCCTTACCCCGCCCGCCTTGCGCGATGGTCAGCCCGGCCTTGGCCATGTTCAGCCCAACGCCATGCTCCCCCACGATCAGTGCCGTGTCGCGCCCACCATTGCCCTGCATGTTCGCCCAATCGTCGCTGGCACTGATCACGAATACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACAT  >  NC_002947/5601483‑5601776                                                                                                                                                     |                                                                                                                                                  tCAGCAAGCTGTCGCCCGAGCCGCCCTTGATGAAAACTCCGCGTTGACCGCCGCTGGCGATGATGTCCTTACCCCGCCCGCCTTGCGCGATGGTCAGCCCGGCCTTGGCCATGTTCAGCCCAACGCCATGCTCCCCCACGATCAGTGCC                                                                                                                                                   >  7:249153‑M1/1‑149 (MQ=255)     cAAGCTGTCGCCCGAGCCGCCCTTGATGAAAACTCCGCGTTGACCGCCGCTGGCGATGATGTCCTTACCCCGCCCGCCTTGCGCGATGGTCAGCCCGGCCTTGGCCATGTTCAGCCCAACGCCATGCTCCCCCACGATCAGTGCCgtgt                                                                                                                                               >  3:268235/1‑149 (MQ=255)     cAAGCTGTCGCCCGAGCCGCCCTTGATGAAAACTCCGCGTTGACCGCCGCTGGCGATGATGTCCTTACCCCGCCCGCCTTGCGCGATGGTCAGCCCGGCCTTGGCCATGTTCAGCCCAACGCCATGCTCCCCCACGATCAGTGCCgtgt                                                                                                                                               >  7:164463/1‑149 (MQ=255)      aaGCTGTCGCCCGAGCCGCCCTTGATGAAAACTCCGCGTTGACCGCCGCTGGCGATGATGTCCTTACCCCGCCCGCCTTGCGCGATGGTCAGCCCGGCCTTGGCCATGTTCAGCCCAACGCCATGCTCCCCCACGATCAGTGCCGTGTc                                                                                                                                              <  8:164463/149‑1 (MQ=255)       aGCTGTCGCCCGAGCCGCCCTTGATGAAAACTCCGCGTTGACCGCCGCTGGCGATGATGTCCTTACCCCGCCCGCCTTGCGCGATGGTCAGCCCGGCCTTGGCCATGTTCAGCCCAACGCCATGCTCCCCCACGATCAGTGCCGTGTcg                                                                                                                                             <  5:238708/149‑1 (MQ=255)            tCGCCCGAGCCGCCCTTGATGAAAACTCCGCGTTGACCGCCGCTGGCGATGATGTCCTTACCCCGCCCGCCTTGCGCGATGGTCAGCCCGGCCTTGGCCATGTTCAGCCCAACGCCATGCTCCCCCACGATCAGTGCCGTGTCGCGccc                                                                                                                                        <  4:186975/149‑1 (MQ=255)              gCCCGAGCCGCCCTTGATGAAAACTCCGCGTTGACCGCCGCTGGCGATGATGTCCTTACCCCGCCCGCCTTGCGCGATGGTCAGCCCGGCCTTGGCCATGTTCAGCCCAACGCCATGCTCCCCCACGATCAGTGCCGTGTCGCGCccac                                                                                                                                      >  5:238999/1‑149 (MQ=255)                 cGAGCCGCCCTTGATGAAAACTCCGCGTTGACCGCCGCTGGCGATGATGTCCTTACCCCGCCCGCCTTGCGCGATGGTCAGCCCGGCCTTGGCCATGTTCAGCCCAACGCCATGCTCCCCCACGATCAGTGCCGTGTCGCGCCCACCAt                                                                                                                                   <  7:235492/149‑1 (MQ=255)                  gAGCCGCCCTTGATGAAAACTCCGCGTTGACCGCCGCTGGCGATGATGTCCTTACCCCGCCCGCCTTGCGCGATGGTCAGCCCGGCCTTGGCCATGTTCAGCCCAACGCCATGCTCCCCCACGATCAGTGCCGTGTCGCGCCCACCAtt                                                                                                                                  <  4:62095/149‑1 (MQ=255)                     ccgccCTTGATGAAAACTCCGCGTTGACCGCCGCTGGCGATGATGTCCTTACCCCGCCCGCCTTGCGCGATGGTCAGCCCGGCCTTGGCCATGTTCAGCCCAACGCCATGCTCCCCCACGATCAGTGCCGTGTCGCGCCCACCATTGcc                                                                                                                               >  2:83062/1‑149 (MQ=255)                          cTTGATGAAAACTCCGCGTTGACCGCCGCTGGCGATGATGTCCTTACCCCGCCCGCCTTGCGCGATGGTCAGCCCGGCCTTGGCCATGTTCAGCCCAACGCCATGCTCCCCCACGATCAGTGCCGTGTCGCGCCCACCATTGCCCTGCa                                                                                                                          >  6:277585/1‑149 (MQ=255)                           ttGATGAAAACTCCGCGTTGACCGCCGCTGGCGATGATGTCCTTACCCCGCCCGCCTTGCGCGATGGTCAGCCCGGCCTTGGCCATGTTCAGCCCAACGCCATGCTCCCCCACGATCAGTGCCGTGTCGCGCCCACCATTGCCCTGCAt                                                                                                                         >  8:261868/1‑149 (MQ=255)                                 aaaaCTCCGCGTTGACCGCCGCTGGCGATGATGTCCTTACCCCGCCCGCCTTGCGCGATGGTCAGCCCGGCCTTGGCCATGTTCAGCCCAACGCCATGCTCCCCCACGATCAGTGCCGTGTCGCGCCCACCATTGCCCTGCATGTTCGc                                                                                                                   >  5:169067/1‑149 (MQ=255)                                  aaaCTCCGCGTTGACCGCCGCTGGCGATGATGTCCTTACCCCGCCCGCCTTGCGCGATGGTCAGCCCGGCCTTGGCCATGTTCAGCCCAACGCCATGCTCCCCCACGATCAGTGCCGTGTCGCGCCCACCATTGCCCTGCATGTTCGcc                                                                                                                  <  8:249153/149‑1 (MQ=255)                                      tCCGCGTTGACCGCCGCTGGCGATGATGTCCTTACCCCGCCCGCCTTGCGCGATGGTCAGCCCGGCCTTGGCCATGTTCAGCCCAACGCCATGCTCCCCCACGATCAGTGCCagcgacgattgggcgaacatgcagggcaatggtgggc                                                                                                              >  4:112706‑M1/1‑112 (MQ=255)                                         gcgTTGACCGCCGCTGGCGATGATGTCCTTACCCCGCCCGCCTTGCGCGATGGTCAGCCCGTCCTTGGCCATGTTCAGCCCAACTCCATGCTCCCCCACGATCAGTGCCGTGTCGCGCCCACCATTGCCCTGCATGTTCGCCCAAtcgt                                                                                                           <  6:164869/149‑1 (MQ=255)                                              gACCGCCGCTGGCGATGATGTCCTTACCCCGCCCGCCTTGCGCGATGGTCAGCCCGGCCTTGGCCATGTTCAGCCCAACGCCATGCTCCCCCACGATCAGTGCCGTGTCGCGCCCACCATTGCCCTGCATGTTCGCCCAATCGTCGCt                                                                                                       <  4:107051/148‑1 (MQ=255)                                                 cgccgcTGGCGATGATGTCCTTACCCCGCCCGCCTTGCGCGATGGTCAGCCCGGCCTTGGCCATGTTCAGCCCAACGCCATGCTCCCCCACGATCAGTGCCGTGTCGCGCCCACCATTGCCCTGCATGTTCGCCCAATCGTCGCTggca                                                                                                   >  4:234831‑M2/1‑145 (MQ=255)                                                  gccgcTGGCGATGATGGCCTTACCCCGCACTCCTTGCGCGATGGGCAGCCCGGCCTTGGCCATGTTCAGCCCAACGCCATGCTCCCCCACGATCAGTGCCGTGTCGCGCCCACCATTGCCCTGCATGTTCGCCCAATCGTCGCTggcac                                                                                                  >  8:151798‑M2/1‑144 (MQ=255)                                                      cTGGCGATGATGTCCTTACCCCGCCCGCCTTGCGCGATGGTCAGCCCGGCCTTGGCCATGTTCAGCCCAACGCCATGCTCCCCCACGATCAGTGCCGTGTCGCGCCCACCATTGCCCTGCATGTTCGCCCAATCGTCGCTggcac                                                                                                  >  6:118335‑M2/1‑140 (MQ=255)                                                      cTGGCGATGATGTCCTTACCCCGCCCGCCTTGCGCGATGGTCAGCCCGGCCTTGGCCATGTTCAGCCCAACGCCATGCTCCCCCACGATCAGTGCCGTGTCGCGCCCACCATTGCCCTGCATGTTCGCCCAATCGTCGCTggcac                                                                                                  <  5:118335‑M2/145‑6 (MQ=255)                                                         gCGATGATGTCCTTACCCCGCCCGCCTTGCGCGATGGTCAGCCCGGCCTTGGCCATGTTCAGCCCAACGCCATGCTCCCCCACGATCAGTGCCagcgacgattgggcgaacatgcagggcaatggtgggcgcgacacggcactgatcgt                                                                                           <  7:98891‑M1/149‑57 (MQ=255)                                                         gCGATGATGTCCTTACCCCGCCCGCCTTGCGCGATGGTCAGCCCGGCCTTGGCCATGTTCAGCCCAACGCCATGCTCCCCCACGATCAGTGCCagcgacgattgggcgaacatgcagggcaatggtgggcgcgacacggcactgatcg                                                                                            <  4:292022‑M1/148‑56 (MQ=255)                                                         gCGATGATGTCCTTACCCCGCCCGCCTTGCGCGATGGTCAGCCCGGCCTTGGCCATGTTCAGCCCAACGCCATGCTCCCCCACGATCAGTGCCGTGTCGCGCCCACCATTGCCCTGCATGTTCGCCCAATCGTCGCTGGCACTGATCAc                                                                                           <  4:259520/149‑1 (MQ=255)                                                         gCGATGATGTCCTTACCCCGCCCGCCTTGCGCGATGGTCAGCCCGGCCTTGGCCATGTTCAGCCCAACGCCATGCTCCCCCACGATCAGTGCCGTGTCGCGCCCACCATTGCCCTGCATGTTCGCCCAATCGTCGCTGGCACTGATCAc                                                                                           <  1:224882/149‑1 (MQ=255)                                                          cGATGATGTCCTTACCCCGCCCGCCTTGCGCGATGGTCAGCCCGGCCTTGGCCATGTTCAGCCCAACGCCATGCTCCCCCACGATCAGTGCCGTGTCGCGCCCACCATTGCCCTGCATGTTCGCCCAATCGTCGCTGGCACTGATCACg                                                                                          <  1:83062/149‑1 (MQ=255)                                                           gatgatGTCCTTACCCCGCCCGCCTTGCGCGATGGTCAGCCCGGCCTTGGCCATGTTCAGCCCAACGCCATGCTCCCCCACGATCAGTGCCagcgacgattgggcgaacatgcagggcaatggtgggcgcgacacggcactgatcgtgg                                                                                         <  3:112706‑M1/149‑59 (MQ=255)                                                           gatgatGTCCTTACCCCGCCCGCCTTGCGCGATGGTCAGCCCGGCCTTGGCCATGTTCAGCCCAACGCCATGCTCCCCCACGATCAGTGCCGTGTCGCGCCCACCATTGCCCTGCATGTTCGCCCAATCGTCGCTGGCACTTATCACg                                                                                          >  2:327295/1‑148 (MQ=255)                                                                tgTCCTTACCCCGCCCGCCTTGCGCGATGGTCAGCCCGGCCTTGGCCATGTTCAGCCCAACGCCATGCTCCCCCACGATCAGTGCCGTGTCGCGCCCACCATTGCCCTGCATGTTCGCCCAATCGTCGCTGGCACTGATCACGAATACa                                                                                    <  5:81580/149‑1 (MQ=255)                                                                  tCCTTACCCCGCCCGCCTTGCGCGATGGTCAGCCCGGCCTTGGCCATGTTCAGCCCAACGCCATGCTCCCCCACGATCACTGCCGTGTCGCGCCCACCATTGCCCTGCATGTTCGCCCAATCGTCGCTGGCACTGATCACGAATACAtc                                                                                  <  4:268235/149‑1 (MQ=255)                                                                    cTTACCCCGCCCGCCTTGCGCGATGGTCAGCCCGGCCTTGGCCATGTTCAGCCCAACGCCATGCTCCCCCACGATCAGTGCCGTGTCGCGCCCACCATTGCCCTGCATGTTCGCCCAATCGTCGCTGGCACTGATCACGAATACAtcgt                                                                                >  8:76335/1‑149 (MQ=255)                                                                         cccgcccgCCTTGCGCGATGGTCAGCCCGGCCTTGGCCATGTTCAGCCGCACGCCATGCTCCCCCACGATCAGTGCCGTGTCGCGCCCACCAGTGCCCTGCATGTTCGCCCCATCGTCGCTGGCAGGCAGCACGAATACATCGAcgccc                                                                           >  3:199331/1‑149 (MQ=255)                                                                          ccgcccgcCTTGCGCGATGGTCAGCCCGGCCTTGGCCATGTTCAGCCCAACGCCATGCTCCCCCACGATCAGTGCCGTGTCGCGCCCACCATTGCCCTGCATGTTCGCCCAATCGTCGCTGGCACTGATCACGAATACATCGTcgccag                                                                          >  5:313965/1‑147 (MQ=255)                                                                                        cgctGGTGAGCCCGGCCTTGGCCAAGTTCAGCCCCACGCCATGCTCCCCCACGATCAGTGCCGTGTCGCGCCCACCATTTCCCTGCATGTTCGCCCAATCGTCGCTGGCACTGATCACGAATACTTCGTCGCCCGCCCCGCCGTGATAg                                                            <  6:238999/146‑1 (MQ=255)                                                                                         gATGGTCAGCCCGGCCTTGGCCATGTTCAGCCCAACGCCATGCTCCCCCACGATCAGTGCCGTGTCGCGCCCACCATTGCCCTGCATGTTCGCCCAATCGTCGCTGGCACTGATCACGAATACATCGTCGCCCGCCCCGCCGTGATAgg                                                           >  7:270250/1‑149 (MQ=255)                                                                                             gTCAGCCCGGCCTTGGCCATGTTCAGCCCAACGCCATGCTCCCCCACGATCAGTGCCGTGTCGCGCCCACCATTGCCCTGCATGTTCGCCCAATCGTCGCTGGCACTGATCACGAATACATCGTCGCCCGCCCCGCCGTGATAGGTGtt                                                       >  8:119769/1‑149 (MQ=255)                                                                                             gTCAGCCCGGCCTTGGCCATGTTCAGCCCAACGCCATGCTCCCCCACGATCAGTGCCGTGTCGCGCCCACCATTGCCCTGCATGTTCGCCCAATCGTCGCTGGCACTGATCACGAATACATCGTCGCCCGCCCCGCCGTGATAGGTGtt                                                       <  5:277585/149‑1 (MQ=255)                                                                                              tccgCCCGGCCGTGGCCATGTTCAGCCGAACGCCCTGCTCCCCCACGATCAGTGCCGTGGCGCGCCCACCATTGCCCTGCATGTTCGCCCAATCGTCGCTGGCACTGATCACGAATACATCGTCGCCCGCCCCGCCGTGATAGGTGTtg                                                      <  4:199331/146‑1 (MQ=255)                                                                                                  cccGGCCTTGGCCATGTTCAGCCCAACGCCATGCTCCCCCACGATCAGTGCCGTGTCGCGCCCACCATTGCCCTGCATGTTCGCCCAATCGTCGCTGGCACTGATCACGAATACATCGTCGCCCGCCccgccg                                                                  <  7:245504/133‑1 (MQ=255)                                                                                                  cccGGCCTTGGCCATGTTCAGCCCAACGCCATGCTCCCCCACGATCAGTGCCGTGTCGCGCCCACCATTGCCCTGCATGTTCGCCCAATCGTCGCTGGCACTGATCACGAATACATCGTCGCCCGCCccgccg                                                                  >  8:245504/1‑133 (MQ=255)                                                                                                       ccTTGGCCATGTTCAGCCCAACGCCATGCTCCCCCACGATCAGTGCCGTGTCGCGCCCACCATTGCCCTGCATGTTCGCCCAATCGTCGCTGGCACTGATCACGAATACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTgccaccg                                             >  8:73603/1‑149 (MQ=255)                                                                                                                        ccAACGCCATGCTCCCCCACGATCAGTGCCGTGTCGCGCCCACCATTGCCCTGCATGTTCGCCCAATCGTCGCTGGCACTGATCACGAATACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACt                            >  6:146648/1‑149 (MQ=255)                                                                                                                                aTGCTCCCCCACGATCAGTGCCGTGTCGCGCCCACCATTGCCCTGCATGTTCGCCCAATCGTCGCTGGCACTGATCACGAATACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTg                    <  7:73603/149‑1 (MQ=255)                                                                                                                                  gCTCCCCCACGATCAGTGCCGTGTCGCGCCCACCATTGCCCTGCATGTTCGCCCAATCGTCGCTGGCACTGATCACGAATACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTggg                  >  8:161176/1‑149 (MQ=255)                                                                                                                                        ccACGATCAGTGCCGTGTCGCGCCCACCATTGCCCTGCATGTTCGCCCAATCGTCGCTggcactgatcgtgggggagcatggcgttgggctgaacatggccaaggccgggctgaccatcgcgcaaggcgggcggggtaaggacatcatc            <  3:112706‑M2/149‑92 (MQ=255)                                                                                                                                         cACGATCAGTGCCGTGTCGCGCCCACCATTGCCCTGCATGTTCGCCCAATCGTCGCTGGCACTGATCACGAATACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCgg           >  7:292736/1‑149 (MQ=255)                                                                                                                                          aCGATCAGTGCCGTGTCGCGCCCACCATTGCCCTGCATGTTCGCCCAATCGTCGCTggcactgatcgtgggggagcatggcgttgggctgaacatggccaaggccgggctgaccatcgcgcaaggcgggcggggtaaggacatcatcgc          <  7:98891‑M2/149‑94 (MQ=255)                                                                                                                                           cGATCAGTGCCGTGTCGCGCCCACCATTGCCCTGCATGTTCGCCCAATCGTCGCTggcactgatcgtgggggagcatggcgttgggctgaacatggccaaggccgggctgaccatcgcgcaaggcgggcggggtaaggacatcatcgc          <  4:292022‑M2/148‑94 (MQ=255)                                                                                                                                             aTCAGTGCCGTGTCGCGCCCACCATTGCCCTGCATGTTCGCCCAATCGTCGCTGGCACTGATCACGAATACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCaa       <  7:161176/149‑1 (MQ=255)                                                                                                                                               cAGTGCCGTGTCGCGCCCACCATTGCCCTGCATGTTCGCCCAATCGTCGCTGGCACTGATCACGAATACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGa     <  7:76335/149‑1 (MQ=255)                                                                                                                                                aGTGCCGTGTCGCGCCCACCATTGCCCTGCATGTTCGCCCAATCGTCGCTGGCACTGATCACGAATACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGa     <  7:261868/148‑1 (MQ=255)                                                                                                                                                 gTGCCagcgacgattgggcgaacatgcagggcaatggtgggcgcgacacggcactgatcgtgggggagcatggcgttgggctgaacatggccaaggccgggctgcccatcgcgcaaggagtgcggggtaaggccatcatcgccagcggc   >  8:151798‑M1/1‑5 (MQ=255)                                                                                                                                                 gTGCCagcgacgattgggcgaacatgcagggcaatggtgggcgcgacacggcactgatcgtgggggagcatggcgttgggctgaacatggccaaggccgggctgaccatcgcgcaaggcgggcggggtaaggacatcatcgccag       >  6:118335‑M1/1‑5 (MQ=255)                                                                                                                                                 gTGCCagcgacgattgggcgaacatgcagggcaatggtgggcgcgacacggcactgatcgtgggggagcatggcgttgggctgaacatggccaaggccgggctgaccatcgcgcaaggcgggcggggtaaggacatcatcgccag       <  5:118335‑M1/145‑141 (MQ=255)                                                                                                                                                 gTGCCGTGTCGCGCCCACCATTGCCCTGCATGTTCGCCCAATCGTCGCTGGCACTGATCACGAATACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGAca   >  3:103076/1‑149 (MQ=255)                                                                                                                                                  tGCCagcgacgattgggcgaacatgcagggcaatggtgggcgcgacacggcactgatcgtgggggagcatggcgttgggctgaacatggccaaggccgggctgaccatcgcgcaaggcgggcggggtaaggacatcatcgccagcggcg  >  4:234831‑M1/1‑4 (MQ=255)                                                                                                                                                    ccGTGTCGCGCCCACCATTGCCCTGCATGTTCGCCCAATCGTCGCTGGCACTGATCACGAATACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGcc                         <  6:81963/124‑1 (MQ=255)                                                                                                                                                    ccGTGTCGCGCCCACCATTGCCCTGCATGTTCGCCCAATCGTCGCTGGCACTGATCACGAATACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGcc                         >  5:81963/1‑124 (MQ=255)                                                                                                                                                     |                                                                                                                                                  TCAGCAAGCTGTCGCCCGAGCCGCCCTTGATGAAAACTCCGCGTTGACCGCCGCTGGCGATGATGTCCTTACCCCGCCCGCCTTGCGCGATGGTCAGCCCGGCCTTGGCCATGTTCAGCCCAACGCCATGCTCCCCCACGATCAGTGCCGTGTCGCGCCCACCATTGCCCTGCATGTTCGCCCAATCGTCGCTGGCACTGATCACGAATACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACAT  >  NC_002947/5601483‑5601776

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.