New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	4284658 =	42 (0.930)	4 (0.090)	4/240	NT	8.8%	coding (19/849 nt)	hbd	3‑hydroxybutyryl‑CoA dehydrogenase
	?	NC_002947	4284718 =	42 (0.960)					coding (79/849 nt)	hbd	3‑hydroxybutyryl‑CoA dehydrogenase

GCGATCATCGCCACCAAGGCCATCCATGAGCTGCAGCGCATCCAGGGCCGCTACGCCCTGGCCACGATGTGTATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAGCGCCAGGTCAGC  >  NC_002947/4284510‑4284804                                                                                                                                                     |                                                                                                                                                   gCGATCATCGCCACCAAGGCCATCCATGAGCTGCAGCGCATCCAGGGCCGCTACGCCCTGGCCACGATGTGTATCGGCGGTGGCCAAGTCATCGCCGTCGTGTTCGAGCGCGTCGGAGGGAGGCTGACACATGAGCATTGAAAAGATCGc                                                                                                                                                   >  2:455531/1‑150 (MQ=255)             accaAGGTCATCGATGAGTTGCAGCGCAACCAGGGCCGCTCCGCACTGGCCACGATGTCTAGCGGGGGAGACCAAGGCATCGCCGTCGTGTTCGTGCGCGTCTCAGGGAGGCTGAGACATGACACTTGAAAGGATCGCGTGGATcgccgc                                                                                                                                       <  2:524721/150‑1 (MQ=255)             accaAGGCCATCCATGAGCTGCAGCGCATCCAGGGCCGCTACGCCCTGGCCACGATGTGTATCGGCGGTGGCCAAGGCCTCGCCGTCGTGTTCGACCGCGTCCCAAGGAGTCCGACACATGAGCATTGAAAAGATCGCCGTGATCGgcgc                                                                                                                                       <  1:209841/150‑1 (MQ=255)                  ggCCATCCATGAGCTGCAGCGCATCCAGGGCCGCTACGCCCTGGCCACGATGTGTATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCa                                                                                                                                  >  1:175867/1‑150 (MQ=255)                              gCTGCAGCGCATCCAGGGCCGCTACGCCCTGGCCACGATGTGTATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACg                                                                                                                      <  1:412030/150‑1 (MQ=255)                                        aTCCAGGGCCGCTACGCCCTGGCCACGATGTGTATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCa                                                                                                            <  2:549750/150‑1 (MQ=255)                                        aTCCAGGGCCGCTACGCCCTGGCCACGATGTGTATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCa                                                                                                            >  1:726370/1‑150 (MQ=255)                                           cAGGGCCGCTACGCCCTGGCCACGATGTGTATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGACCTCACCATGGGCAACGGCATTGCCCAGgt                                                                                                         <  1:822903/150‑1 (MQ=255)                                            aGGGCCGCTACGCCCTGGCCACGATGTGTATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGCAGGCTGACACATGAGCATTGAAAAGATCGCCGTGACCGGCGCGGGCAGCAGGGGCAACGGCATTGCCCAGgtg                                                                                                        >  1:681466/1‑150 (MQ=255)                                               gCCGCTACGCCCTGGCCACGAAGTGTATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTgc                                                                                                     <  2:395653/150‑1 (MQ=255)                                                   cTACGCCCTGGCCACGATGTGTATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCa                                                                                                 >  2:156008/1‑150 (MQ=255)                                                       gCCCTGGCCACGATGTGTATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAACAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGc                                                                                             >  2:1068841/1‑150 (MQ=255)                                                        gcccttgctgacctggcgctcgaggttcttgctcaaggtggccacaccgcgctcgagcgctgcgtcggaaacatccaccagcagcacctggtaGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGcc                                                                                            <  1:765854‑M1/57‑1 (MQ=255)                                                         cccttgctgacctggcgctcgaggttcttgctcaaggtggcaacaacgagatcgagcgctgcgtcggaaacagccaccagcagcacctggtaGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCg                                                                                           <  1:626617‑M1/58‑1 (MQ=255)                                                              ccACGATGTGTATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTAc                                                                                      >  1:1128285/1‑150 (MQ=255)                                                                tgacctggcgctcgaggttcttgctcaaggtggccacaccgcgctcgagcgctgcgtcggaaacatccaccagcagcacctggtaGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCa                                                                                    >  2:626617‑M1/86‑150 (MQ=255)                                                                   aTGTGTATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGAGCCATTGCCGGCTACCAGGt                                                                                 <  2:685236/150‑1 (MQ=255)                                                                     gtgtATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGtgc                                                                               >  1:361510/1‑150 (MQ=255)                                                                     gtgtATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGtgc                                                                               >  2:905306/1‑150 (MQ=255)                                                                        tATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGtgctgc                                                                            <  1:156008/150‑1 (MQ=255)                                                                           cggcggTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCtgg                                                                         <  2:492275/150‑1 (MQ=255)                                                                            ggcggTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCtggt                                                                        <  2:175867/150‑1 (MQ=255)                                                                            cgaggttcttgctcaaggtggccacaccgcgctcgagcgctgcgtcggaaacatccaccagcagcacctggtaGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCtggt                                                                        >  2:765854‑M1/74‑150 (MQ=255)                                                                                   gCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGtt                                                                 >  1:87959/1‑150 (MQ=255)                                                                                      aaGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGCGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCGGGTGGATGTTTcc                                                              <  2:566647/150‑1 (MQ=255)                                                                                            tCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCCGGCCTTGCTGGTGGATGTTTCCGACGCa                                                        <  1:1068841/150‑1 (MQ=255)                                                                                            ccGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCa                                                        <  2:902085/149‑1 (MQ=255)                                                                                                        tCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAgcgc                                            <  1:692481/150‑1 (MQ=255)                                                                                                         cGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGcgcg                                           <  1:507495/150‑1 (MQ=255)                                                                                                           aGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGgt                                         >  2:1139651/1‑150 (MQ=255)                                                                                                           aGCGCGTCTGAGGGAGGCCGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGGCGCAGCGCTCGAGCGCGgt                                         <  2:70927/150‑1 (MQ=255)                                                                                                              gcgTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGtgtg                                      <  2:957357/150‑1 (MQ=255)                                                                                                              gcgTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGtgtg                                      <  2:312751/150‑1 (MQ=255)                                                                                                                gTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGc                                    >  2:164665/1‑150 (MQ=255)                                                                                                                   tGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCAt                                                                                                                  <  2:438036/69‑1 (MQ=255)                                                                                                                   tGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCAt                                                                                                                  >  1:438036/1‑69 (MQ=255)                                                                                                                        gaggCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGa                            <  2:1128285/150‑1 (MQ=255)                                                                                                                        gaggCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGa                            <  2:1097947/150‑1 (MQ=255)                                                                                                                        gaggCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGa                            <  1:455531/150‑1 (MQ=255)                                                                                                                             tGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAg                       >  2:478813/1‑150 (MQ=255)                                                                                                                               acacATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCCACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGaa                     >  1:30090/1‑150 (MQ=255)                                                                                                                                   aTGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCATCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTc                 <  2:653543/150‑1 (MQ=255)                                                                                                                                       gCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCt                                                   <  1:490172/112‑1 (MQ=255)                                                                                                                                       gCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCt                                                   >  2:490172/1‑112 (MQ=255)                                                                                                                                       gCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAgc             <  2:361510/150‑1 (MQ=255)                                                                                                                                         aTTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAgcgc           >  1:502283/1‑150 (MQ=255)                                                                                                                                             aaaaGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAGCGCCAgg       >  1:344972/1‑150 (MQ=255)                                                                                                                                                  aTCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAGCGCCAGGTCAGc  <  2:726370/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   GCGATCATCGCCACCAAGGCCATCCATGAGCTGCAGCGCATCCAGGGCCGCTACGCCCTGGCCACGATGTGTATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAGCGCCAGGTCAGC  >  NC_002947/4284510‑4284804

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