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

GGTTGTCCATGTAACATCGAGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCGATGTACAATACCATTTTCCGTATCAATAATTCCCGGCAAGGTTATGGCAA  >  NC_000913/1667684‑1668119                                                                                                                                                                                                            |                                                                                                                                                                                                                                         ggTTGTCCATGTAACATCGAGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCcg                                                                                                                                                                                          >  1:67262/1‑251 (MQ=255)     gtccatgttacctccagcaaaTGGATTGATTAAGACGCAGCTGTGCCAGCTAAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTACCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCGGCGTGGAGAAGATGACCATCGGTAATGACGCC‑ACCCCCACGTTGTGATCGATAACACCCAGAATCACACCACGCACCCcgcg                                                                                                                                                                                        >  2:132427/21‑249 (MQ=255)                                  ttAAGACGCAGCTGTGCCAGCTCAAGACTACTGTCCACGCTGGCGATGGTTTAGCGCCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGCGGCCTATTTCCACGAGACTACTGCTGCCTCCGTGTAGCAGATGACCATCGGTAATGACGAC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCACACAAGGCCTCTGCCATTGTc                                                                                                                                                         >  1:134661/1‑251 (MQ=255)                                    aaGACCCAGCTGCGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCAACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTcc                                                                                                                                                        >  2:218902/1‑250 (MQ=255)                                     aGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACAGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCCCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCAt                                                                                                                                                      >  1:184191/1‑251 (MQ=255)                                           aGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTg                                                                                                                                                >  1:215113/1‑251 (MQ=255)                                                        cAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGa                                                                                                                                    >  1:112508/1‑250 (MQ=255)                                                               tcTGTCCACGCTGGGGATTGTTTCGAGGCAGCCCTGATTCCCGCAATAGCAGCGTGTCCCATACGGGTCGACCTGTGTGTCGCCTATTTCCACGAGCCTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑TCCCCCTCGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCTCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATaaa                                                                                                                             <  1:98740/249‑1 (MQ=255)                                                                                   tttCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCACATACGGGTAGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGGAATGACGCC‑GCCCCCACGTTGTGATCGATAAACACCTGAATCACATCGCGAGCCCCGCGAGAGGCACCAAACAAGGCATCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGAAtgatggt                                                                                                        >  1:146759/1‑246 (MQ=255)                                                                                            aGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCt                                                                                               >  2:81889/1‑251 (MQ=255)                                                                                             gCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTc                                                                                              <  1:81889/251‑1 (MQ=255)                                                                                              ccGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCAATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCg                                                                                             <  1:177008/251‑1 (MQ=255)                                                                                                  gATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCg                                                                                          >  2:81436/1‑250 (MQ=255)                                                                                                            aTAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCAtctc                                                                               <  2:33288/251‑1 (MQ=255)                                                                                                                        cccATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTc                                                                   <  1:1343/251‑1 (MQ=255)                                                                                                                             aCGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAg                                                               <  2:215113/250‑1 (MQ=255)                                                                                                                               gggTCGACCGGTGTGTGGCCTATTTCCACGCGACTCCTGCTGTCTGCGTGTAGCAGATCACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCAGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAAc                                                            <  2:184191/251‑1 (MQ=255)                                                                                                                                    gACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCAt                                                       >  1:191933/1‑251 (MQ=255)                                                                                                                                         gtgtgtgGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCGAt                                                  <  2:7074/251‑1 (MQ=255)                                                                                                                                             gtgGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTGACATCCTCGTAGAACGGCATGCGATGTa                                               <  1:166273/250‑1 (MQ=255)                                                                                                                                                            gagaCTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCGGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCGATGTACAATACCATTTTCCGt                               >  1:207249/1‑251 (MQ=255)                                                                                                                                                                                      gATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCGATGTACAATACCATTTTCAGTATCAATAATTCCCGGCAAGGTTATgg     >  2:177526/1‑251 (MQ=255)                                                                                                                                                                                         gACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCGATGTACAATACCATTTTCCGTATCAATAATTCCCGGCAAGGTTATGGCaa  <  1:15095/251‑1 (MQ=255)                                                                                                                                                                                                            |                                                                                                                                                                                                                                         GGTTGTCCATGTAACATCGAGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCGATGTACAATACCATTTTCCGTATCAATAATTCCCGGCAAGGTTATGGCAA  >  NC_000913/1667684‑1668119

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