New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	5175709 =	42 (0.740)	3 (0.060)	3/248	NT	5.8%	coding (512/945 nt)	PP_4551	alpha/beta family hydrolase
	?	NC_002947	5175751 =	57 (1.050)					coding (554/945 nt)	PP_4551	alpha/beta family hydrolase

ACCTGCTGCATCACAGTGGCAGCCTGCACGGGGCGATCCTTAGCGGTTCCAACTACCAACCGGCCGTGCTGTACCGCTTTGCGCGGCTGATTGCCCGGCTGGAGTCCTGGCGCCAGGGGCCACTGGGCAAGAGTGCACTGATCGAATGGTTGTCCTTCGGTTCGTTCAACAACGCCTTCAAGCCCAACCGCACGGCCTTCGACTGGCTCAGCCGCGACCCGGGCGAGGTCGATCAGTACGTGAATGACCCGTTGTGCGGCTTTCGCTGCACCAACCAGTTGTGGCTCGATCTGCT  >  NC_002947/5175562‑5175856                                                                                                                                                    |                                                                                                                                                    aCCTGCTGCATCACAGTGGCAGCCTGCACGGGGCGATCCTTAGCGGTTCCAACTACCAACCGGCCGTGCTGTACCGCTTTGCGCGGCTGATTGCCCGGCTGCAGTCCTGGCGCCAGGGGCCACTGGGCAAGAGTGCACTGATCGAATgg                                                                                                                                                    >  4:18098/1‑149 (MQ=255)  ccTGCTGCATCACAGTGGCAGCCTGCACGGGGCGATCCTTAGCGGTTCCAACTACCAACCGGCCGTGCTGTACCGCTTTGCGCGGCTGATTGCCCGGCTGGAGTCCTGGCGCCAGGGGCCACTGGGCAAGAGTGCACTGATCGAATgg                                                                                                                                                    <  5:139617/148‑1 (MQ=255)   ctgctgCATCACAGTGGCAGCCTGCACGGGGCGATCCTTAGCGGTTCCAACTACCAACCGGCCGTGCTGTACCGCTTTGCGCGGCTGATTGCCCGGCTGGAGTCCTGGCGCCAGGGGCCACTGGGCAAGAGTGCACTGATCGAATGGtt                                                                                                                                                  <  6:92246/149‑1 (MQ=255)         cATCACAGTGGCAGCCTGCACGGGGCGATCCTTAGCGGTTCCAACTACCAACCGGCCGTGCTGTACCGCTTTGCGCGGCTGATTGCCCGGCTGGAGTCCTGGCGCCAGGGGCCACTGGGCAAGAGTGCACTGATCGAATGGTTGTCCtt                                                                                                                                            >  1:9511/1‑149 (MQ=255)         cATCACAGTGGCAGCCTGCACGGGGCGATCCTTAGCGGTTCCAACTACCAACCGGCCGTGCTGTACCGCTTTGCGCGGCTGATTGCCCGGCTGGAGTCCTGGCGCCAGGGGCCACTGGGCAAGAGTGCACTGATCGAATGGTTGTCCtt                                                                                                                                            >  3:342948/1‑149 (MQ=255)          aTCACAGTGGCAGCCTGCACGGGGCGATCCTTAGCGGTTCCAACTACCAACCGGCCGTGCTGTACCGCTTTGCGCGGCTGATTGCCCGGCTGGAGTCCTGGCGCCAGGGGCCACTGGGCAAGAGTGCACTGATCGAATGGTTGTCCTTc                                                                                                                                           <  7:116567/149‑1 (MQ=255)            cacaGTGGCAGCCTGCACGGGGCGATCCTTAGCGGTTCCAACTACCAACCGGCCGTGCTGTACCGCTTTGCGCGGCTGATTGCCCGGCTGGAGTCCTGGCGCCAGGGGCCACTGGGCAAGAGTGCACTGATCGAATGGTTGTCCTTCgg                                                                                                                                         >  3:208138/1‑149 (MQ=255)                          gCACGGGGCGATCCTTAGCGGTTCCAACTACCAACCGGCCGTGCTGTACCGCTTTGCGCGGCTGATTGCCCGGCTGGAGTCCTGGCGCCAGGGGCCACTGGGCAAGAGTGCACTGATCGAATGGTTGTCCTTCGGTTCGTTCAACAACg                                                                                                                           >  8:28071/1‑149 (MQ=255)                                  cGATCCTTAGCGGTTCCAACTACCAACCGGCCGTGCTGTACCGCTTTGCGCGGCTGATTGCCCGGCTGGAGTCCTGGCGCCAGGGGCCACTGGGCAAGAGTGCACTGATCGAATGGTTGTCCTTCGGTTCGTTCAACAACGCCTTCAAg                                                                                                                   >  7:240571/1‑149 (MQ=255)                                    aTCCTTAGCGGTTCCAACTACCAACCGGCCGTGCTGTACCGCTTTGCGCGGCTGATTGCCCGGCTGGAGTCCTGGCGCCAGGGGCCACTGGGCAAGAGTGCACTGATCGAATGGTTGTCCTTCGGTTCGTTCAACAACGCCTTCAAGcc                                                                                                                 <  8:229656/149‑1 (MQ=255)                                     tCCTTAGCGGTTCCAACTACCAACCGGCCGTGCTGTACCGCTTTGCGCGGCTGATTGCCCGGCTGGAGTCCTGGCGCCAGGGGCCACTGGGCAAGAGTGCACTGATCGAATGGTTGTCCTTCGGTTCGTTCAACAACGCCTTCAAGccc                                                                                                                <  7:28071/149‑1 (MQ=255)                                           gCGGTTCCAACTACCAACCGGCCGTGCTGTACCGCTTTGCGCGGCTGATTGCCCGGCTGGAGTCCTGGCGCCAGGGGCCACTGGGCAAGAGTGCACTGATCGAATGGTTGTCCTTCGGTTCGTTCAACAACGCCTTCAAGCCCAACCGc                                                                                                          >  6:208262/1‑149 (MQ=255)                                               atcgagccacaactggttggtgcagcgaaagccgcacaacgggtcattcacgtactgatcgacctcgcccgggtcgcggctgagccagtcgaaggccgtgcGGTTGTCCTTCGGTTCGTTCAACAACGCCTTCAAGCCCAACCGCACgg                                                                                                      >  8:79345‑M1/102‑149 (MQ=255)                                                    aCTACCCACCGGCCGTGCTGTACCGCTTTGCGCGGCTGATTGCCCGGCTGGAGTCCTGGCGCCAGGGGCCACTGGGCAAGAGTGCACTGATCGAATGGTTGTCCTTCGGTTCGTTCAACAACGCCTTCAAGCCCAACCGCACGGCCTTc                                                                                                 <  5:84958/149‑1 (MQ=255)                                                      tACCAACCGGCCGTGCTGTACCGCTTTGCGCGGCTGATTGCCCGGCTGGAGTCCTGGCGCCAGGGGCCACTGGGCAAGAGTGCACTGATCGAATGGTTGTCCTTCGGTTCGTTCAACAACGCCTTCAAGCCCAACCGCACg                                                                                                       >  6:305731/1‑141 (MQ=255)                                                      tACCAACCGGCCGTGCTGTACCGCTTTGCGCGGCTGATTGCCCGGCTGGAGTCCTGGCGCCAGGGGCCACTGGGCAAGAGTGCACTGATCGAATGGTTGTCCTTCGGTTCGTTCAACAACGCCTTCAAGCCCAACCGCACg                                                                                                       <  5:305731/141‑1 (MQ=255)                                                       aCCAACCGGCCGTGCTGTACCGCTTTGCGCGGCTGATTGCCCGGCTGGAGTCCGGGCGCCACGGGCCCCTGGGCAAGAGTTCACAGATACAATGGTTGTCCTTAGGTTCGTTGAACCACGCCGTCAAGCCCCACCGCACGGccataga                                                                                               >  1:94114/1‑143 (MQ=255)                                                        ccAACCGGCCGTGCTGTACCGCTTTGCGCGGCTGATTGCCCGGCTGGAGTCCTGGCGCCAGGGCCAACTGGGCAAGAGTGCACTGATCGAATGGTTGTCCTTCGGTTAGTTCAACAACGCCTTCAAGCCCAACCGCACGGCCTTCGaat                                                                                             >  6:14452/1‑147 (MQ=255)                                                             cGGCCGTGCTGTACCGCTTTGCGCTGCTGATTGCCCGGATGGAGTCCTGGCGCGAGGGGCCACTGGGCAAGAGTGCACTGATCGAATGGTTGTCCTTCGGTTCGTTCAACAACTCCTTCAAGCCCAACCGCACGGCCTTCGACTGGCTc                                                                                        <  5:14452/149‑1 (MQ=255)                                                              ggCCGTGCTGTACCGCTTTGCGCGGCTGATTGCCCGGCTGGAGTCCTGGCGCCAGGGGCCACTGGGCAAGAGTGCACTGATCGAATGGGTGTCCTTCGGTTCGTTCAACAACGCCTTCAAGCCCAACCGCACGGCCTTCGACTGGCTCa                                                                                       <  2:317859/149‑1 (MQ=255)                                                                 cGTGCTGTACCGCTTTGCGCGGCTGATTGCCCGGCTGGAGTCCTGGCGCCAGGGGCCACTGGGCAAGAGTGCACTAATCGAATGGTTGTCCTTCGGTTCGTTCAACAACGCCTTCAAGCCCAACCGCACGGCCTTCGACTGGCTCAGcc                                                                                    >  1:71936/1‑149 (MQ=255)                                                                     cTGTACCGCTTTGCGCGGCTGATTGCCCGGCTGCAGTCCTGGCGCCAGGGGCCACTGGGCAAGAGTGCACTGATCGAATGGTTGTCCTTCGGTTCGTTCAACAACGCCTTCAAGCCCAACCGCACGGCCTTCGACTGGCTCAGCCGCGa                                                                                <  3:18098/149‑1 (MQ=255)                                                                       gTACCGCTTTGCGCGGCTGATTGCCCGGCTGGAGTCCTGGCGCCAGGGGCCACTGGGCAAGAGTGCACTGCTCGAATGGTTGTCCTTCGGTTCGTTAAACAACGCCTTCAAGCCCAACCGCACGGCCTTCGACTGGCTCAGCCCCCAcc                                                                              >  3:19065/1‑149 (MQ=255)                                                                        tACCGCTTTGCGCGGCTGATTGCCCGGCTGGAGTCCTGGCGCCAGGGGCCACTGGGCAAGAGTGCACTGATCGAATGGTTGTCCTTCGGTTCGTTCAACAACGCCTTCAAGCCCAACCGCACGGCCTTCGACTGGCTCAGCCGCGAccc                                                                             >  7:207382/1‑149 (MQ=255)                                                                         aCCGCTTTGCGCGGCTGATTGCCCGGCTGGAGTCCTGGCGCCAGGGGCCACTGGGCAAGAGTGCACTGATCGAATGGTTGTCCTTCGGTTCGTCCAACAACGCCTTCAAGCCCAACCGCACGGCCTTCGACTGGCTCAGCCGCGACCCg                                                                            >  8:308222/1‑149 (MQ=255)                                                                           cGCTTTGCGCGGCTGATTGCCCGGCTGGAGTCCTGGCGCCAGGGGCCACTGGGCAAGAGTGCACTGATCGAATGGTTGTCCTTCGGTTCGTTCAACAACGCCTTCAAGCCCAACCGCACGGCCTTCGACTGGCTCAGCCGCGACCCggg                                                                          <  2:94114/149‑1 (MQ=255)                                                                              ccgcacaacgggtcattcacgtactgatcgacctcgcccgggtcgcggctgagccagtcgaaggccgtgcGGTTGTCCTTCGGTTCGTTCAACAACGCCTTCAAGCCCAACCGCACg                                                                                                       <  4:97560‑M1/47‑1 (MQ=255)                                                                              ccgcacaacgggtcattcacgtactgatcgacctcgcccgggtcgcggctgagccagtcgaaggccgtgcGGTTGTCCTTCGGTTCGTTCAACAACGCCTTCAAGCCCAACCGCACg                                                                                                       >  3:97560‑M1/71‑117 (MQ=255)                                                                                  cgcgGCTGATTGCCCGGCTGGAGTCCTGGCGCCAGGGGCCACTGGGCAAGAGTGCACTGATCGAATGGTTGTCCTTCGGTTCGTTCAACAACGCCTTCAAGCCCAACCGCACGGCCTTCGACTGGCTCAGCCGCGACCCGGGCGAGGTc                                                                   >  6:344301/1‑149 (MQ=255)                                                                                       cTGATTGCCCGGCTGGAGTCCTGGCGCCAGGGGCCACTGGGCAAGAGTGCACTGATCGAATGGTTGTCCTTCGGTTCGTCCAACAACGCCTTCAAGCCCAACCGCACGGCCTTCGACTGGCTCAGCCGCGACCCGGGCGAGGTCGATCa                                                              <  7:308222/149‑1 (MQ=255)                                                                                             gCCCGGCTGGAGTCCTGGCGCCAGGGGCCACTGGGCAAGAGTGCACTGATCGAATGGTTGTCCTTCGGTTCGTTCAACAACGCCTTCAAGCCCAACCGCACGGCCTTCGACTGGCTCAGCCGCGACCCGGGCGAGGTCGATCAGTACGt                                                        >  2:87957/1‑149 (MQ=255)                                                                                              cccGGCTGGAGTCCTGGCGCCAGGGGCCACTGGGCAAGAGTGCACTGATCGAATGGTTGTCCTTCGGTTCGTTCAACAACGCCTTCAAGCCCAACCGCACGGCCTTCGa                                                                                               <  5:26263/109‑1 (MQ=255)                                                                                              cccGGCTGGAGTCCTGGCGCCAGGGGCCACTGGGCAAGAGTGCACTGATCGAATGGTTGTCCTTCGGTTCGTTCAACAACGCCTTCAAGCCCAACCGCACGGCCTTCGa                                                                                               >  6:26263/1‑109 (MQ=255)                                                                                               ctggCTGGAGTCCTGGCGCCAGGGGCCACTGGGCAAGAGTGCACTGATCGAATGGTTGTCCTTCGGTTCGTTCAACAACGCCTTCAAGCCCAACCGCACGGCCTTCGACTGGCTCAGCCGCGACCCGGGc                                                                         <  1:66616/128‑1 (MQ=255)                                                                                               ctggCTGGAGTCCTGGCGCCAGGGGCCACTGGGCAAGAGTGCACTGATCGAATGGTTGTCCTTCGGTTCGTTCAACAACGCCTTCAAGCCCAACCGCACGGCCTTCGACTGGCTCAGCCGCGACCCGGGc                                                                         >  2:66616/3‑130 (MQ=255)                                                                                                 ggcgtgAGGCCTGGCGCCAGGGGCCACTGGGCAAGAGTGCACTGATCGAATGGTTGTCCTTCGGTTCGTTCAACAACGCCTTCAAGCCCAACCGCACGGCCTTCGACTGGCTCAGCCGCGACCCGGGCGAGGTCGATCAGTACGTGAAt                                                    <  8:240571/144‑1 (MQ=255)                                                                                                   cTGGAGTCCTGGCGCCAGGGGCCACTGGGCAAGAGTGCACTGATCGAATGGTTGTCCTTCGGTTCGTTCAACAACGCCTTCAAGCCCAACCGCACGGCCTTCGACTGGCTCAGCCGCGACCCGGGCGAGGTCGATCAGTACGTGAATGa                                                  >  4:53975/1‑149 (MQ=255)                                                                                                        ccagtGGCCCCTGGCGCCACTGGGCAAGAGTGCACTGATCGAATGGTTGTCCTTCGGTTCGTTCAACAACGCCTTCAAGCCCAACCGCACGGCCTTCGACTGGCTCAGCCGCGACCCGGGCGAGGTCGATCAGTACGTGAATGACCCGt                                             >  6:343828/5‑149 (MQ=255)                                                                                                            tGGCGCCAGGGGCCACTGGGCAAGAGTGCACTGATCGAATGGTTGTCCTTCGGTTCGTTCAACAACGCCTTCAAGCCCAACCGCACGGCCTTCGACTGGCTCAGCCGCGACCCGGGCGAGGTCGATCAGTACGTGAATGACCCGTtgtg                                         >  8:89006/1‑149 (MQ=255)                                                                                                                 ccAGGGGCCACTGGGCAAGAGTGCACTGATCGAATGGTTGTCCTTCGGTTCGTTCAACAACGCCTTCAAGCCCAACCGCACGGCCTTCGACTGGCTCAGCCGCGACCCGGGCGAGGTCGATCAGTACGTGAATGACCCGTTGTGCGGCt                                    >  4:112237/1‑149 (MQ=255)                                                                                                                      ggCCACTGGGCAAGAGTGCACTGATCGAATGGTTGTCCTTCGGTTCGTTCAACAACGCCTTCAAGCCCAACCGCACGGCCTTCGACTGGCTCAGCCGCGACCCGGGCGAGGTCGATCAGt                                                            >  8:104120/1‑120 (MQ=255)                                                                                                                      ggCCACTGGGCAAGAGTGCACTGATCGAATGGTTGTCCTTCGGTTCGTTCAACAACGCCTTCAAGCCCAACCGCACGGCCTTCGACTGGCTCAGCCGCGACCCGGGCGAGGTCGATCAGt                                                            <  7:104120/120‑1 (MQ=255)                                                                                                                        ccACTGGGCAAGAGTGCACTGATCGAATGGTTGTCCTTCGGTTCGTTCAACAACGCCTTCAAGCCCAACCGCACGGCCTTCGACTGGCTCAGCCGCGACCCGGGCGAGGTCGATCAGTACGTGAATGACCCGTTGTGCGGCTTTCGCTg                             <  5:208262/149‑1 (MQ=255)                                                                                                                               gCAAGAGTGCACTGATCGAATGGTTGTCCTTCGGTTCGTTCAACAACGCCTTCAAGCCCAACCGCACGGCCTTCGACTGGCTCAGCCGCGACCCGGGCGAGGTCGATCAGTACGTGAATGACCCGTTGTGCGGCTTTCGCTGCaccaac                      >  6:131047/1‑149 (MQ=255)                                                                                                                                 aaGAGTGCACTGATCGAATGGTTGTCCTTCGGTTCGTTCAACAACGCCTTCAAGCCCAACCGCACGGCCTTCGACTGGCTCAGCCGCGACCCGGGCGAGGTCGATCAGTACGTGAATGACCCGTTGTGCGGCTTTCGCTGCaccaacca                    <  3:53975/149‑1 (MQ=255)                                                                                                                                      tGCACTGATCGAATGGTTGTCCTTCGGTTCGTTCAACAACGCCTTCAAGCCCAACCGCACGGCCTTCGACTGGCTCAGCCGCGACCCGGGCGAGGTCGATCAGTACGTGAATGACCCGTTGTGCGGCTTTCGCTGCACCCACCagtagt               >  4:144893/1‑149 (MQ=255)                                                                                                                                        cACTAATCGAATGGTTGTCCTTCGGTTCGTTCAACAACGCCTTCAAGCCCAACCGCACGGCCTTCGACTGGCTCAGCCGCGACCCGGGCGAGGTCGATCAGTACGTGAATGACCCGTTGTGCGGCTTTCGCTGCACCAACCAGTTGTgg             <  2:71936/149‑1 (MQ=255)                                                                                                                                          cTGATCGAATGGTTGTCCTTCGGTTCGTTCAACAACGCCTTCAAGCCCAACCGCACGGCCTTCGACTGGCTCAGCCGCGACCCGGGCGAGGTCGATCAGTACGTGAATGACCCGTTGTGCGGCTTTCGCTGCACCAACCAGTTGTGGCt           <  7:89006/149‑1 (MQ=255)                                                                                                                                                 aaTGGTTGTCCTTCGGTTCGTTCAACAACGCCTTCAAGCCCAACCGCACGGCCTTCGACTGGCTCAGCCGCGACCCGGGCGAGGTCGATCAGTACGTGAATGACCCGTTGTGCGGCTTTCGCTGCACCAACCAGTTGTGGCTCGATctg    >  6:222004/1‑149 (MQ=255)                                                                                                                                                  aTGGTTGTCCTTCGGTTCGTTCAACAACGCCTTCAAGCCCAACCGCACGGCCTTCGACTGGCTCAGCCGCCACCCGGGCGAGGTCGATCAGTACGTGAATGACCCGTTGTGCGGCTTTCGCTTCACCAACCACTTGTGGATCGATctgc   >  6:202706/1‑149 (MQ=255)                                                                                                                                                   tGGTTGTCATTCGGTTCGTTCAACAACGCCTTCAAGCCCAACCGCACGGCCTTCGACTGGCTCAGCCGCGACCCGGGCGAGGTCGATCAGTACGTGAATGACCCGTTGTGCGGCTTTCGCTGCACCAACCAGTTGTGGCTCGATctgct  >  1:40798/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                    ACCTGCTGCATCACAGTGGCAGCCTGCACGGGGCGATCCTTAGCGGTTCCAACTACCAACCGGCCGTGCTGTACCGCTTTGCGCGGCTGATTGCCCGGCTGGAGTCCTGGCGCCAGGGGCCACTGGGCAAGAGTGCACTGATCGAATGGTTGTCCTTCGGTTCGTTCAACAACGCCTTCAAGCCCAACCGCACGGCCTTCGACTGGCTCAGCCGCGACCCGGGCGAGGTCGATCAGTACGTGAATGACCCGTTGTGCGGCTTTCGCTGCACCAACCAGTTGTGGCTCGATCTGCT  >  NC_002947/5175562‑5175856

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.