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

GCAATGCCGCCTGACACAATGAGTCCACGGTTAACGGTTGTCCATGTAACATCGAGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCGATGTACAATACCATTTTCCGTATCAATAATTCCCGGCAAGGTTA  >  NC_000913/1667649‑1668113                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                   gcAATGCCGCCTGACACAATGAGTCCACGGTTAACGGTTGTCCATGTAACATCGAGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTtgt                                                                                                                                                                                                                       <  1:167439/251‑1 (MQ=255)      gccgccTGACACAATGAGTCCACGGTTAACGGTTGTCCATGTAACATCGAGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATc                                                                                                                                                                                                                   <  1:309294/250‑1 (MQ=255)                         ccACGGTTAACGGTTGTCCATGTAACATCGAGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATcaca                                                                                                                                                                                                <  1:343723/250‑1 (MQ=255)                            cGGTTAGCGGTTGTCCATGTAACATCGCGCTCCTGGCTTGTTTAAGACGCAGCTGTGCCAGCTCAAGAATCCTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATACCCGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTACCCGATGACAATCGGTAATGAAGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATcgc                                                                                                                                                                                            >  2:34735/1‑251 (MQ=255)                                tAACGGTTGTCCATGTAACATCGAGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCgcgcgc                                                                                                                                                                                        >  2:306678/1‑251 (MQ=255)                                          ccATGTAACATCGAGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCTCCTGAATCACATCGCGCGCCCCGCGTGAg                                                                                                                                                                              <  2:185356/251‑1 (MQ=255)                                                     cGAGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGTCCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACaa                                                                                                                                                                   <  2:185627/251‑1 (MQ=255)                                                                            gCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATgcgc                                                                                                                                            >  1:332873/1‑251 (MQ=255)                                                                                    gCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATc                                                                                                                                     <  1:109404/250‑1 (MQ=255)                                                                                         cTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTg                                                                                                                               >  1:349008/1‑251 (MQ=255)                                                                                                    tGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCg                                                                                                                    >  1:60351/1‑251 (MQ=255)                                                                                                      tCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTGCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGAGCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCgg                                                                                                                   <  1:64480/250‑1 (MQ=255)                                                                                                                                 ccGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCg                                                                                       >  2:143479/1‑251 (MQ=255)                                                                                                                                     gATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCg                                                                                    <  2:9092/250‑1 (MQ=255)                                                                                                                                              aaTAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCAtct                                                                          >  1:354579/1‑251 (MQ=255)                                                                                                                                                           cccATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCTCGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCGCCAGCTCCTGGCCTAGCGGCAGCGCTTTTACATCCt                                                              <  2:227151/250‑1 (MQ=255)                                                                                                                                                           cccATACGGGTCGACCTGTGTGTGACCTATTTCCACGAGGCTACTGCTGCCTGCGTGTAGCAGATGACCATGGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCt                                                              <  1:346965/250‑1 (MQ=255)                                                                                                                                                                         ccTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGc                                               <  1:1056/251‑1 (MQ=255)                                                                                                                                                                          cTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGTACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCg                                              >  2:263307/1‑251 (MQ=255)                                                                                                                                                                            gtgtgtgGCCTATTTCCACGAGGCTACTGCTGCCTGCGTGGAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGGGAGGCACCAAACAAGGGCTCTGCCATCGTCCATGCTCTGATATCATGCTGAAGATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCGTCTCTTGTACATCCTCGTAGAACGGCATGCGa                                             <  2:320418/250‑1 (MQ=255)                                                                                                                                                                            gtgtgtgGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCGAt                                            >  1:241447/1‑251 (MQ=255)                                                                                                                                                                                    ccTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCGATGTACAAt                                     >  1:131759/1‑250 (MQ=255)                                                                                                                                                                                     cTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCGATGTACAATAc                                   >  1:165365/1‑251 (MQ=255)                                                                                                                                                                                                     actGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCGATGTACAATACCATTTGCCGTATCaat                   <  1:263307/251‑1 (MQ=255)                                                                                                                                                                                                           gcctgcGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTATAACGGCATGCGATGTACAATACCATTTTCCGTATCAATAATTc              <  1:161894/250‑1 (MQ=255)                                                                                                                                                                                                                cGTGTAGCAGATGACCATCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCGATGTACAATACCATTTTCCGTATCAATAATTCCCGGCa        <  2:165365/251‑1 (MQ=255)                                                                                                                                                                                                                      gCAGATGACCAGCGGTAATGACGCC‑GCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCGATGTACAATACCATTTTCCGTATCAATAATTCCCGGCAAGGTTa  <  2:333965/251‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                   GCAATGCCGCCTGACACAATGAGTCCACGGTTAACGGTTGTCCATGTAACATCGAGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCGATGTACAATACCATTTTCCGTATCAATAATTCCCGGCAAGGTTA  >  NC_000913/1667649‑1668113

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