New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3995709 =	49 (0.870)	3 (0.060)	3/250	NT	6.4%	coding (127/708 nt)	yigA	conserved protein YigA
	?	NC_000913	3995732 =	42 (0.790)					coding (150/708 nt)	yigA	conserved protein YigA

GTAGAGTAGGCGATAAAACAGGCCTTCGTTGGCGATAGCCTGTTCCATCAACAGCGCCATGTTCTCTTCCAGAACATGAATATGATTACGTGCGCGGGCCATGTGCCACTCGACCAACGAAACGGTGCCGCGTACCGGATGCGGCAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3995855‑3995709 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcggcac‑CGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCA  >  NC_000913/3995732‑3995871                                                                                                                                                                                                                                                                                                  GTAGCGTAGGCGACAAACCAGGCCTTCGTTGGCGATAGCCTGTTCCATCAACCGCGCCATGTTCTCTGCCAGAACATGACTATGATTACGTGCGCGGGCCATGTGCCACTCGACCAACGAAACGGTGCCGCGTACCGGATGCG                                                                                                                                                   <  2:18210/143‑1                                                                                                         GCCACTCGACCAACGAAACGGTGCCGCGTACCGGATGCGGCAC‑CGTTTCGTTGGTCGAGTGGCACATCGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAA                                          >  1:719435/1‑143                                                                                                                GACCAACGAAACGGTGCCGCGTACCGGATGCGGCAC‑CGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATCGCGCGGTTGAGTGAACAGGCTGTACCCAACGAGGGCCTGc                                   >  1:201970/1‑142                                                                                                                   CAACGAAACGGTGCCGCGTACCGGATGCGGCACACGTATCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTT                                 >  2:503223/1‑143                                                                                                                                             GCGGCAC‑CGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGCCCGGTTTTATCGCCCACTCTGCCTGCAGAGCcgc      >  2:977978/1‑140                                                                                                                                               GGCAC‑CGTTTCGGTGTGCGCGTTCCCCGTAGCCCTCGCAGGTACTCATGTGCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCT    <  2:1065929/143‑1                                                                                                                                                 CAC‑CGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCCGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTCCTCTACCTGCAGCCCAGagcaa  >  2:614171/1‑138                                                                                                                                                                                                                                                                                                  GTAGAGTAGGCGATAAAACAGGCCTTCGTTGGCGATAGCCTGTTCCATCAACAGCGCCATGTTCTCTTCCAGAACATGAATATGATTACGTGCGCGGGCCATGTGCCACTCGACCAACGAAACGGTGCCGCGTACCGGATGCGGCAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3995855‑3995709 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcggcac‑CGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCA  >  NC_000913/3995732‑3995871

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