New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 4354950	39 (0.860)	4 (0.090)	4/242	NT	8.9%	coding (106/681 nt)	modB	molybdate ABC transporter permease ModB
	?	NC_002947	= 4354978	44 (1.000)					coding (134/681 nt)	modB	molybdate ABC transporter permease ModB

CGAAAGCCGCAGCGGTGATCAAATCCTATGGTTATGAACGCTGATGCCACTCGATGCCAGTGACCTGGGCGCCATCTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGG  >  NC_002947/4354802‑4355095                                                                                                                                                     |                                                                                                                                                  cGAAAGCCGCAGCGGTGATCAAATCCTATGGTTATGAACGCTGATGCCACTCGATGCCAGTGACCTGGGCGCCATCTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCt                                                                                                                                                  >  1:59278/1‑150 (MQ=255)    aaGCCGCAGCGGTGATCAAATCCTATGGTTATGAACGCTGATGCCACTCGATGCCAGTGACCTGGGCGCCATCTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTcgc                                                                                                                                               <  2:749807/150‑1 (MQ=255)     aGCCGCAGCGGTGATCAAATCCTATGGTTATGAACGCTGATGCCACTCGATGCCAGTGACCTGGGCGCCATCTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTcgcg                                                                                                                                              <  2:932391/150‑1 (MQ=255)                  aTCAAATCCTATGGTTATGAACGCTGATGCCACTCTATGCCAGTGACCTGGGCGCCATCTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCtt                                                                                                                                 <  2:522367/150‑1 (MQ=255)                     aaaTCCTATGGTTATGAACGCTGATGCCACTCGATGCCAGTGACCTGGGCGCCATCTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGc                                                                                                                              >  1:793584/1‑150 (MQ=255)                     aaaTCCTATGGTTATGAACGCTGATGCCACTCGATGCCAGTGACCTGGGCGCCATCTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCCCGCACGCGCTCTTGGc                                                                                                                              >  2:63384/1‑150 (MQ=255)                       aTCCTATGGTTATGAACGCTGATGCCACTCGATGCCAGTGACCTGGGCGCCATCTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCccgcgcagccaagagcgcgtgcg                                                                                                                            >  1:406226‑M1/1‑127 (MQ=255)                          cTATGGTTATGAACGCTGATGACACTCGATGCCAGTGACCTGGGCGCCATCTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTgcgc                                                                                                                         >  2:1039463/1‑150 (MQ=255)                           tATGGTTATGAACGCTGATGCCACTCGATGCCAGTGACCTGGGCGCCATCTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGAAcgcg                                                                                                                                      <  1:307003/136‑1 (MQ=255)                           tATGGTTATGAACGCTGATGCCACTCGATGCCAGTGACCTGGGCGCCATCTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGAAcgcg                                                                                                                                      >  2:307003/1‑136 (MQ=255)                                    gAACGCTGATGCCACTCGATGCCAGTGACCTGGGCGCCATCTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCccgcgcagccaagagcgcgtgcgcgcgagccaccag                                                                                                               <  2:406226‑M1/150‑37 (MQ=255)                                    gAACGCTGATGCCACTCGATGCCAGTGACCTGGGCGCCATCTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCggtgg                                                                                                               <  1:716180/150‑1 (MQ=255)                                     aaCGCTGATGCCACTCGATGCCAGTGACCTGGGCGCCATCTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCggtggg                                                                                                              >  2:328930/1‑150 (MQ=255)                                          tGATGCCACTCGATGCCAGTGACCTGGGCGCCATCTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCgg                                                                                                         >  2:962718/1‑150 (MQ=255)                                                cACTCGATGCCAGTGACCTGGGCGCCATCTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCccgcgcagccaagagcgcgtgcgcgcgagccacca                                                                                                                <  1:375275‑M1/137‑36 (MQ=255)                                                cACTCGATGCCAGTGACCTGGGCGCCATCTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCccgcgcagccaagagcgcgtgcgcgcgagccacca                                                                                                                >  2:375275‑M1/1‑102 (MQ=255)                                                cACTCGATGCCAGTGACCTGGGCGCCATCTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGtggtgg                                                                                                   <  2:214641/150‑1 (MQ=255)                                                   tCGATGCCAGTGACCTGGGCGCCATCTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCAc                                                                                                <  2:793584/150‑1 (MQ=255)                                                          cAGTGACCTGGGCGCCATCTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCt                                                                                         >  2:1048510/1‑150 (MQ=255)                                                           aGTGACCTGGGCGCCATCTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACACGTCACCTGGTGGGGCGCACGCTCGCGCGCTTGGCTGCGCGGGCCGGTGGGAGCGGTGGTGTCGCTGCCGCTg                                                                                        >  1:142600/1‑150 (MQ=255)                                                                   gggCGCCATCTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGcc                                                                                >  2:1159252/1‑150 (MQ=255)                                                                     gcgcCATCTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCAc                                                                              <  2:655570/150‑1 (MQ=255)                                                                           tCTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACgg                                                                        >  2:4766/1‑150 (MQ=255)                                                                                cTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTCCCACCGACGGTGATc                                                                   <  1:340477/150‑1 (MQ=255)                                                                                  gACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGTTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCgg                                                                 <  2:453202/150‑1 (MQ=255)                                                                                     cGTAAAACTGGCCAGCCTGACCCCCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGttt                                                              <  2:162486/150‑1 (MQ=255)                                                                                         aaaCTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTAc                                                          <  1:343715/150‑1 (MQ=255)                                                                                                      tGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGccc                                             >  1:993470/1‑150 (MQ=255)                                                                                                          caccCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTCgg                                         >  2:246423/1‑150 (MQ=255)                                                                                                             ccTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGccc                                      <  1:826989/150‑1 (MQ=255)                                                                                                             ccTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGCGCGGTAGTTGCCATGCCCACACTTCGGCGACGGACGGTGTGCGGTTTTTACCTGTTGATCGCCCTTGGccc                                      <  1:328930/150‑1 (MQ=255)                                                                                                                    ctgctgATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACggc                               >  2:757800/1‑150 (MQ=255)                                                                                                                      gctgATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTTGCCCGCACggctg                             >  2:862241/1‑150 (MQ=255)                                                                                                                           tcatcGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGTTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCCCCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTg                        <  1:267361/150‑1 (MQ=255)                                                                                                                           tcatcGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTg                        <  1:193060/150‑1 (MQ=255)                                                                                                                           tcatcGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTg                        <  2:680801/150‑1 (MQ=255)                                                                                                                           tcatcGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTg                        <  1:4766/150‑1 (MQ=255)                                                                                                                             atcGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTggg                      >  1:1094552/1‑149 (MQ=255)                                                                                                                              tcGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTc                     >  2:982345/1‑150 (MQ=255)                                                                                                                                  cacaCCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGc                 <  1:962718/150‑1 (MQ=255)                                                                                                                                     acCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACggctggct                          >  1:939636/1‑138 (MQ=255)                                                                                                                                     acCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACggctggct                          <  2:939636/138‑1 (MQ=255)                                                                                                                                     acCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCAc              >  2:638071/1‑150 (MQ=255)                                                                                                                                       ccGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCAccc            >  1:117909/1‑150 (MQ=255)                                                                                                                                        ggTCGCCTGGTGGCGCGCGAGCAAGAGCTACAGCCTGCGAGAGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTGCTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCa           <  2:132517/149‑1 (MQ=255)                                                                                                                                              cTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGgccaccaggcgacgggtgtgccgatgatcagcaggatcagggtggtcaggctggccagttttacggtcagccagatggcgcccaggtcactggcatcgagtggcatcagcgttc     <  2:406226‑M2/150‑115 (MQ=255)                                                                                                                                               tggtggCTCGCGCGCACGCGCTCTTGGCTGCGCGGgccaccaggcgacgggtgtgccgatgatcagcaggatcagggtggtcaggctggccagttttacggtcagccagatggcgcccaggtcactggcatcgagtg                 >  2:375275‑M2/1‑35 (MQ=255)                                                                                                                                               tggtggCTCGCGCGCACGCGCTCTTGGCTGCGCGGgccaccaggcgacgggtgtgccgatgatcagcaggatcagggtggtcaggctggccagttttacggtcagccagatggcgcccaggtcactggcatcgagtg                 <  1:375275‑M2/137‑103 (MQ=255)                                                                                                                                               tggtggCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTa    <  1:1048510/150‑1 (MQ=255)                                                                                                                                                 gtggCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCGGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGGTCGGTTTTTCCCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTaag  >  2:441679/1‑148 (MQ=255)                                                                                                                                                 gtggCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAgg  >  2:654988/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                  CGAAAGCCGCAGCGGTGATCAAATCCTATGGTTATGAACGCTGATGCCACTCGATGCCAGTGACCTGGGCGCCATCTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAGG  >  NC_002947/4354802‑4355095

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.