New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	6063024 =	47 (1.590)	4 (0.140)	4/236	NT	8.1%	intergenic (‑81/‑144)	rubA/ubiC	rubredoxin‑1/chorismate pyruvate‑lyase
	?	NC_002947	6063066 =	46 (1.620)					intergenic (‑123/‑102)	rubA/ubiC	rubredoxin‑1/chorismate pyruvate‑lyase

AGACAATACATTGCCACTTTTTCATCAGGTCTATTTCCTCGGTTCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGCAAGCAGTTGGGGCAATCATGCTAAGCTCGCCGCCTCTCTGGCTGTAACAAAGCATACCGACGTGTCGTACGAATCCCCGCAAGCAGCCGCTGTCGCGTGG  >  NC_002947/6062919‑6063206                                                                                                                                                    |                                                                                                                                             aGACAATACATTGCCACTTTTTCATCAGGTCTATTTCCTCGGTTCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCaa                                                                                                                                            >  1:562052/1‑150 (MQ=255)    cAATACATTGCCACTTTTTCATCAGGTCTATTTCCTCGGTTCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGc                                                                                                                                         >  2:13520/1‑150 (MQ=255)       tACACTGCCACTTTTTCATCAGGTCTATTTCCTCGGTTCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTg                                                                                                                                      <  1:435705/150‑1 (MQ=255)          aTTGCCACTTTTTCATCAGGTCTATTTCCTCGGTTCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTgggg                                                                                                                                   <  1:13520/150‑1 (MQ=255)          aTTGCCACTTTTTCATCAGGTCTATTTCCTCGGTTCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTgggg                                                                                                                                   <  2:529705/150‑1 (MQ=255)              ccACTTTTTCATCAGGTCTATTTCCTCGGTTCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCtt                                                                                                                               >  2:471545/1‑150 (MQ=255)              ccACTTTTTCATCAGGTCTATTTCCTCGGTTCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCtt                                                                                                                               <  1:501046/150‑1 (MQ=255)              ccACTTTTTCATCAGGTCTATTTCCTCGGTTCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCtt                                                                                                                               <  2:225332/150‑1 (MQ=255)                     ttCATCAGGTCTATTTCCTCGGTTCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCg                                                                                                                        <  2:596208/150‑1 (MQ=255)                            ggTCTATTTCCTCGGTTCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTgc                                                                                                                 >  2:605297/1‑150 (MQ=255)                                       tCGGTTCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGc                                                                                                      <  2:336158/150‑1 (MQ=255)                                         ggTTCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGCaa                                                                                                    <  1:117534/150‑1 (MQ=255)                                           ttCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCg                                                                                                                        <  2:592015/128‑1 (MQ=255)                                           ttCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCg                                                                                                                        >  1:592015/1‑128 (MQ=255)                                                  tGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGCAAGCAGTTggg                                                                                           <  1:137112/150‑1 (MQ=255)                                                   ggggCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGCAAGCAGTTgggg                                                                                          <  1:735081/150‑1 (MQ=255)                                                    gggCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGCAAGCAGTTGGGGc                                                                                         >  1:533762/1‑150 (MQ=255)                                                      gCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGAAAGAAGTAGGTGCaa                                                                                       >  2:454541/1‑150 (MQ=255)                                                                ctgctgATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTa                                                                                                         >  1:445149/1‑122 (MQ=255)                                                                ctgctgATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTa                                                                                                         <  2:445149/122‑1 (MQ=255)                                                                    tgATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCa                                                                                                            <  1:469622/115‑1 (MQ=255)                                                                    tgATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCa                                                                                                            >  2:469622/1‑115 (MQ=255)                                                                          ccTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTggg                                                                                                                                    >  1:487204/1‑85 (MQ=255)                                                                          ccTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTggg                                                                                                                                    <  2:487204/85‑1 (MQ=255)                                                                            tctTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTgggg                                                                                                                                   <  1:673963/84‑1 (MQ=255)                                                                            tctTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTgggg                                                                                                                                   <  1:673992/84‑1 (MQ=255)                                                                            tctTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTgggg                                                                                                                                   >  2:673963/1‑84 (MQ=255)                                                                            tctTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTgggg                                                                                                                                   >  2:673992/1‑84 (MQ=255)                                                                            tctTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGCAAGCAGTTGGGGCAATCATGCTAAGCTCGCCGCctct                                                                 >  1:518132/1‑150 (MQ=255)                                                                             ctTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTg                                                                                                                             >  2:558010/1‑89 (MQ=255)                                                                             ctTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTg                                                                                                                             <  1:558010/89‑1 (MQ=255)                                                                                cgcgGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGCAAGCAGTTGGGGCAATCATGCTAAGCTCGCCGCCTCTCTgg                                                             <  2:518132/150‑1 (MQ=255)                                                                                     gCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGCACGCCTCGGAGGCAATCATGCTAAGCTCGCCGCCTCTCTGGCTGTa                                                        >  1:407191/1‑150 (MQ=255)                                                                                          ccGCCCCTACAGTGACCGCGCGGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGCAAGCAGTTGGGGCAATCATGCTAAGCTCGCCGCCTCTCTGGCTGTAACaaa                                                   <  2:695034/150‑1 (MQ=255)                                                                                                cTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGCAAGCAGTTGGGGCAATCATGCTAAGCTCGCCGCCTCTCTGGCTGTAACAAAGCATAc                                             >  2:96242/1‑150 (MQ=255)                                                                                                cTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCATGCCGCCAGTATGCAAGCAGTTGGGGCAATCATGCTAAGCTCGCCGCCTCTCTGGCTGTAACAAAGCATAc                                             >  1:13422/1‑150 (MQ=255)                                                                                                 tACAGTGACCGCGCAGGCCTGTCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGCAAGCAGTTGGGGCAATCATGCTAAGCTCGCCGCCTCTCTGGCTGTAACAAAGCATAcc                                            >  1:63817/1‑150 (MQ=255)                                                                                                    aGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATACGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGCAAGCAGTTGGGGCAATCATGCTAAGCTCGCCGCCTCTCTGGCTGTAACAAAGCATACCGAc                                         >  1:297858/1‑150 (MQ=255)                                                                                                      atctcgcgctggccctgtgcgcgggcaagcccgccaggcctgcgcgCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGCAAGCAGTTGGGGCAATCATGCTAAGCTCGCCGCCTCTCTGGCTGTAACAAAGCATACCGACgt                                       <  2:219394‑M2/104‑1 (MQ=255)                                                                                                           gcgctggccctgtgcgcgggcaagcaagccaggcctgcgcgCAATGCTTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGCAAGCAGTTGGGGCAATCATGCTAAGCTCGCCGCCTCTCTGGCTGTAACAAAGCATACcgacgtgagat                                  >  1:219394‑M2/42‑146 (MQ=255)                                                                                                            cgctggccctgtgcgcgggcaagcccgccaggcctgcgcgCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGCAAGCAGTTGGGGCAATCATGCTAAGCTCGCCGCCTCTCTGGCTGTAACAAAGCATAcc                                            <  1:596554‑M2/99‑1 (MQ=255)                                                                                                            cgctggccctgtgcgcgggcaagcccgccaggcctgcgcgCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGCAAGCAGTTGGGGCAATCATGCTAAGCTCGCCGCCTCTCTGGCTGTAACAAAGCATAcc                                            >  2:596554‑M2/41‑139 (MQ=255)                                                                                                                        gggCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGCAAGCAGTTGGGGCAATCATGCTAAGCTCGCCGCCTCTCTGGCTGTAACAAAGCATACCGACGTGTCGTACGAATCCCCGCa                     >  1:793280/1‑150 (MQ=255)                                                                                                                         ggCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGCAAGCAGTTGGGGCAATCATGCTAAGCTCGCCGCCTCTCTGGCTGTAACAAAGCATACCGACGTGTCGTACGAATCCCCGCaa                    <  1:226928/150‑1 (MQ=255)                                                                                                                            ttGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGCAAGCAGTTGGGGCAATCATGCTAAGCTCGCCGCCTCTCTGGCTGTAACAAAGCATACCGACGTGTCGTACGAATCCCCGCAagca                 >  1:60900/1‑150 (MQ=255)                                                                                                                             tGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGCAAGCAGTTGGGGCAATCATGCTAAGCTCGCCGCCTCTCTGGCTGTAACAAAGCATACCGACGTGTCGTACGAATCCCCGCAagcag                <  2:533762/150‑1 (MQ=255)                                                                                                                               cccGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGCAAGCAGTTGGTTCAATCATGCTAAGCTCGCCGCCTCTCTGGCTGTAACAAAGCATACCGACGTGTCGTACGAATCCCCGCAAGCAGcc              <  2:63817/150‑1 (MQ=255)                                                                                                                                 cgcgCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTATCGGCCCTGCCGCCAGTATGCAAGCAGTTGGGGCAATCATGCTAAGCTCGCCGCCTCTCTGGCTGTAACAAAGCATACCGACGTGTCGTACGAATCCCCGCAAGCAGCCGc            >  2:481245/1‑150 (MQ=255)                                                                                                                                  gcgcACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGCAAGCAGTTGGGGCAATCATGCTAAGCTCGCCGCCTCTCTGGCTGTAACAAAGCATACCGACGTGTCGTACGAATCCCCGCAAGCAGCCGCt           >  1:401940/1‑150 (MQ=255)                                                                                                                                    gcACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGCAAGCAGTTGGGGCAATCATGCTAAGCTCGCCGCCTCTCTGGCTGTAACAAAGCATACCGACGTGTCGTACGAATCCCCGCAAGCAGCCGCTGt         <  1:96242/150‑1 (MQ=255)                                                                                                                                         gggCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGCAAGCAGTTGGGGCAATCATGCTAAGCTCACCGCCTCTCTGGCAATCACAATGCATACCGACGTGTCGTACGAATCCCCGCAAGCAGCCGCTGTCGCGt    <  1:454541/150‑1 (MQ=255)                                                                                                                                           gCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGCAAGCAGTTGGGGCAATCATGCTAAGCTCGCCGCCTCTCTGGCTGTAACAAAGCATACCGACGTGTCGTACGAATCCCCGCAAGCAGCCGCTGTCGCGTgg  <  1:471545/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                             AGACAATACATTGCCACTTTTTCATCAGGTCTATTTCCTCGGTTCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGCAAGCAGTTGGGGCAATCATGCTAAGCTCGCCGCCTCTCTGGCTGTAACAAAGCATACCGACGTGTCGTACGAATCCCCGCAAGCAGCCGCTGTCGCGTGG  >  NC_002947/6062919‑6063206

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