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	18 [0]	[0] 18	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)	18 (1.350)	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

GTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATC  >  NZ_CP009273/2399383‑2399541                                                                                  |                                                                              gTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagt                                                                       >  2:153279‑M1/1‑82 (MQ=255)      ggCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgcca                                                                  <  1:153279‑M1/90‑14 (MQ=255)        cTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatg                                                                >  1:109645‑M1/1‑75 (MQ=255)          catcatCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCGCTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatgga                                                              >  1:367683‑M1/1‑73 (MQ=255)             catcTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaata                                                           <  2:109645‑M1/90‑21 (MQ=255)               tctcAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatc                                                         <  1:233035‑M1/90‑23 (MQ=255)               tctcAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatc                                                         <  2:339486‑M1/90‑23 (MQ=255)                             gCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcg                                           <  2:71359‑M1/90‑37 (MQ=255)                             gCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcg                                           <  1:209909‑M1/90‑37 (MQ=255)                                  cACGCCAAGACCGGCTTTCGCTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgcc                                      <  2:367683‑M1/90‑42 (MQ=255)                                                   tCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaa                     <  2:283779‑M1/90‑59 (MQ=255)                                                    cACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaac                    <  2:272799‑M1/90‑60 (MQ=255)                                                    cACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaac                    >  2:288137‑M1/1‑31 (MQ=255)                                                         ccGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggct               <  2:196505‑M1/90‑65 (MQ=255)                                                           gCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccatcaccatatcgcgaaacgggctag             >  2:205658‑M1/1‑24 (MQ=255)                                                              cGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggctaggat          <  2:20439‑M1/90‑70 (MQ=255)                                                                    gCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggctaggatcatcca    >  1:357235‑M1/1‑15 (MQ=255)                                                                      cGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggctaggatcatccagc  <  1:288137‑M1/90‑78 (MQ=255)                                                                                  |                                                                              GTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATC  >  NZ_CP009273/2399383‑2399541

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑