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	37 [0]	[0] 37	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)	37 (0.970)	31/238	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

TCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAGGGATCGGCTCGCCTTTTTGCAATATGTATTCCGCTGCGCAG  >  NZ_CP009273/2399339‑2399599                                                                                                                              |                                                                                                                                        tCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgcca                                                                                                                            <  2:131667‑M1/139‑14 (MQ=255)         cggcggTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatat                                                                                                                    >  2:705637‑M1/1‑118 (MQ=255)         cggcggTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatat                                                                                                                    <  1:395726‑M1/139‑22 (MQ=255)            cggTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctg                                                                                                                 <  2:457390‑M1/139‑25 (MQ=255)                          cTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtc                                                                                                   <  1:694116‑M1/139‑39 (MQ=255)                          cTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagc                                                                                                      <  1:373776‑M1/136‑36 (MQ=255)                            ccACGCAGCACACGCAGGCCCGGGCCCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTGTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgc                                                                                                 <  2:692691‑M1/139‑41 (MQ=255)                                       acGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatat                                                                                      <  1:693186‑M1/139‑52 (MQ=255)                                            ggTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcacaatatcgcg                                                                                  >  2:62864‑M1/1‑83 (MQ=255)                                             gTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaa                                                                                <  2:461521‑M1/139‑58 (MQ=255)                                              tCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaa                                                                               <  2:31206‑M1/139‑59 (MQ=255)                                              tCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaa                                                                               >  1:353566‑M1/1‑81 (MQ=255)                                                 gggCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgg                                                                            >  2:504802‑M1/1‑78 (MQ=255)                                                       atcatcTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggctag                                                                       >  1:158957‑M1/1‑72 (MQ=255)                                                       atcatcTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggctag                                                                       >  1:608203‑M1/1‑72 (MQ=255)                                                            ctcAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagt                                                                                                                                 <  2:207393‑M1/75‑9 (MQ=255)                                                            ctcAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagt                                                                                                                                 >  1:207393‑M1/1‑67 (MQ=255)                                                               aaCCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggctaggatcatcc                                                               <  1:348234‑M1/138‑75 (MQ=255)                                                                aCCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggctaggatcatccag                                                             <  2:757514‑M1/139‑77 (MQ=255)                                                                 ccGGCCTTGCCGTCACGCCAAGACCGGCTTTCTCTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggctaggatcatccag                                                             >  1:670763‑M1/1‑62 (MQ=255)                                                                 ccGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggctaggatcatccag                                                             <  2:670763‑M1/138‑77 (MQ=255)                                                                   ggCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggctaggatcatccagcaa                                                          <  2:642383‑M1/139‑80 (MQ=255)                                                                    gCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggctaggatcatccagcaa                                                          >  1:62156‑M1/1‑59 (MQ=255)                                                                    gCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaa                                                                               <  2:230120‑M1/117‑59 (MQ=255)                                                                    gCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaa                                                                               >  1:230120‑M1/1‑59 (MQ=255)                                                                               aCGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggctaggatcatccagcaatacaagagggat                                              >  2:759038‑M1/1‑48 (MQ=255)                                                                                      gACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggctaggatcatccagcaatacaagagggatcggctcg                                       >  1:335209‑M1/1‑41 (MQ=255)                                                                                                 aCTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggctaggatcatccagcaatacaagagggatcggctcgcctttttgca                             <  1:421229‑M1/138‑109 (MQ=255)                                                                                                  cTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggctaggatcatccagcaatacaagagggatcggctcgcctttttgcaat                           <  1:759038‑M1/139‑111 (MQ=255)                                                                                                  cTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggctaggatcatccagcaatacaagagggatcggctcgcctttttgcaat                           <  1:705637‑M1/139‑111 (MQ=255)                                                                                                             aCGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggctaggatcatccagcaatacaagagggatcggctcgcctttttgcaatatgtattccgc                >  2:580173‑M1/1‑18 (MQ=255)                                                                                                                 ccGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggctaggatcatccagcaatacaagagggatcggctcgcctttttgcaatatgtattccgctgcg            <  2:335209‑M1/139‑126 (MQ=255)                                                                                                                 ccGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggctaggatcatccagcaatacaagagggatcggctcgcctttttgcaatatgtattccgctgcg            <  2:608203‑M1/139‑126 (MQ=255)                                                                                                                       gCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggctaggatcatccagcaatacaagagggatcggctcgcctttttgcaatatgtattccgctgcgcagtac      <  2:158957‑M1/139‑132 (MQ=255)                                                                                                                         gTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggctaggatcatccagcaatacaagagggatcggctcgcctttttgcaatatgtattccgctgcgcagtac      >  1:31024‑M1/1‑6 (MQ=255)                                                                                                                          tCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggctaggatcatccagcaatacaagagggatcggctcgcctttttgcaatatgtattccgctgcgcagtaccag   >  1:52067‑M1/1‑5 (MQ=255)                                                                                                                           cGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggctaggatcatccagcaatacaagagggatcggctcgcctttttgcaatatgtattccgctgcgcagtaccagt  >  1:707813‑M1/1‑4 (MQ=255)                                                                                                                              |                                                                                                                                        TCAGGTAGCGGCGGTAAGCCATCTACTCCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAGGGATCGGCTCGCCTTTTTGCAATATGTATTCCGCTGCGCAG  >  NZ_CP009273/2399339‑2399599

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