New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 4172981	48 (1.060)	3 (0.070)	3/242	NT	5.8%	coding (571/825 nt)	PP_3671	aldo/keto reductase family oxidoreductase
	?	NC_002947	= 4173006	50 (1.140)					coding (596/825 nt)	PP_3671	aldo/keto reductase family oxidoreductase

CTGGGCGCGCCGTTGTCGGCAATTCAGGTGGAATACCATGTTCTGCTCAGCCAGCAACCGCTACTCGACTTCGCCCGCGGCCATGACCTGTTGCTCACCGCCTACACCCCGCTGGCACGCGGCCAGGCAGCGGCGCAGCCGGTGATCCAGGCGATCGCCCGCAAGCATGGCGTGCTGCCCAGCCAGGTGGCGCTGAAATGGCTGCTGGACCAGGACGGCGTGGCGGCGATCCCCAAAGCCAGCAGCCGCGAAAACCAGCTGGCCAACCTGGCCGCGTTGACCGTGCCGCTG  >  NC_002947/4172858‑4173148                                                                                                                                                     |                                                                                                                                               cTGGGCGCGCCGTTGTCGGCAATTCAGGTGGAATACCATGTTCTGCTCAGCCAGCAACCGCTACTCGACTTCGCCCGCGGCCATGACCTGTTGCTCACCGCCTACACCCCGCTGGCACGCGGCCAGGCAGCGGCGCAGCCGGTGATCCAg                                                                                                                                               <  2:307833‑M2/150‑2 (MQ=255)  tGGGCGCGCCGTTGTCGGCAATTCAGGTGGAATACCATGTTCTGCTCAGCCAGCAACCGCTACTCGACTTCGCCCGCGGCCATGACCTGTTGCTCACCGCCTACACCCCGCTGGCACTCGGCCAGGCAGCGGCGCAGCCGGTGATCCAgg                                                                                                                                              >  2:350583‑M2/1‑148 (MQ=255)               gTCGGCAATTCAGGTGGAATACCATGTTCTGCTCAGCCAGCAACCGCTACTCGACTTCGCCCGCGGCCATGACCTGTTGCTCACCGCCTACACCCCGCTGGCACGCGGCCAGGCAGCGGCGCAGCCGGTGATCCAGGCGATCGCCCGCaa                                                                                                                                 <  1:951630/150‑1 (MQ=255)                tCGGCAATTCAGGTGGAATACCATGTTCTGCTCAGCCAGCAACCGCTACTCGACTTCGCCCGCGGCCATGACCTGTTGCTCACCGCCTACACCCCGCTGGCACGCGGCCAGGCAGCGGCGCAGCCGGTGATCCAGGCGATCGCCCGCAAg                                                                                                                                >  1:1004660/1‑150 (MQ=255)                  ggCAATTCAGGTGGAATACCATGTTCTGCTCAGCCAGCAACCGCTACTCGACTTCGCCCGCGGCCATGACCTGTTGCTCACCGCCTACACCCCGCTGGCACGCGGCCAGGCAGCGGCGCAGCCGGTGATCCAGGCGATCGCCCGCAAGCa                                                                                                                              >  1:682457/1‑150 (MQ=255)                  ggCAATTCAGGTGGAATACCATGTTCTGCTCAGCCAGCAACCGCTACTCGACTTCGCCCGCGGCCATGACCTGTTGCTCACCGCCTACA‑CCCGCTGGCACGCGGCCAGGCAGCGGCGCAGCCGGTGATCCAGGCGATCGCCAGCAAGCAt                                                                                                                             >  2:570876/1‑150 (MQ=255)                     aaTTCAGGTGGAATACCATGTTCTGCTCAGCCAGCAACCGCTACTCGACTTCGCCCGCGGCCATGACCTGTTGCTCACCGCCTACACCCCGCTGGCACGCGGCCAGGCAGCGGCGCAGCCGGTGATCCAGGCGATCGCCCGCAAGCATgg                                                                                                                           <  1:998988/150‑1 (MQ=255)                      aTTCAGGTGGAATACCATGTTCTGCTCAGCCAGCAACCGCTACTCGACTTCGCCCGCGGCCATGACCTGTTGCTCACCGCCTACACCCCGCTGGCACGCGGCCAGGCAGCGGCGCAGCCGGTGATCCAGGCGATCGCCCGCAAGCATGGc                                                                                                                          <  1:307833/150‑1 (MQ=255)                      aTTCAGGTGGAATACCATGTTCTGCTCAGCCAGCAACCGCTACTCGACTTCGCCCGCGGCCATGACCTGTTGCTCACCGCCTACACCCCGCTGGCACGCGGCCAGGCAGCGGCGCAGCCGGTGATCCAGGCGATCGCCCGCAAGCATGGc                                                                                                                          >  2:48721/1‑150 (MQ=255)                                                  gccagcAACCGCTACTCGACTTCGCCCGCGGCCATGACCTGTTGCTCACCGCCTACACCCCGCTGGCACGCGGCCAGGCAGCGGCGCAGCCGGTGATCCAGGCGATCGCCCGCAAGCATGGAGTGCTGCCCAGCCAGGTGGCGCTGAAAt                                                                                              <  1:705877/150‑1 (MQ=255)                                                    cagcAACCGCTACTCGACTTCGCCCGCGGCCATGACCTGTTGCTCACCGCCTACACCCCGCTGGCACGCGGCCAGGCAGCGGCGCAGCCGGTGATCCAGGCGATCGCCCGCAAGCATGGCGTGCTGCCCAGCCAGGTGGCGCTGAAATgg                                                                                            >  2:883558/1‑150 (MQ=255)                                                     agcAACCGCTACTCGACTTCGCCCGCGGCCATGACCTGTTGCTCACCGCCTACACCCCGCTGGCACGCGGCCAGGCAGCGGCGCAGCCGGTGATCCAGGCGATCGCCCGCAAGCATGGCGTGCTGCCCAGCCAGGTGGCGCTGAAATGgc                                                                                           <  2:51332/150‑1 (MQ=255)                                                      gcAACCGCTACTCGACTTCGCCCGCGGCCATGACCTGTTGCTCACCGCCTACACCCCGCTGGCACGCGGCCAGGCAGCGGCGCAGCCGGTGATCCAGGCGATCGCCCGCAAGCATGGCGTGCTGCCCAGCCAGGTGGCGCTGAAATGgct                                                                                          >  1:302202/1‑150 (MQ=255)                                                      gcAACCGCTACTCGACTTCGCCCGCGGCCATGACCTGTTGCTCACCGCCTACACCCCGCTGGCACGCGGCCAGGCAGCGGCGCAGCCGGTGATCCAGGCGATCGCCCGCAAGCATGGCGTGCTGCCCAGCCAGGTGGCGCTGAAATGgct                                                                                          <  2:1004660/150‑1 (MQ=255)                                                                cTCGACTTCGCCCGCGGCCATGACCTGTTGCTCACCGCCTACACCCCGCTGGCACGCGGCCAGGCAGCGGCGCAGCCGGTGATCCAGGCGATCGCCCGCAAGCATGGCGTGCTGCCCAGCCAGGGGTCGCTGAAATGGCTGCTGGACCAg                                                                                >  1:772/1‑150 (MQ=255)                                                                           ccGCGGCCATGACCTGTTGCTCACCGCCTACACCCCGCTGGCAGGCGTCCAGCCAGCCGCGCAGCCGGTGATCCAGGCGATCGCCCGCAAGCATGGCGAGCTGCCCAGCCAGGTGGCGCAGAAATGGCTGCTAGACCAGGACGGCGTggc                                                                     >  2:233722/1‑150 (MQ=255)                                                                           ccGCGGCCATGACCTGTTGCTCACCGCCTACACCCCGCTGGCACGCGGCCAGGCAGCGACGCAGCCGGTGATCCAGGCGATCGCCCGCAAGCATGGCGTGCTGCCCAGCCAGGTGGCGCTGAAATGGCTGCTGGACCAGGACGGCGTggc                                                                     >  2:851983/1‑150 (MQ=255)                                                                                 ccATGACCTGTTGCTCACCGCCTACACCCCGCTGGCACCCGGCCAGGCAGCGGCGCAGCCGGTGATCCAGGCGATCGCCCGCAAGCATGGCGTGCTGCCCAGCCAGGTGGCGCTGAAATGGCTGCTGGACCAGGACGGCGTGGCGGCGAt                                                               >  1:936371/1‑150 (MQ=255)                                                                                  cATGACCTGTTGCTCACCGCCTACACCCCGCTGGCACGCGGCCAGGCAGCGGCGCAGCCGGTGATCCAGGCGATCGCCCGCAAGCATGGCGTGCTGCCCAGCCAGGTGGCGCTGAAATGGCTGCTGGACCAGGACGGCGTGGCGGCGATc                                                              >  2:1102109/1‑150 (MQ=255)                                                                                      aCCTGTTGCTCACCGCCTACACCCCGCTGGCACGCGGCCAGGCAGCGGCGCAGCCGGTGATCCAGGCGATCGCCCGCAAGCATGGCGTGCTGCCCAGCCAGGTGGCGCTGAAATGGCTGCTGGACCAGGACGGCGTGGCGGCGATCCCCa                                                          >  2:1073242/1‑150 (MQ=255)                                                                                          gTTGCTCACCGCCTACACCCCGCTGGCACGCGGCCAGGCAGCGGCGCAGCCGGTGATCCAGGCGATCGCCCGCAAGCATGGCGTGCTGCCCAGCCAGGTGGCGCTGAAATGGCTGCTGGACCAGGACGGCGTGGCGGCGATc                                                              >  1:246504/1‑142 (MQ=255)                                                                                          gTTGCTCACCGCCTACACCCCGCTGGCACGCGGCCAGGCAGCGGCGCAGCCGGTGATCCAGGCGATCGCCCGCAAGCATGGCGTGCTGCCCAGCCAGGTGGCGCTGAAATGGCTGCTGGACCAGGACGGCGTGGCGGCGATc                                                              <  2:246504/142‑1 (MQ=255)                                                                                           ttGCTCACCGCCTACA‑CCCGCTGGCACGCGGCCAGGCAGCCGAGCTGCTTGTGATCCATGCTATCGCCCGCAAGCATGGCGTGCTGCCCAGCCAGGTGGCGCTGAAATGGCTGCTGGACCAGGACGGCGTGGCGGCGATCCCCAAAGCca                                                    <  1:570876/150‑1 (MQ=255)                                                                                            tGCTCACCGCCTACACCCCGCTGGCACGCGGCCAGGCAGCGGCGCAGCCGGTGATCCAGGCGATCGCCCGCAAGCATGGCGTGCTGCCCAGCCAGGTGGCGCTGAAATGGCTGCTGGACCAGGACGGCGTGGCGGCGATCCCCAAAGCca                                                    <  1:1102109/150‑1 (MQ=255)                                                                                               tCACCGCCTACACCCCGCTGGCACGCGGCCAGGCAGCGGCGCAGCCGGTGATCCAGGCGATCGCCCGCAAGCATGGCGtgct                                                                                                                     <  2:756908/82‑1 (MQ=255)                                                                                               tCACCGCCTACACCCCGCTGGCACGCGGCCAGGCAGCGGCGCAGCCGGTGATCCAGGCGATCGCCCGCAAGCATGGCGtgct                                                                                                                     >  1:756908/1‑82 (MQ=255)                                                                                               tCACCGCCTACACCCCGCTGGCACGCGGCCAGGCAGCGGCGCAGCCGGTGATCCAGGCGATCGCCCGCAAGCATGGCGTGCTGCCCAGCCAGGTGGCGCTGAAATGGCTGCTGGACCAGGACGGCGTGGCGGCGATCCCCAAAGCcagca                                                 >  1:724946/1‑150 (MQ=255)                                                                                                 aCCGCCTACACCCCGCTGGCACGCGGCCAGGCAGCGGCGCAGCCGGTGATCCAGGCGATCGCCCGCAAGCATGGCGTGCTGCCCAGCCAGGTGGCGCTGAAATGGCTGCTGGACCAGGACGGCGTGGCGGCGATCCCCAAAGCCagcagc                                               >  1:206378/1‑150 (MQ=255)                                                                                                  ccGCCTACACCCCGCTGGCACGCGGCCAGGCAGCGGCGCAGCCGGTGATCCAGGCGATCGCCCGCAAGCATGGCGTGCTGCCCAGCCAGGTGGCGCTGAAATgg                                                                                            >  2:1126585/1‑104 (MQ=255)                                                                                                  ccGCCTACACCCCGCTGGCACGCGGCCAGGCAGCGGCGCAGCCGGTGATCCAGGCGATCGCCCGCAAGCATGGCGTGCTGCCCAGCCAGGTGGCGCTGAAATgg                                                                                            <  1:1126585/104‑1 (MQ=255)                                                                                                      cTACACCCCGCTGGCACGCGGCCAGGCAGCGGCGCAGCCGGTGATCCAGGCGATCGCCCGCAAGCATGGCGTGCTGCCCAGCCAGGTGGCGCTGAAATGGCTGCTGGACCa                                                                                 >  1:383740/1‑111 (MQ=255)                                                                                                      cTACACCCCGCTGGCACGCGGCCAGGCAGCGGCGCAGCCGGTGATCCAGGCGATCGCCCGCAAGCATGGCGTGCTGCCCAGCCAGGTGGCGCTGAAATGGCTGCTGGACCa                                                                                 <  2:383740/111‑1 (MQ=255)                                                                                                        acacCCCGCTGGCACGCGGCCAGGCAGCGGCGCAGCCGGTGATCCAGGCGATCGCCCGCAAGCATGGCGTGCTGCccagccag                                                                                                           >  2:1048583/1‑83 (MQ=255)                                                                                                        acacCCCGCTGGCACGCGGCCAGGCAGCGGCGCAGCCGGTGATCCAGGCGATCGCCCGCAAGCATGGCGTGCTGCccagccag                                                                                                           <  1:1048583/83‑1 (MQ=255)                                                                                                        acacCCCGCTGGCACGCGGCCAGGCAGCGGCGCAGCCGGTGATCCAGGCGATCGCCCGCAAGCATGGCGTGCTGCCCAGCCAGGTGGCGCTGAAATGGCTGCTGGACCAGGACGGCGTggcg                                                                    >  1:533624/1‑122 (MQ=255)                                                                                                        acacCCCGCTGGCACGCGGCCAGGCAGCGGCGCAGCCGGTGATCCAGGCGATCGCCCGCAAGCATGGCGTGCTGCCCAGCCAGGTGGCGCTGAAATGGCTGCTGGACCAGGACGGCGTggcg                                                                    <  2:533624/122‑1 (MQ=255)                                                                                                            cccGCTGGCACGCGGCCAGGCAGCGGCGCAGCCGGTGATCCAGGCGATCGCCCGCAAGCATGGCGTGCTGCGCAGCCAGGTGGCGCTGAAATGGCTGCTGGACCAGGACGGCGTGGCGGCGATCCCAAAAGCCAGCAGCCGCGAAAACCa                                    >  1:79469/1‑150 (MQ=255)                                                                                                            cccGCTGGCACGCGGCCAGGCAGCGGCGCAGCCGGTGATCCAGGCGATCGCCCGCAAGCATGGCGTGCTGCCCAGCCAGGTGGCGCTGAAATGgctgct                                                                                       >  1:753990/1‑99 (MQ=255)                                                                                                            cccGCTGGCACGCGGCCAGGCAGCGGCGCAGCCGGTGATCCAGGCGATCGCCCGCAAGCATGGCGTGCTGCCCAGCCAGGTGGCGCTGAAATGgctgct                                                                                       <  2:753990/99‑1 (MQ=255)                                                                                                            cccGCTGGCACGCGGCCAGGCAGCGGCGCAGCCGGTGATCCAGGCGATCGCCCGCAAGCATGGCGTGCTGCCCAGCCAGGTGGCGCTGAAATGGCTGCTGGACCAGGACGGCGTGGCGGCGATccc                                                            <  1:1117181/126‑1 (MQ=255)                                                                                                            cccGCTGGCACGCGGCCAGGCAGCGGCGCAGCCGGTGATCCAGGCGATCGCCCGCAAGCATGGCGTGCTGCCCAGCCAGGTGGCGCTGAAATGGCTGCTGGACCAGGACGGCGTGGCGGCGATccc                                                            >  2:1117181/1‑126 (MQ=255)                                                                                                            cccGCTGGCACGCGGCCAGGCAGCGGCGCAGCCGGTGATCCAGGCGATCGCCCGCAAGCATGGCGTGCTGCCCAGCCAGGTGGCGCTGAAATGGCTGCTGGACCAGGACGGCGTGGCGGCGATCCCCAAAGCca                                                    >  1:618508/1‑134 (MQ=255)                                                                                                            cccGCTGGCACGCGGCCAGGCAGCGGCGCAGCCGGTGATCCAGGCGATCGCCCGCAAGCATGGCGTGCTGCCCAGCCAGGTGGCGCTGAAATGGCTGCTGGACCAGGACGGCGTGGCGGCGATCCCCAAAGCca                                                    <  2:618508/134‑1 (MQ=255)                                                                                                            cccGCTGGCACGCGGCCAGGCAGCGGCGCAGCCGGTGATCCAGGCGATCGCCCGCAAGCATGGCGTGCTGCCCAGCCAGGTGGCGCTGAAATGGCTGCTGGACCAGGACGGCGTGGCGGCGATCCCCAAAGCCAGCAGCCGCGAAAACCa                                    <  2:682457/150‑1 (MQ=255)                                                                                                                    cACGCGGCCAGGCAGCGGCGCAGCCGGTGATCCAGGCGATCGCCCGCAAGCATGGCGTGCTGCCCAGCCAGGTGGCGCTGAAATGGCTGCTGGACCAGGACGGCGTGGCGGCGATCCCAAAAGCCAGCAGCCCCTGAAACCAGCTGGCCa                            >  2:560418/1‑150 (MQ=255)                                                                                                                     aCGCGGCCAGGCAGCGGCGCAGCCGGTGATCCAggccgcgtgccagcggggtgtaggcggtgagcaacaggtcatggccgcgggcgaagtcgagtagcggttgctggctgagcagaacatggtattccacctgaattgccgacaacggcg                           <  1:859850‑M2/150‑118 (MQ=255)                                                                                                                       gcgGCCAGGCAGCGGCGCAGCCGGTGATCCAggccgcgtgccagcggggtgtaggcggtgagcaacaggtcatgg                                                                                                    >  1:1001152‑M2/1‑31 (MQ=255)                                                                                                                       gcgGCCAGGCAGCGGCGCAGCCGGTGATCCAggccgcgtgccagcggggtgtaggcggtgagcaacaggtcatgg                                                                                                    <  2:1001152‑M2/75‑45 (MQ=255)                                                                                                                         ggCCAGGCAGCGGCGCAGCCGGTGATCCAGGCGATCGCCCGCAAGCATGGCGTGCTGCCCAGCCAGGTGGCGCTGAAATGGCTGCTGGACCAGGACGGCGTGGCGGCGATCCCCAAAGCCAGCAGCCGCGAAAACCAGCTGGCCAACCTg                       <  2:377670/150‑1 (MQ=255)                                                                                                                         ggCCAGGCAGCGGCGCAGCCGGTGATCCAGGCGATCGCCCGCAAGCATGGCGTGCTGCCCAGCCAGGTGGCGCTGAAATGGCTGCTGGACCAGGACGGCGTGGCGGCGATCCCCAAAGCCAGCAGCCGCGAAAACCAGCTGGCCAACCTg                       <  2:759153/150‑1 (MQ=255)                                                                                                                            cAGGCAGCGGCGCAGCCGGTGATCCAGGCGATCGCCCGCAAGCATGGCGTGCTGCCCAGCCAGGTGGCGCTGAAATGGCTGCTGGACCAGGACGGCGTGGCGGCGATCCCCAAAGCCAGCAGCCGCGAAAACCAGCTGGCCAACCTGGcc                    >  2:758407/1‑150 (MQ=255)                                                                                                                                cAGCGGCGCAGCCGGTGATCCAGGCGATTTCCCGCAAGCATGGCGTGCTGCCCAGCCAGGTGGCGCTGAAATGGCTGCTGGACCAGGACGGCGTGGCGGCGATCCCCAAAGCCAGCAGCCGCGAAAACCAGCTGGCCAACCTGGCCGCGt                <  1:691/150‑1 (MQ=255)                                                                                                                                       ggaGCCGGTGATCCAGGCGATCGCCCGCAAGCATGGCGTGCTGCCCAGCCAGGTGGCGCTGAAATGGCTGCTGGACCAGGACGGCGTGGCGGCGATCCCCAAAGCCAGCAGCCGCGAAAACCAGCTGGCCAACCTGGCCGCGTTGACCGt         <  1:773094/148‑1 (MQ=255)                                                                                                                                            cGGTGATCCAGGCGATCGCCCGCAAGCATGGCGTGCTGCCCAGCCAGGTGGCGCTGAAATGGCTGCTGGACCAGGACGGCGTGGCGGCGATCCCCAAAGCCAGCAGCCGCGAAAACCAGCTGGCCAACCTGGCCGCGTTGACCGTGCCGc    <  1:851983/150‑1 (MQ=255)                                                                                                                                              gTGATCCAGGCGATCGCCCGCAAGCATGGCGTGCTGCCCAGCCAGGTGGCGCTGAAATGGCTGCTGGACCAGGACGGCGTGGCGGCGATCCCCAAAGCCAGCAGCCGCGAAAACCAGCTGGCCAACCTGGCCGCGTTGACCGTGCCGCTg  >  2:320462/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                               CTGGGCGCGCCGTTGTCGGCAATTCAGGTGGAATACCATGTTCTGCTCAGCCAGCAACCGCTACTCGACTTCGCCCGCGGCCATGACCTGTTGCTCACCGCCTACACCCCGCTGGCACGCGGCCAGGCAGCGGCGCAGCCGGTGATCCAGGCGATCGCCCGCAAGCATGGCGTGCTGCCCAGCCAGGTGGCGCTGAAATGGCTGCTGGACCAGGACGGCGTGGCGGCGATCCCCAAAGCCAGCAGCCGCGAAAACCAGCTGGCCAACCTGGCCGCGTTGACCGTGCCGCTG  >  NC_002947/4172858‑4173148

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