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

GAACGCTGATGCCACTCGATGCCAGTGACCTGGGCGCCATCTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCA  >  NC_002947/4354837‑4355121                                                                                                                                              |                                                                                                                                                 gAACGCTGATGCCACTCGATGCCAGTGACCTGGGCGCCATCTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCggtgg                                                                                                                                          <  1:716180/150‑1 (MQ=255)  aaCGCTGATGCCACTCGATGCCAGTGACCTGGGCGCCATCTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCggtggg                                                                                                                                         >  2:328930/1‑150 (MQ=255)       tGATGCCACTCGATGCCAGTGACCTGGGCGCCATCTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCgg                                                                                                                                    >  2:962718/1‑150 (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:4766/150‑1 (MQ=255)                                                                                        tcatcGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTg                                                   <  2:680801/150‑1 (MQ=255)                                                                                        tcatcGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTg                                                   <  1:193060/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)                                                                                                              gtggCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAg‑g                            >  2:654988/1‑150 (MQ=255)                                                                                                                   tggcgagcACGAGCTCTTGGCTGCCTGGGCCGGTGGGTTCGTGGGTGGCACTGCCGCCGGAACTTCCACCGACGGTGATCGGTTTTTAGATGTTGATGGCCCTTGGCTCGCGCGGGTGGCTTGGTCAGGCCACCCAGGCGCTAG‑TCCTgg                       >  2:355278/3‑150 (MQ=255)                                                                                                                   tCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTgg                       >  2:46979/1‑150 (MQ=255)                                                                                                                     gcgcgcACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGC‑CTAGCGCCTgggt                     <  1:757800/150‑2 (MQ=255)                                                                                                                       gcgcACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGC‑CTAGCGCCTgggtgg                   <  1:982345/150‑4 (MQ=255)                                                                                                                        cgcACGCGCTCTTGGCTGCGCGGgccaccaggcgacgggtgtgccgatgatcagcaggatcagggtggtcaggctggccagttttacggtcagccagatggcgcccaggtcactggcatcgagtggcatcagcgttcataaccataggat                   >  1:406226‑M2/1‑23 (MQ=255)                                                                                                                           aCGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGtgtg               >  1:420555/1‑150 (MQ=255)                                                                                                                                  ctTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCa        <  1:1039463/150‑1 (MQ=255)                                                                                                                                   ttGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAg       >  1:638228/1‑150 (MQ=255)                                                                                                                                   ttGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAg       <  1:46979/150‑1 (MQ=255)                                                                                                                                      gCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCtt    >  1:757613/1‑150 (MQ=255)                                                                                                                                        tGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTca  >  1:805128/1‑150 (MQ=255)                                                                                                                                              |                                                                                                                                                 GAACGCTGATGCCACTCGATGCCAGTGACCTGGGCGCCATCTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCA  >  NC_002947/4354837‑4355121

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