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

ATTTGCTGGCGCTGGGCGTTACCTGGCTGGAGTGCATTAACGGCGGGGATCATCCTCGGTTCCTGGTGGGCCTACTGCGAACTCGGCTGGGGCGGCTGGTGGTTCTGGGATCCGGTGGAAAACGCCTCTTTATTACCCTGGCTTTCTGCCACTGCGCTGCTGCACAGTTTGTCCCTGACACGCCAGCGGGGGATTTTTTGCCACTGGTCGCTGTTGCTGGCGATAGTTACTCTGATGCTGTCGCTGCTGGGCACCTTAATTGTCCGTTCTGGCATTCTGGTTTCGGTTCATGCGTTCGCGCTGGATAACGTCCGCGCCGTGCCGTTGTTCAGCCTGTTTGCACTGATTAGCCTTGCGTCTCTGGCTCTGTATGGCTGGCGAGCGCGGGACGGTGGCCCGGCGGTGCATTTTTCGGGGTTATCGCGGGAAATGTTA  >  CP009273/4282004‑4282438                                                                                                                                                                                                                |                                                                                                                                                                                                                                    tcattactgtccagcgcgaggacactgtcgatcattgcgcgagtctgggcggtcatatgctggctggcctggtaataaagcaccgtcgccagcgccgcatacagttcagcgttctcaccacgtagttgcaacgcctgacggtacgccagcaacgaattgctgtaatcgttttgccacagatagnactcgcccagtaacgcccactgtTCGCTGTTGCTGGCGATAGTTACTCTGATGCTGTCGCTGCTGGCCACCTTAATTGTCCGTTCTGGCAt                                                                                                                                                                  <  1:234269‑M1/68‑1 (MQ=255)            cagcgcgagggctttgtcgatcattgcgcgagtctgggcggtcatatgctggctggcctggtaataaagcaccgtcgccagcgccgcatacagttcagcgttctcaccacgtagttgcaacgcctgacggtacgccagcaacgaattgctgtaatcgttttgccacagatagtactcgcccagtaacgcccactgtTCGCTg                                                                                                                                                                                                                                <  1:140414‑M1/6‑1 (MQ=255)                      ccTGGCTGGAGTGCATTAACGGCGGGGATCATCCTCGGTTCCTGGTGGGCCTACTGCGAACTCGGCTGGGGCGGCTGGTGGTTCTGGGATCCGGTGGAAAACGCCTCTTTATTACCCTGGCTTTCTGCCACTGCGCTGCTGCACAGTTTGTCCCTGACACGCCAGCGGGGGATTTTTTGCCACTGGTCGCTGTTGCTGGCGATAGTTACTCTGATGCTGTCGCTGCTGGGCACCTTAATTGTCCGTTCTGGCATTCTGGTTTCGGTTCATGCGt                                                                                                                                              <  1:100339/274‑1 (MQ=255)                            tagatgattgggggagtctgggcgatcttctgctggctggccttgtaataacgcgccggcgtcagctccgcatacagttcatctttctcaccacgtagttgcaacgcctgacggtacgccagcaacgaattgctgtaatcgttttgccacagatagtactcgcccagtaacgcccactgtTCGCTGTTGCTGGCGATAGTTACTCTGATGCTGTCGCTGCTGGGCACCTTAATTATCCGTTCAGGCAATCTGGTTTCGGTTCCTGCGTTCgcgc                                                                                                                                        >  2:234262‑M1/181‑274 (MQ=255)                                                  tcatcCTCGGNTCCTGGTGGGTATACTGCGAACTCGGCTGGGGCGGCTGGTGGTTCTGGGATCCGGTGGAAAACGCCTCTTTATTACCCTGGCTTTCTGCCACTGCGCTGCTGCACAGTTTGTCCCTGACACGCCAGCGGGGGATTTTTTGCCACTGGTCGCTGTTGCt                                                                                                                                                                                                                           <  1:303606/169‑1 (MQ=255)                                                                                                 ttttggtTTTGGGCTCCGGTTGAAATCGCCTCTTTATTACCCTTGCTTTCTGCCATGGCGCTGCTGCACAGTTTGTCCCTTACACGCCAGCGGGGGATTTTTTGTCACTGGTCGCTGTTGCTGGCGATAGTTACTCTGATGCTGTCGCTGCTGGGCACCTTAATTGTCCGTTCTGGCATTCTGGTTTCGGTTCATGCGTTCGCGCTGGATAACGTCCGCGCCGTGCCGTTGTTCAGCCTGTTTGCACTGATTAGCCTTGCGTCTCTGGCTCTGTa                                                                  <  2:92839/272‑1 (MQ=255)                                                                                                       agttcagcgttctcaccacgtagttgcaacgcctgacggtacgccagcaacgaattgctgtaatcgttttgccacagatagtactcgcccagtaacgcccactgtTCGCTGTTGCTGGCGATAGTTACTCTGATGCTGTCGCTGCTGGGCACCTTAATTGTCCGTTCTGGCATTCTGGTTTCGGTTCATGCGTTCGCGCTGGATAACGTCCGCGCCGTGCCGTTGTTCAGCCTGTTTGCACTGATTAGCCTTGCGTCTCTGGCTCTGTAtggctg                                                            <  1:135071‑M1/170‑1 (MQ=255)                                                                                                                    tGGAAAACGCCTCTTTATTACCCTGGCTTTCTGCCACTGCGCTGCTGCACAGTTTGTCCCTGACACGCCAGCGGGGGATTTTTTGCCACTGGTCGCTGTTGCTGGCGATAGTTACTCTGATGCTGTCGCTGCTGGGCACCTTAATTGTCCGTTCTGGCATTCTGGTTTCGGTTCATGCGTTCGCGCTGGATAACGTCCGCGCCGTGCCGTTGTTCAGCCTGTTTGCa                                                                                               <  1:65404/227‑1 (MQ=255)                                                                                                                                             gCTTTCTGCCACTGCGCTGCTGCACAGTTTGTCCCTGACACGCCAGCGGGGGATTTTTTGCCACTGGTCGCTGTTGCTGGCGATAGTTACTCTGATGCTGTCGCTGCTGGGCACCTTAATTGTc                                                                                                                                                                             >  2:160554/1‑124 (MQ=255)                                                                                                                                                 tCTGCCACTGCGCTGCTGCACAGTTTGTCCCTGACACGCCAGCGGGGGATTTTTTGCCACTGGTCGCTGTTGCTGGCGATAGTTACTCTGATGCTGTCGCTGCTGGGCACCTTAATTGTCCGtt                                                                                                                                                                         <  1:160558/124‑1 (MQ=255)                                                                                                                                                  cTGCCACTGCGCTGCTGCACAGTTTGTCCCTGACACGCCAGCGGGGGATTTTTTGCCACTGGTCGCTGTTGCTGGCGATAGTTACTCTGATGCTGTCGCTGCTGGGCACCTTAATTGTCCGTTCTGGCATTCTGGTTTCGGTTCATGCGTTCGCGCTGGATAACGTCCGCgccg                                                                                                                      >  2:304069/1‑174 (MQ=255)                                                                                                                                                      cACTGCGCTGCTGCACAGTTTGTCCCTGACACGCCAGCGGGGGATTTTTTGCCACTGGTCGCTGTTGCTGGCGATAGTTACTCTGATGCTGTCGCTGCTGGGCACCTTAATTGTCCGTTCTGGCATTCTGGTTTCGGTTCATGCGTTCGCGCTGGATAACGTCCGCgccgtgcc                                                                                                                  <  1:304080/174‑1 (MQ=255)                                                                                                                                                        cTGCGCTGCTGCACTGTTTGTCCCTGACTCGCCCGCGTGGGCCTTTTTGCCACTGTTTGCTGTTGCTGGCGATAGTTCCTCTGACGCTGTCGCTGCTGGGCACCTTAATTTTCCGTTCTGGCATTCTGGTTTCGTTTCATGCGTTCTCGCTGGATaa                                                                                                                                 >  2:75887/1‑157 (MQ=255)                                                                                                                                                           cgcTGCTGCACAGTTTGTCCCTGACACGCCAGCGGGGGATTTTTTGCCACTGGTCGCTGtt                                                                                                                                                                                                                              >  1:105684/1‑61 (MQ=255)                                                                                                                                                            gctgctGCACAGTTTGTCCCTGACACGCCAGCGGGGGATTTTTTGCCACTGGTCGCTGTTGCTGGCGATAGTTACTCTGATGCTGTCGCTGCTGGGCACCTTAATTGTCCGTTCTGGCATTCTGGTTTCGGTTCATGCGTTCGCGCTGGATAACGTc                                                                                                                             <  1:75889/157‑1 (MQ=255)                                                                                                                                                               gctgcACAGTTTGTCCCTGACACGCCAGCGGGGGATTTTTTGCCACTGGTCGCTGTTGCTg                                                                                                                                                                                                                          <  2:105682/61‑1 (MQ=255)                                                                                                                                                                                      gCCAGCGGGGGATTTTTTGCCACTGGTCGCTGTTGCTGGCGATAGTTACTCTGATGCTGTCGCTGCTGGGCACCTTAATTGTCCGTTCTGGCATTCTGGTTTCGGTTCATGCGTTCGCGCTGGATAACGTCCGCGCCGTGCCgttgt                                                                                                             >  1:238285/1‑147 (MQ=255)                                                                                                                                                                                          gCGGGGGATTTTTTGCCACTGGTCGCTGTTGCTGGCGATAGTTACTCTGATGCTGTCGCTGCTGGGCACCTTAATTGTCCGTTCTGGCATTCTGGTTTCGGTTCATGCGTTCGCGCTGGATAACGTCCGCGCCGTGCCGTTGTTCAg                                                                                                         <  2:238278/147‑1 (MQ=255)                                                                                                                                                                                                            cTGGTCGCTGTTGCTGGCGATAGTTACTCTGATGCTGTCGCTGCTGGGCACCTTAATTGTCCGTTCTGGCATTCTg                                                                                                                                                              >  2:226045/1‑76 (MQ=255)                                                                                                                                                                                                             tGGTCGCTGTTGCTGGCGATAGTTACTCTGATGCTGTCGCTGCTGGGCACCTTAATTGTCCGTTCTGGCATTCTGGTTTCGGTTCATGCGTTCGCGCTGGATAACGTCCGCGCCGTGCCGTTGTTCAGCCTGTTTGCACTGATTAGCCTTGCGTCTCTGGCTCTGTATGGCTGGCGAGCGCGGGACGGCGGCCCGGCGGTGCATTTTTCGGGGTTATCGCGGGAAATGTTa  >  1:210594/1‑231 (MQ=255)                                                                                                                                                                                                                TCGCTGTTGCTGGCGATAGTTACTCTGATGCTGTCGCTGCTGGGCACCTTAATTGTCCGTTCTGGCATTCTGGttt                                                                                                                                                          >  1:226051‑M1/1‑76 (MQ=255)                                                                                                                                                                                                                |                                                                                                                                                                                                                                    ATTTGCTGGCGCTGGGCGTTACCTGGCTGGAGTGCATTAACGGCGGGGATCATCCTCGGTTCCTGGTGGGCCTACTGCGAACTCGGCTGGGGCGGCTGGTGGTTCTGGGATCCGGTGGAAAACGCCTCTTTATTACCCTGGCTTTCTGCCACTGCGCTGCTGCACAGTTTGTCCCTGACACGCCAGCGGGGGATTTTTTGCCACTGGTCGCTGTTGCTGGCGATAGTTACTCTGATGCTGTCGCTGCTGGGCACCTTAATTGTCCGTTCTGGCATTCTGGTTTCGGTTCATGCGTTCGCGCTGGATAACGTCCGCGCCGTGCCGTTGTTCAGCCTGTTTGCACTGATTAGCCTTGCGTCTCTGGCTCTGTATGGCTGGCGAGCGCGGGACGGTGGCCCGGCGGTGCATTTTTCGGGGTTATCGCGGGAAATGTTA  >  CP009273/4282004‑4282438

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