New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	5739225 =	53 (0.890)	4 (0.070)	4/234	NT	7.4%	coding (99/1365 nt)	hutF	formiminoglutamate deiminase
	?	NC_002947	5739261 =	49 (0.840)					coding (135/1365 nt)	hutF	formiminoglutamate deiminase

AACAGGATTCCACCGATGTCAGCCTACTTCGCCGAACGCGCCCTGCTACCTACGGGTTGGGCCAGCCATGTGCGTATCGAGGTCGCCAGCGATGGCCGTGTGGCCAGCATCGAGCCTGGCGCTTCGGCCGAAGACGCCGAGCGCCTGGCCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCACTCGCATGCCTTCCAGCGCGCCATGGCAGGCCTGGCGGAAGTCGCCGGCAACCCCAACGACAGTTTCTGGACCTGGCGAGACCTGATGTACCGCCTGGTCGGGCAGA  >  NC_002947/5739112‑5739400                                                                                                                                                      |                                                                                                                                            aaCAGGATTCCACCGATGTCAGCCTACTTCGCCGAACGCGCCCTGCTACCTACGGGTTGGGCCAGCCATGTGCGTATCGAGGTCGCCAGCGATGGCCGTGTGGCCAGCATCGAGCCTGGCGCTTCGGCCGAAGACGCCGAGCGCCTggcc                                                                                                                                             <  1:607898/150‑1 (MQ=255)      gATTCCACCGATGTCAGCCTACTTCGCCGAACGCGCCCTGCTACCTACGGGTTGGGCCAGCCATGTGCGTATCGAGGTCGCCAGCGATGGCCGTGTGGCCAGCATCGAGCCTGGCGCTTCGGCCGAAGACGCCGAGCGCCTGGCCGGccc                                                                                                                                        >  2:1237144/1‑150 (MQ=255)            aCCGATGTCAGCCTACTTCGCCGAACGCGCCCTGCTACCTACGGGTTGGGCCAGCCATGTGCGTATCGAGGTCGCCAGCGATGGCCGTGTGGCCAGCATCGAGCCTGGCGCTTCGGCCGAAGACGCCGAGCGCCTGGCCGGCCCgctgct                                                                                                                                  >  1:690120/1‑150 (MQ=255)             ccGATGTCAGCCTACTTCGCCGAACGCGCCCTGCTACCTACGGGTTGGGCCAGCCATGTGCGTATCGAGGTCGCCAGCGATGGCCGTGTGGCCAGCATCGAGCCTGGCGCTTCGGCCGAAGACGCCGAGCGCCTGGCCGGCCCGctgctg                                                                                                                                 >  1:1306443/1‑150 (MQ=255)                   tgtcgttggggttgccggcgacttccgccaggcctgccatggcgcgctggaaggcatgcgagtgcaggttgggcatgccgggcagcagcgggccgGCCTGGCGCTTCGGCCGAAGACGCCGAGCGCCTggcc                                                                                                                                             >  2:143996‑M1/96‑132 (MQ=255)                   tgtcgttggggttgccggcgacttccgccaggcctgccatggcgcgctggaaggcatgcgagtgcaggttgggcatgccgggcagcagcgggccgGCCTGGCGCTTCGGCCGAAGACGCCGAGCGCCTggcc                                                                                                                                             <  1:143996‑M1/37‑1 (MQ=255)                       ccTACTTCGCCGAACGCGCCCTGCTACCTACGGGTTGGGCCAGCCATGTGCGTATCGAGGTCGCCAGCGATGGCCGTGTGGCCAGCATCGAGCCTGGCGCTTCGGCCGAAGACGCCGAGCGCCTGGCCGGCCCGCTGCTGCCCGGCATGc                                                                                                                       <  2:1436523/150‑1 (MQ=255)                        cTACTTCGCCGAACGCGCCCTGCTACCTACGGGTTGGGCCAGCCATGTGCGTATCGAGGTCGCCAGCGATGGCCGTGTGGCCAGCATCGAGCCTGGCGCTTCGGCCGAAGACGCCGAGCGCCTGGCCGGCCCGCTGCTGCCCGGCATGcc                                                                                                                      <  1:685978/150‑1 (MQ=255)                            ttCGCCGAACGCGCCCTGCTACCTACGGGTTGGGCCAGCCATGTGCGTATCGAGGTCGCCAGCGATGGCCGTGTGGCCAGCATCGAGCCTGGCGCTTCGGCCGAAGACGCCGAGCGCCTGGCCGGCCCGCTGCTGCCCGGCATGCCCAAc                                                                                                                  >  2:1628602/1‑150 (MQ=255)                            ttCGCCGAACGCGCCCTGCTACCTACGGGTTGGGCCAGCCATGTGCGTATCGAGGTCGCCAGCGATGGCCGTGTGGCCAGCATCGAGCCTGGCGCTTCGGCCGAAGACGCCGAGCGCCTGGCCGGCCCGCTGCTGCCCGGAATGCCCAAc                                                                                                                  >  2:1132259/1‑150 (MQ=255)                                  gAACGCGCCCTGCTACCTACGGGTTGGGCCAGCCATGTGCGTATCGAGGTCGCCAGCGATGGCCGTGTGGCCAGCATCGAGCCTGGCGCTTCGGCCGAAGACGCCGAGCGCCTGGCCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCAc                                                                                                            >  2:847238/1‑150 (MQ=255)                                  gAACGCGCCCTGCTACCTACGGGTTGGGCCAGCCATGTGCGTATCGAGGTCGCCAGCGATGGCCGTGTGGCCAGCATCGAGCCTGGCGCTTCGGCCGAAGACGCCGAGCGCCTGGCCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCAc                                                                                                            <  1:573260/150‑1 (MQ=255)                                         cccTGCTACCTACGGGTTGGGCCAGCCATGTGCGTATCGAGGTCGCCAGCGATGGCCGTGTGGCCAGCATCGAGCCTCGCGCTTCGGCCGAAGACGCCGAGCGCCTGGCCGGCCCGCTGCGGCCCGGCATGCCCAACCTGCACTCGCATg                                                                                                     <  1:622278/150‑1 (MQ=255)                                          ccTGCTACCTACGGGTTGGGCCAGCCATGTGCGTATCGAGGTCGCCAGCGATGGCCGTGTGGCCAGCATCGAGCCTGGCGCTTCGGCCGAAGACGCCGAGCGCCTGGCCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCACTCGCATGc                                                                                                    <  1:1237144/150‑1 (MQ=255)                                                acctacGGGTTGGGCCAGCCATGTGCGTATCGAGGTCGCCAGCGATGGCCGTGTGGCCAGCATCGAGCCTGGCGCTTCGGCCGAAGACGCCGAGCGCCTGGCCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCACTCGCATGCCTTCCa                                                                                              <  1:638214/150‑1 (MQ=255)                                                acctacGGGTTGGGCCAGCCATGTGCGTATCGAGGTCGCCAGCGATGGCCGTGTGGCCAGCATCGAGCCTGGCGCTTCGGCCGAAGACGCCGAGCGCCTGGCCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCACTCGCATGCCTTCCa                                                                                              <  2:310631/150‑1 (MQ=255)                                                 cttacGGGTTGGGCCAGCCATGTGCGTATCGAGGTCGCCAGCGATGGCCGTGTGGCCAGCATCGAGCCTGGCGCTTCGGCCGAAGACGCCGAGCGCCTGGCCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCACTCGCATGCCTTCCAg                                                                                             >  2:408607/3‑150 (MQ=255)                                                  ctacGGGTTGGGCCAGCCATGTGCGTATCGAGGTCGCCAGCGATGGCCGTGTGGCCAGCATCGAGCCTGGCGCTTCGGCCGAAGACGCCGAGCGCCTGGCCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCACTCGCATGCCTTCCAgc                                                                                            >  1:1129220/1‑150 (MQ=255)                                                   tacGGGTTGGGCCAGCCATGTGCGTATCGAGGTCGCCAGCGATGGCCGTGTGGCCAGCATCGAGCCTGGCGCTTCGGCCGAAGACGCCGAGCGCCTGGCCGGCCCgc                                                                                                                                      >  2:1207800/1‑107 (MQ=255)                                                   tacGGGTTGGGCCAGCCATGTGCGTATCGAGGTCGCCAGCGATGGCCGTGTGGCCAGCATCGAGCCTGGCGCTTCGGCCGAAGACGCCGAGCGCCTGGCCGGCCCgc                                                                                                                                      <  1:1207800/107‑1 (MQ=255)                                                        gTTGGGCCAGCCATGTGCGTATCGAGGTCGCCAGCGATGGCCGTGTGGCCAGCATCGAGCCTGGCGCTTCGGCCGAAGACGCCGAGCGCCTGGCCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCACTCGCATGCCTTCCAGCGCGCCa                                                                                      >  2:729949/1‑150 (MQ=255)                                                         ttGGGCCAGCCATGTGCGTATCGAGGTCGCCAGCGATGGCCGTGTGGCCAGCATCGAGCCTGGCGCTTCGGCCGAAGACGCCGAGCGCCTGGCCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCACTCGCATGCCTTCCAGCGCGCCAt                                                                                     <  2:1306443/150‑1 (MQ=255)                                                            ggCCAGCCATGTGCGTATCGAGGTCGCCAGCGATGGCCGTGTGGCCAGCATCGAGCCTGGCGCTTCGGCCGAAGACGCCGAGCGCCTGGCCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCACTCGCATGCCTTCCAGCGCGCCCTGGc                                                                                  >  2:328790/1‑150 (MQ=255)                                                             gccagccaTGTGCGTATCGAGGTCGCCAGCGATGGCCGTGTGGCCAGCATCGAGCCTGGCGCTTCGGCCGAAGACGCCGAGCGCCTGGCCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCACTCGCATGCCTTCCAGCGCGCCATGGCa                                                                                 <  1:408607/150‑1 (MQ=255)                                                              ccagccaTGTGCGTATCGAGGTCGCCAGCGATGGCCGTGTGGCCAGCATCGAGCCTGGCGCTTCGGCCGAAGACGCCGAGCGCCTGGCCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCACTCGCATGCCTTCCAGCGCGCCATGGCAg                                                                                <  2:1129220/150‑1 (MQ=255)                                                               cagccaTGTGCGTATCGAGGTCGCCAGCGATGGCCGTGTGGCCAGCATCGAGCCTGGCGCTTCGGCCGAAGACGCCGAGCGCCTGGCCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCACTCGCATGCCTTCCAGCGCGCCATGGCAgg                                                                               <  1:1628602/150‑1 (MQ=255)                                                                  ccaTGTGCGTATCGAGGTCGCCAGCGATGGCCGTGTGGCCAGCATCGAGCCTGGCGCTTCGGCCGAAGACGCCGAGCGCCTGGCCGGCCCGCTGCTGCCCGGCATGCCCCCGCTGCACTCGCATGCCTTCCAGCGCGCCATGGCAGGCCt                                                                            >  1:1320015/1‑150 (MQ=255)                                                                    aTGTGCGTATCGAGGTCGCCAGCGATGGCCGTGTGGCCAGCATCGAGCCTGGCGCTTCGGCCGAAGACGCCGAGCGCCTGGCCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCACTCGCATGCCTTCCAGCGCGCCATGGCAGGCCTgg                                                                          <  1:1324816/150‑1 (MQ=255)                                                                        gCGTATCGAGGTCGCCAGCGATGGCCGTGTGGCCAGCATCGAGCCTGGCGCTTCGGCCGAAGACGCCGAGCGCCTGGCCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCACTCGCATGcc                                                                                                   >  2:1156341/1‑121 (MQ=255)                                                                        gCGTATCGAGGTCGCCAGCGATGGCCGTGTGGCCAGCATCGAGCCTGGCGCTTCGGCCGAAGACGCCGAGCGCCTGGCCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCACTCGCATGcc                                                                                                   <  1:1156341/121‑1 (MQ=255)                                                                         cGTATCGAGGTCGCCAGCGATGGCCGTGTGGCCAGCATCGAGCCTGGCGCTTCGGCCGAAGACGCCGAGCGCCTGGCCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCACTCGCATGCCTTCCAGCGCGCCATGGCAGGCCTGGCGGaa                                                                     <  2:690120/150‑1 (MQ=255)                                                                         cGTATCGAGGTCGCCAGCGATGGCCGTGTGGCCAGCATCGAGCCTGGCGCTTCGGCCGAAGACGCCGAGCGCCTGGCCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCACTCGCATGCCTTCCAGCGCGCCATGGCAGGCCTGGCGGaa                                                                     <  1:1065889/150‑1 (MQ=255)                                                                            aTCGAGGTCGCAAGCGATGGCCGTGTGGCCAGAATCGAGCCTGGCGCTTCGGCCGAAGACGCCGAGCCCCTGGCCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCACTCGCATGCCTTCCAGCGCGCCATGGCAGGCCTGGCGGAAGTc                                                                  <  2:54864/150‑1 (MQ=255)                                                                                          cGATGGCCGTGTGGCCAGCATCGAGCCTTGCGCTTCGGCCGAAGACGCCGAGCGCCTGGCCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCACTCGCATGCCTTCCAGCGCGCCATGGCAGGCCTGGCGGAAGTCGCCGGCAACCCCaa                                                    <  1:757047/150‑1 (MQ=255)                                                                                                 cGTGTGGCCAGCATCGAGCCTGGCGCTTCGGCCGAAGACGCCGAGCGCCTGGCCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCACTCGCATGCCTTCCAGCGCGCCATGGCCGGCCTGGCGGAAGTCGCCGGCAACCCCAACGACAGt                                             >  1:9001/1‑150 (MQ=255)                                                                                                  gtgtGGCCAGCATCGAGCCTGGCGCTTCGGCCGAAGACGCCGAGCGCCTGGCCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCACTCGCATGCCTTCCAGCGCGCCATGGCAGGCCTGGCGGAAGTCGCCGGCAACCCCAACGACAGtt                                            >  1:1126927/1‑150 (MQ=255)                                                                                                     tgGCCAGCATCGAGCCTGGCGCTTCGGCCGAAGACGCCGAGCGCCTGGCCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCACTCGCATGCCTTCCAGCGCGCCATGGCAGGCCTGGCGGAAGTCGCCGGCAACCCCAACGACAGTTTCt                                         <  2:1126927/150‑1 (MQ=255)                                                                                                         cAGCATCGAGCCTGGCGCTTCGGCCGAAGACGCCGAGCGCCTGGTCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCa                                                                                                             >  2:1253763/1‑78 (MQ=255)                                                                                                         cAGCATCGAGCCTGGCGCTTCGGCCGAAGACGCCGAGCGCCTGGTCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCa                                                                                                             <  1:1253763/78‑1 (MQ=255)                                                                                                          aGCATCGAGCCTGGCGCTTCGGCCGAAGACGCCGAGCGCCTGGCCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCACTCGCATGCCTTCCAGCGCGCCATGGCAGGCCTGGCGGAAGTCGCCGGCAACCCCAACGACAGTTTCTGGAcc                                    <  1:729949/150‑1 (MQ=255)                                                                                                             aTCGAGCCTGGCGCTTCGGCCGAAGACGCCGAGCGCCTGGCCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCACTCGCATGCCTTCCAGCGCGCCATGGCAGGCCTGGCGGAAACCGCCGGCAACCCCAACGACAGTTTCTGGACCTgg                                 >  2:744708/1‑150 (MQ=255)                                                                                                              tCGAGCCTGGCGCTTCGGCCGAAGACGCCGAGCGCCTGGCCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCACTCGCATGCCTTCCAGCGCGCCATGGCAGGCCTGGCGGAAGTCGCCGGCAACCCCAACGACAGTTTCTGGACCTGGc                                >  1:594043/1‑150 (MQ=255)                                                                                                                gAGCCTGGCGCTTCGCCCGAAGACGCCGAGCGCCTGGCCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCACTCGCATGCCTTCCAGCGCGCCATGGCAGGCCTGACGGAAGTCGCCGGCAACCCCAACGACAGTTTCTGGACCTGGCga                              >  2:1185934/1‑150 (MQ=255)                                                                                                                 ggccaggcgctcggcgtcttcggccgaagcgccaggcCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCACTCGCATGCCTTCCAGCGCGCCATGGCAGGCCTGGCGGAAGTCGCCGGCAACCCCAACGACa                                               >  2:143996‑M2/38‑132 (MQ=255)                                                                                                                 ggccaggcgctcggcgtcttcggccgaagcgccaggcCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCACTCGCATGCCTTCCAGCGCGCCATGGCAGGCCTGGCGGAAGTCGCCGGCAACCCCAACGACa                                               <  1:143996‑M2/95‑1 (MQ=255)                                                                                                                   ccaggcgctcggcgtcttcggccgaagcgccaggcCGGCCCGCTGCTGCCCGGCATGc                                                                                                                       >  2:1225990‑M2/36‑58 (MQ=255)                                                                                                                   ccaggcgctcggcgtcttcggccgaagcgccaggcCGGCCCGCTGCTGCCCGGCATGc                                                                                                                       <  1:1225990‑M2/23‑1 (MQ=255)                                                                                                                    CTGGCGCTTCGGCCGAAGACGCCGAGCGCCTGGCCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCACTCGCATGCCTTCCAGCGCGCCATGGCAGGCCTGGCGGAAGTCGCCGGCAACCCCAACGACAGTTTCTGGACCTGGCGAGAcc                          <  2:1117447‑M1/150‑1 (MQ=255)                                                                                                                                ccGAAGACGCCGAGCGCCTGGCCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCACTCGCATGCCTTCCAGCGCGCCATGGCAGGCCTGGCGGAAGTCGCCGGCAACCCCAACGACATTTTCTGGACCTTGCGAGACCTGATGTACCGcc              <  1:328790/150‑1 (MQ=255)                                                                                                                                 cGAAGACGCCGAGCGCCTGGCCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCACTCGCATGCCTTCCAGCGCGCCATGGCAGGCCTGGCGGAAGTCGCCGGCAACCCCAACGACAGTTTCTGGACCTGGCGAGACCTGATGTACCGCCt             >  2:845520/1‑150 (MQ=255)                                                                                                                                  gAAGACGCCGAGCGCCTGGCCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCACTCGCATGCCTTCCAGCGCGCCATGGCAGGCCTGGCGGAAGTCGCCGGCAACCCCAACGACAGTTTCTGGACCTGGCGAGACCTGATGTACCGCCTg            >  2:994946/1‑150 (MQ=255)                                                                                                                                  gAAGACGCCGAGCGCCTGGCCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCACTCGCATGCCTTCCAGCGCGCCATGGCAGGCCTGGCGGAAGTCGCCGGCAACCCCAACGACAGTTTCTGGACCTGGCGAGACCTGATGTACCGCCTg            >  1:1181674/1‑150 (MQ=255)                                                                                                                                   aaGACGCCGAGCGCCTGGCCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCACTCGCATGCCTTCCAGCGCGCCATGGCAGGCCTGGc                                                                         <  2:447515/88‑1 (MQ=255)                                                                                                                                   aaGACGCCGAGCGCCTGGCCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCACTCGCATGCCTTCCAGCGCGCCATGGCAGGCCTGGc                                                                         >  1:447515/1‑88 (MQ=255)                                                                                                                                       cGCCGAGCGCCTGGCCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCACTCGCATGCCTTCCAGCGCGCCATGGCAGGCCTGGCGGAAGTCGCCGGCAACCCCAACGACAGTTTCTGGACCTGGCGAGACCTGATGTACCGCCTGGTCgg       <  2:594043/150‑1 (MQ=255)                                                                                                                                          cGAGCGCCTGGCCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCACTCGCATGCCTTCCAGCGCGCCATGGCAGGCCTGGCGGAAGTCGCCGGCAACCCCAACGACAGTTTCTGGACCTGGCGAGACCTGATGTACCGCCTGGTCGGGCa    <  1:994946/150‑1 (MQ=255)                                                                                                                                            aGCGCCTGGCCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCACTCGCATGCCTTCCAGCGCGCCATGGCAGGCCTGGCGGAAGTCGCCGGCAACCCCAACGACAGTTTCTGGACCTGGCGAGACCTGATGTACCGCCTGGTCGGGCAGa  <  1:847238/150‑1 (MQ=255)                                                                                                                                                ccTGGCCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCACTCGCATGCCTTCCAGCGCGCCATGGCAGGCCTGGCGGAAGTCGCCGGCAACCCCAACGACAGt                                             <  2:558299/103‑1 (MQ=255)                                                                                                                                                ccTGGCCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCACTCGCATGCCTTCCAGCGCGCCATGGCAGGCCTGGCGGAAGTCGCCGGCAACCCCAACGACAGt                                             >  1:558299/1‑103 (MQ=255)                                                                                                                                                      |                                                                                                                                            AACAGGATTCCACCGATGTCAGCCTACTTCGCCGAACGCGCCCTGCTACCTACGGGTTGGGCCAGCCATGTGCGTATCGAGGTCGCCAGCGATGGCCGTGTGGCCAGCATCGAGCCTGGCGCTTCGGCCGAAGACGCCGAGCGCCTGGCCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCACTCGCATGCCTTCCAGCGCGCCATGGCAGGCCTGGCGGAAGTCGCCGGCAACCCCAACGACAGTTTCTGGACCTGGCGAGACCTGATGTACCGCCTGGTCGGGCAGA  >  NC_002947/5739112‑5739400

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