New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4292753 =	52 (0.990)	5 (0.100)	5/270	6.4	9.3%	coding (1242/1659 nt)	nrfE	heme lyase (NrfEFG) for insertion of heme into c552, subunit NrfE
	?	NC_000913	4292776 =	48 (0.950)					coding (1265/1659 nt)	nrfE	heme lyase (NrfEFG) for insertion of heme into c552, subunit NrfE
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CTCTCCGTTGGCGCGCCGTATTTTAACCGCGCGACGTTACCGTTTGGTCTGTTGATGCTGGTGGTGATTGTGCTGGCGACGTTTGTCTCTGGCAAACGCGTGCAGCTTCCGGCGCTGGTAGCTCATGCGGGCGTGCTGTTATTTGCCGCTGGGGTCGTGGTCTCCAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGTGCCTCGATCTACAAGCCAAAGGCAATTACACCAGCGAAAAAGCGATAG  >  NC_000913/4292607‑4292901                                                                                                                                                   |                                                                                                                                                     ctctCCGTTGGCGCGCCGTATTTTAACCGCGCGCCGTTACCGTTTGGTCTGTTGATGCTGGTGGTGATTGTGCTGGCGACGTTTGTCTCTGGCAAACGCGTGCAGCTTCCGGCGCTGGTAGCTCATGCGGGCGTGCTGTTATTTGCCGc                                                                                                                                                    >  2:658613/1‑149 (MQ=255) ctctCCGTTGGCGCGCCGTATTTTAACCGCGCGACGTTACCGTTTGGTCTGTTGATGCTGGTGGTGATTGTGCTGGCGACGTTTGTCTCTGGCAAACGCGTGCAGCTTCCGGCGCTGGTAGCTCATGCGGGCGTGCTGTTATTTGCCGc                                                                                                                                                    >  1:477215/1‑149 (MQ=255) ctctCCGTTGGCGCGCCGTATTTTAACCGCGCGACGTTACCGTTTGGTCTGTTGATGCTGGTGGTGATTGTGCTGGCGACGTTTGTCTCTGGCAAACGCGTGCAGCTTCCGGCGCTGGTAGCTCATGCGGGCGTGCTGTTATTTGCCGc                                                                                                                                                    >  1:386957/1‑149 (MQ=255)      cGTTGGCGCGCCGTATTTTAACCGCGCGACGTTACCGTTTGGTCTGTTGATGCTGGTGGTGATTGTGCTGGCGACGTTTGTCTCTGGCAAACGCGTGCAGCTTCCGGCGCTGGTAGCTCATGCGGGCGTGCTGTTATTTGCCGCTgggg                                                                                                                                               <  1:500283/149‑1 (MQ=255)           gcgcgcCGTATTTTAACCGCGCGACGTTACCGTTTGGTCTGTTGATGCTGGTGGTGATTGTGCTGGCGACGTTTGTCTCTGGCAAACGCGTGCAGCTTCCGGCGCTGGTAGCTCATGCGGGCGTGCTGTTATTTGCCGCTGGGGTCGTg                                                                                                                                          <  2:886611/149‑1 (MQ=255)           gcgcgcCGTATTTTAACCGCGCGACGTTACCGTTTGGTCTGTTGATGCTGGAGGTGATTGTGCTGGCGACGTTTGTCTCTGGCAAACGCGTGCAGCTTCCGGCGCTGGTAGCTCATGCGGGCGTGCTGTTAGTTGCCGCTGGGGTCgcg                                                                                                                                          >  1:800033/1‑147 (MQ=255)               gcCGTATTTTAACCGCGCGACGTTACCGTTTGGTCTGTTGATGCTGGTGGTGATTGTGCTGGCGACGTTTGTCTCTGGCAAACGCGTGCAGCTTCCGGCGCTGGTAGCTCATGCGGGCGTGCTGTTATTTGCCGCTGGGGTCGTGGtct                                                                                                                                      <  2:477215/149‑1 (MQ=255)                ccGTATTTTAACCGCGCGACGTTACCGTTTGGTCTGTTGATGCTGGTGGTGATTGTGCTGGCGACGTTTGTCTCTGGCAAACGCGTGCAGCTTCCGGCGCTGGTAGCTCATGCGGGCGTGCTGTTATTTGCCGCTGGGGTCGTGGtctc                                                                                                                                     >  2:577242/1‑149 (MQ=255)                 cGTATTTTAACCGCGCGACGTTACCGTTTGGTCTGTTGATGCTGGTGGTGATTGTGCTGGCGACGTTTGTCTCTGGCAAACGCGTGCAGCTTCCGGCGCTGGTAGCTCATGCGGGCGTGCTGTTATTTGCCGCTGGGGTCGTGGTCTcc                                                                                                                                    <  1:863965/149‑1 (MQ=255)                    aTTTTAACCGCGCGACGTTACCGTTTGGTCTGTTGATGCTGGTGGTGATTGTGCTGGCGACGTTTGTCTCTGGCAAACGCGTGCAGCTTCCGGCGCTGGTAGCTCATGCGGGCGTGCTGTTATTTGCCGCTGGGGTCGTGGTCTCCAGc                                                                                                                                 <  1:910385/149‑1 (MQ=255)                            cgcgcgACGTTACCGTTTGGTCTGTTGATGCTGGTGGTGATTGTGCTGGCGACGTTTGTCTCTGGCAAACGCGTGCAGCTTCCGGCGCTGGTAGCTTATTCGGGCGTGCTGTTGTTTGCCGCTGGGTGCGTGGTCTCCAACGTCAGCCg                                                                                                                         >  1:425275/1‑149 (MQ=255)                               gcgACGTTACCGTTTGGTCTGTTGATGCTGGTGGCGATTTTTTTGTCGACGTTTGTCTCTGGCAAACGCGTGCAGCTTCCGGCGCTGGTAGCTCATGCGGGCGTGCTGTTATTTGCCGCTGGGGTCGTGGTCTCCAGCGTCAGCCGTCa                                                                                                                      >  1:220569/1‑149 (MQ=255)                                 gACGTTAACGTTTGGTCTGTTGATGCTGGTGGGGATTGTGCTGGCGACGTTTGTCTCTGGCAAACACGTGCAGCTTCCGGCGCTGGTAGCTCATGCGGGCGTGCTGTCGTTTGCCGCTGGGGACGAGGTCGCCAGCGTCAGCCGTCAgg                                                                                                                    >  2:513320/1‑149 (MQ=255)                                  aCGTTACCGTTTGGTCTGTTGATGCTGGTGGTGATTGTGCTGGCGACGTTTGTCTCTGGCAAACGCGTGCAGCTTCCGGCGCTGGTAGCTCATGCGGGCGTGCTGTTATTTGCCGCTGGGGTCGTGGTCTCCAGCGTCAGCCGTCAGga                                                                                                                   <  2:386957/149‑1 (MQ=255)                                                tCTGTTGATGCTGGTGGTGATTGTGCTGGCGACGTTTGTCTCTGGCAAACGCGTGCAGCTTCCGGCGCTGGTAGCTCATGCGGGCGTGCTGTTATTTGCCGCTGGGGTCGTGGTCTCCAGCGTCAGCCGTCAGGAGATCAGCCTGAAtt                                                                                                     >  1:799835/1‑149 (MQ=255)                                                   gTTGATGCTGGTGGTGATTGTGCTGGCGACGTTTGTCTCTGGCAAACGCGTGCAGCTTCCGGCGCTGGTAGCTCATGCGGGCGTGCTGTTATTTGCCGCTGGGGTCGTGGTCTCCAGCGTCAGCCGTCAGGAGATCAGCCTGACTTTAc                                                                                                  >  1:759300/1‑149 (MQ=255)                                                    ttGATGCTGGTGGTGATTGTGCTGGCGACGTTTGTCTCTGGCAAACGCGTGCAGCTTCCGGCGCTGGTAGCTCATGCGGGCGTGCTGTTATTTGCCGCTGGGGTCGTGGTCTCCAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACa                                                                                                 <  1:49845/149‑1 (MQ=255)                                                    ttGATGCTGGTGGTGAACGTATTGGTTAAGTTTGTCTCTGGCAAACGCGTGCAGCTTCCGGCGCTGGTAGCTCATGCGGGCGTGCTGTTATTTGCCGCTGGGGTCGTGGTCTCCAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACa                                                                                                 >  2:293048/1‑149 (MQ=255)                                                        tGCGGGTGGGGATTGTGCTGGCGACGGTTGTCTCTGGCAAACGCGGGCAGCTTCCGGCGCGGGTAGCTCATGCGGGCGTGCTGTTATTTGCCGCTGGGGTCGTGGTCTCCAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCg                                                                                             <  2:254981/149‑1 (MQ=255)                                                         gCTGGTGGTGATTGTGCTGGCGACGTTTGTCTCTGGCAAACGCGTGCAGCTTCCGGCGCTGGTAGCTCATGCGGGCGTGCTGTTATTTGCCGCTGGGGTCGTGGTCTCCAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCgg                                                                                            <  2:125613/149‑1 (MQ=255)                                                           tggtggTGATTGTGCTGGCGACGTTTGTCTCTGGCAAACGCGTGCAGCTTCCGGCGCTGGTAGCTCATGCGGGCGTGCTGTTATTTGCCGCTGGGGTCGTGGTCTCCAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGt                                                                                          >  1:21207/1‑149 (MQ=255)                                                              tggtgATTGTGCTGGCGACGTTTGTCTCTGGCAAACGCGTGCAGCTTCCGGCGCTGGTAGCTCATGCGGGCGTGCTGTTATTTGCCGCTGGGGTCGTGGTCTCCAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTcag                                                                                       >  2:483714/1‑149 (MQ=255)                                                               ggtgATTGTGCTGGCGACGTTTGTCTCTGGCAACCGCGTGAAGGTTCCGCCTTTGGTAGCTCATGCGGGCGTGCTGTTATTTGCCGCTGGGGTCGTGGTCTCCAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTcagc                                                                                      <  2:557801/149‑1 (MQ=255)                                                                   aTTGTGCTGGCGACGTTTGTCTCTGGCAAACGCGTGCAGCTTCCGGCGCTGGTAGCTCATGCGGGCGTGCTGTTATTTGCCGCTGGGGTCGTGGTCTCCAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGt                                                                                  <  1:83297/149‑1 (MQ=255)                                                                            gCGACGTTTGTCTCTGGCAAACGCGTGCAGCTTCCGGCGCTGGTAGCTCATGCGGGCGTGCTGTTATTTGCCGCTGGGGTCGTGGTCTCCAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGc                                                                         <  2:760800/149‑1 (MQ=255)                                                                            gCGACGTTTGTCTCTGGCAAACGCGTGCAGCTTCCGGCGCTGGTAGCTCATGCGGGCGTGCTGTTATTTGCCGCTGGGGTCGTGGTCTCCAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGc                                                                         >  1:32968/1‑149 (MQ=255)                                                                             ctaCGTTTGCCTCTGGCAAACGCGTGCAGCTTCCGTCGCTGGTAGCTCAGGCTGTCGTGCTGTTATTTGCCGCTGGGGTCGTGGTCTCCATCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCa                                                                        <  1:444857/147‑1 (MQ=255)                                                                              gACGTTTGTCTCTGGCAAACGCGTGCAGCTTCCGGCGCTGGTAGCTCATGCGGGCGTGCTGTTATTTGCCGCTGGGGTCGTGGTCTCCAGCGTCAGAAGTCAGTAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAg                                                                       >  1:670988/1‑149 (MQ=255)                                                                                      tctcTGGCAAACGCGTGCAGCTTCCGGCGCTGGTAGCTCATGCGGGCGTGCTGTTATTTGCCGCTGGGGTCGTGGTCTCCAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACAcc                                                               >  1:540246/1‑149 (MQ=255)                                                                                      tctcTGGCAAACGCGTGCAGCTTCCGGCGCTGGTAGCTCATGCGGGCGTGCTGTTATTTGCCGCTGGGGTCGTGGTCTCCAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACAcc                                                               >  1:509776/1‑149 (MQ=255)                                                                                       ctctGGCAAACGCGTGCAGCTTCCGGCGCTGGTAGCTCATGCGGGCGTGCTGTTATTTGCCGCTGGGGTCGTGGTCTCCAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCt                                                              >  1:881798/1‑149 (MQ=255)                                                                                       ctctGGCAAACGCGTGCAGCTTCCGGCGCTGGTAGCTCATGCGGGCGTGCTGTTATTTGCCGCTGGGGTCGTGGTCTCCAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCt                                                              <  1:883330/149‑1 (MQ=255)                                                                                        tctGGCAAACGCGTGCAGCTTCCGGCGCTGGTAGCTCATGCGGGCGTGCTGTTATTTGCCGCTGGGGTCGTGGTCTCCAGCGTCAGCCGTCAGGAGATCAGCCTGAATGTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCtt                                                             >  2:470119/1‑149 (MQ=255)                                                                                         ctGGCAAGCGCGTGCAGCTTCCGGCGCTGGTAGCTCATGCGGGCGTGCTGTTATTTGCCGCTGGGGTCGTGGTCTCCAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTc                                                            <  2:16396/149‑1 (MQ=255)                                                                                            gCAAACGCGTGCAGCTTCCGGCGCTGGTAGCTCATGCGGGCGTGCTGTTATTTGCCGCTGGGGTCGTGGTCTCCAGCGTCAGCCGTCAGGAGATCAGCCt                                                                                                          <  1:591930/100‑1 (MQ=255)                                                                                            gCAAACGCGTGCAGCTTCCGGCGCTGGTAGCTCATGCGGGCGTGCTGTTATTTGCCGCTGGGGTCGTGGTCTCCAGCGTCAGCCGTCAGGAGATCAGCCt                                                                                                          >  2:591930/1‑100 (MQ=255)                                                                                               aaCGCGTGCAGCTTCCGGCGCTGGTAGCTCATGCGGGCGTGCTGTTATTTGCCGCTGGGGTCGTGGTCTCCAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGtttt                                                      <  2:540246/149‑1 (MQ=255)                                                                                                    gTGCAGCTTCCGGCGCTGGTAGCTCATGCGGGCGTGCTGTTATTTGCCGCTGGGGTCGTGGTCTCCAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGTg                                                 >  1:116708/1‑149 (MQ=255)                                                                                                      gtaGCGTCGGGCGTTGGTAGCTGTTGCGGGTCTGGTGTTATTTGCTGCTGGGGTCGTGGTTTCTAGCGTCTGCTGTCAGGGGATTAGCCTGAATTTACAGCCGGGTCAGCAGGTGTCGCTGGCAGGATACACCTTCCGTTTTGAGTGcc                                               <  2:257482/147‑1 (MQ=255)                                                                                                            tCCGGCGCTGGTAGCTCATGCGGGCGTGCTGTTATTTGCCGCTGGGGTCGTg                                                                                                                                          >  1:358229/1‑52 (MQ=255)                                                                                                            tCCGGCGCTGGTAGCTCATGCGGGCGTGCTGTTATTTGCCGCTGGGGTCGTg                                                                                                                                          <  2:358229/52‑1 (MQ=255)                                                                                                               ggCGCTGGTAGCTCATGCGGGCGTGCTGTTATTTGCCGCTGGGGTCGTGGTCTCCAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGTGCCTCGATCTAc                                      <  1:658613/149‑1 (MQ=255)                                                                                                                gcgcTGGTAGCTCATGCGGGCGTGCTGTTATTTGCCGCTGGGGTCGTGGTCTCCAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGTGCCTCGATCTACa                                     >  1:699700/1‑149 (MQ=255)                                                                                                                      gTAGCTCATGCGGGCGTGCTGTTATTTGCCGCTGGGGTCGTGGTCTCCAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCCGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGGGCCTCGATCTACAAGCCaa                               >  1:803545/1‑149 (MQ=255)                                                                                                                        aGCTCATGCGGGCGTGCTGTTATTTGCCGCTGGGGTCGTGGTCTCCAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGTGCCTCGATCTACAAGCCAAAg                             >  1:765962/1‑149 (MQ=255)                                                                                                                        aGCTCATGCGGGCGTGCTGTTATTTGCCGCTGGGGTCGTGGTCTCCAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGTGCCTCGATCTACAAGCCAAAg                             >  2:590414/1‑149 (MQ=255)                                                                                                                         gCTCATGCGGGCGTGCTGTTATTTGCCGCTGGGGTCGTGGTCTCCAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGTGCCTCGATCTACAAGCCAAAgg                            <  2:509776/149‑1 (MQ=255)                                                                                                                           tCATGCGGGCGTGCTGTTATTTGCCGCTGGGGTCGTGGTCTCCAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGTGCCTCGATCTACAAGCCAAAGGCa                          >  1:340680/1‑149 (MQ=255)                                                                                                                           tCATGCGGGCGTGCTGTTATTTGCCGCTGGGGTCGTGGTCTCCAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGCCAGGATACACCTTCCGTTTTGATGGCCTCGATCTAAAACCCACAgcgt                          >  2:461457/1‑146 (MQ=255)                                                                                                                                   gCGTGCTGTTATTTGCCGCTGGGGTCGTGGTCTCCAGCGTCAGCCGTCAGGAGCTCAGCCTGAATTTCCAGCCGGGTCAGCAGGTGACGCTGGCAGGATCCACCCTCCGTTTTGAGCGCATCGATCTCCAAGCACAAGGCCATTACAcc                  >  1:884025/1‑149 (MQ=255)                                                                                                                                           ttATTTGCCGCTGGGGTCGTGGTCTCCAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGTGCCTCGATCTACAAGCCAAAGGCAATTACACCAGCGaaaa          >  1:30273/1‑149 (MQ=255)                                                                                                                                             aTTTGCCGCTGGGGTCGTGGTCTCCAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGTGCCTCGATCTACAAGCCAAAGGCAATTACACCAGCGAAAAAg        <  2:30273/149‑1 (MQ=255)                                                                                                                                                   cGCTGGGGTCGTGGTCTCCAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGTGCCTCGATCTACAAGCCAAAGGCAATTACACCAGCGAAAAAGCgctcg  >  2:482253/1‑145 (MQ=255)                                                                                                                                                   cGCTGGGGTCGTGGTCTCCAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGTGCCTCGATCTACAAGCCAAAGGCAATTACACCAGCGAAAAAGCGATAg  <  1:502427/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                     CTCTCCGTTGGCGCGCCGTATTTTAACCGCGCGACGTTACCGTTTGGTCTGTTGATGCTGGTGGTGATTGTGCTGGCGACGTTTGTCTCTGGCAAACGCGTGCAGCTTCCGGCGCTGGTAGCTCATGCGGGCGTGCTGTTATTTGCCGCTGGGGTCGTGGTCTCCAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGTGCCTCGATCTACAAGCCAAAGGCAATTACACCAGCGAAAAAGCGATAG  >  NC_000913/4292607‑4292901

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