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.

CCACTATCGCTTTTTCGCTGGTGTAATTGCCTTTGGCTTGTAGATCGAGGCACTCAAAACGGAAGGTGTATCCTGCCAGCGTCACCTGCTGACCCGGCTGTAAATTCAGGCTGATCTCCTGACGGCTGACGCTGGAGACCACGACCCCAGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4292904‑4292753
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccagcgTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGTGCCTCGATCTACAAGCCAAAGGCAATTACACCAGCGAAAAAGCGATAGTGGCACTGTTTGACCATC  >  NC_000913/4292776‑4292919
                                                                                                                                                                                                                                                                                                        
CCACTATCGCTTTTTCGCTGGTGTAATTGCCTTTGGCTTGTAGATCGAGGCACTCAAAACGGAAGGTGTATCCTGCCAGCGTCACCTGCTGACCCGGCTGTAAATTCAGGCTGATCTCCTGACGGCTGACGCTGGAGACCACGACCCCA                                                                                                                                                     >  2:21207/1‑149
CCACTATCGCTTTTTCGCTGGTGTAATCCGCCTAGACTTGTAGATCGAGGCACTCAAAACGGAAGGTGTATCCTGCCAGCGTCACCTGCTGACCCGGCTGTAAATTCAGGCTGATCTCCTGACGGCTGACGCTGGAGACCACGCCCCCA                                                                                                                                                     >  1:455092/1‑149
 CACTATTGCTTTTTCGCTGGTGTAATTGCCTTTGGCTTGTAGATCGAGGCACTCAAAACGGAAGGTGTATCCTGCCAGCGTCACCTGCTGACCCGGCTGTAAATTCAGGCTGATCTCCTGACGGCTGACGCTGGAGACCACGACCCCAG                                                                                                                                                    <  2:456112/149‑1
   cgagCGCTTTTTCGCTGGTGTAATTGCCTTTGGCTTGTAGATCGAGGCACTCAAAACGGAAGGTGTATCCTGCCAGCGTCACCTGCTGACCCGGCTGTAAATTCAGGCTGATCTCCTGACGGCTGACGCTGGAGACCACGACCCCAGCG                                                                                                                                                  <  2:482253/145‑1
   CTATCGCTTTTTCGCTGGTGTAATTGCCTTTGGCTTGTAGATCGAGGCACTCAAAACGGAAGGTGTATCCTGCCAGCGTCACCTGCTGACCCGGCTGTAAATTCAGGCTGATCTCCTGACGGCTGACGCTGGAGACCACGACCCCAGCG                                                                                                                                                  >  1:502427/1‑149
                                                                                                       ATTCAGGCTGATCTCCTGACGGCTGACGCTGGAGACCACGACCCCAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGTGCCTCGATCTACAAGCCAAAGG                                              >  1:196736/1‑149
                                                                                                                            gcgGACACTGGAGACCACGACCCCAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGTGCCTCGATCTACAAGCCTAAGGCAATTACACCAGCGAAAAAGC                         <  1:831128/146‑1
                                                                                                                             CTGACGCTGGAGACCACGACCCCAGCGGCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGTGCCTCGATCTACAAGCCAAAGGCAATTACACCAGCGAAAAAGCG                        >  1:764401/1‑149
                                                                                                                                   CTGGAGACCACGACCCCAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGTGCCTCGATCTACAAGCCAAAG                                               <  1:116274/120‑1
                                                                                                                                   CTGGAGACCACGACCCCAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGTGCCTCGATCTACAAGCCAAAG                                               >  2:116274/1‑120
                                                                                                                                                   CAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGTGCCTCGATCTACAAGCCAAAGGCAATTACACCAGCGAAAAAGCGATAGTGGCACTGTTTGACCATC  >  2:474266/1‑149
                                                                                                                                                                                                                                                                                                        
CCACTATCGCTTTTTCGCTGGTGTAATTGCCTTTGGCTTGTAGATCGAGGCACTCAAAACGGAAGGTGTATCCTGCCAGCGTCACCTGCTGACCCGGCTGTAAATTCAGGCTGATCTCCTGACGGCTGACGCTGGAGACCACGACCCCAGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4292904‑4292753
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccagcgTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGTGCCTCGATCTACAAGCCAAAGGCAATTACACCAGCGAAAAAGCGATAGTGGCACTGTTTGACCATC  >  NC_000913/4292776‑4292919

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