Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,520,526	C→T	R89C (CGC→TGC)	mcbR →	DNA‑binding transcriptional dual regulator McbR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,520,526	0	C	T	100.0%	105.6 / NA	32	R89C (CGC→TGC)	mcbR	DNA‑binding transcriptional dual regulator McbR
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/16);  total (16/16)

AATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAACCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTGAGATGATTGAGCAACTGTGGGTCAGGATGGGGCCTGGTTTACATTATCTCTAT  >  NC_000913/1520279‑1520771                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      aaTGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAACCTGGCGCACGCCTGATTACTAAAACTCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCTGCt                                                                                                                                                                                                                                                    >  1:365971/1‑251 (MQ=255)   tGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAACCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCTGCTAc                                                                                                                                                                                                                                                  >  1:528051/1‑251 (MQ=255)         aTGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTAGTGGCGCATTAAAACCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCTGCTACGAACTg                                                                                                                                                                                                                                            <  2:23562/251‑1 (MQ=255)             cAGTTTGACTCTGCAGGTTGAGAACGGCCTGAAACACCAGCTTAGTATTGGCGCATTAAAACCTGGCGCACGCCTGATTACGAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAGTGCGCTTTCTGTCGCACCTGCACAAGCAGTTACAGTACCGGAAGTGGGTAAACGTCGATGAGACGAAAGCAATCTGGTCTGCGACGAACTGGAAt                                                                                                                                                                                                                                        <  2:60624/251‑1 (MQ=255)                  tGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAACCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCCCCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCTGCTACGAACTGGAATTAATg                                                                                                                                                                                                                                   >  2:437757/1‑251 (MQ=255)                  tGACGCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAACCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCTGCTACGAACTGGAATTAAt                                                                                                                                                                                                                                    <  2:332145/250‑1 (MQ=255)                         gCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAACCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCTGCTACGAACTGGAATTAATGGCAGtt                                                                                                                                                                                                                             <  1:587915/250‑1 (MQ=255)                         gCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAACCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCTGCTACGAACTGGAATTAATGGCAGTTg                                                                                                                                                                                                                            <  1:183544/251‑1 (MQ=255)                              ttGAGAACGACCTGAAACATCAGTTTAGTATTGGCGCATTAAAACCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTCGACGAAATCAATCGGATCTGCTACGAACTGGAATTACTGGCAGTTGCTCTg                                                                                                                                                                                                                       >  1:189733/1‑251 (MQ=255)                                      gACCTGAAACATCAGCTTAGTATTGGCGCATTAAAACCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCTGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTg                                                                                                                                                                                                                >  1:244786/1‑250 (MQ=255)                                                                          ccTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCTGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAgg                                                                                                                                                                            >  1:64623/1‑250 (MQ=255)                                                                             ggCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCTGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAAtt                                                                                                                                                                        <  2:208822/251‑1 (MQ=255)                                                                                 cACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCTGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTc                                                                                                                                                                    <  1:518918/251‑1 (MQ=255)                                                                                        gATTACTAAAAATCTGGCGGAGCAGTTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCTGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAAt                                                                                                                                                             >  1:368062/1‑251 (MQ=255)                                                                                            aCTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTAAAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCTGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGACTTGCTCGGGAAATTAc                                                                                                                                                          >  1:22129/1‑250 (MQ=255)                                                                                              tAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCTGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAAtttaagc                                                                                                                                                       >  1:324694/1‑246 (MQ=255)                                                                                                      tGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCTGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCa                                                                                                                                                >  1:251001/1‑250 (MQ=255)                                                                                                                                                 attattACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGACGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCTGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAACCAGGCTAt                                                                                                    >  2:235491/1‑251 (MQ=255)                                                                                                                                                   tattaCGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCTGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAACTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTAtt                                                                                                   <  1:311671/250‑1 (MQ=255)                                                                                                                                                    attaCGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCTGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCGCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTAttt                                                                                                  <  2:388448/250‑1 (MQ=255)                                                                                                                                                              gTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCTGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGAt                                                                                       >  2:10066/1‑251 (MQ=255)                                                                                                                                                                      gAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCTGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTatcatc                                                                               >  1:120942/1‑251 (MQ=255)                                                                                                                                                                      gAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCTGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCGATTTCGCTTAGCGATTTatcatc                                                                               <  2:368062/251‑1 (MQ=255)                                                                                                                                                                                  tGTCGCCCCTGCAAAAGCATCTACAGTTCCGGAAGTGGGGAAACGCCAATTAGCCGAAATCAATCGGATCTGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGATCAAATATGc                                                                   >  2:105895/1‑251 (MQ=255)                                                                                                                                                                                          cTGCACACGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCTGCTACGAACTGGAATTAATGGCAGTTGCTATGGATGTTGAAAACATCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGGCCTGGAACAAATAATTAATGTAAACAGGCTCTTTCGCTTAGCGCTTTatcat                                                                                >  2:496432/1‑230 (MQ=255)                                                                                                                                                                                          cTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCTGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTatcat                                                                                <  1:496432/230‑1 (MQ=255)                                                                                                                                                                                                aaGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCTGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCGCAAATATGCTCATCCTGTGTgag                                                     <  2:64623/251‑1 (MQ=255)                                                                                                                                                                                                      ttACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCTGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATGTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTTTGAgatgat                                                <  2:194575/250‑1 (MQ=255)                                                                                                                                                                                                           gTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCTGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTGAGATGATTGAGCa                                          <  1:373078/251‑1 (MQ=255)                                                                                                                                                                                                           gTTCCGGAAGCGGGGAAACGTCAATTAGACGAAATCAATCGGATCTGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCACACCCCGCAAGACCTTGCGGAAATCCAGGAATTGCTCGAGAAATTAAAGCAGGCGAAAGCAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTACAGATTTATCATCGCTCAATTATGGCCATCCTGTGTGAGATGATTGAgg                                           >  2:392640/1‑249 (MQ=255)                                                                                                                                                                                                                                 aaTTAGACGAAATCAATCGGATCTGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTCAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTGAGATGATTGAGCAACTGTGGGTCAGGATGGGGCCt                    <  2:419787/251‑1 (MQ=255)                                                                                                                                                                                                                                                   ggATCTGCGACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGTAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTGAGATGATTGAGCAACTGTGGTTCAGGATGGGGTCTGGTTTACATTATCTCTAt  <  1:75596/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      AATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAACCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTGAGATGATTGAGCAACTGTGGGTCAGGATGGGGCCTGGTTTACATTATCTCTAT  >  NC_000913/1520279‑1520771

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