New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 5196440	44 (1.390)	8 (0.270)	7/246	2.5	15.7%	coding (123/690 nt)	ybgJ	allophanate hydrolase subunit 1
	?	NC_002947	= 5196465	45 (1.540)					coding (98/690 nt)	ybgJ	allophanate hydrolase subunit 1
Rejected: Frequency below/above cutoff threshold.

GTGTCGGCGGCCGCCACTGCCTGTGTCAGGCTGCAAGCCGTCCAGCGCCTGAGCGATTCGCGCCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCACTTCGATTCGCGGCTTCATCACACCTCCAGCCCATCGAGGGCCTGGCGA  >  NC_002947/5196299‑5196589                                                                                                                                              |                                                                                                                                                      gtgtCGGCGGCCGCCACTGCCTGTGTCAGGCTGCAAGCCGTCCAGCGCCTGAGCGATTCGCGCCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGccagccgg                                                                                                                                               >  1:602511‑M1/1‑142 (MQ=255)       gcggcCGCCACTGCCTGTGTCAGGCTGCAAGCCGTCCAGCGCCTGAGCGATTCGCGCCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCaa                                                                                                                                         >  1:448123/1‑150 (MQ=255)       gcggcCGCCACTGCCTGTGTCAGGCTGCAAGCCGTCCAGCGCCTGAGCGATTCGCGCCCTGGCCTCGCCTGCCGCCACGACAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCaa                                                                                                                                         <  2:135373/150‑1 (MQ=255)             gccACTGCCTGTGTCAGGCTGCAAGCCGTCCAGCGCCTGAGCGATTCGCGCCCTGGCCTCGCCTGGCGCCAGGTCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGccagccggcgcttgagcgcg                                                                                                                                   >  1:698149‑M1/1‑130 (MQ=255)               cACTGCCTGTGTCAGGCTGCAAGCCGTCCAGCGCCTGAGCGATTCGCGCCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGc                                                                                                                                 >  2:693110/1‑150 (MQ=255)                    ccTGTGTCAGGCTGCAAGCCGTCCAGCGCCTGAGCGATTCGCGCCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGaacgccggcgcttgagcgcggcgttcg                                                                                                                            <  2:141482‑M1/150‑28 (MQ=255)                      tgtgTCAGGCTGCAAGCCGTCCAGCGCCTGAGCGATTCGCGCCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGccagccggcgcttgagcgcggcgttcggc                                                                                                                          <  2:602511‑M1/150‑30 (MQ=255)                       gtgtCAGGCTGCAAGCCGTCCAGCGCCTGAGCGATTCGCGCCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGccagccggcgcttgagcgcggcgttcggcg                                                                                                                         <  2:698149‑M1/150‑31 (MQ=255)                                 gCAAGCCGTCCAGCGCCTGAGCTATTCGCGCCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATg                                                                                                               <  2:448123/150‑1 (MQ=255)                                    aGCCGTCCAGCGCCTGAGCGATTCGCGCCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCa                                                                                                            <  2:12302/150‑1 (MQ=255)                                         tCCAGCGCCTGAGCGATTCGCGCCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAGATGCTCGCCGccagccggcgcttgagcgcggcgttcggcgagcatttgctggacctgg                                                                                                       <  1:634085‑M1/150‑49 (MQ=255)                                           catcgcCTGTGCGATTCCCGCCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCGCGCCGAACGCCTCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGAc                                                                                                     <  1:264913/147‑1 (MQ=255)                                           cAGCGCCTGAGCGATTCGCGCCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGAc                                                                                                     <  2:143520/150‑1 (MQ=255)                                           cAGCGCCTGAGCGATTCGCGCCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGAc                                                                                                     <  1:80391/150‑1 (MQ=255)                                                ccTGAGCGATTCGCGCCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCAt                                                                                                <  1:62928/150‑1 (MQ=255)                                                    aGCGATTCGCGCCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCATGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGAt                                                                                            >  1:656871/1‑150 (MQ=255)                                                      cGATTCGCGCCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATAc                                                                                          >  1:377759/1‑150 (MQ=255)                                                         ttagcgcCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCAGCGCTCAAGCGCCGACCGGCGGCGAGCATCCATGGCAGGTTGACCTCATCGATACggg                                                                                       >  2:175145/4‑149 (MQ=255)                                                          tcgaggccTGGCCCCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATCAGGCGGAAGCAGGTCCAGCAAATGCTCGCTGAACGCCGCGCTCCAGCGCCGGCTGGACGCGAGCATCCTTGGCATGTTGACCTCATCGATACGGTc                                                                                      <  2:335982/144‑1 (MQ=255)                                                          tCGCGCCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTc                                                                                      <  2:76915/150‑1 (MQ=255)                                                            gcgcCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGAc                                                                                                     >  1:571318/1‑133 (MQ=255)                                                            gcgcCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGAc                                                                                                     <  2:571318/133‑1 (MQ=255)                                                                ccTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAg                                                                                >  1:13385/1‑150 (MQ=255)                                                                ccTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAg                                                                                >  1:245698/1‑150 (MQ=255)                                                                 cTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAgg                                                                               <  2:656871/150‑1 (MQ=255)                                                                    gCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCg                                                                             <  1:741893/149‑1 (MQ=255)                                                                    gCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCg                                                                             >  2:741893/1‑149 (MQ=255)                                                                     ccTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTa                                                                           <  1:693110/150‑1 (MQ=255)                                                                          ccTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATc                                                                      <  2:494988/150‑1 (MQ=255)                                                                           cTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGccagccggcgcttgagcgcggcgttcggcgagcatttgctggacctgg                                                                                                       >  2:613533‑M1/1‑68 (MQ=255)                                                                           cTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGccagccggcgcttgagcgcggcgttcggcgagcatttgctggacctgg                                                                                                       <  1:613533‑M1/116‑49 (MQ=255)                                                                                  aaGGGCCGGCTGCACCCTCAGGGTAGTAGAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTc                                                              <  2:251001/149‑1 (MQ=255)                                                                                     ggCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGAt                                                           >  1:682943/1‑150 (MQ=255)                                                                                     ggCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGAt                                                           >  1:622111/1‑150 (MQ=255)                                                                                                        gtagtaTAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCCCTTCGATTcgc                                        <  2:796669/150‑1 (MQ=255)                                                                                                           gtaTAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCACTTCGATTcgc                                        <  2:273431/147‑1 (MQ=255)                                                                                                           gtaTAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCACTTCGATTcgc                                        >  1:273431/1‑147 (MQ=255)                                                                                                              taAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCACTTCGATTCGCGGCTtc                                  >  1:631131/1‑150 (MQ=255)                                                                                                                  aCGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCACTTCGATTCGCGGCTtcatca                              >  1:753942/1‑150 (MQ=255)                                                                                                                        ccAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGcggcgagcatttgctggacctggtgccgtcttatactacgctgatggtgcagtctgccctggcgccag                                                          <  1:613533‑M2/116‑69 (MQ=255)                                                                                                                        ccAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGcggcgagcatttgctggacctggtgccgtcttatactacgctgatggtgcagtctgccctggcgccag                                                          >  2:613533‑M2/1‑48 (MQ=255)                                                                                                                        ccAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGcggcgagcatctgctggacctggtgccgtcttatactacgctgatggtgcagtctgccctggcgccaggcgaggccagggcgcgaatcgctcaggcgctgga                        <  1:634085‑M2/150‑103 (MQ=255)                                                                                                                          aGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCACTTCGATTCGCGGCTTCATCACACCTCCa                      <  2:245698/150‑1 (MQ=255)                                                                                                                                aGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCACTTCGATTCGCGGCTTCATCACACCTCCAGCCCAt                <  2:377759/150‑1 (MQ=255)                                                                                                                                         tCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAg                                                                                <  2:99697/77‑1 (MQ=255)                                                                                                                                         tCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAg                                                                                >  1:99697/1‑77 (MQ=255)                                                                                                                                         tCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGACTGTCGATGGCGACCACTTCGATTCGCGGCTTCATCACACCTCCAGCCCGTCGGGGGCCt       >  2:3690/1‑150 (MQ=255)                                                                                                                                          cGCCGAACGCCGCGCTCAAGCGCCGGCTGGcggcgagcatttgctggacctggtgccgtcttatactacgctgatggtgcagtctgccctggcgccaggcgaggccagggcgcgaatcgctcaggcgctggacggcttgcagcctgacac      <  2:698149‑M2/150‑121 (MQ=255)                                                                                                                                           gCCGAACGCCGCGCTCAAGCGCCGGCTGGcggcgagcatttgctggacctggtgccgtcttatactacgctgatggtgcagtctgccctggcgccaggcgaggccagggcgcgaatcgctcaggcgctggacggcttgcagcctgacaca     <  2:602511‑M2/150‑122 (MQ=255)                                                                                                                                           gCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATcc                                                                                                                  >  2:214240/1‑41 (MQ=255)                                                                                                                                           gCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATcc                                                                                                                  <  1:214240/41‑1 (MQ=255)                                                                                                                                           gCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGc                                                                              <  2:523310/77‑1 (MQ=255)                                                                                                                                           gCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGc                                                                              >  1:523310/1‑77 (MQ=255)                                                                                                                                           gCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTc                                                              >  2:264568/1‑93 (MQ=255)                                                                                                                                           gCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTc                                                              <  1:264568/93‑1 (MQ=255)                                                                                                                                            ccGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCACTTCGATTCGCGGATTCATCACACCTCCAGCCCATCGAGGGCCTGGc    <  2:13385/150‑1 (MQ=255)                                                                                                                                             cGAACGCCGCGCTCAAGCGCCGGCGTTcggcgagcatttgctggacctggtgccgtcttatactacgctgatggtgcagtctgccctggcgccaggcgaggccagggcgcgaatcgctcaggcgctggacggcttgcagcctgacacagg   <  2:141482‑M2/150‑124 (MQ=255)                                                                                                                                              gAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCACTTCGATTCGCGGCTTCATCACACCTCCAGCCCATCGAGGGCCTGGCGa  >  2:609497/1‑150 (MQ=255)                                                                                                                                              |                                                                                                                                                      GTGTCGGCGGCCGCCACTGCCTGTGTCAGGCTGCAAGCCGTCCAGCGCCTGAGCGATTCGCGCCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCACTTCGATTCGCGGCTTCATCACACCTCCAGCCCATCGAGGGCCTGGCGA  >  NC_002947/5196299‑5196589

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