New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 6063028	47 (1.590)	4 (0.140)	4/236	NT	8.1%	intergenic (‑85/‑140)	rubA/ubiC	rubredoxin‑1/chorismate pyruvate‑lyase
	?	NC_002947	= 6063060	46 (1.620)					intergenic (‑117/‑108)	rubA/ubiC	rubredoxin‑1/chorismate pyruvate‑lyase

CCACAGACAATACATTGCCACTTTTTCATCAGGTCTATTTCCTCGGTTCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGCAAGCAGTTGGGGCAATCATGCTAAGCTCGCCGCCTCTCTGGCTGTAACAAAGCATACCGACGTGTCGTACGAATCCCCGCAAGCAGCCGCTGTCGCGTGG  >  NC_002947/6062915‑6063206                                                                                                                                                  |                                                                                                                                                   ccACAGACAATACATTGCCACTTTTTCATCAGGTCTATTTCCTCGGTTCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAgcgc                                                                                                                                                <  2:659600‑M2/150‑5 (MQ=255) ccACAGACAATACATTGCCACTTTTTCATCAGGTCTATTTCCTCGGTTCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAgcgc                                                                                                                                                >  2:22771‑M2/1‑146 (MQ=255) atctcacgtcggtatgctttgttacagccagagaggcggcgagcttagcatgattgccccaactgcttgcatactggcggcagggccgctacaaagcccccaagcattgCGCGCAGGCCTGGCTTGCTTGCCCGCGCACAGGGCCAgcgc                                                                                                                                                >  1:219394‑M1/110‑150 (MQ=255)     aGACAATACATTGCCACTTTTTCATCAGGTCTATTTCCTCGGTTCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCaa                                                                                                                                            >  1:562052/1‑150 (MQ=255)      acgtcggtatgctttgttacagccagagaggcggcgagcttagcatgattgccccaactgcttgcatactggcggcagggccgctacaaagcccccacgcattgCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGgccagcgcgagat                                                                                                                                           <  2:219394‑M1/46‑5 (MQ=255)        cAATACATTGCCACTTTTTCATCAGGTCTATTTCCTCGGTTCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGc                                                                                                                                         >  2:13520/1‑150 (MQ=255)           tACACTGCCACTTTTTCATCAGGTCTATTTCCTCGGTTCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTg                                                                                                                                      <  1:435705/150‑1 (MQ=255)           ggtatgctttgttacagccagagaggcggcgagcttagcatgattgccccaactgcttgcatactggcggcagggccgctacaaagcccccacgcattgCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAgcg                                                                                                                                                 >  2:596554‑M1/100‑139 (MQ=255)           ggtatgctttgttacagccagagaggcggcgagcttagcatgattgccccaactgcttgcatactggcggcagggccgctacaaagcccccacgcattgCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAgcg                                                                                                                                                 <  1:596554‑M1/40‑1 (MQ=255)              aTTGCCACTTTTTCATCAGGTCTATTTCCTCGGTTCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTgggg                                                                                                                                   <  2:529705/150‑1 (MQ=255)              aTTGCCACTTTTTCATCAGGTCTATTTCCTCGGTTCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTgggg                                                                                                                                   <  1:13520/150‑1 (MQ=255)                  ccACTTTTTCATCAGGTCTATTTCCTCGGTTCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCtt                                                                                                                               <  2:225332/150‑1 (MQ=255)                  ccACTTTTTCATCAGGTCTATTTCCTCGGTTCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCtt                                                                                                                               >  2:471545/1‑150 (MQ=255)                  ccACTTTTTCATCAGGTCTATTTCCTCGGTTCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCtt                                                                                                                               <  1:501046/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                                                                                                                        >  1:592015/1‑128 (MQ=255)                                               ttCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCg                                                                                                                        <  2:592015/128‑1 (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                                                                                                         <  2:445149/122‑1 (MQ=255)                                                                    ctgctgATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTa                                                                                                         >  1:445149/1‑122 (MQ=255)                                                                        tgATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCa                                                                                                            >  2:469622/1‑115 (MQ=255)                                                                        tgATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCa                                                                                                            <  1:469622/115‑1 (MQ=255)                                                                            ggCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAgcgcg                                                                                                                                               >  2:82211‑M2/1‑71 (MQ=255)                                                                            ggCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAgcgcg                                                                                                                                               <  1:82211‑M2/76‑6 (MQ=255)                                                                              ccTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTggg                                                                                                                                    <  2:487204/85‑1 (MQ=255)                                                                              ccTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTggg                                                                                                                                    >  1:487204/1‑85 (MQ=255)                                                                                tctTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTgggg                                                                                                                                   >  2:673963/1‑84 (MQ=255)                                                                                tctTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTgggg                                                                                                                                   >  2:673992/1‑84 (MQ=255)                                                                                tctTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTgggg                                                                                                                                   <  1:673992/84‑1 (MQ=255)                                                                                tctTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTgggg                                                                                                                                   <  1:673963/84‑1 (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)                                                                                                            aCCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAgcgcggtcactgtaggag                                                                                                                                  <  2:728123‑M2/57‑19 (MQ=255)                                                                                                            aCCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAgcgcggtcactgtaggag                                                                                                                                  >  1:728123‑M2/1‑39 (MQ=255)                                                                                                             ccGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAgcgcggtcactgtaggagcggcttgcccgcgaagaggccatcagcagcaacacaagccccagcctgaaccgaggaaatagacctgatgaaaaagtggcaa                                                >  2:32977‑M2/1‑38 (MQ=255)                                                                                                             ccGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAgcgcggtcactgtaggagcggcttgcccgcgaagaggccatcagcagcaacacaagccccagcctgaaccgaggaaatagacctgatgaaaaagtggcaa                                                <  1:32977‑M2/138‑101 (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)                                                                                                                                                  |                                                                                                                                                   CCACAGACAATACATTGCCACTTTTTCATCAGGTCTATTTCCTCGGTTCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGCAAGCAGTTGGGGCAATCATGCTAAGCTCGCCGCCTCTCTGGCTGTAACAAAGCATACCGACGTGTCGTACGAATCCCCGCAAGCAGCCGCTGTCGCGTGG  >  NC_002947/6062915‑6063206

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.