Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,667,887	Δ1 bp	coding (678/1221 nt)	mlc ←	DNA‑binding transcriptional repressor Mlc

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,667,887	0	C	.	100.0%	139.5 / NA	31	coding (678/1221 nt)	mlc	DNA‑binding transcriptional repressor Mlc
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (13/18);  total (13/18)

AATGCCGCCTGACACAATGAGTCCACGGTTAACGGTTGTCCATGTAACATCGAGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCGATGTACAATACCATTTTCCGTATCAATAATTCCCGGCAAG  >  NC_000913/1667651‑1668109                                                                                                                                                                                                                                             |                                                                                                                                                                                                                               aaTGCCGCCTGACACAATGAGTCCACGGTTAACGGTTGTCCATGTAACATCGAGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTtgtg                                                                                                                                                                                                                  >  1:72825/1‑250 (MQ=255)            acacAATGAGTCCACGGTTAACGGTTGTCCATGTAACATCGAGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTAGTTCCACGAGACTACTGCTGCCTGCGGGTAGCAGATGACCATCGGTAATGACGCC‑GCCGCCACGTTGTGATCGATAaccac                                                                                                                                                                                                      <  2:166542/251‑1 (MQ=255)                           ggTTAACGGTTGTCCATGTACCATCGAGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATcgcg                                                                                                                                                                                       >  1:120835/1‑251 (MQ=255)                               aaCGGTTGTCCATGTAACATCGAGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGcc                                                                                                                                                                                   >  1:331902/1‑251 (MQ=255)                                  ggTTGTCCTTGTAACATGGAGCTCATGGCTTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGGCATTCCCGCAATAACCGCGTTTCCCATACGGGTCGCCCTGTGTGTGGCCTATTTCCACGCGACTACTGCTGCCTGCGTGTAGCAGCTGACCATCGGTATTGACGCC‑GCCCCCACGTTGTGATCGAAAACCACCTCAATCACATCGCGCGccacg                                                                                                                                                                                >  2:165987/1‑248 (MQ=255)                                        ccATGTAACATCGAGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGa                                                                                                                                                                           <  1:158182/250‑1 (MQ=255)                                                            ggATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTg                                                                                                                                                       <  2:322776/250‑1 (MQ=255)                                                            ggATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGc                                                                                                                                                      <  2:115658/251‑1 (MQ=255)                                                                 gattAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATc                                                                                                                                                  <  1:169885/250‑1 (MQ=255)                                                                          gCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGCGGCACCAAACAAGGCCTCTGCCATCGTCCATgcgc                                                                                                                                        >  2:372431/1‑251 (MQ=255)                                                                              atgtgACAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTAGAGGCAGCCGTGATTCCCGCAATAACAGCGGTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCCCGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCCAATCGCGCGCCCCGCGGGCGGCACCAAACAAGGCCTCTTCCATCGTCCCTGCGCTGat                                                                                                                                    >  2:299419/2‑251 (MQ=255)                                                                                              aTACTGACCAAGCTGGCGATGGTTTCGAGGCAGCAGTGATTCCCGCAATAACAGCGTTTCCCATTCGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCGGCTTGTAGCAGACGACCATCAGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCGCCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATTCGCTGATATCATGCTGAATATaa                                                                                                                    <  2:227449/251‑1 (MQ=255)                                                                                                aCTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAAc                                                                                                                  <  2:50999/251‑1 (MQ=255)                                                                                                   gTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCGATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACAGCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTGTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCgg                                                                                                               <  1:121291/251‑1 (MQ=255)                                                                                                     ccACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGa                                                                                                              <  2:40118/250‑1 (MQ=255)                                                                                                                  ggTTTCGAGGCAGCCGTGCTTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGCGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGTTATAATGCTGAATATAAACCGTAACGCCGGTAtgctg                                                                                                >  2:326282/1‑251 (MQ=255)                                                                                                                            cAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACCGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGcc                                                                                      >  1:376302/1‑251 (MQ=255)                                                                                                                              gCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGATACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTc                                                                                    <  2:325296/251‑1 (MQ=255)                                                                                                                               ccGTGATTCCCGCAATAACAGAGTTTCCCATACGGGTCGACATGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGCTGTGATCGATAACAACCTGAATAACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCTTCTGCCATCGTCCATGCGATGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGAGCCTcc                                                                                   >  2:57927/1‑250 (MQ=255)                                                                                                                                                     gTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACAt                                                             <  2:281224/251‑1 (MQ=255)                                                                                                                                                       ttCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATcc                                                           >  1:203720/1‑251 (MQ=255)                                                                                                                                                              aCGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAg                                                     <  1:187501/250‑1 (MQ=255)                                                                                                                                                              aCGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCGCTTTTACATCCTCGTAGa                                                    <  1:322760/251‑1 (MQ=255)                                                                                                                                                                     gACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGCCCACCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCAt                                             <  2:72825/251‑1 (MQ=255)                                                                                                                                                                          gtgtgtgGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCGa                                         <  1:372431/250‑1 (MQ=255)                                                                                                                                                                                 gCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCGATGTACaa                                  >  2:117483/1‑250 (MQ=255)                                                                                                                                                                                          cACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCATCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCGATGTACAATACCATTTTc                        >  1:271111/1‑251 (MQ=255)                                                                                                                                                                                                   actGCTGCTTGCGTGTAGCAGATTACCATCGGTAATGACGCC‑GCCCCCCCGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGTAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTTGTAGAACGGCATGCGATGTACAATACCATTTTCCGTATCa                 <  2:116446/249‑1 (MQ=255)                                                                                                                                                                                                         gcctgcGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCGATGTACAATACCATTTTCCGTATCAATAAtt           <  2:372234/249‑1 (MQ=255)                                                                                                                                                                                                               gtgtAGAAGATGACCATCGGTAATGTCGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGACACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTCTTACATCCTCGTAGAACGGCATGCGATGAACAATACCTTTTTCCGTATCAATAATTccggtca    >  2:318742/1‑245 (MQ=255)                                                                                                                                                                                                                 gtAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCGATGTACAATACCATTTTCCGTATCAATAATTCCCGGCAAg  <  2:120835/250‑1 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                               AATGCCGCCTGACACAATGAGTCCACGGTTAACGGTTGTCCATGTAACATCGAGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCGATGTACAATACCATTTTCCGTATCAATAATTCCCGGCAAG  >  NC_000913/1667651‑1668109

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