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	26 [0]	[0] 26	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)	26 (0.960)	19/170	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

GCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAA  >  NZ_CP009273/2399389‑2399555                                                                                    |                                                                                    ggtccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAg                                                                               >  2:180592‑M2/84‑90 (MQ=255)  gtccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGt                                                                              >  1:251603‑M2/83‑90 (MQ=255)   tccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTc                                                                             <  2:222245‑M2/9‑1 (MQ=255)   tccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTc                                                                             <  1:449673‑M2/9‑1 (MQ=255)      gggctcatcgtcacaaccggccttgccgtcacgcgaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGcc                                                                          <  1:124561‑M2/12‑1 (MQ=255)      gggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGcc                                                                          >  2:24434‑M2/79‑90 (MQ=255)              cacctcgacccgccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAata                                                                  <  1:23509‑M2/20‑1 (MQ=255)                tctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATc                                                                <  1:42748‑M2/22‑1 (MQ=255)                     accggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCtt                                                           <  1:137287‑M2/27‑1 (MQ=255)                        ggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTttat                                                        >  1:359962‑M2/61‑90 (MQ=255)                             tgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGc                                                   >  2:247833‑M2/56‑90 (MQ=255)                                 gtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCg                                               <  1:247833‑M2/39‑1 (MQ=255)                                               ggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATcg                                 >  1:497910‑M2/38‑90 (MQ=255)                                                     cactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAAc                           <  2:494697‑M2/59‑1 (MQ=255)                                                     cactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAAc                           <  2:473224‑M2/59‑1 (MQ=255)                                                     cactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAAc                           <  2:472064‑M2/59‑1 (MQ=255)                                                     cactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAAc                           <  2:429946‑M2/59‑1 (MQ=255)                                                     cactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAAc                           <  2:140639‑M2/59‑1 (MQ=255)                                                     cactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAAc                           <  1:700034‑M2/59‑1 (MQ=255)                                                             cacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAgg                   <  2:251603‑M2/67‑1 (MQ=255)                                                                    ggccggaagcgtcgggAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATcc            <  2:273083‑M2/74‑1 (MQ=255)                                                                    ggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATcc            >  2:140712‑M2/17‑90 (MQ=255)                                                                       cggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGc         <  2:191047‑M2/77‑1 (MQ=255)                                                                              gtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGa                  <  2:20745‑M2/68‑1 (MQ=255)                                                                              gtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGa                  >  1:20745‑M2/7‑74 (MQ=255)                                                                              gtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACaa  >  1:407854‑M2/7‑90 (MQ=255)                                                                                    |                                                                                    GCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAA  >  NZ_CP009273/2399389‑2399555

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