New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	CP009273	4282211 =	14 (0.810)	3 (0.180)	3/288	NT	11.7%	coding (772/1659 nt)	nrfE	heme lyase (NrfEFG) for insertion of heme into c552, subunit NrfE
	?	CP009273	4283692 =	32 (1.950)					coding (222/597 nt)	nrfG	heme lyase (NrfEFG) for insertion of heme into c552, subunit NrfG

TGTTGTTACTGCTGATGGCACTGATCCTCTGGCGAGTGAGGGCGAAGCGATGAAACAGCCCAAAATACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGGCGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGATGCTGCTGGCTTCCGATGCGTTTATGCAGGCGAACTACGCGCAAGCCAT  >  CP009273/4283422‑4283961                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                              catttttACTGCTGGTGGCACTGGTCTGCTGGCGAGTGAGGCCGAAGCGATGAAACAGCCCAAAATACCGGGGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCANTTATCTGTTAAGACCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGt                                                                                                                                                                                                                                                                           <  1:287419/270‑1 (MQ=255)               aagaacatgtgagtctGGGGAGGGAGGGTGGAGCGATGGAACGGCCCAAAATACCGGGGAAATTGCTTACAACACGCACTATATTGTGGGTTTTTCTCTGGGTCGTCAGTTATCTGGTAAGCCCAACATGGCAGGCGGTACGTGCGGAGTATCAGCGTAAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTggg                                                                                                                                                                                                                                                              <  2:280801/261‑1 (MQ=255)                tctcccTGATCCTCTGGCGAGTGAGGGCGAACCGATTAAACAGCCCAAAAATCCGGTGAAAATGCTGACAACCCTCACGATTTTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGGc                                                                                                                                                                                                                                                             <  1:86667/269‑1 (MQ=255)                             ctGGCGAGTGAGGGCGAAGCGATGAAACAGCCCAAAATACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTCTCTGTTAAGCCCAAAATGGCAGTCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAACCAGCGAACATCAGACGATACTGGACGCTTACTATCtgtg                                                                                                                                                                                                                                               >  1:134004/1‑275 (MQ=255)                                  gAGTGAGGGCGAAGCGATGAAACAGCCCAAAATACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGGCGAGTACTATCTGTGGCaaa                                                                                                                                                                                                                                          >  1:137312/1‑275 (MQ=255)                                                                            aaTGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATCAACTCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGGCGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTg                                                                                                                                                                                                 >  1:60228/1‑274 (MQ=255)                                                                                                tatTGATGTTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGGCGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCt                                                                                                                                                                            <  1:170676/275‑1 (MQ=255)                                                                                                     cagccatacagagccagagacgcaaggctaatcagtgcaaacaggctgaacaacggcacggcgcggacgttatccagcgcgaacgcatgaaccgaaaccagaatgccagaacggacaattaaggtgcccagcagcgacagcatcagagtaactatcgccagcaacagcgaACAGTGGGCGTTACTGGGCGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACt                                                                                                                                                                       <  1:135071‑M2/105‑1 (MQ=255)                                                                                                                                      aGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCg                                                                                                                                                                                                                                                                      >  1:74677/1‑147 (MQ=255)                                                                                                                                        cccAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATACGCGCTAATCCACAAAACAGCGAACAGTGGACGTTACTGGGCGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCaa                                                                                                                                     >  1:54926/1‑273 (MQ=255)                                                                                                                                          cAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTAc                                                                                                                                                                                                                                                                  <  2:74675/147‑1 (MQ=255)                                                                                                                                             aagggCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGACAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGGCGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGgccagc                                                                                                                               <  1:198799/272‑1 (MQ=255)                                                                                                                                              agggaagtcGGTACGTGCGGATCATCGGCGTCGGCGCGATCCGTTACATAAGTGTGGCAGCCAGCGAAACCCGTAAGCGCATCTTGAGGAATGGCAGGATAAAATCCGCGCTAATCTACAAAACAGGGAACAGTGGGCGTTAGTGGGCGACTACTATCTGTGGCAGAACGATTACAGGAATTCGTGGCTCGCGTACCGTTAGTCGTTGCAACAACGTGGTGAGAACGCTGAACAGAATGCGGCGCTGGCGACGGTGCTTTATTACCAGgccagc                                                                                                                               <  2:333320/266‑1 (MQ=255)                                                                                                                                                         tACGTGCGGACTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGGCGAGGACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTAGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGATGGCGACGGTGCTTTATTACCAGGCCAGCCAGCAAATGAcc                                                                                                                   >  1:333332/1‑275 (MQ=255)                                                                                                                                                                              agcgCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGGCGAGTACTATCTGTGGCaaaa                                                                                                                                                                                                                                         >  2:258523/1‑136 (MQ=255)                                                                                                                                                                                 gcgcgaacgcaggaaccgaaaccagattgcctgaacggataattaaggtgcccagcagcgacagcatcagagtaactatcgccagcaacagcgaACAGTGGGCGTTACTGGGCGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAAAGATGAACTGTATGCGGAGCTGACGCCGGCGCGTTATTACAAGGCCAGCCAGCAGAAGATCGCCCAGACTCCCCcaatcatcta                                                                                            >  2:234262‑M2/95‑272 (MQ=255)                                                                                                                                                                                  cgATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGGCGAGTACTATCTGTGGCAAAACGAt                                                                                                                                                                                                                                     <  1:258531/136‑1 (MQ=255)                                                                                                                                                                                    aTCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGGCGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACg                                                                                                                                                                                         >  2:193303/1‑178 (MQ=255)                                                                                                                                                                                        gCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGGCGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACgtggt                                                                                                                                                                                     <  1:193307/178‑1 (MQ=255)                                                                                                                                                                                                   ttGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGGCGTGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGa                                                                          <  1:180586/274‑1 (MQ=255)                                                                                                                                                                                                      ccagccagCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGGCGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACgtg                                                                                                                                                                                       >  1:274154/1‑162 (MQ=255)                                                                                                                                                                                                          ccagcAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGGCGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTga                                                                                                                                                                                   <  2:274143/162‑1 (MQ=255)                                                                                                                                                                                                          ccagcAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGGCGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGc                                                                                                                                                        >  2:146417/1‑189 (MQ=255)                                                                                                                                                                                                           atgccagaacggacaattaaggtggccagcagcgacagcatcagagtaactatcgccagcaacagcgaACAGTGGGCGTTACTGGGCGAGTNCTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAGTGTCCTCGCGCTGGACAGTAATGa                                                                 <  1:234269‑M2/207‑1 (MQ=255)                                                                                                                                                                                                              cAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGGCGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACg                                                                                                                                                    <  1:146420/189‑1 (MQ=255)                                                                                                                                                                                                                                gatCAGGCATGGCAGGAAAAAATCCCCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGGCGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGGTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGATgctgct                                            <  2:137310/273‑1 (MQ=255)                                                                                                                                                                                                                                          tGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGGCGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTAcc                                                                                                                                                                                                            >  1:191023/1‑105 (MQ=255)                                                                                                                                                                                                                                              ggATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGGCGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCa                                                                                                                                                                                                        <  2:191019/105‑1 (MQ=255)                                                                                                                                                                                                                                               gATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGGCGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGa                                                                                                                                                      >  2:320525/1‑154 (MQ=255)                                                                                                                                                                                                                                                   aaaTCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGGCGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGt                                                                                                                                                  <  1:320537/154‑1 (MQ=255)                                                                                                                                                                                                                                                    aaTCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGGCGAGTa                                                                                                                                                                                                                                                        >  2:55819/1‑51 (MQ=255)                                                                                                                                                                                                                                                    aaTCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGGCGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGc                                                                                                                                                                                                     >  1:32491/1‑102 (MQ=255)                                                                                                                                                                                                                                                        cgcgCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGGCGAGTACTAt                                                                                                                                                                                                                                                    <  1:55821/51‑1 (MQ=255)                                                                                                                                                                                                                                                        cgcgCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGGCGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTg                                                                                                                                                                                                 <  2:32491/102‑1 (MQ=255)                                                                                                                                                                                                                                                          cgcTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGGCGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCTCTg                                                                                                                                                          >  2:24367/1‑139 (MQ=255)                                                                                                                                                                                                                                                              aaTCCACAAAACAGCGAACAGTGGGCGTTACTGGGCGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGa                                                                                                                                                      <  1:24367/139‑1 (MQ=255)                                                                                                                                                                                                                                                                     aaaaCAGCGAACAGTGGGCGTTACTGGGCGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTcgc                                                                                >  2:140411/1‑202 (MQ=255)                                                                                                                                                                                                                                                                         cagcgaACAGTGGGCGTTACTGGGCGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTg                                                                            <  1:140414‑M2/196‑1 (MQ=255)                                                                                                                                                                                                                                                                           gagcaCAGGGGGCGTGACTGTCCGAGTACTATCTGTGGCCAAACGAGTACAGCAAGTCGTTGCTGGCGTACGGTCCGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTGATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGATGCTGCTGGCTTCCGATGCGTTTATGCAGGCGAACTACGCGCAAGCCAt  <  2:60226/270‑1 (MQ=255)                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                              TGTTGTTACTGCTGATGGCACTGATCCTCTGGCGAGTGAGGGCGAAGCGATGAAACAGCCCAAAATACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGGCGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGATGCTGCTGGCTTCCGATGCGTTTATGCAGGCGAACTACGCGCAAGCCAT  >  CP009273/4283422‑4283961

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