New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	3680591 =	55 (1.110)	3 (0.060)	3/250	NT	5.5%	coding (345/468 nt)	crtZ	beta‑carotene hydroxylase
	?	NC_002947	3680632 =	50 (1.040)					coding (386/468 nt)	crtZ	beta‑carotene hydroxylase

GCGCCTTTGCAGTGGGTGGGCGGTGGTGTGGCAGGCTATGGAGCACTGTATGTACTGGCCCACGACGGTTTCTTTCACCGGCACTGGCCGCGCAAGCCGCGGCCGGTCAACCGC‑TACCTGAAACGCTTGCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACC  >  NC_002947/3680460‑3680738                                                                                                                                     |                                                                                                                                                    gcgcCTTTGCAGTGGGTGGGCGGTGGTGTGGCAGGCTATGGAGCACTGTATGTACTGGCCCACGACGGTTTCTTTCACCGGCACTGGCCGCGCAAGCCGCGGCCGGTCAACCGC‑TACCTGAAACGCTTGCACCGCGCGCACCGCTTGc                                                                                                                                     <  2:287183/148‑1 (MQ=255) gcgcCTTTGCAGTGGGTGGGCGGTGGTGTGGCAGGCTATGGAGCACTGTATGTACTGGCCCACGACGGTTTCTTTCACCGGCACTGGCCGCGCAAGCCGCGGCCGGTCAACCGC‑TACCTGAAACGCTTGCACCGCGCGCACCGCTTGCa                                                                                                                                    >  1:14498/1‑149 (MQ=255) gcgcCTTTGCAGTGGGTGGGCGGTGGTGTGGCAGGCTATGGAGCACTGTATGTACTGGCCCACGACGGTTTCTTTCACCGGCACTGGCCGCGCAAGCCGCGGCCGGTCAACCGC‑TACCTGAAACGCTTGCACCGCGCGCACCGCTTGCa                                                                                                                                    >  8:278317/1‑149 (MQ=255)             tgggtgggCGGTGGTGTGGCAGGCTATGGAGCACTGTATGTACTGGCCCACGACGGTTTCTTTCACCGGCACTGGCCGCGCAAGCCGCGGCCGGTCAACCGC‑TACCTGAAACGCTTGCACCGCGCGCACCGCTTGCACCATGCGGTGaa                                                                                                                        <  2:53432/149‑1 (MQ=255)                     cGGTGGTGTGGCAGGCTATGGAGCACTGTATGTACTGGCCCACGACGGTTTCTTTCACCGGCACTGGCCGCGCAAGCCGCGGCCGGTCAACCGC‑TACCTGAAACGCTTGCACCGCGCGCACCGCTTGCCCCATGCGGTGAAGGGGCGCa                                                                                                                >  8:169192/1‑149 (MQ=255)                       gtggtgTGGCAGGCTATGGAGCACTGTATGTACTGGCCCACGACGGTTTCTTTCACCGGCACTGGCCGCGCAAGCCGCGGCCGGTCAACCGC‑TACCTGAAACGCTTGCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACg                                                                                                              <  4:27649/149‑1 (MQ=255)                        tggtgTGGCAGGCTATGGAGCACTGTATGTACTGGCCCACGACGGTTTCTTTCACCGGCACTGGCCGCGCAAGCCGCGGCCGGTCAACCGC‑TACCTGAAACGCTTGCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACgg                                                                                                             >  7:124840/1‑149 (MQ=255)                         aggtgccgggcgtgccagggatgcgtcacgattggctgcgcctgctgcgcaattgctgcttcaacaccttcagcggcggcgcatagaagaacccgaacgacacgctccCCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACgg                                                                                                             >  4:283606‑M1/109‑149 (MQ=255)                           tgtgGCAGGCTATGGAGCACTGTCTGTACTGGCCCACGACGGTTTCTTTCACCGGCACTGGCCGCGCAAGCCGCGGCCGGTCAACCGCTTCCCTGAAACGCTTGCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACgggg                                                                                                           <  1:90822/149‑1 (MQ=255)                               gCAGGCTATGGAGCACTGTATGTACTCGCCCACGACGGTTTCTTTCACCGGCACTGGCCGCGCAAGCCGCGGCCGGTCAACCGC‑TACCTGAAACGCTTGCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCgt                                                                                                      <  7:278317/149‑1 (MQ=255)                                     tATGGAGCACTGTATGTACTGGCCCACGACGGTTTCTTTCACCGGCACTGGCCGCGCAAGCCGCGGCCGGTCAACCGC‑TACCTGAAACGCTTGCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGtt                                                                                                >  3:76257/1‑149 (MQ=255)                                      aTGGAGCACTGTATGTACTGGCCCACGACGGTTTCTTTCACCGGCACTGGCCGCGCAAGCCGCGGCCGGTCAACCGC‑TACCTGAAACGCTTGCACCgcgcgc                                                                                                                                              >  6:225215/1‑102 (MQ=255)                                      aTGGAGCACTGTATGTACTGGCCCACGACGGTTTCTTTCACCGGCACTGGCCGCGCAAGCCGCGGCCGGTCAACCGC‑TACCTGAAACGCTTGCACCgcgcgc                                                                                                                                              <  5:225215/102‑1 (MQ=255)                                      aTGGAGCACTGTATGTACTGGCCCACGACGGTTTCTTTCACCGGCACTGGCCGCGCAAGCCGCGGCCGGTCAACCGC‑TACCTGAAACGCTTGCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTc                                                                                               >  2:117552/1‑149 (MQ=255)                                              cTGTATGTACTGGCCCACGACGGTTTCTTTCACCGGCACTGGCCGCGCAAGCCGCGGCCGGTCAACCGC‑TACCTGAAACGCTTGCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGttctt                                                                                       <  1:117552/149‑1 (MQ=255)                                               tgtatgtaCTGGCCCACGACGGTTTCTTTCACCGGCACTGGCCGCGCAAGCCGCGGCCGGTCCACCGC‑TACCTGAAACGCTTGCACCGCGCGCCCCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGttcttc                                                                                      >  8:66153/1‑149 (MQ=255)                                               tgtatgtaCTGGCCCACGACGGTTTCTTTCACCGGCACTGGCCGCGCAAGCCGCGGCCGGTCAACCGC‑TACCTGAAACGCTTGCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGgc                                                                                                                   >  6:266009/1‑120 (MQ=255)                                               tgtatgtaCTGGCCCACGACGGTTTCTTTCACCGGCACTGGCCGCGCAAGCCGCGGCCGGTCAACCGC‑TACCTGAAACGCTTGCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGgc                                                                                                                   <  5:266009/120‑1 (MQ=255)                                               tgtatgtaCTGGCCCACGACGGTTTCTTTCACCGGCACTGGCCGCGCAAGCCGCGGCCGGTCAACCGC‑TACCTGAAACGCTTGCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCgg                                                                                             >  5:38483/1‑142 (MQ=255)                                               tgtatgtaCTGGCCCACGACGGTTTCTTTCACCGGCACTGGCCGCGCAAGCCGCGGCCGGTCAACCGC‑TACCTGAAACGCTTGCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCgg                                                                                             <  6:38483/142‑1 (MQ=255)                                                 tatgtaCTGGCCCACGACGGTTTCTTTCACCGGCACTGGCCGCGCAAGCCGCGGCCGGTCAACCGC‑TACCTGAAACGCTTGCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTtcttct                                                                                     <  8:124840/148‑1 (MQ=255)                                                  atgtaCTGGCCCACGACGGTTTCTTTCACCGGCACTGGCCGCGCAAGCCGCGGCCGGTCAACCGC‑TACCTGAAACGCTTGCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTAt                                                                                   >  3:15238/1‑149 (MQ=255)                                                   tgtaCTGGCCCACGACGGTTTCTTTCACCGGCACTGGCCGCGCAAGCCGCGGCCGGTCAACCGC‑TACCTGAAACGCTTGCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTAt                                                                                   <  7:182020/148‑1 (MQ=255)                                                     taCTGGCCCACGACGGTTTCTTTCACCGGCACTGGCCGCGCAAGCCGCGGCCGGTCAACCGC‑TACCTGAAACGCTTGCACCGCGCGCACCGCTTGCACCATGCGGTGAAggg                                                                                                                     >  5:287994/1‑112 (MQ=255)                                                     taCTGGCCCACGACGGTTTCTTTCACCGGCACTGGCCGCGCAAGCCGCGGCCGGTCAACCGC‑TACCTGAAACGCTTGCACCGCGCGCACCGCTTGCACCATGCGGTGAAggg                                                                                                                     <  6:287994/112‑1 (MQ=255)                                                        tGGCCCACGACGGTTTCTTTCACCGGCACTGGCCGCGCAAGCCGCGGCCGGTCAACCGC‑TACCTGAAACGCTTGCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGcgccg                                                                             >  5:61792/1‑149 (MQ=255)                                                           cccACGACGGTTTCTTTCACCGGCACTGGCCGCGCAAGCCGCGGCCGGTCAACCGC‑TACCTGAAACGCTTGCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGcgccgccg                                                                          >  4:246010/1‑149 (MQ=255)                                                             cACGACGGTTTCTTTCACCGGCACTGGCCGCGCAAGCCGCGGCCGGTCAACCGC‑TACCTGAAACGCTTGCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGgc                                                                                                                   <  7:28975/106‑1 (MQ=255)                                                             cACGACGGTTTCTTTCACCGGCACTGGCCGCGCAAGCCGCGGCCGGTCAACCGC‑TACCTGAAACGCTTGCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGgc                                                                                                                   >  8:28975/1‑106 (MQ=255)                                                                   ggTTTCTTTCACCGGCACTGGCCGCGCAAGCCGCGGCCGGTCAACCGC‑TACCTGAAACGCTTGCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGgt                                                                  <  4:76257/149‑1 (MQ=255)                                                                    gTTTCTTTCACCGGCACTGGCCGCGCAAGCCGCGGCCGGTCAACCGC‑TACCTGAAACGCTAGCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGgtg                                                                 >  5:309377/1‑149 (MQ=255)                                                                     tttctttcACCGGCACTGGCCGCGCAAGCCGCGGCCGGTCAACCGC‑TACCTGAAACGCTTGCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATgcgc                                                                               >  5:36882/1‑134 (MQ=255)                                                                     tttctttcACCGGCACTGGCCGCGCAAGCCGCGGCCGGTCAACCGC‑TACCTGAAACGCTTGCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATgcgc                                                                               <  6:36882/134‑1 (MQ=255)                                                                     tttctttcACCGGCACTGGCCGCGCAAGCCGCGGCCGGTCAACCGC‑TACCTGAAACGCTTGCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGaa                                                                     >  6:318922/1‑144 (MQ=255)                                                                     tttctttcACCGGCACTGGCCGCGCAAGCCGCGGCCGGTCAACCGC‑TACCTGAAACGCTTGCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGaa                                                                     <  5:318922/144‑1 (MQ=255)                                                                      ttctttcACCGGCACTGGCCGCGCAAGCCGCGGCCGGTCAACCGC‑TACCTGAAACGCTTGCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGtt                                                               >  3:56493/1‑149 (MQ=255)                                                                       tctttcACCGGCACTGGCCGCGCAAGCCGCGGCCGGTCAACCGC‑TACCTGAAACGCTTGCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGtt                                                               <  6:309377/148‑1 (MQ=255)                                                                         tttcACCGGCACTGGCCGCGCAAGCCGCGGCCGGTCAACCGC‑TACCTGAAACGCTTGCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTg                                                              >  7:274577/1‑147 (MQ=255)                                                                             aCCGGCACTGGCCGCGCAAGCCGCGGCCGGTCAACCGC‑TACCTGAAACGCTTGCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAagca                                                         <  4:15238/148‑1 (MQ=255)                                                                                 tgcttcaacaccttcagcggcggcgcatagaagaacccgaacgacacgctccCCGCGCGCACCGCTTGCACCATGCGGTGAAggg                                                                                                                     <  8:187549‑M1/33‑1 (MQ=255)                                                                                 tgcttcaacaccttcagcggcggcgcatagaagaacccgaacgacacgctccCCGCGCGCACCGCTTGCACCATGCGGTGAAggg                                                                                                                     >  7:187549‑M1/53‑85 (MQ=255)                                                                                           cgcAAGCCGCGGCCGGTCAACCGC‑TACCTGAAACGCTTGCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGcagcag                                          <  1:53454/149‑1 (MQ=255)                                                                                           cgcAAGCCGCGGCCGGTCAACCGC‑TACCTGAAACGCTTGCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGcagcag                                          <  8:274577/149‑1 (MQ=255)                                                                                           cgcAAGCCGCGGCCGGTCAACCGC‑TACCTGAAACGCTTGCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGcagcag                                          <  3:271984/149‑1 (MQ=255)                                                                                                 ccGCGGCCGGTCAACCGC‑TACCTGAAACGCTTGCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAg                                    <  2:14498/149‑1 (MQ=255)                                                                                                       ccGGTCAACCGC‑TACCTGAAACGCTTGCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATc                              <  3:246010/149‑1 (MQ=255)                                                                                                           tCAACCGC‑TACCTGAAACGCTTGCACTGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCAGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGa                          >  1:168639/1‑149 (MQ=255)                                                                                                           tCAACCGC‑TACCTGAAACGCTTGCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGa                          <  4:56493/149‑1 (MQ=255)                                                                                                                 gC‑TACCTGAAACGCTTGCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGAc                         >  7:210142/1‑144 (MQ=255)                                                                                                                 gC‑TACCTGAAACGCTTGCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGAc                         <  8:210142/144‑1 (MQ=255)                                                                                                                    tACCTGAAACGCTTGCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCgca                                              <  5:168625/121‑1 (MQ=255)                                                                                                                    tACCTGAAACGCTTGCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCgca                                              >  6:168625/1‑121 (MQ=255)                                                                                                                         gAAACGCTTGCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCa             >  6:36599/1‑149 (MQ=255)                                                                                                                          aaaCGCTTGCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCGGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGc              <  7:69599/147‑1 (MQ=255)                                                                                                                            aCGCTTGCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGtt                                                               >  7:272817/1‑96 (MQ=255)                                                                                                                            aCGCTTGCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGtt                                                               <  8:272817/96‑1 (MQ=255)                                                                                                                              gCTTGCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGTGCAGCCAATCGTGACGCATCCCTGGCACGccc        >  2:23656/1‑149 (MQ=255)                                                                                                                              gCTTGCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGccc        >  6:228217/1‑149 (MQ=255)                                                                                                                                    aCCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCAcc  >  4:154048/1‑149 (MQ=255)                                                                                                                                     |                                                                                                                                                    GCGCCTTTGCAGTGGGTGGGCGGTGGTGTGGCAGGCTATGGAGCACTGTATGTACTGGCCCACGACGGTTTCTTTCACCGGCACTGGCCGCGCAAGCCGCGGCCGGTCAACCGC‑TACCTGAAACGCTTGCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACC  >  NC_002947/3680460‑3680738

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.