Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,668,069	G→A	R166C (CGC→TGC)	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,668,069	0	G	A	100.0%	143.6 / NA	43	R166C (CGC→TGC)	mlc	DNA‑binding transcriptional repressor Mlc
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/24);  total (19/24)

CCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCGATGTACAATACCATTTTCCGTATCAATAATTCCCGGCAAGGTTATGGCAATCGAAGTTAGACGCTCAAGTTTTTTCTGGTGGCGGATAAAAAACTGATCGATATGGGAAATAATACGATCCAGCAATGGCA  >  NC_000913/1667933‑1668200                                                                                                                                         |                                                                                                                                    ccGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCAATg                                                                                                                                  >  1:31190/1‑140 (MQ=255) ccGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCAATg                                                                                                                                  <  2:848906/140‑1 (MQ=255)   gcgTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGGATCTCTTTTACATCCTCGTAGAACGGCATGCAATGTa                                                                                                                                <  2:877392/140‑1 (MQ=255)     gTGAGGCACCAAACAAGGCATCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCAATGTACa                                                                                                                              <  2:172988/140‑1 (MQ=255)       gAGGAACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCAATGTACAAt                                                                                                                            <  2:881234/140‑1 (MQ=255)           cACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCAATGTACAATACCa                                                                                                                        <  1:254131/140‑1 (MQ=255)             ccAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCAATGTACAATACCAtt                                                                                                                      >  1:106623/1‑140 (MQ=255)              cAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGAAGATCGACATGCAGTGTACAATACCAttt                                                                                                                     >  1:479101/1‑140 (MQ=255)                   aaGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCAATGTACAATACCATTTTCCGt                                                                                                                <  1:1009012/140‑1 (MQ=255)                   aaGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCAATGTACAATACCATTTTCCGt                                                                                                                <  2:669895/140‑1 (MQ=255)                                 cGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCAATGTACAATACCATTTTCCGTATCAATAATTCCCg                                                                                                  >  1:977331/1‑140 (MQ=255)                                  gTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGGAATGTACAATACCAATTTCCGTATCAATAATTCCCgg                                                                                                 >  1:439431/1‑140 (MQ=255)                                   tCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCAATGTACAATACCATTTTCCGTATCAATAATTCCCGGc                                                                                                <  2:772852/140‑1 (MQ=255)                                   tCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCAATGTACAATACCATTTTCCGTATCAATAATTCCCGGc                                                                                                <  2:767896/140‑1 (MQ=255)                                   tCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCAATGTACAATACCATTTTCCGTATCAATAATTCCCGGc                                                                                                <  1:204710/140‑1 (MQ=255)                                   tCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCAATGTACAATACCATTTTCCGTATCAATAATTCCCGGc                                                                                                <  1:171143/140‑1 (MQ=255)                                       tGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCAATGTACAATACCATTTTCCGTATCAATAATTCCCGGCAAgg                                                                                            >  2:616237/1‑140 (MQ=255)                                         cgcTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCAATGTACAATACCATTTTCCGTATCAATAATTCCCGGCAAGGtt                                                                                          <  1:616237/140‑1 (MQ=255)                                              atatCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCAATGTACAATACCATTTTCCGTATCAATAATTCCCGGCAAGGTTATGGc                                                                                     >  1:1009633/1‑140 (MQ=255)                                                   aTGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCAATGTACAATACCATTTTCCGTATCAATAATTCCCGGCAAGGTTATGGCAATCg                                                                                >  1:244114/1‑140 (MQ=255)                                                      cTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCAATGTACAATACCATTTTCCGTATCAATAATTCCCGGCAAGGTTATGGCAATCGAAg                                                                             >  1:696307/1‑140 (MQ=255)                                                                   ggAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCAATGTACAATACCATTTTCCGTATCAATAATTCCCGGCAAGGTTATGGCAATCGAAGTTAGACGCTCAAg                                                                >  1:674824/1‑140 (MQ=255)                                                                     aaCGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCAATGTTCAATACCATTTTCCGTATCAATAATTCCCGGCAAGGTTATGGCAATCGAAGTTAGACGCTCAAGtt                                                              >  1:981061/1‑140 (MQ=255)                                                                     aaCGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCAATGTACAATACCATTTTCCGTATCAATAATTCCCGGCAAGGTTATGGCAATCGAAGTTAGACGCt                                                                    <  2:836559/134‑1 (MQ=255)                                                                     aaCGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCAATGTACAATACCATTTTCCGTATCAATAATTCCCGGCAAGGTTATGGCAATCGAAGTTAGACGCt                                                                    >  1:836559/1‑134 (MQ=255)                                                                     aaCGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCAATGTACAATACCATTTTCCGTATCAATAATTCCCGGCAAGGTTATGGCAATCGAAGTTAGACGCTCAAGtt                                                              >  2:41775/1‑140 (MQ=255)                                                                                    ctCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCAATGTACAATACCATTTTCCGTATCAATAATTCCCGGCAAGGTTATGGCAATCGAAGATAGCCGATCAAGTTTTTTCTGGTGGCGGa                                               >  1:458377/1‑140 (MQ=255)                                                                                       cAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCAATGTACAATACCATTTTCCGTATCAATAATTCCCGGCAAGGTTATGGCAATCGAAGTTAGACGGTCAAGTTTTCTCTGGTGGCGGATaa                                            >  1:355791/1‑140 (MQ=255)                                                                                              tCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCAACGTACAATACCATTTTCCGTATCAATAATTCCCGGCAAGGTTATGGCAATCGAAGTTAGACGCTCAAGTTTTTTCTGGTGGCGGATAAAAAACTg                                     <  2:696307/140‑1 (MQ=255)                                                                                                gCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCAATGTACAATACCATTTTCCGTATCAATAATTCCCGGCAAGGTTATGGCAATCGAAGTTAGACGCTCAAGTTTTTTCTGGTGGCGGATAAAAAACtaat                                   <  2:244114/140‑4 (MQ=255)                                                                                                gCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCAATGTACAATACCATTTTCCGTATCAATAATTCCCGGCAAGGTTATGGCAATCGAAGTTAGACGCTCAAGTTTTTTCTGGTGGCGGATAAAAAACTGAt                                   <  2:1009633/140‑1 (MQ=255)                                                                                                gCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCAATGTACAATACCATTTTCCGTATCAATAATTCCCGGCAAGGTTATGGCAATCGAAGTTAGACGCTCAAGTTTTTTCTGGTGGCGGATAAAAAACTGAt                                   <  2:31190/140‑1 (MQ=255)                                                                                                     gCGGCATCTCTTTTAGATCCTCATACAACGGCATGCAATGTACAATTCCATTTTCCGTATCAATAATTCCCGGCAAGGTTATGGCAATCGAAGTTAGACGCTCAAGTTTTTTCTGGTGGCGGATAAAAAACTGATCGata                              <  1:473692/140‑1 (MQ=255)                                                                                                       ggCATCTCTTTTACATCCTCGTAGAACGGCATGCAATGTACAATACCATTTTCCGTATCAATAATTCCCGGCAAGGTTATGGCAATCGAAGTTAGACGCTCAAGTTTTTTCTGGTGGCGGATAAAAAACTGATCGATATg                            <  2:977331/140‑1 (MQ=255)                                                                                                       ggCATCTCTTTTACATCCTCGTAGAACGGCATGCAATGTACAATACCATTTTCCGTATCAATAATTCCCGGCAAGGTTATGGCAATCGAAGTTAGACGCTCAAGTTTTTTCTGGTGGCGGATAAAAAACTGATCGATATg                            <  1:178120/140‑1 (MQ=255)                                                                                                       ggCATCTCTTTTACATCCTCGTAGAACGGCATGCAATGTACAATACCATTTTCCGTATCAATAATTCCCGGCAAGGTTATGGCAATCGAAGTTAGACGCTCAAGTTTTTTCTGGTGGCGGATAAAAAACTGATCGATATg                            <  1:980193/140‑1 (MQ=255)                                                                                                          aTCTCTTTTACATCCTCGTAGAACGGCATGCAATGTACAATACCATTTTCCGTATCAATAATTCCCGGCAAGGTTATGGCAATCGAAATTAGACGCTCAAGTTTTTTCTGGTGGCGGATAAAAAACTGATCGATATGGGa                         <  2:604264/140‑1 (MQ=255)                                                                                                            ctctTTTACATCCTCATAGAACGGCATGCCATGTACACTACCATTTTCCGTATCAATAATTCCCGGCAAGGTTATGGCAATCGAAGTTAGACGCTCAAGTTTTTTCTGGTGGCGGATAAAAAACTGATCGATATGGGaaa                       <  1:424571/140‑1 (MQ=255)                                                                                                                       ccTCGTAGAACGGCATGCAATGTACAATACCATTTTCCGTATCAATAATTCCCGGCAAGGTTATGGCAATCGAAGTTAGACGCTCAAGTTTTTTCTGGTGGCGGATAAAAAACTGATCGATATGGGaaa                       >  2:616153/1‑129 (MQ=255)                                                                                                                       ccTCGTAGAACGGCATGCAATGTACAATACCATTTTCCGTATCAATAATTCCCGGCAAGGTTATGGCAATCGAAGTTAGACGCTCAAGTTTTTTCTGGTGGCGGATAAAAAACTGATCGATATGGGaaa                       <  1:616153/129‑1 (MQ=255)                                                                                                                         tCGTAGAACGGCATGCAATGTACAATACCATTTTCCGTATCAATAATTCCCGGCAAGGTTATGGCAATCGAAGTTAGACGCTCAAGTTTTTTCTGGTGGCGGATAAAAAACTGATCGATATGGGAAATAATACGATCCAg          >  1:705918/1‑140 (MQ=255)                                                                                                                          cGTAGAACGGCATGCAATGTACAATACCATTTTCCGTATCAATAATTCCCGGCAAGGTTATGGCAATCGAAGTTAGACGCTCAAGTTTTTTCTGGTGGCGGATAAAAAACTGATCGATATGGGAAATAATACGATCCAGc         >  2:50969/1‑140 (MQ=255)                                                                                                                          cGTAGAACGGCATGCAATGTACAATACCATTTTCCGTATCAATAATTCCCGGCAAGGTTATGGCAATCGAAGTTAGACGCTCAAGTTTTTTCTGGTGGCGGAGAAAAAACTGATCGACATGGGAAATAATACGATCCAGc         >  1:443773/1‑140 (MQ=255)                                                                                                                                 cGGCATGCAATGTACAATACCATTTTCCGTATCAATAATTCCCGGCAAGGTTATGGCAATCGAAGTTAGACGCTCAAGTTTTTTCTGGTGGCGGATAAAAAACTGATCGATATGGGAAATAATACGATCCAGCAATGGCa  >  2:751904/1‑140 (MQ=255)                                                                                                                                 cGGCATGCAATGTACAATACCATTTTCCGTATCAATAATTCCCGGCAAGGTTATGGCAATCGAAGTTAGACGCTCAAGTTTTTTCTGGTGGCGGATAAAAAACTGATCGATATGGGAAATAATACGATCCAGCAATGGCa  <  2:674824/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                    CCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCGATGTACAATACCATTTTCCGTATCAATAATTCCCGGCAAGGTTATGGCAATCGAAGTTAGACGCTCAAGTTTTTTCTGGTGGCGGATAAAAAACTGATCGATATGGGAAATAATACGATCCAGCAATGGCA  >  NC_000913/1667933‑1668200

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