New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 3976261	50 (1.580)	8 (0.260)	5/240	NT	15.6%	coding (651/1014 nt)	PP_3506	magnesium chelatase subunit ChII
	?	NC_002947	= 3976311	38 (1.240)					coding (601/1014 nt)	PP_3506	magnesium chelatase subunit ChII

GCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTCCAGCGCCACGTTCAAGCCGAAACGGTCGAGCAACTGCGGGCGCAGTT  >  NC_002947/3976162‑3976451                                                                                                                                                      |                                                                                                                                             gcgcagcCACACCAGGTCGGCGCGCAGGCCGCCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCGGGTAATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGTGAGCGCTCGCGCAGTTGCgcctgggcctg                                                                                                                                              <  2:720787/150‑1 (MQ=255)  cgcAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAATGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGg                                                                                                                                             >  2:248765‑M2/1‑149 (MQ=255)       ccACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCa                                                                                                                                        <  1:552206/150‑1 (MQ=255)        cacaCCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCTGCGGCCGGGGCCCAg                                                                                                                                       <  2:253544/150‑1 (MQ=255)         acacCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAgt                                                                                                                                      >  2:625137/1‑150 (MQ=255)           acCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCACTGTCGGGCAAGGCCTGGCGCGCCGGGTGGGGGGGCGGGGGGAGTGGCGCCGGGGCCGGGGCCCAgtgt                                                                                                                                    <  2:438064/150‑1 (MQ=255)           acCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAgtgt                                                                                                                                    <  1:212483/150‑1 (MQ=255)                        gcAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCt                                                                                                                       >  1:161215/1‑150 (MQ=255)                          aGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGa                                                                                                                     >  1:492894/1‑150 (MQ=255)                               gTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTc                                                                                                                >  1:772260/1‑150 (MQ=255)                                   aCGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTg                                                                                                            <  1:241458/150‑1 (MQ=255)                                   aCGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTg                                                                                                            >  1:369315/1‑150 (MQ=255)                                     gCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTc                                                                                                          <  2:840869/150‑1 (MQ=255)                                        ggcggcGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGaa                                                                                                       <  2:846218/150‑1 (MQ=255)                                        ggcggcGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGaa                                                                                                       <  2:369315/150‑1 (MQ=255)                                        ggcggcGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGaa                                                                                                       <  1:514833/150‑1 (MQ=255)                                        ggcggcGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGaa                                                                                                       <  2:689529/150‑1 (MQ=255)                                         gcggcgTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAg                                                                                                      <  2:630689/150‑1 (MQ=255)                                         gcggcgTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAg                                                                                                      <  2:161215/150‑1 (MQ=255)                                               tAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGa                                                                                                <  1:647707/150‑1 (MQ=255)                                                    gcgcTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGAcggcg                                                                                           <  1:390525/150‑1 (MQ=255)                                                      gcTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACggcggc                                                                                         >  1:130560/1‑150 (MQ=255)                                                       cTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGgcccaggcccagg                                                                                                                                      <  1:328058/104‑2 (MQ=255)                                                       cTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGgcccaggcccagg                                                                                                                                      >  2:328058/1‑103 (MQ=255)                                                              aTCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGTATCGCTGGCGCAGTTGCGACCGAGCCTTAGCCGACCAAACGCAGATCGCCTGAGCGACATTCTCCAGGGCCAGACGGCGGCGAATGATc                                                                                 <  2:109435/150‑1 (MQ=255)                                                                ccAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATctg                                                                               <  1:792593/150‑1 (MQ=255)                                                                  aGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATctgct                                                                             >  2:392234/1‑150 (MQ=255)                                                                     cAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGcc                                                                          <  2:130560/150‑1 (MQ=255)                                                                               gTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCggg                                                                <  2:209095/150‑1 (MQ=255)                                                                                  aTCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTc                                                             >  2:111287/1‑150 (MQ=255)                                                                                      aGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGggccgacattgccctggatgaccaggcgcttgcctggattaccgagcgctgctacgccgccggcgtggacggcctgcgcgccgac                                                         <  2:30430‑M2/150‑86 (MQ=255)                                                                                           aaTGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGcc                                                    >  2:489311/1‑150 (MQ=255)                                                                                           aaTGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGcc                                                    >  2:539966/1‑150 (MQ=255)                                                                                              gTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTc                                                 <  2:2662/150‑1 (MQ=255)                                                                                              gTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTc                                                 <  1:481727/150‑1 (MQ=255)                                                                                                  gCCAAGGCCTAGCGCGCAGCATGGCAGCGCTCGAGCAGTTGCGAATGCGAATGGACCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTCCAgc                                             >  2:334863/1‑150 (MQ=255)                                                                                                   ccAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGggccgacattgccctggatgaccaggcgcttgcctggattaccgagcgctgctacgccgccggcgtggacggcctgcgcgccgacctggtgtggctgc                                            <  2:614721‑M2/150‑99 (MQ=255)                                                                                                   ccAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGggccgacattgccctggatgaccaggcgcttgcctggattaccgagcgctgctacgccgccggcgtggacggcctgcgcgccgacctggtgtggctgc                                            >  1:30430‑M2/1‑52 (MQ=255)                                                                                                   ccAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGggccgacattgccctggatgaccaggcgcttgcctggattaccgagcgctgctacgccgccggcgtggacggcctgcgcgccgacctggtgtggctgc                                            <  2:555901‑M2/150‑99 (MQ=255)                                                                                                   ccAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGggccgacattgccctggatgaccaggcgcttgcctggattaccgagcgctgctacgccgccggcgtggacggcctgcgcgccgacctggggtggcttc                                            <  2:208050‑M2/150‑99 (MQ=255)                                                                                                   ccAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGgccgacattgccctggatgaccaggcgcttgcctggattaccgagcgctgctacgccgccggcgtggacggcctgcgcgccgacctggtgtggctgcg                                            <  2:465255‑M2/150‑99 (MQ=255)                                                                                                    cAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCg                                                                  <  2:509842/127‑1 (MQ=255)                                                                                                    cAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCg                                                                  >  1:509842/1‑127 (MQ=255)                                                                                                       ggCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTCCAGCGCCAc                                        <  2:659823/150‑1 (MQ=255)                                                                                                           tGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTCCAGCGCCACGTTc                                    <  1:111287/150‑1 (MQ=255)                                                                                                                          cAGCGCTCGCGCAGTTGCGCCTGGGCCTGgccgacattgccctggatgaccaggcgcttgcctggattaccgagcgctgctacgccgccggcgtggacggcctgcgcgccgacctggtgtggctgcgcgctgcgcgcgcccactgtgcct                     >  1:465255‑M2/1‑29 (MQ=255)                                                                                                                           aGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTCCAGCGCCACGTTCAAGCCGAAACGGTCGa                    <  1:351680/150‑1 (MQ=255)                                                                                                                                tCGCGCAGTTGCGCCTGGGCCTGggccgacattgccctggatgaccaggcgcttgcctggattaccgagcgctgctacgccgccggcgtggacggcctgcgcgccgacctggtgtggctgcgcgctgcgcgcgcccactgtgcctggcgt               >  1:614721‑M2/1‑23 (MQ=255)                                                                                                                                     cAGCTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTCCAGCCCCACGTTCAAGCCGAAACGGTCGAGCAACTGCgg          <  2:693037/150‑1 (MQ=255)                                                                                                                                            gcctgggcctggGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTTCAGCGACATGTTCAAGCCGAAACGTGCGAGCAACTGCGGGCGCAGt   >  2:649397/1‑150 (MQ=255)                                                                                                                                             actgggcccaggcCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTCCAGCGCCACGTTCAAGCCGAAACGGTCGAGCAACTGCGGGCGCAGtt  >  2:197263/2‑150 (MQ=255)                                                                                                                                              ctgggcctgggcCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCaa                                                       >  2:650840/1‑96 (MQ=255)                                                                                                                                              ctgggcctgggcCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCaa                                                       <  1:650840/96‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                             GCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTCCAGCGCCACGTTCAAGCCGAAACGGTCGAGCAACTGCGGGCGCAGTT  >  NC_002947/3976162‑3976451

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.