New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4226628 =72 (0.970)3 (0.040) 3/206 7.3 4.2% coding (2801/3684 nt) metH homocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent
?NC_000913 4226644 = 69 (0.970)coding (2817/3684 nt) metH homocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TGGCGTCTTTAAACAGCCGCTGCGCCTCAACGCCCACCACTTCATCTTCCAGAATGCGCGGATACTTCCCGGCCAGCGACCAGGTCATAAAGAACGGTGTCCAGTCGATGTAATTACGCAGCGTTTCGATGCTGGCTTCGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4226768‑4226628
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tcgaAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACAT  >  NC_000913/4226644‑4226775
                                                                                                                                                                                                                                                                                 
TGGCGTCTTTAAACAGCCGCTGCGCCTCAACGCCCACCACTTCATCTTCCAGAATGCGCGGATACTTCCCGGCCAGCGACCAGGTCATAAAGAACGGTGTCCAGTCGATGTAATTACGCAGCGTTTCGATGCTGGCTTC                                                                                                                                        <  2:717589/139‑1
                                                                            CGACCAGGTCATAAAGAACGGTGTCCAGTCGATGTAATTACGCAGCGTTTCGATGCTGGCTTCG                                                                                                                                       <  1:1740750/64‑1
                                                                            CGACCAGGTCATAAAGAACGGTGTCCAGTCGATGTAATTACGCAGCGTTTCGATGCTGGCTTCG                                                                                                                                       >  2:1740750/1‑64
                                                                                        AAAGAACGGTGTCCAGTCGATGTAATTACGCAGCGTTTCGATGCTGGCTT                                                                                                                                         <  1:1684660/50‑1
                                                                                        AAAGAACGGTGTCCAGTCGATGTAATTACGCAGCGTTTCGATGCTGGCTT                                                                                                                                         >  2:1684660/1‑50
                                                                                                                      CAGCGTTTCGATGCTGGCTTCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAG                                                                         <  1:399061/84‑1
                                                                                                                      CAGCGTTTCGATGCTGGCTTCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAG                                                                         >  2:399061/1‑84
                                                                                                                                      GCTTCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACAT  <  1:353860/139‑1
                                                                                                                                                                                                                                                                                 
TGGCGTCTTTAAACAGCCGCTGCGCCTCAACGCCCACCACTTCATCTTCCAGAATGCGCGGATACTTCCCGGCCAGCGACCAGGTCATAAAGAACGGTGTCCAGTCGATGTAATTACGCAGCGTTTCGATGCTGGCTTCGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4226768‑4226628
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tcgaAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACAT  >  NC_000913/4226644‑4226775

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