New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4226590 =48 (1.410)5 (0.150) 5/256 3.5 9.7% coding (2763/3684 nt) metH homocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent
?NC_000913 4226621 = 48 (1.480)coding (2794/3684 nt) metH homocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TCAACGCCCACCACTTCATCTTCCAGAATGCGCGGATACTTCCCGGCCAGCGACCAGGTCATAAAGAACGGTGTCCAGTCGATGTAATTACGCAGCGTTTCGATGCTGGCTTCGACTTCCTGCACGCCGAGACGGTGCGCCACCGGCGGCGTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4226742‑4226590
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cggcgtgCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAA  >  NC_000913/4226621‑4226761
                                                                                                                                                                                                                                                                                                      
tccACGCCCACCACTTCATCTTCCAGAATGCGCGGATACTTCCCGGCCAGCGACCAGGTCATAAAGAACGGTGTCCAGTCGATGTAATTACGCAGCGTTTCGATGCTGGCTTCGACTTCCTGCACGCCGAGACGGTGCGCCACCGGCGG                                                                                                                                                   >  8:102842/4‑149
TCAACGCCCACCACTTCATCTTCCAGAATGCGCGGATACTTCCCGGCCAGCGACCAGGTCATAAAGAACGGTGTCCAGTCGATGTAATTACGCAGCGTTTCGATGCTGGCTTCGACTTCCTGCACGCCGAGACGGTGCGCCACCGGCGG                                                                                                                                                   >  8:101020/1‑149
  AACGCCCACCACTTCATCTTCCAGAATGCGCGGATACTTCCCGGCCAGCGACCAGGTGATAAAGAACGGTGTCCAGTCGATGTAAGTACGCAGCGTTTCGATGCTGGCTTCGACTTCCTCCACGCCGAGAAGGTGCGCCACAGGCGGCG                                                                                                                                                 >  4:86100/1‑149
                                                                                                                    TTCCTGCACGCCGAGACGGTGCGCCACCGGCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGG                               >  5:80033/1‑149
                                                                                                                        TGCACGCCGAGACGGTGCGCCACCGGCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAA                                                            <  5:126062/116‑1
                                                                                                                        TGCACGCCGAGACGGTGCGCCACCGGCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAA                                                            >  6:126062/1‑116
                                                                                                                           ACGCCGAGACGGTGCGCCACCGGCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTCTCCGCGCATTCTGGAAGATGAAGTGGTGGGCGT                        >  4:73102/1‑149
                                                                                                                            CGCCGAGACGGTGCGCCACCGGCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTT                       <  6:80033/149‑1
                                                                                                                                                   GGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAA  <  2:65249/147‑1
                                                                                                                                                                                                                                                                                                      
TCAACGCCCACCACTTCATCTTCCAGAATGCGCGGATACTTCCCGGCCAGCGACCAGGTCATAAAGAACGGTGTCCAGTCGATGTAATTACGCAGCGTTTCGATGCTGGCTTCGACTTCCTGCACGCCGAGACGGTGCGCCACCGGCGGCGTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4226742‑4226590
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cggcgtgCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAA  >  NC_000913/4226621‑4226761

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