Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC JC	NC_000913	1,757,082	IS5 (–) +4 bp	26.2%	coding (1385‑1388/1413 nt)	pykF →	pyruvate kinase I

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1397238	NA (NA)	10 (0.110)	10/266	NT	28.8%	noncoding (1195/1195 nt)	IS5	repeat region
	?	NC_000913	= 1757085	24 (0.270)					coding (1388/1413 nt)	pykF	pyruvate kinase I
*	?	NC_000913	1757082 =	24 (0.270)	7 (0.080)	6/266	NT	22.0%	coding (1385/1413 nt)	pykF	pyruvate kinase I
	?	NC_000913	= 2102943	NA (NA)					noncoding (1/1195 nt)	IS5	repeat region

CGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/1397106‑1397238 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TTAGTAGTGCCGCTCGGTACCAGTGCACCAGAAACCATAACTACAACGTCACCTTTGTGTGCCAGACCGCTCTGCAGAG  <  NC_000913/1757085‑1757007                                                                                                                                                                                                                      CGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAGTAGTGCCGCTCG                                                                 <  3:90968/149‑1        AATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAGTAGTGCCGCTCGGTACCA                                                           <  3:118536/148‑1                AGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAGTAGTGCCGCTCGGTACCAGTGCAC                                                     >  3:148705/1‑146                        TCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAGTAGTGCCGCTCGGTACCAGTGCACCAGAAACCAT                                           >  5:371381/1‑148                           AAATAGGCTGATTCAAGGCATTTACGGGAGAGAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAAATTCCCCGCTTATTCGCACCTTCCTTAGTAGTGCCGCTCGGTACCAGAGCACCAGAAACCATAAC                                        >  5:387451/1‑148                                         AAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAGTAGTGCCGCTCGGTACCAGTGCACCAGAAACCATAACTACAACGTCACCTTT                         <  1:183543/149‑1                                           GGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAGTAGTGCCGCTCGGTACCAGTGCACCAGAAACCATAACTACAACGTCACCTTTGT                       <  4:148705/149‑1                                               TTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAGTAGTGCCGCTCGGTACCAGTGCACCAGAAACCATAACTACAACGTCACCTTTGTGTGC                   >  4:266223/1‑149                                                              ATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAGTAGTGCCGCTCGGTACCAGTGCACCAGAAACCATAACTACAACGTCACCTTTGTGTGCCAGACCGCTCTGCAG    <  3:266223/149‑1                                                                 GGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAGTAGTGCCGCTCGGTACCAGTGCACCAGAAACCATAACTACAACGTCACCTTTGTGTGCCAGACCGCTCTGCAGAG  <  6:296754/148‑1                                                                                                                                                                                                                      CGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/1397106‑1397238 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TTAGTAGTGCCGCTCGGTACCAGTGCACCAGAAACCATAACTACAACGTCACCTTTGTGTGCCAGACCGCTCTGCAGAG  <  NC_000913/1757085‑1757007

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.