New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4587017 =	94 (0.950)	3 (0.030)	3/270	NT	3.3%	coding (69/915 nt)	mrr	methylated adenine and cytosine restriction protein
	?	NC_000913	4587057 =	86 (0.890)					coding (109/915 nt)	mrr	methylated adenine and cytosine restriction protein
Rejected: Frequency below/above cutoff threshold.

CTTTGCGACAACCCGGCACGTTTTAAACGGTCATGCGCCCAGCCTGCACGATTTTTATAAACAAGTTGTCCGCTGGTAATGACTTTCGCTCGCTGGCTGTCATCCAGTCCTAATGCATCCGCGGCAGCCTCATGAACATCACGCGCGGCTGCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4587169‑4587017 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tgcaTTAGGACTGGATGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGCGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGG  >  NC_000913/4587057‑4587205                                                                                                                                                                                                                                                                                                                CTTTGCGACAACCCGGCACGTTTTAAACGGTCATGCGCCCAGCCTGCACGATTTTTATAAACAAGTTGTCCGCTGGTAATGACTTTCGCTCGCTGGCTGTCATCCAGTCCTAATGCATCCGCGGCAGCCTCATGAACATCACGCGCGGCT                                                                                                                                                          <  2:920930/150‑1  TTTGCGACAACCCGGCACGTTTTAAACGGTCATGCGCCCAGCCTGCACGATTTTTATAAACAAGTTGTCCGCTGGTAATGACTTTCGCTCGCTGGCTGTCATCCAGTCCTAATGCATCCGCGGCAGCCTCATGAACATCACGCGCGGCTG                                                                                                                                                         <  2:574440/150‑1   TTGCGACAACCCGGCACGTTTTAAACGGTCATGCGCCCAGCCTGCACGATTTTTATAAACAAGTTGTCCGCTGGTAATGACTTTCGCTCGCTGGCTGTCATCCAGTCCTAATGCATCCGCGGCAGCCTCATGAACATCACGCGCGGCTGC                                                                                                                                                        <  1:968409/150‑1                                                                                                            TCCTAATGCATCCGCGGCAGCCTCATGAACATCACGCGCGGCTGCATTAGGACTGGATGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGCGCATGACCGTTTAAAACGTGCCGGGTT                                               >  2:544698/1‑150                                                                                                             CCTAATGCATCCGCGGCAGCCTCATGAACATCACGCGCGGCTGCATTAGGACTGGATGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGCGCATGACCGTTTAAAACGTGCCGGGTTG                                              >  2:607798/1‑150                                                                                                              CTAATGCATCCGCGGCAGCCTCATGAACATCACGCGCGGCTGCATTAGGACTGGATGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGCGCATGACCGTTTAAAACGTGCCGGGTTGT                                             >  1:1409350/1‑150                                                                                                                                                      TGCATTAGGACTGGGTGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGCGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGC     <  1:1217401/150‑1                                                                                                                                                        CATTAGGACTGGATGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGCGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGG   >  1:625075/1‑150                                                                                                                                                        CATTAGGACTGGATGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGCGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGG   >  2:593345/1‑150                                                                                                                                                         ATTAGGACTGGATGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGCGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGG  <  2:1233302/150‑1                                                                                                                                                         ATTAGGACTGGATGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGCGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGG  <  2:1409728/150‑1                                                                                                                                                                                                                                                                                                                CTTTGCGACAACCCGGCACGTTTTAAACGGTCATGCGCCCAGCCTGCACGATTTTTATAAACAAGTTGTCCGCTGGTAATGACTTTCGCTCGCTGGCTGTCATCCAGTCCTAATGCATCCGCGGCAGCCTCATGAACATCACGCGCGGCTGCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4587169‑4587017 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tgcaTTAGGACTGGATGACAGCCAGCGAGCGAAAGTCATTACCAGCGGACAACTTGTTTATAAAAATCGTGCAGGCTGGGCGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGG  >  NC_000913/4587057‑4587205

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