New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4227233 =	19 (0.940)	3 (0.160)	3/280	2.1	20.1%	coding (3406/3684 nt)	metH	homocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent
	?	NC_000913	4227352 =	6 (0.310)					coding (3525/3684 nt)	metH	homocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent

GATCCGCGAAAACTACCAGGGCATCCGTCCGGCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTGGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGTACTTCAGCCACCCGGACAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTA  >  NC_000913/4227202‑4227498                                                                                                                                                       |                                                                                                                                                   gATCCGCGAAAACTACCAGGGCATCCGTCCGGCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTGGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCa                                                                                                                                                    >  2:259511/1‑151 (MQ=255)   tCCGCGAAAACTACCAGGGCATCCGTCCGGCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTGGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATc                                                                                                                                                  <  2:192143/151‑1 (MQ=255)             cTACCAGGGCATCCGTCCGGCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTGGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGt                                                                                                                                        >  1:229076/1‑151 (MQ=255)              taccaacccgaaaccgatgcaccgggccacatggcgaaagattctgtgagtttcatgccagtgtgtttttccacttccagcagctcccagatggtggctttttccgtatgttccgggcaggccggatagcccggtgcATCGGTTTCggg                                                                                                                                         <  2:264566‑M2/12‑1 (MQ=255)                    cccgaaaccgatGCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTGGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGTa                                                                                                                                   <  2:264566‑M1/137‑1 (MQ=255)                                                ccTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTGGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGTACTTCAGCCACCCGGACAGCAAGTACTACGc                                                                                                     <  2:264472/151‑1 (MQ=255)                                                              aCGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTGGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGTACTTCAGCCACCCGGACAGCAAGTACTACGCTGTAGCACAAATTc                                                                                       <  1:131732/151‑1 (MQ=255)                                                                           tgtgtttttccactgccagcagctcccagatggtgtctttttccgtatgttccgggcaggccggatagcccggtgcAGCGGTTTCGGGTTGGTACTTCAGCCACCCGGACAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGt                                                                          <  2:42906‑M2/75‑1 (MQ=255)                                                                                                   cccagatggtggctttttccgtatgttccgggcaggccggatagcccggtgcATCGGTTTCGGGTTGGTACTTCAGCCACCCGGACAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAgg                                                  >  2:38880‑M2/53‑151 (MQ=255)                                                                                                                                       ccATGTGGCCCGGTGCATCGGTTTCGGGTTGGTACTTCAGCCCCCCGGACAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCAc               <  1:14604/150‑1 (MQ=255)                                                                                                                                                  ggtgcTTCGTTTTCGGTTTGTTACTTCACCCACCCGAACAGCAAGTACTCCGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAAATTAGGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGACACCGAATCTgggg    <  1:133569‑M2/145‑1 (MQ=255)                                                                                                                                                   gtgcATCGGTTTCGGGTTGGTACTTCAGCCACCCGGACAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTa  >  2:77152‑M2/5‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                   GATCCGCGAAAACTACCAGGGCATCCGTCCGGCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTGGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGTACTTCAGCCACCCGGACAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTA  >  NC_000913/4227202‑4227498

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