New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 3387771	42 (0.710)	3 (0.050)	3/242	NT	6.9%	coding (552/795 nt)	PP_2989	short‑chain dehydrogenase/reductase family oxidoreductase
	?	NC_002947	= 3387797	42 (0.760)					coding (578/795 nt)	PP_2989	short‑chain dehydrogenase/reductase family oxidoreductase

GCCCGAGTTCGGCATGTCCATCTACGGCGCTACCAAAGCCTTCGTGCTGTTTCTGTCACAGGGGCTGAGTCAGGAGCTGTCACCCAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGTGATGGAGGTCGGTGAGCTGGTCGATGCTGCACTTGTCGGTTTCGACCGTCGCGAGTTGGTAACCATCCCGCCACTTCATGTTGCCGAGCGCTGGGATGCACTGGA  >  NC_002947/3387654‑3387935                                                                                                                                                |                                                                                                                                           gcCCGAGTTCGGCATGTCCATCTACGGCGCTACCAAAGCCTTCGTGCTGTTTCTGTCACAGGGGCTGAGTCAGGAGCTGTCACCCAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGtgccg                                                                                                                                       >  8:85697‑M2/1‑144 (MQ=255) gcCCGAGTTCGGCATGTCCATCTACGGCGCTACCAAAGCCTTCGTGCTGTTTCTGTCACAGGGGCTGAGTCAGGAGCTGTCACCCAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGtgccg                                                                                                                                       <  7:41213‑M2/149‑6 (MQ=255) gcCCGAGTTCGGCATGTCCATCTACGGCGCTACCAAAGCCTTCGTGCTGTTTCTGTCACAGGGGCTGAGTCAGGAGCTGTCACCCAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCAGCACCGAGATCTGGGCGCGtgccg                                                                                                                                       <  5:183443‑M2/149‑6 (MQ=255)   cggttttcggaaTGTCCATCTACGGCGCTTCCAAAGCCTTCGTGCTGTTTCTGTGACAGGGGCTGAGTCAGGAGCTGTCACCCAAGGGCGTGTACGTGCAATCTGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGTCGCGtgccg                                                                                                                                       >  6:192922‑M2/6‑142 (MQ=255)   ccGAGTTCGGCATGTCCATCTACGGCGCTACCAAAGCCTTCGTGCTGTTTCTGTCACAGGGGCTGAGTCAGGAGCTGTCACCCAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGtgccgga                                                                                                                                     <  5:192922‑M2/149‑8 (MQ=255)      aGTTCGGCATGTCCATCTACGGCGCTACCAAAGCCTTCGTGCTGTTTCTGTCACAGGGGCTGAGTCAGGAGCTGTCACCCAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGTGCGCGTGCCGGCAtt                                                                                                                                  <  6:183443/149‑1 (MQ=255)           ggCATGTCCATCTACGGCGCTACCAAAGCCTTCGTGCTGTTTCTGTCACAGGGGCTGAGTCAGGAGCTGTCACCCAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGt                                                                                                                             >  4:361995/1‑149 (MQ=255)            gCATGTCCATCTACGGCGCTACCAAAGCCTTCGTGCTGTTTCTGTCACAGGGGCTGAGTCAGGAGCTGTCACCCAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGa                                                                                                                                <  4:42005/145‑1 (MQ=255)            gCATGTCCATCTACGGCGCTACCAAAGCCTTCGTGCTGTTTCTGTCACAGGGGCTGAGTCAGGAGCTGTCACCCAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGa                                                                                                                                >  3:42005/1‑145 (MQ=255)              aTGTCCATCTACGGCGCTACCAAAGCCTTCGTGCTGTTTCTGTCACAGGGGCTGAGTCAGGAGCTGTCACCCAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCaa                                                                                                                          <  5:359060/149‑1 (MQ=255)              aTGTCCATCTACGGCGCTACCAAAGCCTTCGTGCTGTTTCTGTCACAGGGGCTGAGTCAGGAGCTGTCACCCAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCaa                                                                                                                          >  6:290404/1‑149 (MQ=255)                               tACCAAAGCCTTCGTGCTGTTTCTGTCACAGGGGCTGAGTCAGGAGCTGTCACCCAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAaca                                                                                                                        >  3:1114/1‑134 (MQ=255)                               tACCAAAGCCTTCGTGCTGTTTCTGTCACAGGGGCTGAGTCAGGAGCTGTCACCCAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAaca                                                                                                                        <  4:1114/134‑1 (MQ=255)                                 ccAAAGCCTTCGTGCTGTTTCTGTCACAGGGGCTGAGTCAGGAGCTGTCACCCAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGTGATg                                                                                                       >  4:278911/1‑149 (MQ=255)                                  cAAAGCCTTCGTGCTGTTTCTGTCACAGGGGCTGAGTCAGGAGCTGTCACCCAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGTGATgg                                                                                                      <  3:278911/149‑1 (MQ=255)                                   aaaGCCTTCGTGCTGTTTCTGTCACAGGGGCTGAGTCAGGAGCTGTCACCCAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGTGATg                                                                                                       >  1:222709/1‑147 (MQ=255)                                   aaaGCCTTCGTGCTGTTTCTGTCACAGGGGCTGAGTCAGGAGCTGTCACCCAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTACCGGCATTGACg                                                                                                                              <  6:58397/124‑1 (MQ=255)                                   aaaGCCTTCGTGCTGTTTCTGTCACAGGGGCTGAGTCAGGAGCTGTCACCCAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTACCGGCATTGACg                                                                                                                              >  5:58397/1‑124 (MQ=255)                                    aaGCCTTCGTGCTGTTTCTGTCACAGGGGCTGAGTCAGGAGCTGTCACCCAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGTGATGGAg                                                                                                    <  8:41213/149‑1 (MQ=255)                                     aGCCTTCGTGCTGTTTCTGTCACAGGGGCTGAGTCAGGAGCTGTCACCCAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGTGATGGAgg                                                                                                   >  2:99717/1‑149 (MQ=255)                                        cTTCGTGCTGTTTCTGTCACAGGGGCTGAGTCAGGAGCTGTCACCCAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGTGATGGAGGTCg                                                                                                <  3:361995/149‑1 (MQ=255)                                             tGCTGTTTCTGTCACAGGGGCTGAGTCAGGAGCTGTCACCCAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGTGATGGAGGTCGGTGAg                                                                                           <  1:99717/149‑1 (MQ=255)                                                  tttCTGTCACAGGGGCTGAGTCAGGAGCTGTCACCCAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGTGATGGAGGTCGGTGAGCTg                                                                                        <  6:66824/147‑1 (MQ=255)                                                        tCACAGGGGCTGAGTCAGGAGCTGTCACCCAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGTGATGGAGGTCGGTGAGCTGGt                                                                                      >  2:211191/1‑143 (MQ=255)                                                        tCACAGGGGCTGAGTCAGGAGCTGTCACCCAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGTGATGGAGGTCGGTGAGCTGGt                                                                                      <  1:211191/143‑1 (MQ=255)                                                        tCACAGGGGCTGAGTCAGGAGCTGTCACCCAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGTGATGGAGGTCGGTGAGCTGGTCGAt                                                                                  <  5:290404/147‑1 (MQ=255)                                                          acaGGGGCTGAGTCAGGAGCTGTCACCCAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGTGATGGAGGTCGGTGAg                                                                                           >  1:328017/1‑136 (MQ=255)                                                          acaGGGGCTGAGTCAGGAGCTGTCACCCAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGTGATGGAGGTCGGTGAg                                                                                           <  2:328017/136‑1 (MQ=255)                                                                   gAGTCAGGAGCTGTCACCCAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGTGATGGAGGTCGGTGAGCTGGTc                                                                                     <  3:154948/133‑1 (MQ=255)                                                                   gAGTCAGGAGCTGTCACCCAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGTGATGGAGGTCGGTGAGCTGGTc                                                                                     >  4:154948/1‑133 (MQ=255)                                                                      tCAGGAGCTGTCACCCAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGTGATGGAGGTCGGTGAGCTGGTCGAtgc                                                                                >  3:246615/1‑135 (MQ=255)                                                                      tCAGGAGCTGTCACCCAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGTGATGGAGGTCGGTGAGCTGGTCGAtgc                                                                                <  4:246615/135‑1 (MQ=255)                                                                      tCAGGAGCTGTCACCCAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGTGATGGAGGTCGGTGAGCTGGTCGATGCTGCACTTGTCg                                                                     <  2:9066/146‑1 (MQ=255)                                                                      tCAGGAGCTGTCACCCAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGTGATGGAGGTCGGTGAGCTGGTCGATGCTGCACTTGTCg                                                                     >  1:9066/1‑146 (MQ=255)                                                                             cTGTCACCCAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGTGATGGAGGTCGGTGAGCTGGTCGATGCTGCACTTGTCGGTTTCGACCg                                                           >  5:127242/1‑149 (MQ=255)                                                                               gTCACCCAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGTGATGGAGGTCGGTGAGCTGGTCGATGCTGCACTTGTCGGTTTCGACCGTc                                                         >  2:126389/1‑149 (MQ=255)                                                                                tCACCCAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGTGATGGAGGTCGGTGAGCTGGTCGATGCTGCACTTGTCGGTTTCGACCGTcg                                                        >  7:18540/1‑149 (MQ=255)                                                                                 cACCCAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGTGATGGAGGTCGGTGAGCTGGTCGATGCTGCACTTGTCGGTTTCGACCGTcgc                                                       >  7:130263/1‑149 (MQ=255)                                                                                   cccAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGTGATGGAGGTCGGTGAGCTGGTCGATGCTGCACTTGTCGGTTTCGACCGTCGCGa                                                     <  6:127242/149‑1 (MQ=255)                                                                                    ccAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGTGATGGAGGTCGGTGAGCTGGTCGATGCTGCACTTGTCGGTTTCGACCGTCGCGAg                                                    >  5:253844/1‑149 (MQ=255)                                                                                           cGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGTGATGGAGGTCGGTGAGCTGGTCGATGCTGCACTTGTCGGTTTCGACCGTCGCGAGTTGGTaa                                             <  8:130263/149‑1 (MQ=255)                                                                                               tACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGt                                                                                                           >  3:134446/1‑83 (MQ=255)                                                                                               tACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGt                                                                                                           <  4:134446/83‑1 (MQ=255)                                                                                                    gaaaGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGTGATGGAGGTCGGTGAGCTGGTCGATGCTGCACTTGTCGGTTTCGACCGTCGCGAGTTGGTAACCATCCCGc                                    >  7:176462/3‑149 (MQ=255)                                                                                                    gCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGtgccggcagcaccgcttgcacgtacacgcccttgggtgacagctcctgactcagcccctgtgacagaaacagcacgaa                                                              <  7:343663‑M2/123‑79 (MQ=255)                                                                                                    gCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGtgccggcagcaccgcttgcacgtacacgcccttgggtgacagctcctgactcagcccctgtgacagaaacagcacgaa                                                              >  8:343663‑M2/1‑45 (MQ=255)                                                                                                           gTGCTGCCGGCAGCGACCCGCACAGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGTGGTGATGGAGGTCGGTGAGCTGGTCGATGCTGCACTTGTCGGTTTCGACCGTCGCGAGTTGGTAACCATCCCGCCACTTCa                             <  8:231628/149‑1 (MQ=255)                                                                                                                  cGGCAGCGACCCGCACCGAGATCTGGGCGCGtgccggcagcaccgcttgcacgtacacgcccttgggtgacagctcctgactcagcccctgtgacagaaacagcacgaaggctttggtagcgccgtagatggacatgccgaactcgggc                      <  3:126245‑M2/149‑119 (MQ=255)                                                                                                                       gCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGTGATGGAGGTCGGTGAGCTGGTCGATGCTGCACTTGTCGGTTTCGACCGTCGCGAGTTGGTAACCATCCCGCCACTTCATGTTGCCGAgcg                 >  2:1408/1‑149 (MQ=255)                                                                                                                          aCCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGTGATGGAGGTCGGTGAGCTGGTCGATGCTGCACTTGTCGGTTTCGACCGTCGCGAGTTGGTAACCATCCCGCCACTTCATGTTGCCGAGCGCTg              <  6:253844/149‑1 (MQ=255)                                                                                                                                           ggCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGTGATGGAGGTCGGTGAGCTGGTCGATGCTGCACTTGTCGGTTTCGACCGTCGCGAGTTGGTAACCATCCCGCCACTTCATGTTGCCGAGCGCTgg             <  7:154452/133‑1 (MQ=255)                                                                                                                                           ggCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGTGATGGAGGTCGGTGAGCTGGTCGATGCTGCACTTGTCGGTTTCGACCGTCGCGAGTTGGTAACCATCCCGCCACTTCATGTTGCCGAGCGCTgg             >  8:154452/1‑133 (MQ=255)                                                                                                                                           ggCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGTGATGGAGGTCGGTGAGCTGGTCGATGCTGCACTTGTCGGTTTCGACCGTCGCGAGTTGGTAACCATCCCGCCACTTCATGTTGCCGAGCGCTGGGATGCACTGGa  <  1:1408/144‑1 (MQ=255)                                                                                                                                                |                                                                                                                                           GCCCGAGTTCGGCATGTCCATCTACGGCGCTACCAAAGCCTTCGTGCTGTTTCTGTCACAGGGGCTGAGTCAGGAGCTGTCACCCAAGGGCGTGTACGTGCAAGCGGTGCTGCCGGCAGCGACCCGCACCGAGATCTGGGCGCGTGCCGGCATTGACGTCAACACACTCCCGGAGGTGATGGAGGTCGGTGAGCTGGTCGATGCTGCACTTGTCGGTTTCGACCGTCGCGAGTTGGTAACCATCCCGCCACTTCATGTTGCCGAGCGCTGGGATGCACTGGA  >  NC_002947/3387654‑3387935

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.