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] 32	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)	31 (1.420)	18/170	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

GTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAGGGA  >  NZ_CP009273/2399383‑2399561                                                                                          |                                                                                          cacgcaggtccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgaga                                                                                           <  1:17218‑M2/1‑1 (MQ=255)        gtccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGt                                                                                    >  2:278072‑M2/83‑90 (MQ=255)         tccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTc                                                                                   <  1:607368‑M2/9‑1 (MQ=255)            gggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGcc                                                                                <  2:570378‑M2/12‑1 (MQ=255)               ctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATg                                                                             <  2:338358‑M2/15‑1 (MQ=255)               ctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATg                                                                             <  2:264312‑M2/15‑1 (MQ=255)               ctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATg                                                                             <  1:48632‑M2/15‑1 (MQ=255)               ctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATg                                                                             <  1:102578‑M2/15‑1 (MQ=255)                 catcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGa                                                                           >  2:397627‑M2/74‑90 (MQ=255)                      tctcaaccggccttgccgtcaccccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATc                                                                      <  1:254001‑M2/22‑1 (MQ=255)                        tcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTg                                                                    <  2:334807‑M2/24‑1 (MQ=255)                           accggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCtt                                                                 >  2:447013‑M2/64‑90 (MQ=255)                              ggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTttat                                                              >  1:99365‑M2/61‑90 (MQ=255)                               gccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTttatt                                                             >  2:562985‑M2/60‑90 (MQ=255)                               gccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTttatt                                                             >  2:6856‑M2/60‑90 (MQ=255)                                       gtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCg                                                     <  2:368282‑M2/39‑1 (MQ=255)                                         cacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGcc                                                   <  2:449982‑M2/41‑1 (MQ=255)                                         cacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGcc                                                   <  2:509271‑M2/41‑1 (MQ=255)                                         cacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGcc                                                   <  1:321371‑M2/41‑1 (MQ=255)                                         cacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGcc                                                   <  2:596636‑M2/41‑1 (MQ=255)                                                          tcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGaaa                                  <  1:562985‑M2/58‑1 (MQ=255)                                                           cactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAAc                                 <  2:194725‑M2/59‑1 (MQ=255)                                                           cactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAAc                                 <  2:232012‑M2/59‑1 (MQ=255)                                                           cactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAAc                                 <  2:170041‑M2/59‑1 (MQ=255)                                                           cactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAAc                                 <  1:527597‑M2/59‑1 (MQ=255)                                                           cactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAAc                                 <  1:505334‑M2/59‑1 (MQ=255)                                                           cactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAAc                                 <  1:159746‑M2/59‑1 (MQ=255)                                                           cactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAAc                                 <  1:132176‑M2/59‑1 (MQ=255)                                                                        acggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGatcat                    <  2:439070‑M2/72‑1 (MQ=255)                                                                             cggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGc               <  1:541787‑M2/77‑1 (MQ=255)                                                                                  gcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATAc          <  1:327798‑M2/82‑1 (MQ=255)                                                                                          aaGCCAGTCGCCATGGAACATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAGGGa  <  1:154744/90‑1 (MQ=255)                                                                                          |                                                                                          GTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAGGGA  >  NZ_CP009273/2399383‑2399561

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