Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,662,649	T→G	H58P (CAC→CCC)	mlc ←	DNA‑binding transcriptional repressor Mlc

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,662,649	0	T	G	100.0%	59.7 / NA	20	H58P (CAC→CCC)	mlc	DNA‑binding transcriptional repressor Mlc
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (8/12);  total (8/12)

GCTTCAGTTTCAACCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGTGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGACGGGAAAGATCGATACGC  >  NZ_CP009273/1662571‑1662726                                                                               |                                                                              gCTTCAGTTTCAACCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGc                                                                    >  1:131143/1‑90 (MQ=255)  cTTCAGTTTCAACCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCa                                                                   <  2:528031/90‑1 (MQ=255)    tCAGTTTCAACCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCAtc                                                                 >  1:561900/1‑90 (MQ=255)    tCAGTTTCAACCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCAtc                                                                 >  1:5247/1‑90 (MQ=255)          tCAACCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACgg                                                           <  2:425873/90‑1 (MQ=255)          gcAACCACCAGCCCCACCGCCGGCCGGCCACGTTTCCCCGCTTCTGTGGTTTCTAGCTCTTGCACCTGGGGTGGTTATAGCATCCCACgg                                                           >  1:290591/2‑90 (MQ=255)            aaCCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGAc                                                         <  2:190073/90‑1 (MQ=255)                      cccACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTg                                               <  1:442043/90‑1 (MQ=255)                       ccACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGa                                              >  1:208723/1‑90 (MQ=255)                         aCCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATa                                            >  2:59621/1‑90 (MQ=255)                           cgccgGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACt                                          <  1:118846/90‑1 (MQ=255)                               ggACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACt                                          <  2:392446/86‑1 (MQ=255)                               ggACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACt                                          >  1:392446/1‑86 (MQ=255)                                        cGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGt                             <  2:201608/90‑1 (MQ=255)                                             ccccGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTgcgc                        <  1:270455/90‑1 (MQ=255)                                                  cttcttTGATTTCCAGCTCTTGCACCGGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGAc                   <  1:426773/90‑1 (MQ=255)                                                      tttGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGACGGGa               >  2:90352/1‑90 (MQ=255)                                                          aTTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGACGGGAAAGa           <  1:422598/90‑1 (MQ=255)                                                                   tctTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGACGGGAAAGATCGATAcgc  <  1:572320/90‑1 (MQ=255)                                                                   tctTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGACGGGAAAGATCGATAcgc  <  1:111599/90‑1 (MQ=255)                                                                               |                                                                              GCTTCAGTTTCAACCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGTGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGACGGGAAAGATCGATACGC  >  NZ_CP009273/1662571‑1662726

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