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	51 [0]	[0] 51	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)	51 (1.210)	26/166	0.3	100%	coding (657/939 nt)	lrhA	transcriptional regulator LrhA
	?	NZ_CP009273	2399472 =	0 (0.000)					coding (649/939 nt)	lrhA	transcriptional regulator LrhA

GTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAA  >  NZ_CP009273/2399383‑2399550                                                                                          |                                                                               cacgcaggtccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgaga                                                                                >  2:771714‑M2/90‑90 (MQ=255) cacgcaggtccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgaga                                                                                >  2:754123‑M2/90‑90 (MQ=255) cacgcaggtccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgaga                                                                                >  1:824692‑M2/90‑90 (MQ=255)   cgcaggtccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAg                                                                              >  2:702216‑M2/88‑90 (MQ=255)   cgcaggtccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAg                                                                              >  2:368045‑M2/88‑90 (MQ=255)         tccgggctcatcatctcaacctgccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTc                                                                        <  1:649935‑M2/9‑1 (MQ=255)            gggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGcc                                                                     <  2:960717‑M2/12‑1 (MQ=255)             ggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCa                                                                    <  2:555975‑M2/13‑1 (MQ=255)              gctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCAt                                                                   <  1:30236‑M2/14‑1 (MQ=255)               ctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATg                                                                  <  1:911185‑M2/15‑1 (MQ=255)                 catcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTc                                                                        <  1:1232804‑M2/9‑1 (MQ=255)                 catcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTc                                                                        >  2:1232804‑M2/74‑82 (MQ=255)                 catcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGa                                                                >  1:93148‑M2/74‑90 (MQ=255)                 catcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGa                                                                >  2:687904‑M2/74‑90 (MQ=255)                      tctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATc                                                           <  2:755369‑M2/22‑1 (MQ=255)                           accggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCtt                                                      >  2:1058489‑M2/64‑90 (MQ=255)                               gccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTttatt                                                  >  2:163412‑M2/60‑90 (MQ=255)                                   tgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGc                                              >  1:1228298‑M2/56‑90 (MQ=255)                                    gccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCg                                             <  1:1058489‑M2/36‑1 (MQ=255)                                    gccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCg                                             <  1:1015224‑M2/36‑1 (MQ=255)                                         cacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGcc                                        >  1:410015‑M2/50‑90 (MQ=255)                                         cacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGcc                                        <  1:1137573‑M2/41‑1 (MQ=255)                                         cacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGcc                                        >  1:114008‑M2/50‑90 (MQ=255)                                         cacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGcc                                        <  2:410015‑M2/41‑1 (MQ=255)                                              caagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCAc                                   <  2:723393‑M2/46‑1 (MQ=255)                                                     ggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAg                                                                          <  2:159309‑M2/7‑1 (MQ=255)                                                     ggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAg                                                                          >  1:159309‑M2/38‑44 (MQ=255)                                                     ggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCg                                             <  1:129618‑M2/36‑1 (MQ=255)                                                     ggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCg                                             >  2:129618‑M2/38‑73 (MQ=255)                                                          tcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGaaa                       <  2:824692‑M2/58‑1 (MQ=255)                                                          tcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGaaa                       <  1:346009‑M2/58‑1 (MQ=255)                                                           cactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAAc                      <  2:1208984‑M2/59‑1 (MQ=255)                                                           cactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAAc                      <  1:58947‑M2/59‑1 (MQ=255)                                                           cactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAAc                      <  2:528476‑M2/59‑1 (MQ=255)                                                           cactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAAc                      <  1:163412‑M2/59‑1 (MQ=255)                                                           cactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAAc                      <  2:254869‑M2/59‑1 (MQ=255)                                                           cactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAAc                      <  1:962923‑M2/59‑1 (MQ=255)                                                               gccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATc                                                           <  2:346933‑M2/22‑1 (MQ=255)                                                               gccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATc                                                           >  1:346933‑M2/28‑49 (MQ=255)                                                               gccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGc                  <  2:473749‑M2/63‑1 (MQ=255)                                                                ccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCt                 >  1:341551‑M2/27‑90 (MQ=255)                                                                ccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCt                 <  1:189610‑M2/64‑1 (MQ=255)                                                                ccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCt                 <  2:1228298‑M2/64‑1 (MQ=255)                                                                             cggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGc    <  2:383557‑M2/77‑1 (MQ=255)                                                                             cggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGc    <  1:368045‑M2/77‑1 (MQ=255)                                                                             cggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGc    <  1:771714‑M2/77‑1 (MQ=255)                                                                             cggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGc    <  1:687904‑M2/77‑1 (MQ=255)                                                                             cggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGc    <  1:565320‑M2/77‑1 (MQ=255)                                                                             cggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGc    <  1:534312‑M2/77‑1 (MQ=255)                                                                             cggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGc    <  1:398825‑M2/77‑1 (MQ=255)                                                                               gaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCaa  >  2:272925‑M2/12‑90 (MQ=255)                                                                                          |                                                                               GTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAA  >  NZ_CP009273/2399383‑2399550

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