Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NZ_CP009273	2,399,465	Δ7 bp	coding (650‑656/939 nt)	lrhA ←	transcriptional regulator LrhA

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NZ_CP009273	2399465	2399471	7	20 [0]	[0] 20	lrhA	transcriptional regulator LrhA

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	= 2399464	0 (0.000)	20 (1.580)	16/168	0.1	100%	coding (657/939 nt)	lrhA	transcriptional regulator LrhA
	?	NZ_CP009273	2399472 =	0 (0.000)					coding (649/939 nt)	lrhA	transcriptional regulator LrhA

GTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NZ_CP009273/2399383‑2399464 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑AAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGAT  >  NZ_CP009273/2399472‑2399540                                                                                                                                                         GTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGAAGCCAGT                                                               >  1:244899/1‑90  TCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGAAGCCAGTC                                                              <  1:222244/90‑1  TCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGAAGCCAGTC                                                              <  1:363495/90‑1       GCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGAAGCCAGTCGCCAT                                                         <  1:4436/90‑1        CTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGAAGCCAGTCGCCATG                                                        <  2:244899/90‑1             CATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGAAGCCAGTCGCCATGGAATA                                                   <  2:258610/90‑1                       GGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGAAGCCAGTCGCCATGGAATATC                                                 <  1:150668/82‑1                       GGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGAAGCCAGTCGCCATGGAATATC                                                 >  2:150668/1‑82                        GCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGAAGCCAGTCGCCATGGAATATCTGCTTTATT                                        <  2:177816/90‑1                        GCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGAAGCCAGTCGCCATGGAATATCCGCTTTATT                                        >  2:212273/1‑90                             GCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCG                                   <  1:295739/90‑1                             GCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCG                                   <  1:117840/90‑1                                  CACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCC                              >  1:260144/1‑90                                  CACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCC                              <  1:19582/90‑1                                  CACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCC                              <  2:41370/90‑1                                    CGCCAGGACCGGCTTTCACTGCCGCACGACCGGCCGGAAGCGTCGAGAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAG                            <  2:362897/90‑1                                     GCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGC                           <  2:284889/90‑1                                                    CACTGCCGCACGAACGGCCGGAAGCGTCGAGAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAAC            <  2:310739/90‑1                                                    CACTGCCGCACGAACGGCCGGAAGCGTCGAGAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAAC            <  2:42210/90‑1                                                              CGAACGGCCGGAAGCGTCGAGAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGAT  <  2:234598/90‑1                                                                                                                                                         GTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NZ_CP009273/2399383‑2399464 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑AAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGAT  >  NZ_CP009273/2399472‑2399540

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