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	31 [0]	[0] 30	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)	30 (1.010)	18/164	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

ACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATC  >  NZ_CP009273/2399375‑2399541                                                                                          |                                                                              acacGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAg                                                                               <  1:833773/90‑1 (MQ=255)        ggTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccag                                                                        >  2:87520‑M1/1‑83 (MQ=255)         gTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagt                                                                       >  1:74536‑M1/1‑82 (MQ=255)            cGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgc                                                                    <  2:70915‑M1/90‑12 (MQ=255)            cGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgc                                                                    >  1:155109‑M1/1‑79 (MQ=255)                  catcatCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatgga                                                              >  1:573318‑M1/1‑73 (MQ=255)                  catcatCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatgga                                                              >  2:572080‑M1/1‑73 (MQ=255)                            aCCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctt                                                    >  1:688566‑M1/1‑63 (MQ=255)                               ggCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttat                                                 >  1:680231‑M1/1‑60 (MQ=255)                               ggCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttat                                                 >  2:584574‑M1/1‑60 (MQ=255)                                gCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttatt                                                <  2:74536‑M1/90‑32 (MQ=255)                                     gCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcg                                           <  1:615133‑M1/90‑37 (MQ=255)                                          cgcgcCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgcc                                      >  2:835731‑M1/3‑49 (MQ=255)                                          cACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgcc                                      <  2:680231‑M1/90‑42 (MQ=255)                                          cACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgcc                                      >  2:590033‑M1/1‑49 (MQ=255)                                          cACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgcc                                      <  1:572080‑M1/90‑42 (MQ=255)                                          cACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgcc                                      <  2:464324‑M1/90‑42 (MQ=255)                                          cACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcg                                           >  2:464073‑M1/1‑49 (MQ=255)                                          cACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcg                                           <  1:464073‑M1/85‑37 (MQ=255)                                               cAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcac                                 <  2:837878‑M1/90‑47 (MQ=255)                                                           tCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaa                     <  2:596278‑M1/90‑59 (MQ=255)                                                            cACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaac                    <  2:809783‑M1/90‑60 (MQ=255)                                                            cACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaac                    <  2:688566‑M1/90‑60 (MQ=255)                                                            cACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcg                                           <  2:89907‑M1/67‑37 (MQ=255)                                                            cACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcg                                           >  1:89907‑M1/1‑31 (MQ=255)                                                                gCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggc                >  1:693256‑M1/1‑27 (MQ=255)                                                                            gCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggctaggatcatcca    <  2:816313‑M1/90‑76 (MQ=255)                                                                              cGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggctaggatcatccagc  <  2:573318‑M1/90‑78 (MQ=255)                                                                              cGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggctaggatcatccagc  <  1:795226‑M1/90‑78 (MQ=255)                                                                              cGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggctaggatcatccagc  <  1:567463‑M1/90‑78 (MQ=255)                                                                              cGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggctaggatcatccagc  <  1:253956‑M1/90‑78 (MQ=255)                                                                                          |                                                                              ACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATC  >  NZ_CP009273/2399375‑2399541

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