New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	5213292 =	68 (1.150)	3 (0.050)	3/238	NT	5.3%	coding (420/1182 nt)	PP_4594	cystathionine gamma‑synthase
	?	NC_002947	5213361 =	44 (0.800)					coding (351/1182 nt)	PP_4594	cystathionine gamma‑synthase

ACACCGCCTTGATGTCGACCAGGTGCAGCAACGGGTTGGTCGGCGACTCGATCCAGATCAGCTGGGTTTCC‑GGCTTGATCGCAGCCGCTACGGCTTCAATGTCGTTGAGATCGACATAGGTGGTGGTCAGCCCCGAGGTGCGGCTGCGGTAATCTTCCATGATGCGGAACGTGCCACCGTACACGCCGTTCATCACCACCACGTGGGCATCCTTGGGCAGCAGCTCAAGCACGGTTGCGGTGGCATTGACTCCCGAGGCGCAGGCCACAGCGCCCACGCCTTCTTCCAAGGCAGCGA  >  NC_002947/5213213‑5213509                                                                                                                                                      |                                                                                                                                                     acacCGCCTTGATGTCGACCAGGTGCAGCAACGGGTTGGTCGGCGACTCGATCCAGATCAGCTGGGTTTCC‑GGCTTGATCGCAGCCGCTACGGCTTCAATGTCGTTGAGATCGACATAGGTGGTGGTCAGCCCCGAGGTGCGGCTGCgg                                                                                                                                                      >  2:181370/1‑149 (MQ=255)   gagctgctgcccaaggatgcccacgtggtggtgatgaacggcgtgtacggtggcacgttccgcatcatggaagattacCGCAGCCGCTACGGCTTCAATGTCGTTGAGATCGACATAGGTGGTGGTCAGCCCCGAGGTGCGGCTGCgg                                                                                                                                                      >  4:235977‑M1/79‑148 (MQ=255)        cTTGATGTCGACCAGGTGCAGCAACGGGTTGGTCGGCGACTCGATCCAGATCAGCTGGGTTTCC‑GGCTTGATCGCAGCCGCTACGGCTTCAATGTCGTTGAGATCGACATAGGTGGTGGTCAGCCCCGAGGTGCGGCTGCGGTAATCtt                                                                                                                                               >  5:5593/1‑149 (MQ=255)           gATGTCGACCAGGTGCAGCAACGGGTTGGTCGGCGACTCGATCCAGATCAGCTGGGTTTCC‑GGCTTGATCGCAGCCGCTACGGCTTCAATGTCGTTGAGATCGACATAGGTGGTGGTCAGCCCCGAGGTGCGGCTGCGGTAATCTTcc                                                                                                                                             <  3:150129/148‑1 (MQ=255)           gATGTCGACCAGGTGCAGCAACGGGTTGGTCGGCGACTCGATCCAGATCAGCTGGGTTTCC‑GGCTTGATCGCAGCCGCTACGGCTTCAATGTCGTTGAGATCGACATAGGTGGTGGTCAGCCCCGAGGTGCGGCTGCGGTAATCTTCCa                                                                                                                                            >  2:261964/1‑149 (MQ=255)               tCGACCAGGTGCAGCAACGGGTTGGTCGGCGACTCGATCCAGATCAGCTGGGTTTCC‑GGCTTGATCGCAGCCGCTACGGCTTCAATGTCGTTGAGATCGACATAGGTGGTGGTCAGCCCCGAGGTGCGGCTGCGGTAATCTTCCatgat                                                                                                                                        >  6:146487/1‑149 (MQ=255)                 gACCAGGTGCAGCAACGGGTTGGTCGGCGACTCGATCCAGATCAGCTGGGTTTCC‑GGCTTGATCGCAGCCGCTACGGCTTCAATGTCGTTGAGATCGACATAGGTGGTGGTCTGCCCCGAGGTGCGGCTGCGGTAATCTTCCATGATGc                                                                                                                                      >  8:55354/1‑149 (MQ=255)                      tgTGCAGCAACGGGTTGGTCGGCGACTCGATCCAGATCAGCTGGGTTTCC‑GGCTTGATCGCAGCCGCTACGGCTTCAATGTCGTTGAGATCGACATAGGTGGTGGTCAGCCCCGAGGTGCGGCTGCGGTAATCTTCCATGATGCGGAAc                                                                                                                                 <  6:5593/148‑1 (MQ=255)                       gTGCAGCAACGGGTTGGTCGGCGACTCGATCCAGATCAGCTGGGTTTCC‑GGCTTGATCGCAGCCGCTACGGCTTCAATGTCGTTGAGATCGACATAGGTGGTGGTCAGCCCCGAGGTGCGGCTGCGGTAATCTTCCATGATGCGGAACg                                                                                                                                >  4:356647/1‑149 (MQ=255)                            gcaACGGGTTGGTCGGCGACTCGATCCAGATCAGCTGGGTTTCC‑GGCTTGATCGCAGCCGCTACCGCTTCAATGTCGTTGAGATCGACATAGGAGGTGGTCAGCCCCGAGCTGCGGCTGCGGTAATCTTCCATGATGCGGAACGTGCCa                                                                                                                           >  4:43203/1‑149 (MQ=255)                                  ggTTGGTCGGCGACTCGATCCAGATCAGCTGGGTTTCC‑GGCTTGATCGCAGCCGCTACGGCTTCAATGTCGTTGAGATCGACATAGGTGGTGGTCAGCCCCGAGGTGCGGCTGCGGTAATCTTCCATGATGCGGAACGTGCCACCGTac                                                                                                                     <  1:221279/149‑1 (MQ=255)                                        tCGGCGACTCGATCCAGATCAGCTGGGTTTCCGGGCTTGATCGCAGCCGCTACGGCTTCAATGTCGTTGAGATCGACATAGGTGGTGGTCAGCCCCGAGGTGCGGCTGCGGTAATCTTCCATGATGCGGAACGTGCAACCGTACACGcc                                                                                                                >  1:237911/1‑149 (MQ=255)                                         cGGCGACTCGATCCAGATCAGCTGGGTTTCC‑GGCTTGATCGCAGCCGCTACGGCTTCAATGTCGTTGAGATCGACATAGGGGGTGGTCAGCCCCGAGGTGCGGCTGCGGTAATCTTCCATGATGCGGAACGTGCCACCGTACACGCCGt                                                                                                              <  4:298337/149‑1 (MQ=255)                                          ggCGACTCGATCCAGATCAGCTGGGTTTCC‑GGCTTGATCGCAGCCGCTACGGCTTCAATGTCGTTGAGATCGACATAGGTGGTGGTCAGCCCCGAGGTGCGGCTGCGGTAATCTTCCATGATGCGGAACGTGCCACCGTACACGCCGt                                                                                                              <  1:261964/148‑1 (MQ=255)                                           gCGACTCGATCCAGATCAGCTGGGTTTCC‑GGCTTGATCGCAGCCGCTACGGCTTCAATGTCGTTGAGATCGACATAGGTGGTGGTCAGCCCCGAGGTGCGGCTGCGGTAATCTTCCATGATGCGGAACGTGCCACCGTACACGCCGTtc                                                                                                            <  8:5260/149‑1 (MQ=255)                                             gACTCGCTCCAGATCAGCTGGGTTTCC‑GGCTTGATCGCAGCCGCTACGGCTTCAATGTCGTTGAGATCGACATAGGTGGTGGTCAGCCCCGAGGTGCGGCTGCGGTAATCTTCCATGATGCGGAACGTGCCACCGTACACGCCGTtca                                                                                                           >  6:34462/1‑148 (MQ=255)                                             gACTCGATCCAGATCAGCTGGGTTTCC‑GGCTTGATCGCAGCCGCTACGGCTTCAATGTCGTTGAGATCGACATAGGTGGTGGTCAGCCCCGAGGTGCGGCTGCGGTAATCTTCCATGATGCGGAACGTGCCACCGTACACGCCGTtcat                                                                                                          <  5:146487/149‑1 (MQ=255)                                              aCTCGATCCAGATCAGCTGGGTTTCC‑GGCTTGATCGCAGCCGCTACGGCTTCAATGTCGTTGAGATCGACATAGGTGGTGGTCAGCCCCGAGGTGCGGCTGCGGTAATCTTCCATGATGCGGAACGTGCCACCGTACACGCCGTtcatc                                                                                                         <  5:34462/149‑1 (MQ=255)                                              aCTCGATCCAGATCAGCTGGGTTTCC‑GGCTTGATCGCAGCCGCTACGGCTTCAATGTCGTTGAGATCGACATAGGTGGTGGTCAGCCCCGAGGTGCGGCTGCGGTAATCTTCCATGATGCGGAACGTGCCACCGTACACGCCGTtcatc                                                                                                         <  3:356647/149‑1 (MQ=255)                                                tCGATCCAGATCAGCTGGGTTTCC‑GGCTTGATCGCAGCCGCTACGGCTTCAATGTCGTTGAGATCGACATAGGTGGTGGTCAGCCCCGAGGTGCGGCTGCGGTAATCTTCCATGATGCGGAACGTGCCACCGTACACGCCGTTCATcac                                                                                                       >  3:46575/1‑149 (MQ=255)                                                   aTCCAGATCAGCTGGGTTTCC‑GGCTTGATCGCAGCCGCTACGGCTTCAATGTCGTTGAGATCGACATAGGTGGTGGTCAGCCCCGAGGTGCGGCTGCGGTAATCTTCCATGATGCGGAACGTGCCACCGTACACGCCGTTCATcaccac                                                                                                    >  7:362039/1‑149 (MQ=255)                                                      cAGATCAGCTGGGTTTCC‑GGCTTGATCGCAGCCGCTACGGCTTCAATGTCGTTGAGATCGACATAGGTGGTGGTCAGCCCCGAGGTGCGGCTGCGGTAATCTTCCATGATGCGGAACGTGCCACCGTACACGCCGTTCATcac                                                                                                       >  7:362988/1‑143 (MQ=255)                                                      cAGATCAGCTGGGTTTCC‑GGCTTGATCGCAGCCGATACGGCTTCAATGTCGTTGAGATCGACATAGGTGGTGGTCAGCCCCGAGGTGCGGCTGCGGTAATCTTCCATGATGCGGAACGTGCCACCGTACACGCCGTTCATcac                                                                                                       <  8:362988/143‑1 (MQ=255)                                                          tCAGCTGGGTTTCC‑GGCTTGATCGCAGCCGCTACGGCTTCAATGTCGTTGAGATCGACATAGGTGGTGGTCAGCCCCGAGGTGCGGCTGCGGTAATCTTCCATGATGCGGAACGTGCCACCGTACACGCCGTTCATCACCACCACGTgg                                                                                             <  7:141688/149‑1 (MQ=255)                                                                 ggTTTCC‑GGCTTGATCGCAGCCGCTACGGCTTCAATGTCGTTGAGATCGACATAGGTGGTGGTCAGCCCCGAGGTGCGGCTGCGGTAATCTTCCATGATGCGGAACGTGCCACCGTACACGCCGTTCATCACCACCACGTGGGCATc                                                                                        <  4:46575/147‑1 (MQ=255)                                                                   tttCC‑GGCTTGATCGCAGCCGCTACGGCTTCAATGTCGTTGAGATCGACATAGGTGGTGGTCAGCCCCGAGGTGCGGCTGCGGTAATCTTCCATGATGCGGAACGTGCCACCGTACACGCCGTtc                                                                                                            <  2:236071/125‑1 (MQ=255)                                                                   tttCC‑GGCTTGATCGCAGCCGCTACGGCTTCAATGTCGTTGAGATCGACATAGGTGGTGGTCAGCCCCGAGGTGCGGCTGCGGTAATCTTCCATGATGCGGAACGTGCCACCGTACACGCCGTtc                                                                                                            >  1:236071/1‑125 (MQ=255)                                                                         ggCTTGATCGCAGCCGCTACGGCTTCAATGTCGTTGAGATCGACATAGGTGGTGGTCAGCCCCGAGGTGCGGCTGCGGTAATCTTCCATGATGCGGAACGTGCCACCGTACACGCCGTTCATCACCACCACGTGGGCATCCTTGGgcag                                                                               >  3:204791/1‑149 (MQ=255)                                                                                tCGCAGCCGCTACGGCTTCAATGTCGTTGAGATCGACATAGGTGGTGGTCAGCCCCGAGGTGCGGCTGCGGTAATCTTCCATGATGCGGAACGTGCCa                                                                                                                           >  2:161242/1‑98 (MQ=255)                                                                                ccgcagccgcacctcggggctgaccaccacctatgtcgatctcaacgacattgaagccgtagcggctgcgGTAATCTTCCATGATGCGGAACGTGCCACCGTACACGCCGTTCATCACCACCACGTGGGCATCCTTGGGCAGCAGCTc                                                                         >  4:235977‑M2/71‑148 (MQ=255)                                                                                 cgcagccgcacctcggggctgaccaccacctatgtcgatctcaacgacattgaagccgtagcggctgcgGTAATCTTCCATGATGCGGAACGTGCCACCGTACACGCCGTTCATCACCACCACGTGGGCATCCTTGGGCAGCAGCTCaa                                                                       >  7:18506‑M2/70‑149 (MQ=255)                                                                                  GCAGCCGCTACGGCTTCAATGTCGTTGAGATCGACATAGGTGGTGGTCAGCCCCGAGGTGCGGCTGCGGTAATCTTGCATGATGCGGAACGTGCCACCGTACACGCCGTTCATCACCACGACGTGGGCATCCTTGGGCAGCAGCTCaa                                                                       >  2:237911‑M1/1‑148 (MQ=255)                                                                                              gCTTCAATGTCGTTGAGATCGACATAGGTGGTGGTCAGCCCCGAGGTGCGGCTGCGGTAATCTTCCATGATGCGGAACGTGCCACCGTACACGCCGTtc                                                                                                            >  6:186685/1‑99 (MQ=255)                                                                                                  cAATGTCGTTGAGATCGACATAGGTGGTGGTCAGCCCCGAGGTGCGGCTGCGGTAATCTTCCATGATGCGGAACGTGCCACCGTACACGCCGTTCATCACCACCACGTGGGCATCCTTGGGCAGCAGCTCAAGCACGGTTGCGGTGGCa                                                      >  5:21787/1‑149 (MQ=255)                                                                                                  cAATGTCGTTGAGATCGACATAGGTCGTGGACAGCCACGAGGTGCGGCTGCGGTAATCTTCCATGATGAGGAACGTGCCACAGAACACGACGATCATCACCACCACGTGCGCATCCTTCGGCAGCAGCTCAAGCACGGTTGCGGTGGCa                                                      >  8:96639/1‑149 (MQ=255)                                                                                                    aTGTCGTTGAGATCGACATAGGTGGTGGTCAGCCCCGAGGTGCGGCTGCGGTAATCTTCCATGATGCGGAACGTGCCACCGTACACGCCGTTCATCACCACCACGTGGGCATCCTTGGGCAGCAGCTCAAGCACGGTTGCGGTGGCAtt                                                    >  7:289220/1‑149 (MQ=255)                                                                                                      accaccacctatgtcgatctcaacgacattgaagccgtagcggctgcgGTAATCTTCCATGATGCGGAACGTGCCACCGTACACGCCGTTCATCACCACCACGTGGGCATCCTTGGGCAGCAGCTCAAGCACGGTTGCGGTGGCATTGa                                                  <  8:18506‑M2/101‑1 (MQ=255)                                                                                                       tCGTTGAGATCGACATAGGTGGTGGTCAGCCCCGAGGTGCGGCTGCGGTAATCTTCCATGATGCGGAACGTGCCACCGTACACGCCGTTCATCACCACCACGTGGGCATCCTTGGGCAGCAGCTCAAGCACGGTTGCGGTGGCATTGa                                                  <  8:289220/148‑1 (MQ=255)                                                                                                        cGTTGAGATCGACATAGGTGGTGGTCAGCCCCGAGGTGCGGCTGCGGTAATCTTCCATGATGCGGAACGTGCCACCGTACACGCCGTTCATCACCACCACGTGGGCATCCTTGGGCAGCAGCTCAAGCACGGTTGCGGTGGCATTGAc                                                 <  8:362039/148‑1 (MQ=255)                                                                                                         gTTGAGATCGACATAGGTGGTGGTCAGCCCCGAGGTGCGGCTGCGGTAATCTTCCATGATGCGGAACGTGCCACCGTACACGCCGTTCATCACCACCACGTGGGCATCCTTGGGCAGCAGCTCAAGc                                                                     >  2:231314/1‑127 (MQ=255)                                                                                                         gTTGAGATCGACATAGGTGGTGGTCAGCCCCGAGGTGCGGCTGCGGTAATCTTCCATGATGCGGAACGTGCCACCGTACACGCCGTTCATCACCACCACGTGGGCATCCTTGGGCAGCAGCTCAAGc                                                                     <  1:231314/127‑1 (MQ=255)                                                                                                          ttGAGATCGACATAGGTGGTGGTCAGCCCCGAGGTGCGGCTGCGGTAATCTTCCATGATGCGGAACGTGCCACCGTACACGCCGTTCATCACCACCACGTGGGCATCCTTGGGCAGCAGCTCAAGCACGGTTGCGGTGGCATTGACTcc                                              >  8:21477/1‑149 (MQ=255)                                                                                                             agaTCGACATAGGTGGTGGTCAGCCCCGAGGTGCGGCTGCGGTAATCTTCCATGATGCGGAACGTGCCACCGTACACGCCGTTCATCACCACCACGTGGGCATCCTTGGGCAGCAACTCAAGCACGGTTGCGGTGGCATTGACTCCCGa                                           <  3:43203/149‑1 (MQ=255)                                                                                                                tCGACATAGGTGGTGGGCAGCCCCGAGGTGCGGCTGCGGTAATCTTCCATGATGCGGAACGTGCCACCGTACACGCCGTTCATCACCACCACGTGGGCATGCTTGGGCAGCAGCTCAAGCACGGTTGCGGTGGCATTGACTCCCGAGgc                                        <  1:280933/149‑1 (MQ=255)                                                                                                                    cATAGGTGGTGGTCAGCCCCGAGGTGCGGCTGCGGTAATCTTCCATGATGCGGAACGTGCCACCGTACACGCCGTTCATCACCACCACGTGGGCATCCTTGGGCAGCAGCTCAAGCACGGTTGCGGTGGCATTGACTCCCGAGGCGCAg                                    <  1:181370/149‑1 (MQ=255)                                                                                                                      tAGGTGGTGGTCAGCCCCGAGGTGCGGCTGCGGTAATCTTCCATGATTCGGAACGGGCCACCGTACACGCCGTTCATCACCACCACGTGGGCATCCTTGGGCAGCAGCTCAAGCACGGTTGCGGTGGCATTGACTCCCGAGGCGCAGGc                                  <  7:158549/149‑1 (MQ=255)                                                                                                                      tAGGTGGTGGTCAGCCCCGAGGTGCGGCTGCGGTAATCTTCCATGATGCGGAACGTGCCACCGTACACGCCGTTCATCACCACCACGTGGGCATCCTTGGGCAGCAGCTCCAGCACGGTTGCGGTGGCATTGCCTCCCGAGCCGCCGGc                                  >  4:232588/1‑149 (MQ=255)                                                                                                                              ggtCAGCCCCGAGGTGCGGCTGCGGTAATCTTCCATGATGCGGAACGTGCCACCGTACACGCCGTTCATCACCACCACGTGGGCATCCTTGGGCAGCAGCTCAAGCACGGTTGCGGTGGCATTGACTCCCGAGGCGCAGGCCACAgcgc                          >  7:100947/1‑149 (MQ=255)                                                                                                                               gtCAGCCCCGAGGTGCGGCTGCGGTAATCTTCCATGATGCGGAACGTGCCACCGTACACGCCGTTCATCACCACCACGTGGGCATCCTTGGGCAGCAGCTCAAGCACGGTTGCGGTGGCATTGACTCCCGAGGCGCAg                                    <  8:99805/138‑1 (MQ=255)                                                                                                                               gtCAGCCCCGAGGTGCGGCTGCGGTAATCTTCCATGATGCGGAACGTGCCACCGTACACGCCGTTCATCACCACCACGTGGGCATCCTTGGGCAGCAGCTCAAGCACGGTTGCGGTGGCATTGACTCCCGAGGCGCAg                                    >  7:99805/1‑138 (MQ=255)                                                                                                                                           gTGCGGCTGCGGTAATCTTCCATGATGCGGAACGTGCCACCGTACACGCCGTTCATCACCACCACGTGGGCATCCTTGGGCAGCAGCTCAAGCACGGTTGCGGTGGCATTGACTCCCGAGGCGCAGGCCACAGCGCCCACGCCttcttc             >  1:133297/1‑149 (MQ=255)                                                                                                                                              cggctgcgGTAATCTTCCATGATGCGGAACGTGCCACCGTACACGCCGTTCATCACCACCACGTGGGCATCCTTGGGCAGCAGCTCAAGCACGGTTGCGGTGGCATTGACTCCCGAGGCGCAGGCCACAGCGCCCACGCCTTCTTcc            >  3:178340‑M2/9‑147 (MQ=255)                                                                                                                                                gctgcgGTAATCTTCCATGATGCGGAACGTGCCACCGTACACGCCGTTCATCACCACCACGTGGGCATCCTTGGGCAGCAGCTCAAGCACGGTTGCGGTGGCATTGACTCCCGAGGCGCAGGCCACAGCGCCCACGCCTTCTTCCAAgg        >  5:150482‑M2/7‑149 (MQ=255)                                                                                                                                                     gGTAATCTTCCATGATGCGGAACGTGCCACCGTACACGCCGTTCATCACCACCACGTGGGCATCCTTGGGCAGCAGCTCAAGCACGGTTGCGGTGGCATTGACTCCCGAGGCGCAGGCCACAGCGCCCACGCCTTCTTCCAAGGCAGc    >  1:163436‑M2/2‑148 (MQ=255)                                                                                                                                                      gTAATCTTCCATGATGCGGAACGTGCCACCGTACACGCCGTTCATCACCACCACGTGGGCATCCTTGGGCAGCAGCTCAAGCACGGTTGCGGTGGCATTGACTCCCGAGGCGCAGGCCACAGCGCCCACGCCTTCTTCCAAGGCAGCGa  <  4:178340/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                     ACACCGCCTTGATGTCGACCAGGTGCAGCAACGGGTTGGTCGGCGACTCGATCCAGATCAGCTGGGTTTCC‑GGCTTGATCGCAGCCGCTACGGCTTCAATGTCGTTGAGATCGACATAGGTGGTGGTCAGCCCCGAGGTGCGGCTGCGGTAATCTTCCATGATGCGGAACGTGCCACCGTACACGCCGTTCATCACCACCACGTGGGCATCCTTGGGCAGCAGCTCAAGCACGGTTGCGGTGGCATTGACTCCCGAGGCGCAGGCCACAGCGCCCACGCCTTCTTCCAAGGCAGCGA  >  NC_002947/5213213‑5213509

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.