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

CACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGC  >  NZ_CP009273/2399376‑2399548                                                                                         |                                                                                     cacGCAGGTCCGGGCTCACCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGa                                                                                     >  2:492011‑M1/1‑89 (MQ=255)         tCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtc                                                                             <  1:1030331‑M1/90‑10 (MQ=255)           cGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgc                                                                           <  2:463392‑M1/90‑12 (MQ=255)            gggCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgcc                                                                          <  1:138379‑M1/90‑13 (MQ=255)             ggCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgcca                                                                         >  1:1044957‑M1/1‑77 (MQ=255)             ggCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgcca                                                                         >  1:1075094‑M1/1‑77 (MQ=255)                 catcatCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatgga                                                                     >  2:65811‑M1/1‑73 (MQ=255)                 catcatCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatgga                                                                     >  1:635868‑M1/1‑73 (MQ=255)                    catcTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaata                                                                  <  1:297589‑M1/90‑21 (MQ=255)                      tctcAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatc                                                                <  1:547669‑M1/90‑23 (MQ=255)                      tctcAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatc                                                                <  2:39341‑M1/90‑23 (MQ=255)                           aCCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctt                                                           >  2:855482‑M1/1‑63 (MQ=255)                           aCCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagcca                                                                                <  1:613901‑M1/69‑7 (MQ=255)                           aCCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagcca                                                                                >  2:613901‑M1/1‑63 (MQ=255)                               gCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttatt                                                       <  2:852912‑M1/90‑32 (MQ=255)                               gCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtc                                                                             <  2:718922‑M1/68‑10 (MQ=255)                               gCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtc                                                                             >  1:718922‑M1/1‑59 (MQ=255)                                         cACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgcc                                             <  2:723496‑M1/90‑42 (MQ=255)                                            gCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagc                                          >  1:705682‑M1/1‑46 (MQ=255)                                                     ggCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcg                                 >  1:852653‑M1/1‑37 (MQ=255)                                                     ggCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcg                                 >  2:1009706‑M1/1‑37 (MQ=255)                                                     ggCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcg                                 >  2:113198‑M1/1‑37 (MQ=255)                                                     ggCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcg                                 >  1:673820‑M1/1‑37 (MQ=255)                                                          tCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaa                            <  1:941495‑M1/90‑59 (MQ=255)                                                           cACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaac                           >  2:1080063‑M1/1‑31 (MQ=255)                                                           cACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaac                           <  2:673820‑M1/90‑60 (MQ=255)                                                           cACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaac                           <  2:705682‑M1/90‑60 (MQ=255)                                                           cACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaac                           <  2:872321‑M1/90‑60 (MQ=255)                                                           cACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaac                           <  2:350345‑M1/90‑60 (MQ=255)                                                           cACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaac                           <  1:1009706‑M1/90‑60 (MQ=255)                                                           cACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaac                           <  1:910467‑M1/90‑60 (MQ=255)                                                           cACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaac                           <  1:670725‑M1/90‑60 (MQ=255)                                                                                  gCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggctaggatcatccagcaatac    <  2:430154‑M1/90‑83 (MQ=255)                                                                                    gTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggctaggatcatccagcaatacaa  >  2:472315‑M1/1‑6 (MQ=255)                                                                                     tCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaac                           <  1:970654‑M1/64‑60 (MQ=255)                                                                                     tCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaac                           >  2:970654‑M1/1‑5 (MQ=255)                                                                                         |                                                                                     CACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGC  >  NZ_CP009273/2399376‑2399548

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