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.110)	27/244	0.4	100%	coding (657/939 nt)	lrhA	transcriptional regulator LrhA
	?	NZ_CP009273	2399472 =	0 (0.000)					coding (649/939 nt)	lrhA	transcriptional regulator LrhA

GGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAGGGATCGGCTCGCCTTTTTGCAATATGTATTCCGCTGCGC  >  NZ_CP009273/2399342‑2399597                                                                                                                                   |                                                                                                                              ggtatcaggtagcggcggtaagccatctactccactcagcacacgcaggtccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGt                                                                                                                        <  1:471145‑M2/8‑1 (MQ=255) ggtatcaggtagcggccgcaagccatctactccactcagcacacgcaggtccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGt                                                                                                                        <  2:742927‑M2/8‑1 (MQ=255)             cggcggtaagccatctactccactcagcacacgcagttccgggctcatcatctcaaccggccttgccgtcacgccaataccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAatat                                                                                                           <  1:859604‑M2/21‑1 (MQ=255)             cggcggtaagccatctactccactcagcacacgcaggtccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAatat                                                                                                           >  1:266268‑M2/119‑139 (MQ=255)             cggcggtaagccatctactccactcagcacacgcaggtccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAatat                                                                                                           <  1:542935‑M2/21‑1 (MQ=255)                 ggtaagccatctactccactcagcacacgcaggtccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGc                                                                                                       <  2:20727‑M2/25‑1 (MQ=255)                             actccactcagcacacgcaggtccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGt                                                                                           >  1:170424‑M2/103‑139 (MQ=255)                                 cactcagcacacgcaggtccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGc                                                                                        <  1:158632‑M2/40‑1 (MQ=255)                                    tcagcacgcgcaggtccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAg                                                                                     >  2:350340‑M2/96‑138 (MQ=255)                                    tcagcacacgcaggtccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAg                                                                                     >  1:853389‑M2/96‑138 (MQ=255)                                      agcacacgcaggtccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGt                                                                                                                        <  1:106245‑M2/8‑1 (MQ=255)                                      agcacacgcaggtccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGt                                                                                                                        >  2:106245‑M2/94‑101 (MQ=255)                                      agcacacgcaggtccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCAc                                                                                  >  1:111936‑M2/94‑139 (MQ=255)                                      agcacacgcaggtccgggctcatcatctcaaccggccttgccgtcacgccaagaccggcttgcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCa                                                                                   <  1:234973‑M2/45‑1 (MQ=255)                                             gcaggtccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAatat                                                                                                           <  2:280974‑M2/21‑1 (MQ=255)                                             gcaggtccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAatat                                                                                                           >  1:280974‑M2/87‑107 (MQ=255)                                                ggtccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATcgcg                                                                         >  2:751580‑M2/84‑138 (MQ=255)                                                  tccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacggacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGaaa                                                                      <  1:286619‑M2/58‑1 (MQ=255)                                                  tccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGaaa                                                                      <  2:341351‑M2/58‑1 (MQ=255)                                                  tccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGaaa                                                                      >  2:466451‑M2/82‑139 (MQ=255)                                                  tccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGaaa                                                                      <  2:809843‑M2/58‑1 (MQ=255)                                                  tccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGaaa                                                                      <  1:711231‑M2/58‑1 (MQ=255)                                                  tccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGaaa                                                                      <  1:36685‑M2/58‑1 (MQ=255)                                                  tccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGaaa                                                                      <  1:559264‑M2/58‑1 (MQ=255)                                                  tccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGaaa                                                                      <  1:583876‑M2/58‑1 (MQ=255)                                                           atcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGc                                                                                             <  2:451029‑M2/35‑1 (MQ=255)                                                           atcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGc                                                                                             >  1:451029‑M2/73‑107 (MQ=255)                                                           atcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAgg                                                             >  2:434808‑M2/73‑139 (MQ=255)                                                           atcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAg                                                              >  1:35540‑M2/73‑138 (MQ=255)                                                            tcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGa                                                            >  2:717135‑M2/72‑139 (MQ=255)                                                            tcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGa                                                            >  1:209009‑M2/72‑139 (MQ=255)                                                                   aaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATcc                                                      <  2:266268‑M2/74‑1 (MQ=255)                                                                    accggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAg                                                    <  2:811277‑M2/76‑1 (MQ=255)                                                                    accggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAg                                                    <  2:512777‑M2/76‑1 (MQ=255)                                                                    accggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAg                                                    <  1:420660‑M2/76‑1 (MQ=255)                                                                    accggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAg                                                    <  2:853389‑M2/76‑1 (MQ=255)                                                                    accggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAg                                                    <  1:675698‑M2/76‑1 (MQ=255)                                                                     ccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGc                                                                                             >  2:7029‑M2/63‑97 (MQ=255)                                                                     ccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGc                                                                                             <  1:7029‑M2/35‑1 (MQ=255)                                                                             gccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATAc                                              >  1:671128‑M2/55‑136 (MQ=255)                                                                                        aagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAGGGATCGGCt                                <  1:434808‑M2/96‑1 (MQ=255)                                                                                               gctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAatat                                                                                                           >  2:176235‑M2/37‑57 (MQ=255)                                                                                               gctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAatat                                                                                                           <  1:176235‑M2/21‑1 (MQ=255)                                                                                               gctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAGGGATCGGCTCGCCttt                         >  2:90296‑M2/37‑139 (MQ=255)                                                                                                    cactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAGGGATCGGCTCGCCTTTTTGc                     <  1:751580‑M2/107‑1 (MQ=255)                                                                                                      ctgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAGGGATCGGCTCGCCTTTTTGCAat                  <  2:382517‑M2/110‑1 (MQ=255)                                                                                                         ccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAGGGATCGGCTCGCCTTTTTGCAATATg               >  2:305558‑M2/27‑139 (MQ=255)                                                                                                                     ccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAGGGATCGGCTCGCCTTTTTGCAATATGTATTCCGCTgcg   <  1:466451‑M2/125‑1 (MQ=255)                                                                                                                     ccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAGGGATCGGCTCGCCTTTTTGCAATATGTATTCCGCTgcg   <  1:639932‑M2/125‑1 (MQ=255)                                                                                                                     ccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAGGGATCGGCTCGCCTTTTTGCAATATGTATTCCGCTgcg   <  1:717135‑M2/125‑1 (MQ=255)                                                                                                                      cggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAGGGATCGGCTCGCCTTTTTGCAATATGTATTCCGCTgcgc  <  1:646494‑M2/126‑1 (MQ=255)                                                                                                                                   |                                                                                                                              GGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAGGGATCGGCTCGCCTTTTTGCAATATGTATTCCGCTGCGC  >  NZ_CP009273/2399342‑2399597

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: ‑