New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	5731921 =	49 (1.130)	6 (0.150)	5/240	NT	11.7%	coding (691/1206 nt)	hutI	imidazolonepropionase
	?	NC_002947	5731965 =	44 (1.080)					coding (647/1206 nt)	hutI	imidazolonepropionase

AGCACGGCGACCGTGCCGGCGCTTGCCATGGCCACGGCGTCTTCCTCTGTCATGAACTCCAGATGGTCAGCCGACAGTGCCTGGTAACGCGCGGCCAGGCTTGAGCCGTGCAACGATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTC  >  NC_002947/5731773‑5732065                                                                                                                                                     |                                                                                                                                                 agcaCGGCGACCGTGCCGGCGCTTGCCATGGCCACGGCGTCTTCCTCTGTCATGAACTCCAGATGGTCAGCCGACAGTGCCTGGTAACGCGCGGCCAGGCTTGAGCCGTGCAACGATGACAGTTGTTCGGCATGCAGCTTGACCGGCAgg                                                                                                                                                 <  2:498620/150‑1 (MQ=255)           ccGTGCCGGCGCTTGCCATGGCCACGGCGTCTTCCTCTGTCATGAACTCCAGCTGGTCAGCCGACAGTGCCTGGTAACGCGCGGCCAGGCTTGAGCCGTGCAACGATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTcgc                                                                                                                                       <  2:623241/150‑1 (MQ=255)                 cGGCGCTTGCCATGGCCACGGCGTCTTCCTCTGTCATGAACTCCAGATGGTCAGCCGACAGTGCCTGGTAACGCGCGGCCAGGCTTGAGCCGTGCAACGATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCt                                                                                                                                 >  1:977788/1‑150 (MQ=255)                    cgcTTGCCATGGCCACGGCGTCTTCCTCTGTCATGAACTCCAGATGGTCAGCCGACAGTGCCTGGTAACGCGCGGCCAGGCTTGAGCCGTGCAACGATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTtga                                                                                                                              <  2:880242/150‑1 (MQ=255)                     gcTTGCCATGGCCACGGCGTCTTCCTCTGTCATGAACTCCAGATGGTCAGCCGACAGTGCCTGGTAACGCGCGGCCAGGCTTGAGCCGTGCAACGATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTtgat                                                                                                                             <  2:775293/150‑1 (MQ=255)                     gcTTGCCATGGCCACGGCGTCTTCCTCTGTCATGAACTCCAGATGGTCAGCCGACAGTGCCTGGTAACGCGCGGCCAGGCTTGAGCCGTGCAACGATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTtgat                                                                                                                             <  1:340633/150‑1 (MQ=255)                           cATGGCCACGGCGTCTTCCTCTGTCATGAACTCCAGATGGTCAGCCGACAGTGCCTGGTAACGCGCGGCCAGGCTTGAGCCGTGCAACGATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAg                                                                                                                       >  1:585516/1‑150 (MQ=255)                           cATGGCCACGGCGTCTTCCTCTGTCATGAACTCCAGATGGTCAGCCGACAGTGCCTGGTAACGCGCGGCCAGGCTTGAGCCGTGCAACGATGACAGTTGTTCGGCATGCAGCTTGACCCGCAGGCCCAGTTCGCGCGCCTTGATGAACAg                                                                                                                       >  2:463040/1‑150 (MQ=255)                                    ggtcgggccgacgactacatcgcgcacatctgtgacgagatgttgccggccttggccggcgaaggcctggtggatgcggtggatgccttctgcgaacacctcgcgttctcgccGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGAc                                                                                                              >  1:963065‑M1/114‑150 (MQ=255)                                      tcgggccgacgactacatcgcgcacatctgtgacgagatgttgccggccttggccggcgaaggcctggtggatgcggtggatgccttctgcgaacacctcgcgttctcgccGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCt                                                                                                            >  1:874347‑M1/112‑150 (MQ=255)                                      tcgggccgacgactacatcgcgcacatctgtgacgagatgttgccggccttggccggcgaaggcctggtggatgcggtggatgccttctgcgaacacctcgcgttctcgccGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCt                                                                                                            >  1:874101‑M1/112‑150 (MQ=255)                                      ggTCTTCCTCTGTCATGAACTCCAGATGGTCAGCCGACAGTGCCTGGTAACGCGCGGCCAGGCTTGAGCCGTGCAACGATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCt                                                                                                            <  2:346837/149‑1 (MQ=255)                                       gTCTTCCTCTGTCATGAACTCCAGATGGTCAGCCGACAGTGCCTGGTAACGCGCGGCCAGGCTTGAGCCGTGCAACGATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTg                                                                                                           >  2:924234/1‑150 (MQ=255)                                       cgggccgacgactacatcgcgcacatctgtgacgagatgttgccggccttggccggcgaaggcctggtggatgcggtggatgccttctgcgaacacctcgcgttctcgccGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTg                                                                                                           >  1:251892‑M1/111‑150 (MQ=255)                                        tCTTCCTCTGTCATGAACTCCAGATGGTCAGCCGACAGTGCCTGGTAACGCGCGGCCAGGCTTGAGCCGTGCAACGATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTgg                                                                                                          <  1:302574/150‑1 (MQ=255)                                        tCTTCCTCTGTCATGAACTCCAGATGGTCAGCCGACAGTGCCTGGTAACGCGCGGCCAGGCTTGAGCCGTGCAACGATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTgg                                                                                                          <  2:387858/150‑1 (MQ=255)                                          ttCCTCTGTCATGAACTCCAGATGGTCAGCCGACAGTGCCTGGTAACGCGCGGCCAGGCTTGAGCCGTGCAACGATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGc                                                                                                        >  2:147081/1‑150 (MQ=255)                                             ctctGTCATGAACTCCAGATGGTCAGCCGACAGTGCCTGGTAACGCGCGGCCAGGCTTGAGCCGTGCAACGATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCgg                                                                                                     <  2:907393/150‑1 (MQ=255)                                             ctctGTCATGAACTCCAGATGGTCAGCCGACAGTGCCTGGTAACGCGCGGCCAGGCTTGAGCCGTGCAACGATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCgg                                                                                                     >  1:907393/1‑150 (MQ=255)                                              tctGTCATGAACTCCAGATGGTCAGCCGACAGTGCCTGGTAACGCGCGGCCAGGCTTGAGCCGTGCAACGATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGc                                                                                                    >  1:114178/1‑150 (MQ=255)                                              tctGTCATGAACTCCAGATGGTCAGCCGACAGTGCCTGGTAACGCGCGGCCAGGCTTGAGCCGTGCAACGATGACAGTTGTTCGGCATGCAGCCTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGc                                                                                                    <  2:585516/150‑1 (MQ=255)                                                tGTCATGAACTCCAGATGGTCAGCCGACAGTGCCTGGTAACGCGCGGCCAGGCTTGAGCCGTGCAACGATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGACAGGCAGCTCGACCTGGGCCGGCga                                                                                                  >  1:335199/1‑150 (MQ=255)                                                tGTCATGAACTCCAGATGGTCAGCCGACAGTGCCTGGTAACGCGCGGCCAGGCTTGAGCCGTGCAACGATGACAGTTGTAGATCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGATCTGGGACGGCga                                                                                                  >  1:69899/1‑150 (MQ=255)                                                gactacatcgcgcacatctgtgacgagatgttgccggccttggccggcgaaggcctggtggatgcggtggatgccttctgcgaacacctcgcgttctcgccGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCga                                                                                                  >  2:873733‑M1/102‑150 (MQ=255)                                                     tGAACTCCAGATGGTCAGCCGACAGTGCCTGGTCTGACACGGCCAGGCTTGAGCCGTGCAACGATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAAcg                                                                                             >  1:1112665/1‑150 (MQ=255)                                                         gcgcacatctgtgacgagatgttgccggccttggccggcgaaggcctggtggatgcggtggatgccttctgcgaacacctcgcgttctcgccGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAg                                                                                         >  1:978233‑M1/93‑150 (MQ=255)                                                            cAGATGGTCAGCCGACAGTGCCTGGTAACGCGCGGCCAGGCTTGAGCCGTGCAACGATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGa                                                                                          >  1:358407/1‑146 (MQ=255)                                                            cAGATGGTCAGCCGACAGTGCCTGGTAACGCGCGGCCAGGCTTGAGCCGTGCAACGATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGa                                                                                          <  2:358407/146‑1 (MQ=255)                                                               gtGGTCAGCCGACAGTGCGTGGTGACGCGCGGCCAGGCTTGAGCCGTGCAACGATGCCAGTTGTTCGGCATGCAGCTTGACCGGCCGCTCCAGTTCGCGGGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTc                                                                                   <  2:375387/149‑1 (MQ=255)                                                                     aGCCGACAGTGCCTGGTAACGCGCGGCCAGGATTGAGCAGTTCAACGATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGaa                                                                             <  2:1112665/150‑1 (MQ=255)                                                                      gCCGACAGTGCCTGGTAACGCGCGGCCAGGCTTGAGCCGTGCAACGATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTAGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAg                                                                            >  2:400900/1‑150 (MQ=255)                                                                       ccGACAGTGCCTGGTAACGCGCGGCCAGGCTTGAGCCGTGCAACGATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAg                                                                            >  1:19541/1‑149 (MQ=255)                                                                       ccGACAGTGCCTGGTAACGCGCGGCCAGGCTTGAGCCGTGCAACGATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAg                                                                            <  2:19541/149‑1 (MQ=255)                                                                                  tGGTAACGCGCGGCCAGGCTTGAGCCGTGCAACGATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGgcatccaccgca                                                                <  2:977788/150‑1 (MQ=255)                                                                                    gTAACGCGCGGCCAGGCTTGAGCCGTGCAACGATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGgcatccaccgcatc                                                              >  2:516771/1‑150 (MQ=255)                                                                                     tAACGCGCGGCCAGGCTTGAGCCGTGCAACGATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGa                                                                                          <  2:1040122/121‑1 (MQ=255)                                                                                     tAACGCGCGGCCAGGCTTGAGCCGTGCAACGATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCATTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGa                                                                                          >  1:1040122/1‑121 (MQ=255)                                                                                          cgcgGCCAGGCTTGAGCCGTGCAACGATGACAGTTGTTCGGCATCAAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAg                                                        <  1:674687/150‑1 (MQ=255)                                                                                               ccAGGCTTGAGCCGTGCAACGATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCGCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCtc                                                   <  2:187307/150‑2 (MQ=255)                                                                                                 aGGCTTGAGCCGTGCAACGATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCg                                                 >  2:134997/1‑150 (MQ=255)                                                                                                  ggCTTGAGCCGTGCAACGATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGc                                                <  2:311760/150‑1 (MQ=255)                                                                                                     ttGAGCCGTGCAACGATGACATTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAAGCCTTCGCCgg                                             >  2:100220/1‑150 (MQ=255)                                                                                                              gCAACGATGACATGTGTTCGGCATGCAGCGTGACCGGCAGGCCCAGTTCGCTCGCCTTCATGAACCGGCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGTCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCg                                    <  1:463040/150‑1 (MQ=255)                                                                                                               cAACGATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTTGCCGGCCTTGGCCgg                                   >  2:294933/1‑150 (MQ=255)                                                                                                                    aTGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTTGCCGGCCTTGGCCGGcgaag                              <  1:294933/150‑5 (MQ=255)                                                                                                                       aCAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACAtct                           <  1:126683/150‑1 (MQ=255)                                                                                                                         aGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCaa                                         >  1:612291/1‑134 (MQ=255)                                                                                                                         aGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCaa                                         <  2:612291/134‑1 (MQ=255)                                                                                                                              ttCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTcaca                    >  2:546349/1‑150 (MQ=255)                                                                                                                                 ggCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCAGGCGCGCACATTATGTTTTCGCATAATTCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGAt                 <  1:134997/150‑1 (MQ=255)                                                                                                                                  gCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGaa                                                                             <  2:483362/89‑1 (MQ=255)                                                                                                                                  gCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGaa                                                                             >  1:483362/1‑89 (MQ=255)                                                                                                                                      gCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTgcg            >  2:999346/1‑150 (MQ=255)                                                                                                                                         gCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGa         <  1:1007210/150‑1 (MQ=255)                                                                                                                                              aCCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGAt        >  1:564507/1‑146 (MQ=255)                                                                                                                                              aCCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGAt        <  2:564507/146‑1 (MQ=255)                                                                                                                                                cGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAgtc  <  1:147081/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                 AGCACGGCGACCGTGCCGGCGCTTGCCATGGCCACGGCGTCTTCCTCTGTCATGAACTCCAGATGGTCAGCCGACAGTGCCTGGTAACGCGCGGCCAGGCTTGAGCCGTGCAACGATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTC  >  NC_002947/5731773‑5732065

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.