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	33 [0]	[0] 33	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)	33 (1.280)	21/168	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

GGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAG  >  NZ_CP009273/2399388‑2399556                                                                                     |                                                                                     aggtccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCa                                                                                 >  2:117892‑M2/85‑90 (MQ=255)   gtccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGt                                                                               <  2:715914‑M2/8‑1 (MQ=255)   gtccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGt                                                                               >  1:687803‑M2/83‑90 (MQ=255)       gggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGcc                                                                           <  2:463378‑M2/12‑1 (MQ=255)         gctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCAt                                                                         <  2:752869‑M2/14‑1 (MQ=255)            catcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTc                                                                              >  1:185810‑M2/74‑82 (MQ=255)            catcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTc                                                                              <  2:185810‑M2/9‑1 (MQ=255)            catcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGa                                                                      >  2:681342‑M2/74‑90 (MQ=255)            catcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGa                                                                      >  2:612879‑M2/74‑90 (MQ=255)                     aaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCt                                                             <  2:523688‑M2/26‑1 (MQ=255)                      accggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCtt                                                            >  2:119527‑M2/64‑90 (MQ=255)                          gccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTttatt                                                        <  2:48888‑M2/31‑1 (MQ=255)                               gccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCg                                                   <  2:469529‑M2/36‑1 (MQ=255)                                    cacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGt                                                                               <  1:565862‑M2/8‑1 (MQ=255)                                    cacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGt                                                                               >  2:565862‑M2/50‑57 (MQ=255)                                    cacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGc                                                                            <  1:50419‑M2/11‑1 (MQ=255)                                    cacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGc                                                                            >  2:50419‑M2/50‑60 (MQ=255)                                    cacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGcc                                              >  2:681801‑M2/50‑90 (MQ=255)                                         caagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCAc                                         <  2:561412‑M2/46‑1 (MQ=255)                                                ggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATcg                                  >  2:250036‑M2/38‑90 (MQ=255)                                                ggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATcg                                  >  2:394701‑M2/38‑90 (MQ=255)                                                      cactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAAc                            <  1:125939‑M2/59‑1 (MQ=255)                                                      cactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAAc                            <  2:553489‑M2/59‑1 (MQ=255)                                                      cactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAAc                            <  1:250036‑M2/59‑1 (MQ=255)                                                      cactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAAc                            >  2:315154‑M2/32‑90 (MQ=255)                                                      cactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAAc                            <  1:394701‑M2/59‑1 (MQ=255)                                                      cactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAAc                            <  1:681342‑M2/59‑1 (MQ=255)                                                      cactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAAc                            <  1:612879‑M2/59‑1 (MQ=255)                                                              cacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAgg                    <  1:229063‑M2/67‑1 (MQ=255)                                                                  aacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAAc                            <  2:63673‑M2/59‑1 (MQ=255)                                                                  aacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAAc                            >  1:63673‑M2/20‑78 (MQ=255)                                                                      gccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCa            <  2:668012‑M2/75‑1 (MQ=255)                                                                                tcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAag  >  1:307891‑M2/6‑90 (MQ=255)                                                                                     |                                                                                     GGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAG  >  NZ_CP009273/2399388‑2399556

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