Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,864,572	A→G	Q225R (CAG→CGG)	yidP →	UTRA domain‑containing GntR family putative transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,864,572	0	A	G	100.0%	50.4 / NA	16	Q225R (CAG→CGG)	yidP	UTRA domain‑containing GntR family putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (6/10);  total (6/10)

ACGCCAGTGCTCGCCGATAGACTGCTGGCGCGGCAAATGAAGGTAGCGGAACACACGCCACTGCTGCGGATCACCTCGTTGTCATATAGCGAGAGCGGGGAGTTTTTGAATTATTCAGTGATGTTCAGAAATGCCAGCGAATACCAGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAACATTACCAACGGATAGACAGTGGCGTAAGAG  >  NC_000913/3864427‑3864699                                                                                                                                                  |                                                                                                                                aCGCCAGTGCTCGCCGATAGACTGCTGGCGCGGCAAATGAAGGTAGCGGAACACACGCCACTGCTGCGGATCACCTCGTTGTCATATAGCGAGAGCGGGGAGTTTTTGAATTATTCAGTGATGTTCAGAAATGCCAGCGAATACCGGGTgg                                                                                                                            >  2:20917/1‑151 (MQ=255)  cGCCAGTGCTCGCCGATAGACTGCTGGCGCGGCAAATGAAGGTAGCGGAACACACGCCACTGCTGCGGATCACCTCGTTGTCATATAGCGAGAGCGGGGAGTTTTTGAATTATTCAGTGATGTTCAGAAATGCCAGCGAATACCGGGTGGa                                                                                                                           >  1:271337/1‑151 (MQ=255)                                    aaTGAAGGTAGCGGAACACACGCCACTGCTGCGGATCACCTCGTTGTCATATAGCGAGAGCGGGGAGTTTTTGAATTATTCAGTGATGTTCAGAAATGCCAGCGAATACCGGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGt                                                                                         <  2:311982/151‑1 (MQ=255)                                            tAGCGGAACACACGCCACTGCTGCGGATCACCTCGTTGTCATATAGCGAGAGCGGGGAGTTTTTGAATTATTCAGTGATGTTCAGAAATGCCAGCGAATACCGGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGa                                                                                  <  2:135413/150‑1 (MQ=255)                                                 gAACACACGCCACTGCTGCGGATCACCTCGTTGTCATATAGCGAGAGCGGGGAGTTTTTGAATTATTCAGTGATGTTCAGAAATGCCAGCGAATACCGGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTc                                                                             <  2:111703/150‑1 (MQ=255)                                                 gAACACACGCCACTGCTGCGGATCACCTCGTTGTCATATAGCGAGAGCGGGGAGTTTTTGAATTATTCAGTGATGTTCAGAAATGCCAGCGAATACCGGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTc                                                                             >  2:265981/1‑150 (MQ=255)                                                          ccACTGCTGCGGATCACCTCGTTGTCATATAGCGAGAGCGGGGAGTTTTTGAATTATTCAGTGATGTTCAGAAATGCCAGCGAATACCGGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGc                                                                    <  2:257420/150‑1 (MQ=255)                                                               gctgcGGATCACCTCGTTGTCATATAGCGAGAGCGGGGAGTTTTTGAATTATTCAGTGATGTTCAGAAATGCCAGCGAATACCGGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGc                                                               >  1:155724/1‑150 (MQ=255)                                                                                 gTCATATAGCGAGAGCGGGGAGTTTTTGAATTATTCAGTGATGTTCAGAAATGCCAGCGAATACCGGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATa                                            >  2:178876/1‑151 (MQ=255)                                                                                                   ggAGTTTTTGAATTATTCAGTGATGTTCAGAAATGCCAGCGAATACCGGGTGGAGTACCATTTACGGCGCCTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAACATTAc                           <  2:163441/150‑1 (MQ=255)                                                                                                              attattCAGTGATGTTCAGAAATGCCAGCGAATACCGGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAACATTACCAACGGATAGAc               <  1:60444/151‑1 (MQ=255)                                                                                                                         aTGTTCAGAAATGCCAGCGAATACCGGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAACATTACCAACGGATAGACAGTGGCGTAag    <  1:239911/151‑1 (MQ=255)                                                                                                                           gTTCAGAAATGCCAGCGAATACCGGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAACATTACCAACGGATAGACAGTGGCGTAagag  >  1:156848/1‑151 (MQ=255)                                                                                                                            ttaaGAAATGCCAGCAAATCCCGGGTGGAGTACCATTAACGGCGACTCCACCCGGAAAAGAGTTACCCGATCCTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAACATTACCAACGGATAGGCAGTGGCGTAagag  <  2:165854/147‑1 (MQ=255)                                                                                                                            ttCAGAAATGCCAGCGAATACCGGGTGGAGTCCCATTTACGGCGACCCCACCCGGAAAAGATTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAACATTACCAACGGATAGACAGTGGCGTAagag  <  2:182241/150‑1 (MQ=255)                                                                                                                            tccAGAAATGCCAGCAAATACCGGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGATTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAACATTACCAACGGATAGACAGTGGCGTAagag  <  2:1406/148‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                ACGCCAGTGCTCGCCGATAGACTGCTGGCGCGGCAAATGAAGGTAGCGGAACACACGCCACTGCTGCGGATCACCTCGTTGTCATATAGCGAGAGCGGGGAGTTTTTGAATTATTCAGTGATGTTCAGAAATGCCAGCGAATACCAGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAACATTACCAACGGATAGACAGTGGCGTAAGAG  >  NC_000913/3864427‑3864699

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