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%	61.8 / NA	21	R89C (CGC→TGC)	mcbR	DNA‑binding transcriptional dual regulator McbR
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (13/8);  total (13/8)

GAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAACCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTGAGATGATTGAGCAACTGTGGGTCAGGATGGGGCCTGGTTTACATTATC  >  NC_000913/1520282‑1520766                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 gAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAACCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCTGCTAc                                                                                                                                                                                                                                             <  1:21285/250‑1 (MQ=255)       tGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAACCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCTGCTACGAACTgg                                                                                                                                                                                                                                      >  1:499306/1‑251 (MQ=255)                             gagaACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAACCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCTGCTACGAACTGGAATTAATGGCAGTTGCTCTGGc                                                                                                                                                                                                                <  1:484567/251‑1 (MQ=255)                              agaACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAACCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCTGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCt                                                                                                                                                                                                               >  1:216276/1‑251 (MQ=255)                                    aacTGAAACATCAGCTTAGTATTGGCGCATTAAAACCTGGCGCACGCCTGATTACTAAAAATATGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTAGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCCATAGGGACTGCTACGAACTGGAATTAATGGCAGGTGCCCTGGCTGtaga                                                                                                                                                                                                          >  2:316681/3‑247 (MQ=255)                                      cTGAAACATCAGCTTAGTATTGGCGCATTAAAACCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCTGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGaaa                                                                                                                                                                                                        >  2:198293/1‑250 (MQ=255)                                               cAGCTTATGATTGGCGCTTTAAAACCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTATGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGATGAAATCAAGCGGATCTGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAATCCTCAcccc                                                                                                                                                                                              <  2:557106/251‑1 (MQ=255)                                                       tATTGGCGCATTAAAACCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTAAGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGAATTTACAGGTCCGGAAGTGGGGAAACGTCAATTAGACGAATTCCATCGGATCTGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGTAAAACTCACCCCGAAAGAcc                                                                                                                                                                                      >  1:271000/1‑251 (MQ=255)                                                                                         aCTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCTGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTAca                                                                                                                                                    <  1:316681/251‑1 (MQ=255)                                                                                            aaaaaTCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCTGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTAcag                                                                                                                                                   <  1:472295/249‑1 (MQ=255)                                                                                                        gAGCAATTAAGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCTGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCACGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACATCAGGCGCaaaaaa                                                                                                                                      >  2:221711/1‑250 (MQ=255)                                                                                                                               acCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCTGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCAtt                                                                                                               >  1:56441/1‑250 (MQ=255)                                                                                                                                                                    aaTGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCTGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGAAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCg                                                                         >  1:493908/1‑251 (MQ=255)                                                                                                                                                                      tGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCTGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAAAAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCt                                                                       >  1:116851/1‑251 (MQ=255)                                                                                                                                                                                gTCGCTCCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCTGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGcc                                                             <  1:490669/251‑1 (MQ=255)                                                                                                                                                                                   gCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCTGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCAt                                                          <  1:198293/251‑1 (MQ=255)                                                                                                                                                                                                 aTTTACAGATCCGGAAGTGGGGAAACGTCAGTTAGCCGAAATCAATCGGATCTGCTACGAACTGGAATTAATGGCAGTTGCTCTTGCTGTTGAAAACCTCACCCCGCAAGACCTTGCCGAACTCCAGGAATTTATCGGGAAATTACCGCAGGCGCAAGAGAAGGGCGACCTGGGACAAATCATTAATGTCAAAAGGCTATTTCCCTTAGCGCTTTATCATCGCCCAAATATGCCCATcccgttggagatgt                                            >  2:557001/1‑239 (MQ=255)                                                                                                                                                                                                       aGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCTGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTGAGATGATTGAGc                                      >  1:194768/1‑251 (MQ=255)                                                                                                                                                                                                         ttCCGGAAATGGGGAAACGTCAATTAGACGAAATCAATCGGATCTGCTACGAACTGGAATTAATGGAAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGAAAGACAAGGGCGCCAAGGAACAAATCATTAATGTAAACAGGCTATTACACTTAGGGATTTATCATCGCTAAAAAATTCCCATCCTGTGTGAGATGGTTGAgta                                     >  2:492635/1‑248 (MQ=255)                                                                                                                                                                                                                           gTCAATTAGACGAAATCAATCGGATCTGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCCCCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCCGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATAGCTCAAATATGCCCATTCTGTGTTGGATGATTTGGCAACTGTGGGTCAGGATgggg                  >  2:66451/1‑251 (MQ=255)                                                                                                                                                                                                                                           ccaccGGATCTGCTACGAACTGGAAGTAAGAGCAGTTGCGCTGGCTGTTGCATTCCTCACCCCGCATGAACTTGCATAACTACAGGAATTCCGCGAGAAATTACAGCAGGTGCAAGAAAAGGGCGACATGGAACAAATCATTAATTTAAACAGGCTATTTCGCTTAGCGTTTTATCATCGCTCAAATAGGCCCATCCTTTGTCAGATGATTGAGCATCTGTTGGTCTGGCTGGGGCCTGGTTTACATTAtc  <  2:468233/247‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 GAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAACCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTGAGATGATTGAGCAACTGTGGGTCAGGATGGGGCCTGGTTTACATTATC  >  NC_000913/1520282‑1520766

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