New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_pae_tpiA = 429142745 (1.280)5 (0.150) 5/264 3.9 10.4% intergenic (+79/‑1) nrfD/nrfE formate‑dependent nitrite reductase, membrane subunit/heme lyase (NrfEFG) for insertion of heme into c552, subunit NrfE
?NC_000913_3_pae_tpiA = 4291429 44 (1.350)coding (2/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.

CGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913_3_pae_tpiA/4291272‑4291427
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cggggtcaacACGTTGACCCCGAGACTTAACAACAACGCCAGAAAACCGACTTCAGGAAGGAATACGTCCAGAAACCACCTCTTTTCTGGCTTACGCCACCGTCAACTGTCCGGCGTAGAGAATGAAAAAGCGCAGCATCAACACGCCGACCAG  <  NC_000913_3_pae_tpiA/4291429‑4291286
                                                                                                                                                                                                                                                                                                            
CGTGGGGGGCCGGTGTGGTCGGCGGGTGGATGGTGCGCTTTTTCATTCTCCGCGTCGGACAGTTGACGGTGGCGTAAGCGAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGG                                                                                                                                                         <  2:433138/149‑1
    TGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAATTCTCGGGGTC                                                                                                                                                     >  2:496948/1‑149
                      CGTGTGGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGACCCCGAGA                                                                                                                                   <  1:251219/149‑1
                                             TTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGACCCCGAGACTTAA                                                                                                                              >  1:511058/1‑131
                                             TTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGACCCCGAGACTTAA                                                                                                                              <  2:511058/131‑1
                                                                               CAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGACCCCGAGACTTAACAACAACGCCAGAAAACC                                                                                                            >  1:288750/1‑115
                                                                               CAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTCTTGTTAAGTCTCGGGGTCAACACGTTGACCCCGAGACTTAACAACAACGCCAGAAAACC                                                                                                            <  2:288750/115‑1
                                                                                                                                                       TCAACACGTTGACCCCGAGACTTAACAACAACGCCAGAAAACCGACTTCAGGAAGGAATACGTCCAGAAACCACCTCTTTTCTGGCTTACGCCACCGTCAACTGTCCGGCGTAGAGAATGAAAAAGCGCAGCATCAACACGCCGACCAG  >  1:320994/1‑149
                                                                                                                                                                                                                                                                                                            
CGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913_3_pae_tpiA/4291272‑4291427
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cggggtcaacACGTTGACCCCGAGACTTAACAACAACGCCAGAAAACCGACTTCAGGAAGGAATACGTCCAGAAACCACCTCTTTTCTGGCTTACGCCACCGTCAACTGTCCGGCGTAGAGAATGAAAAAGCGCAGCATCAACACGCCGACCAG  <  NC_000913_3_pae_tpiA/4291429‑4291286

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.