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.

CGCGCGACGTTACCGTTTGGTCTGTTGATGCTGGTGGTGATTGTGCTGGCGACGTTTGTCTCTGGCAAACGCGTGCAGCTTCCGGCGCTGGTAGCTCATGCGGGCGTGCTGTTATTTGCCGCTGGGGTCGTGGTCTCCAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGTGCCTCGATCTACAAGCCAAAGGCAATTACACCAGCGAAAAAGCGATAGTGGCACTGTTTGACCATCAGCAA  >  NC_000913/4292634‑4292924                                                                                                                                               |                                                                                                                                                     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                                                                                                                                 >  2:591930/1‑100 (MQ=255)                                                                 gCAAACGCGTGCAGCTTCCGGCGCTGGTAGCTCATGCGGGCGTGCTGTTATTTGCCGCTGGGGTCGTGGTCTCCAGCGTCAGCCGTCAGGAGATCAGCCt                                                                                                                                 <  1:591930/100‑1 (MQ=255)                                                                    aaCGCGTGCAGCTTCCGGCGCTGGTAGCTCATGCGGGCGTGCTGTTATTTGCCGCTGGGGTCGTGGTCTCCAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGtttt                                                                             <  2:540246/149‑1 (MQ=255)                                                                         gTGCAGCTTCCGGCGCTGGTAGCTCATGCGGGCGTGCTGTTATTTGCCGCTGGGGTCGTGGTCTCCAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGTg                                                                        >  1:116708/1‑149 (MQ=255)                                                                           gtaGCGTCGGGCGTTGGTAGCTGTTGCGGGTCTGGTGTTATTTGCTGCTGGGGTCGTGGTTTCTAGCGTCTGCTGTCAGGGGATTAGCCTGAATTTACAGCCGGGTCAGCAGGTGTCGCTGGCAGGATACACCTTCCGTTTTGAGTGcc                                                                      <  2:257482/147‑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                                                    >  2:590414/1‑149 (MQ=255)                                                                                             aGCTCATGCGGGCGTGCTGTTATTTGCCGCTGGGGTCGTGGTCTCCAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGTGCCTCGATCTACAAGCCAAAg                                                    >  1:765962/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)                                                                                                                          cTGGGGTCGTGGTCTCCAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGTGCCTCGATCTACAAGCCAAAGGCAATTACACCAGCGAAAAAGCAATAGTg                       >  2:456112/1‑149 (MQ=255)                                                                                                                           tGGGGTCGTGGTCTCCAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGTGCCTCGATCTACAAGCCAAAGGCAATTACACCAGCGAAAAAGCGATAGTgg                      <  2:21207/149‑1 (MQ=255)                                                                                                                           tGGGGGCGTGGTCTCCAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGTGCCTCGATCTACAAGTCTAGGCGGATTACACCAGCGAAAAAGCGATAGTgg                      <  1:455092/149‑1 (MQ=255)                                                                                                                                   tGGTCTCCAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTACGTTTTGCGGGCCTAGCTTTACAATCCAAAGGCAATTACACCAGCGAAAAAGCGATAGTGGCACTGttt              >  1:860794/1‑149 (MQ=255)                                                                                                                                      tctcCAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGTGCCTCGATCTACAAGCCAAAGGCAATTGCACCAGCGAAAAAGCGATAGTGGCACTGTTTGAc           >  1:712779/1‑149 (MQ=255)                                                                                                                                      cctcCAGCGTAAGCCGCCAGGAGAGCAGCCTGAATTTACAGCCGGGGCAGCAGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGTGCCTCGATCTACAAGCCAAAGGCAATTACACCAGCGAAAAAGCGATAGTGGCACTGTTTGAc           <  1:293048/148‑1 (MQ=255)                                                                                                                                          cAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGTGCCTCGATCTACAAGCCAAAGGCAATTACACCAGCGAAAAAGCGATAGTGGCACTGTTTGACCATc       >  2:474266/1‑149 (MQ=255)                                                                                                                                               tCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGTGCCTCGATCTACAAGCCAAAGGCAATTACACCAGCGAAAAAGCGATAGTGGCACTGTTTGACCATCAGCaa  >  2:339758/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                                                     CGCGCGACGTTACCGTTTGGTCTGTTGATGCTGGTGGTGATTGTGCTGGCGACGTTTGTCTCTGGCAAACGCGTGCAGCTTCCGGCGCTGGTAGCTCATGCGGGCGTGCTGTTATTTGCCGCTGGGGTCGTGGTCTCCAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGTGCCTCGATCTACAAGCCAAAGGCAATTACACCAGCGAAAAAGCGATAGTGGCACTGTTTGACCATCAGCAA  >  NC_000913/4292634‑4292924

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