Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,123,886	C→T	A191T (GCC→ACC)	cytR ←	DNA‑binding transcriptional repressor CytR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,123,886	0	C	T	100.0%	26.5 / NA	12	A191T (GCC→ACC)	cytR	DNA‑binding transcriptional repressor CytR
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (5/7);  total (5/7)

AGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCCTTGAACCCAGCAACAACATG  >  NC_000913/4123719‑4124079                                                                                                                                                                        |                                                                                                                                                                                                  aGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGTTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATtggatttgcccctatatttccagacatctgttatcacttaan                                                                                                                <  1:63315‑M1/251‑43 (MQ=255)        tgtgGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGTTAGACAGTCAATCCGTGTATGCCCTTGCTCATATAAATAATtggatttgcccctatattgccagacatctgttatcacttaacccatta                                                                                                          <  1:211977‑M1/250‑49 (MQ=255)                  naagCAGCTGCTGCATCGCTTTGCTTACGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGAAGTGACACAGCGGCATCTCTTCGGGACCGGTTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctg                                                                                               >  1:175705‑M1/2‑192 (MQ=255)                           gctgcATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGTTATACAGCCAATCCGTTTATGCCCTTGCTc                                                                                                                                                                     <  1:93673/172‑1 (MQ=255)                                                                gCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGTTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcgga                                                 <  1:102961‑M1/251‑106 (MQ=255)                                                                    cGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGTTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgcc                                             <  1:5280‑M1/251‑110 (MQ=255)                                                                       ncgATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGTTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccatt                                          >  1:25737‑M1/2‑139 (MQ=255)                                                                        cGATGCATTGCGGATCAACAATAATGCCGAAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGTTATACAGCCAATCCGttt                                                                                                                                                                                 >  1:70421/1‑115 (MQ=255)                                                                                    natCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGTTATCCAGCCAATCCGTTTATGCCCTTGCTCCTATAAATAATtggctttgcccctatatttccagacatctgttatcacttaacccattaccagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgcgctcgga                                >  1:63564‑M1/2‑126 (MQ=255)                                                                                                         gACGCAGCGCCTGAACATAGCCTTGCAGGAGGTAGTGACACAGCGGCATCTCTTCGGGACCGGTTATACAGCCCATCCGTTTATGCCCTTGCTCATATAAATAATtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgacgcagatattcccgtggcgcgcgataaaccagcgcactatgcggatgccattcgttatactgctcgaacgcctctgcaaggttctt        >  1:18058‑M1/1‑105 (MQ=255)                                                                                                               gcgcCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGTTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttctttgctgc  <  1:80404‑M1/251‑153 (MQ=255)                                                                                                                                               cacaGCGGCATCTCTTCGGGACCGGTTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATtggatttgcccc                                                                                                                                              <  1:106407‑M1/79‑13 (MQ=255)                                                                                                                                                                        |                                                                                                                                                                                                  AGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCCTTGAACCCAGCAACAACATG  >  NC_000913/4123719‑4124079

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