Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC	NC_002947	3,866,157	Δ3 bp	26.7%	coding (763‑765/1020 nt)	PP_3415 ←	LacI family transcriptional regulator

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 3866156	22 (0.690)	8 (0.250)	8/248	NT	26.7%	coding (766/1020 nt)	PP_3415	LacI family transcriptional regulator
	?	NC_002947	3866160 =	22 (0.690)					coding (762/1020 nt)	PP_3415	LacI family transcriptional regulator

GTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/3866073‑3866156 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ACAGAAGAAGATACCGTCGACATCCGGCGCCTGCTGCATCAGCTGGCTGAACAGCACACTGCCCAGGGCAATCGACGATGGCTCCGGGGCCAATACCTCCAGCTCCGCAGCCT  >  NC_002947/3866160‑3866272                                                                                                                                                                                                       GTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGCCGTCACAGAAGAAGATACCGTCGACATCCGGCGCCTGCTGCATCAGCTGGCTGAACAGCACACTGCCCAG                                                 >  2:315041/1‑150      AACCGACCATCGCCACCTGGCGTGGCACCTCCAACCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCACAGAAGAAGATACCGTCGACATCCGGCGCCTGCTGCATCAGCTGGCTGAACAGCACACTGCCCAGGGCAA                                            >  1:707018/1‑150                                        CCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGCCGTCACAGAAGAAGATACCGTCGACATCCGGCGCCTGCTGCATCAGCTGGCTGAACAGCACACTGCCCAGGGCAATCGACGATGGCTCCTGGGCCAATACCTCCAGCTC          <  1:315041/150‑1                                          CTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCACAGAAGAAGATACCGTCGACATCCGGCGCCTGCTGCATCAGCTGGCTGAACAGCACACTGCCCAGGGCAATCGACGATGGCTCCGG                            <  1:361299/130‑1                                          CTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCACAGAAGAAGATACCGTCGACATCCGGCGCCTGCTGCATCAGCTGGCTGAACAGCACACTGCCCAGGGCAATCGACGATGGCTCCGG                            >  2:361299/1‑130                                           TTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCACAGAAGAAGATACCGTCGACATCCGGCGCCTGCTGCATCAGCTGGCTGAACAGC                                                            <  1:734819/97‑1                                           TTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCACAGAAGAAGATACCGTCGACATCCGGCGCCTGCTGCATCAGCTGGCTGAACAGC                                                            >  2:734819/1‑97                                                GCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCACAGAAGAAGATACCGTCGACATCCGGCGCCTGCTGCATCAGCTGGCTGAACAGCACACTGCCCAGGGCAATCGACGATGGCTCCGGGGCCAATACCTCCAGCTCCGCAGCCT  >  2:403876/1‑150                                                                                                                                                                                                       GTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/3866073‑3866156 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ACAGAAGAAGATACCGTCGACATCCGGCGCCTGCTGCATCAGCTGGCTGAACAGCACACTGCCCAGGGCAATCGACGATGGCTCCGGGGCCAATACCTCCAGCTCCGCAGCCT  >  NC_002947/3866160‑3866272

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