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	21 [0]	[0] 21	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)	21 (1.020)	16/166	0.2	100%	coding (657/939 nt)	lrhA	transcriptional regulator LrhA
	?	NZ_CP009273	2399472 =	0 (0.000)					coding (649/939 nt)	lrhA	transcriptional regulator LrhA

CGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATC  >  NZ_CP009273/2399378‑2399544                                                                                       |                                                                                 cGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaag                                                                               >  1:291530‑M1/1‑87 (MQ=255)      gTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagt                                                                          >  1:588149‑M1/1‑82 (MQ=255)                    tctcAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatc                                                            <  2:291530‑M1/90‑23 (MQ=255)                       cAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgc                                                         <  1:428629‑M1/90‑26 (MQ=255)                                  gCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcg                                              <  2:341028‑M1/90‑37 (MQ=255)                                       cACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgcc                                         <  1:336032‑M1/90‑42 (MQ=255)                                       cACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgcc                                         <  2:79578‑M1/90‑42 (MQ=255)                                               gACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccat                                 <  2:332026‑M1/90‑50 (MQ=255)                                                   ggCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcg                             >  1:56906‑M1/1‑37 (MQ=255)                                                   ggCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcg                             >  2:345131‑M1/1‑37 (MQ=255)                                                         cACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaac                       <  2:334264‑M1/90‑60 (MQ=255)                                                         cACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaac                       <  2:357288‑M1/90‑60 (MQ=255)                                                         cACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaac                       <  2:409224‑M1/90‑60 (MQ=255)                                                         cACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaac                       <  1:336297‑M1/90‑60 (MQ=255)                                                             gCCGCACGAACGGCCGGAAGCGTCGAGaagccagtc                                                                         >  1:574934‑M1/1‑27 (MQ=255)                                                             gCCGCACGAACGGCCGGAAGCGTCGAGaagccagtc                                                                         <  2:574934‑M1/36‑10 (MQ=255)                                                                 cACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggctagg               <  2:413484‑M1/90‑68 (MQ=255)                                                                           cGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggctaggatcatccagc     <  2:83247‑M1/90‑78 (MQ=255)                                                                              aaGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggctaggatcatccagcaat  <  1:235539‑M1/90‑81 (MQ=255)                                                                                   tCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggct                  <  1:484295‑M1/69‑65 (MQ=255)                                                                                   tCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggct                  >  2:484295‑M1/1‑5 (MQ=255)                                                                                       |                                                                                 CGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATC  >  NZ_CP009273/2399378‑2399544

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