New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	5594143 =	61 (1.050)	13 (0.280)	11/208	3.0	20.5%	intergenic (‑39/+106)	PP_4921/thiC	NCS1 family transporter/phosphomethylpyrimidine synthase
	?	NC_002947	5594236 =	53 (1.160)					intergenic (‑132/+13)	PP_4921/thiC	NCS1 family transporter/phosphomethylpyrimidine synthase
Rejected: Coverage evenness skew score above cutoff.

AGGCCGATGCCGAGGGAAAACCACAGTGAAAACAGGTCGCGAGCACCGAAAATGCGCTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACC  >  NC_002947/5593998‑5594289                                                                                                                                                  |                                                                                                                                                   aagCCGATGCCGAGGGAAAACCACAGTGAAAACAGGTCGCGAGCACCGAAAATGCGCTGCTGGATGGGGATTGAATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCt                                                                                                                                                <  1:1343485/148‑1 (MQ=255)        tGCCGAGGGAAAACCACAGTGAAAACAGGTCGCGAGCACCGAAAATGCGCTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAg                                                                                                                                         <  2:412598/150‑1 (MQ=255)                      cacaGTGAAAACAGGTCGCGAGCACCGAAAATGCGCTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCcgcg                                                                                                                           >  2:1183506/1‑150 (MQ=255)                       acaGTGAAAACAGGTCGCGAGCACCGAAAATGCGCTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCtgtg                                                                                                                                             <  1:1018736/131‑1 (MQ=255)                       acaGTGAAAACAGGTCGCGAGCACCGAAAATGCGCTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCtgtg                                                                                                                                             >  2:1018736/1‑131 (MQ=255)                       acaGTGAAAACAGGTCGCGAGCACCGAAAATGCGCTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGa                                                                                                                          >  2:648678/1‑150 (MQ=255)                             aaaaCAGGTCGCGAGCACCGAAAATGCGCTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGGAGGAGTGATGTTAAGTGGCATACAGCACACAGAGAGAGCGGGCTTGCCCGCGAACACCg                                                                                                                    >  2:487161/1‑150 (MQ=255)                                aCAGGTCGCGAGCACCGAAAATGCGCTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCa                                                                                                                 >  1:1103172/1‑150 (MQ=255)                                    gTCGCGAGCACCGAAAATGCGCTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCTTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGc                                                                                                             <  1:487161/150‑1 (MQ=255)                                    gTCGCGAGCACCGAAAATGCGCTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGc                                                                                                             <  2:1023191/150‑1 (MQ=255)                                      cgcgAGCACCGAAAATGCGCTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCg                                                                                                           <  2:1355718/150‑1 (MQ=255)                                        cgAGTACCGAAAATGCGCTGCTGGATGTGTATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGt                                                                                                         <  2:1180625/150‑1 (MQ=255)                                         gAGCACCGAAAATGCGCTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTg                                                                                                        >  2:1251448/1‑150 (MQ=255)                                          aGCACCGAAAATGCGCTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGc                                                                                                       <  2:1103172/150‑1 (MQ=255)                                             aCCGAAAATGCGCTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCAt                                                                                                    >  1:340383/1‑150 (MQ=255)                                                       cgcTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTgg                                                                                                                                            <  2:906629/100‑1 (MQ=255)                                                       cgcTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTgg                                                                                                                                            >  1:906629/1‑100 (MQ=255)                                                          tgctgGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGtgct                                                                                       <  2:472667/150‑1 (MQ=255)                                                          tgctgGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGtgct                                                                                       >  2:1018045/1‑150 (MQ=255)                                                               gATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTgcttgctt                                                                                  >  1:805521/1‑150 (MQ=255)                                                                     gATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCggg                                                                            >  2:1204319/1‑150 (MQ=255)                                                                      aTTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGc                                                                           >  1:707112/1‑150 (MQ=255)                                                                         ggATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAg                                                                        <  2:340383/150‑1 (MQ=255)                                                                          gATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCtt                                                                                                                                  <  2:1410505/91‑1 (MQ=255)                                                                          gATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCtt                                                                                                                                  >  1:1410505/1‑91 (MQ=255)                                                                                   gagagaACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCca                                                              >  2:805993/1‑150 (MQ=255)                                                                                     gagaACTGCCTGGGTGATCTCATGGTCCTGTCCCCAGGAGTGTTGTTAAGTGGCATACACCACCCTGTGGGAGCCGGCTTGCCCGCGAACTCCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTTGCGGGCAAGCCCGCTTCcaca                                                            <  2:380520/150‑1 (MQ=37)                                                                                        aaCTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCggg                                                                            >  1:1253372/1‑131 (MQ=255)                                                                                        aaCTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCggg                                                                            <  2:1253372/131‑1 (MQ=255)                                                                                               tggGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTgg                                                                                                                                            <  2:420535/60‑1 (MQ=255)                                                                                               tggGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTgg                                                                                                                                            >  1:420535/1‑60 (MQ=255)                                                                                               tggGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGtg                                                  <  2:281002/150‑1 (MQ=255)                                                                                                 ggTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCa                                                                                                     <  2:116854/97‑1 (MQ=255)                                                                                                 ggTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCa                                                                                                     >  1:116854/1‑97 (MQ=255)                                                                                                        tcatggtcatgTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGttt                                         >  2:1371338/1‑150 (MQ=255)                                                                                                           tggtcatgTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGGAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCCCAGGGTAAGGTGCTGCGGTTTAGa                                      >  1:49277/1‑150 (MQ=255)                                                                                                              tcatgTCCGCAGGAGCGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCATGCCCGTGAACACCGGTAAAGCCGGTGCCAGGCACCGCAGTGCTTGGTTCGCGGGCAAGCACGCTCCCACAGGGTAAGGCGCTGCGGTTTAGACtt                                   >  1:1073372/1‑150 (MQ=37)                                                                                                                atgTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGTTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTtg                                 <  1:648678/150‑1 (MQ=255)                                                                                                                atgTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTtg                                 <  1:805993/150‑1 (MQ=255)                                                                                                                   tCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTgg                              >  1:1310983/1‑150 (MQ=255)                                                                                                                    ccGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCgg                                                                                                          >  2:346761/1‑73 (MQ=255)                                                                                                                    ccGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCgg                                                                                                          <  1:346761/73‑1 (MQ=255)                                                                                                                     cGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACAGCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAg                                                     >  2:531523/1‑125 (MQ=255)                                                                                                                     cGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACAGCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAg                                                     <  1:531523/125‑1 (MQ=255)                                                                                                                     cGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGtg                                                                                         <  1:369034/89‑1 (MQ=255)                                                                                                                     cGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGtg                                                                                         >  2:369034/1‑89 (MQ=255)                                                                                                                        aGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAg                         >  1:1070832/1‑150 (MQ=255)                                                                                                                               attttAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTaa                                                      >  2:392154/4‑114 (MQ=255)                                                                                                                               aTGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTaa                                                      <  1:392154/114‑1 (MQ=255)                                                                                                                                 gTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTcc                                                                >  1:1143372/1‑102 (MQ=255)                                                                                                                                 gTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTcc                                                                <  2:1143372/102‑1 (MQ=255)                                                                                                                                  ttAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCt               >  1:19200/1‑150 (MQ=255)                                                                                                                                     aGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTcc            <  1:1018045/150‑1 (MQ=255)                                                                                                                                     aGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGAAAACTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTcc            >  1:1114671/1‑150 (MQ=255)                                                                                                                                        ggCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTgg         <  1:1251448/150‑1 (MQ=255)                                                                                                                                        ggCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTgg         <  2:707112/150‑1 (MQ=255)                                                                                                                                         tcttACACCACCATGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCTTGCACCGCAGTGCTTGCTCCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGc        <  2:1073372/147‑1 (MQ=255)                                                                                                                                         gCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGtgct                                                                                       <  2:255600/71‑1 (MQ=255)                                                                                                                                         gCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGtgct                                                                                       >  1:255600/1‑71 (MQ=255)                                                                                                                                         gCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGc        >  1:1176932/1‑150 (MQ=255)                                                                                                                                               accaccCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAAcc  >  1:1365838/1‑150 (MQ=255)                                                                                                                                                  |                                                                                                                                                   AGGCCGATGCCGAGGGAAAACCACAGTGAAAACAGGTCGCGAGCACCGAAAATGCGCTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACC  >  NC_002947/5593998‑5594289

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.