New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4226588 =	52 (0.990)	7 (0.130)	6/280	6.1	11.8%	coding (2761/3684 nt)	metH	homocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent
	?	NC_000913	4226623 =	53 (1.010)					coding (2796/3684 nt)	metH	homocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GCTGGCTTCGACTTCCTGCACGCCGAGACGGTGCGCCACCGGCGGCGTGTA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4226638‑4226588 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGAC  >  NC_000913/4226623‑4226764                                                                                                                                                                                                   GCTGGCTTCGACTTCCTGCACGCCGAGACGGTGCGCCACCGGCGGCGTGTAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGG                                              >  2:102573/1‑149             TTCCTGCACGCCGAGACGGTGCGCCACCGGCGGCGTGTAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGG                                  >  2:5283/1‑149                 TGCACGCCGAGACGGTGCGCCACCGGCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGGGGTGGG                              >  1:742015/1‑149                  GCACGCCGAGACGGTGCGCCACCGGCGGCGTGCAGGAAGTCGAAGCGAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGC                             >  1:573679/1‑149                        CGAGACGGTGCGCCACCGGCGGCGTGTAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAG                       <  1:5283/149‑1                            ACGGTGCGCCACCGGCGGCGTGTAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGC                   <  1:102573/149‑1                                             GCGTGAAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGAC  <  2:290055/149‑1                                                                                                                                                                                                   GCTGGCTTCGACTTCCTGCACGCCGAGACGGTGCGCCACCGGCGGCGTGTA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4226638‑4226588 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGAC  >  NC_000913/4226623‑4226764

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