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%	137.0 / NA	31	coding (678/1221 nt)	mlc	DNA‑binding transcriptional repressor Mlc
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (15/16);  total (15/16)

CGCAATGCCGCCTGACACAATGAGTCCACGGTTAACGGTTGTCCATGTAACATCGAGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCGATGTACAATACCATTTTCCGTATCAATAATTCCCGGCAAGGTTATGGCAATCGAAGTTAGACGCTCAA  >  NC_000913/1667648‑1668137                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                           cgcAATGCCGCCTGACACAATGAGTCCACGGTTAACGGTTGTCCATGTAACATCGAGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTtg                                                                                                                                                                                                                                                <  1:344608/251‑1 (MQ=255)                       ggTCCCCGGTTAACGGTTGTCCATGTAACATCGAGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATca                                                                                                                                                                                                                          <  2:342168/250‑1 (MQ=255)                                                        aGCTCATGGCTTGCTTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAgg                                                                                                                                                                                         >  2:270531/1‑251 (MQ=255)                                                        aGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCCAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGGGAGGCACCAAACAAgg                                                                                                                                                                                         >  2:69259/1‑251 (MQ=255)                                                          cTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGcc                                                                                                                                                                                       >  1:60124/1‑251 (MQ=255)                                                                     attAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACACGGCCTCTGCCATCGt                                                                                                                                                                            >  1:349080/1‑251 (MQ=255)                                                                                             aaGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAat                                                                                                                                                    >  1:93941/1‑251 (MQ=255)                                                                                                   aCTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATaaa                                                                                                                                               <  2:186803/250‑1 (MQ=255)                                                                                                             cTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGGGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCCCCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACgacg                                                                                                                                    >  2:362612/1‑248 (MQ=255)                                                                                                              tgacGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCgg                                                                                                                                   >  1:207273/4‑251 (MQ=255)                                                                                                              tGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCgg                                                                                                                                   >  1:207255/1‑251 (MQ=255)                                                                                                                                 gCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTc                                                                                                                <  2:69951/251‑1 (MQ=255)                                                                                                                                 gCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTc                                                                                                                <  2:50682/251‑1 (MQ=255)                                                                                                                                      gactCCCGCAATAACAGCGTTTCCCATACGGGTCTACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCg                                                                                                            >  2:262041/4‑250 (MQ=255)                                                                                                                                                    cAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCtttt                                                                                             >  2:241248/1‑251 (MQ=255)                                                                                                                                                              cATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGt                                                                                   <  2:246491/251‑1 (MQ=255)                                                                                                                                                              cATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGt                                                                                   <  1:241248/251‑1 (MQ=255)                                                                                                                                                                   gggTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCCCGTTGTGATCGATAACCACCTGAATCACATCGCGCGTCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAAc                                                                              <  2:207273/251‑1 (MQ=255)                                                                                                                                                                   gggTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAAc                                                                              <  2:207255/251‑1 (MQ=255)                                                                                                                                                                    ggTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACg                                                                             <  1:84235/251‑1 (MQ=255)                                                                                                                                                                             gtgtgtgGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCGAt                                                                    <  1:209893/251‑1 (MQ=255)                                                                                                                                                                              tgtgtgGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAg                                                                                                                                                                                                      >  1:407150/1‑120 (MQ=255)                                                                                                                                                                              tgtgtgGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAg                                                                                                                                                                                                      <  2:407150/120‑1 (MQ=255)                                                                                                                                                                                tgtgGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCTCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCGGCTCCAGCGCCTCGCCGAGCGGCATCTCtt                                                                                               <  2:191736/221‑1 (MQ=255)                                                                                                                                                                                tgtgGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCtt                                                                                               >  1:191736/1‑221 (MQ=255)                                                                                                                                                                                        aTTTCCCCGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCGCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGTCATGCGATGTACAATAc                                                           <  1:362612/249‑1 (MQ=255)                                                                                                                                                                                                           tgcctgcGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCGATGTACAATACCATTTTCCGTATCAATAATTc                                      >  2:387476/1‑251 (MQ=255)                                                                                                                                                                                                                         aGATGACCATCGGTAATGACGCC‑GCACCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGACGGTATGCTGCTCCAGCGCCTCGCCGCGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCTATGTAAAATACCATTTTCCGTATCAATAATTCCCGGCAAGGTTATg                        >  2:382070/1‑251 (MQ=255)                                                                                                                                                                                                                             gACCATCGGTAACGACGCC‑GCCCCCACGTTGTCATAGATAACCACCTGAATCACATCGCGCTCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCTATGTACAATACCATTTTCCGTATCAATAATTCCCGGCAAGGTGATGGCaa                    <  2:354695/251‑1 (MQ=255)                                                                                                                                                                                                                               ccATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCGATGTACAATACCATTTTCCGTATCAATAATTCCCGGCAAGGTTATGGCAATc                  <  1:77599/251‑1 (MQ=255)                                                                                                                                                                                                                                        aTGACGCC‑GCCCCCACGTTGCAATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCGATGTACAATACCATTTTCCGTATCAATAATTCCCGGCAAGGTTATGGCAATCGAAGtt            >  1:292928/1‑248 (MQ=255)                                                                                                                                                                                                                                                cgcccCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAATGGCATGCGATGTACAATACCATTTTCCGGATCAATAATTCCCGGCAAGGTTATGGCAATCGAAGTTAGACGCTCaa  <  2:206641/251‑1 (MQ=255)                                                                                                                                                                                                                                                cgcccCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCGATGTACAATACCATGTTCCGTATCAATAATTCCCGGCAAGGTTATGGCAATCGAAGTTAGACGCTCaa  <  1:270531/251‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                           CGCAATGCCGCCTGACACAATGAGTCCACGGTTAACGGTTGTCCATGTAACATCGAGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCGATGTACAATACCATTTTCCGTATCAATAATTCCCGGCAAGGTTATGGCAATCGAAGTTAGACGCTCAA  >  NC_000913/1667648‑1668137

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