New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 4053029225 (1.220)4 (0.020) 3/268 NT 1.8% coding (985/1374 nt) hemN coproporphyrinogen III oxidase, SAM and NAD(P)H dependent, oxygen‑independent
?NC_000913 = 4053039 232 (1.310)coding (995/1374 nt) hemN coproporphyrinogen III oxidase, SAM and NAD(P)H dependent, oxygen‑independent
Rejected: Frequency below/above cutoff threshold.

GCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCGTGCTGCATCGTAACTTCCAGGGCTACACCACTCAGGGCGATACCGATCTGCTGGGGATGGGCGTTTCCGCCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4052878‑4053029
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgccaATCATGCTGATGGCGGAAACGCCCATCCCCAGCAGATCGGTATCGCCCTGAGTGGTGTAGCCCTGGAAGTTACGATGCAGCACGCCTTCACGCTGGGCCACCGCCAGCTCGTCATCCGGACGGGCAAAGTGATCCATACCGATAAACTG  <  NC_000913/4053039‑4052891
                                                                                                                                                                                                                                                                                                             
GCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCGTGCTGCATCGTAACTTCCAGGGCTACACCACTCAGGGCGATACCGATCTGCTGGGGATGGGCGTTTCCGC                                                                                                                                                         >  2:609147/1‑150
                TTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCGTGCTGCATCGTAACTTCCAGGGCTACACCACTCAGGGCGATACCGATCTGCTGGGGATGGGCGTTTCCGCCAATCATGCTGATGGC                                                                                                                                         <  1:583732/150‑1
                TTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCGTGCTGCATCGTAACTTCCAGGGCTACACCACTCAGGGCGATACCGATCTGCTGGGGATGGGCGTTTCCGCCAATCATGCTGATGGC                                                                                                                                         <  1:583800/150‑1
                                                                                               TCCAGGGCTACACCACTCAGGGCGATACCGATCTGCTGGGGATGGGCGTTTCCGCCAATCATGCTGATGGCGGAAAC                                                                                                                                   >  1:2870781/1‑77
                                                                                               TCCAGGGCTACACCACTCAGGGCGATACCGATCTGCTGGGGATGGGCGTTTCCGCCAATCATGCTGATGGCGGAAAC                                                                                                                                   <  2:2870781/77‑1
                                                                                                                                                    GCCAATCATGCTGATGGCGGAAACGCCCATCCCCAGCAGATCGGTATCGCCCTGAGTGGTGTAGCCCTGGAAGTTACGATGCAGCACGCCTTCACGCTGGGCCACCGCCAGCTCGTCATCCGGACGGGCAAAGTGATCCATACCGATAAA     <  1:2321447/150‑1
                                                                                                                                                     CCAATCATGCTGATGGCGGAAACGCCCATCCCCAGCAGATCGGTATCGCCCTGAGTGGTGTAGCCCTGGAAGTTACGATGCAGCACGCCTTCACGCTGGGCCACCGCCAGCTCGTCATCCGGACGGGCAAAGTGATCCATACCGATAAAC    >  2:1281765/1‑150
                                                                                                                                                       AATCATGCTGATGGCGGAAACGCCCATCCCCAGCAGATCGGTATCGCCCTGAGTGGTGTAGCCCTGGAAGTTACGATGCAGCACGCCTTCACGCTGGGCCACCGCCAGCTCGTCATCCGGACGGGCAAAGTGATCCATACCGATAAACTG  >  2:1777812/1‑150
                                                                                                                                                       AATCATGCTGATGGCGGAAACGCCCATCCCCAGCAGATCGGTATCGCCCTGAGTGGTGTAGCCCTGGAAGTTACGATGCAGCACGCCTTCACGCTGGGCCACCGCCAGCTCGTCATCCGGACGGGCAAAGTGATCCATACCGATAAACTG  >  2:829008/1‑150
                                                                                                                                                       AATCATGCTGATGGCGGAAACGCCCATCCCCAGCAGATCGGTATCGCCCTGAGTGGTGTAGCCCTGGAAGTTACGATGCAGCACGCCTTCACGCTGGGCCACCGCCAGCTCGTCATCCGGACGGGCAAAGTGATCCATACCGATAAACTG  >  2:829018/1‑150
                                                                                                                                                                                                                                                                                                             
GCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCGTGCTGCATCGTAACTTCCAGGGCTACACCACTCAGGGCGATACCGATCTGCTGGGGATGGGCGTTTCCGCCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4052878‑4053029
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgccaATCATGCTGATGGCGGAAACGCCCATCCCCAGCAGATCGGTATCGCCCTGAGTGGTGTAGCCCTGGAAGTTACGATGCAGCACGCCTTCACGCTGGGCCACCGCCAGCTCGTCATCCGGACGGGCAAAGTGATCCATACCGATAAACTG  <  NC_000913/4053039‑4052891

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