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	36 [0]	[0] 36	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)	36 (1.130)	26/240	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

TAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAGGGATCGGCTCGCCTTTTTGCAATATGTATTCCGCTGCGCA  >  NZ_CP009273/2399344‑2399598                                                                                                                                 |                                                                                                                               tatcaggtagcggcggtaagccatctactccactcagcacacgcaggtccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGc                                                                                                                      <  2:434549‑M2/11‑1 (MQ=255)           cggcggtaagccatctactccactcagcacacgcaggtccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAatat                                                                                                            <  2:413854‑M2/21‑1 (MQ=255)              cggtaagccatctactccactcagcacacgcaggtccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTg                                                                                                         <  2:413264‑M2/24‑1 (MQ=255)               ggtaagccatctactccactcagcacacgcaggtccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGc                                                                                                        <  1:500472‑M2/25‑1 (MQ=255)                 taagccatctactccactcagcacacgcaggtccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCtt                                                                                                      >  1:582634‑M2/113‑139 (MQ=255)                         ctactccactcagcacacgcaggtccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGc                                                                                              >  2:192883‑M2/105‑139 (MQ=255)                            ctccactcagcacacgcaggtccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGc                                                                                              <  1:251813‑M2/35‑1 (MQ=255)                               cactcagcacacgcaggtccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGc                                                                                         <  2:170946‑M2/40‑1 (MQ=255)                               cactcagcacacgcaggtccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGc                                                                                         >  1:170946‑M2/99‑138 (MQ=255)                               cactcagcacacgcaggtccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGc                                                                                         <  2:450592‑M2/40‑1 (MQ=255)                                              ggtccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATcgcg                                                                          >  1:591192‑M2/84‑138 (MQ=255)                                               gtccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGaa                                                                        <  1:457477‑M2/57‑1 (MQ=255)                                                tccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGaaa                                                                       <  1:553462‑M2/58‑1 (MQ=255)                                                tccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGaaa                                                                       <  2:493091‑M2/58‑1 (MQ=255)                                                tccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGaaa                                                                       <  2:582634‑M2/58‑1 (MQ=255)                                                tccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGaaa                                                                       <  1:617710‑M2/58‑1 (MQ=255)                                                tccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGaaa                                                                       <  2:591192‑M2/58‑1 (MQ=255)                                                tccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGaaa                                                                       <  1:140815‑M2/58‑1 (MQ=255)                                                tccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATAGCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGaaa                                                                       <  2:118451‑M2/58‑1 (MQ=255)                                                   gggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTACGCGTCGCCAGCACCATATCGCGAAACgg                                                                    >  2:274703‑M2/79‑139 (MQ=255)                                                         atcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGaaa                                                                       <  2:189221‑M2/58‑1 (MQ=255)                                                         atcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGaaa                                                                       >  1:189221‑M2/73‑130 (MQ=255)                                                           catctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATcg                                                                            >  1:57007‑M2/71‑123 (MQ=255)                                                           catctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATcg                                                                            <  2:57007‑M2/53‑1 (MQ=255)                                                                  accggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCTAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAg                                                     <  1:211593‑M2/76‑1 (MQ=255)                                                                      gccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCaa                                                  >  2:617324‑M2/60‑138 (MQ=255)                                                                                             gctttcactgccgcatgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACggg                                                                   >  2:274409‑M2/37‑98 (MQ=255)                                                                                             gctttcactgccgcatgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACggg                                                                   <  1:274409‑M2/62‑1 (MQ=255)                                                                                             gctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAGGGATCGGCTCGCCttt                          >  2:108556‑M2/37‑139 (MQ=255)                                                                                                  cactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAGGGATCGGCTCGCCTTTTTGc                      <  1:108556‑M2/107‑1 (MQ=255)                                                                                                       ccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGaaa                                                                       <  1:497349‑M2/58‑1 (MQ=255)                                                                                                       ccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGaaa                                                                       >  2:497349‑M2/27‑84 (MQ=255)                                                                                                        cgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAGGGATCGGCTCGCCTTTTTGCAATATg                >  2:168658‑M2/26‑138 (MQ=255)                                                                                                              aacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAGGGATCGGCTCGCCTTTTTGCAATATGTATTCCg         <  1:283057‑M2/120‑1 (MQ=255)                                                                                                                   ccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAGGGATCGGCTCGCCTTTTTGCAATATGTATTCCGCTgcg    <  2:142460‑M2/125‑1 (MQ=255)                                                                                                                     ggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTATGATCATCCAGCAATACAAGAGGGATCGGCTCGCCTTTTTGCAATATGTATTCCGCTGCGCa  <  1:212988‑M2/127‑1 (MQ=255)                                                                                                                                 |                                                                                                                               TAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAGGGATCGGCTCGCCTTTTTGCAATATGTATTCCGCTGCGCA  >  NZ_CP009273/2399344‑2399598

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