Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,522,608	C→T	R48R (CGG→CGA)	hofM ←	DNA utilization protein HofM

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,522,608	0	C	T	100.0%	25.5 / NA	11	R48R (CGG→CGA)	hofM	DNA utilization protein HofM
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (7/4);  total (7/4)

TCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGCCAGGCCGTTTGCTCCCGCTCACCAAGGGACATCGACGGGCGCGGAAATGACCGCTGTAATGTGCGACTGGCGGGAAACGCCAACATAATGTGATGACGCTGCGGCAGTTCGCGACTCCACGGTAACAACGTTTTAGCCAGCTGCTGCGCATCAACAATCCGCCCATCTTTGATAATGTCGTTCTCCAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCG  >  NC_000913/3522384‑3522822                                                                                                                                                                                                                                 |                                                                                                                                                                                                                       tCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGCCAGGCCGTTTGCTCCCGCTCACCAAGGGACATCGCCGGGCGCGGAAATGACCGCTGTAATGTGCGACTGGCGGGAAACGCCAACATAATGTGATGACGCTGCGGCAGTTCGCGACTCCACGGTAACAACGTTTTAGCCAGCTGCTGCGCATCAACAATTCGCCCATCTTTGATAATGTCGTtc                                                                                                                                                                                                <  1:218846/249‑1 (MQ=255)                                                               naGCCAGGCCGTTTGCTCCCGCTCACCAAGGGACATCGACGGGCGCGGAAATGACCGCTGTAATGTGCGACTGGCGGGAAACGCCCACATAATGTGATGACGCTGCGGCAGTTCGCGACTCCACGGTAACAAAGTTTTAGCCAGCTGCTGCGCATCAACAATTCGCCCATCTTTGATAATGTCGt                                                                                                                                                                                                  >  1:56165/2‑185 (MQ=255)                                                                                      cACCAAGGGACATCGACGGGCGCGGAAATGACCGCTGTAATGTGCGACTGGCGGGAAACGCCACCATAATGTGATGACGCTGAGGCCGTTCGCGACTCCACGGTAACAACGTTTTAGCCCGCTGCTGCGCATCAACAATTCGCCCATCTTTGATAATGTCGTTCTCCAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTtcgc                                                                                                                                                  >  1:439975/1‑210 (MQ=255)                                                                                        ccAAGGGACATCGACGGGCGCGGAAATGACCGCTGTAATGTGCGACTGGCGGCAAACGCCAACATAAAGAGATGACGCTACGGACGTTCGCGACTCCACGGCAACAACGTTTTAGCCAGCTGCAGCGCATCAACAATTCGCCCCTCTTTGATAATTTCGTTCTCCAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTACGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAAtt                                                                                                       >  1:341569/1‑250 (MQ=255)                                                                                                                                                            gtgATGACGCTGCGGCAGTTCGCGACTCCACGGTAACAACGTTTTAGCCAGCTGCTGCGCATCAACAATTCGCACATCTTTGATAATGTCGTTCTCCAGCGGCAACCGc                                                                                                                                                                                 >  1:485380/1‑109 (MQ=255)                                                                                                                                                                    gCTGCGGCAGTTCGCGACTCCACGGTAACAACGTTTTAGCCAGCTGCTGCGCATCAACAATTCGCCCATCTTTGATAATGTCGTTCTCCAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTGGCTGTTGTAAATGCAAACCAATTTTCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCAc                                                     <  1:86587/225‑1 (MQ=255)                                                                                                                                                                                cgcgACTCCACGGTAACAACGTTTTAGCCAGCTGCTGCGCATCAACAATTGGCCCATCTTTGATAATGTCGTTCTCCAGCGGCAACCGCCACCAGAGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACAGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTgt               >  1:185764/1‑250 (MQ=255)                                                                                                                                                                                   naCTCCACGGTAACAACGTTTTAGCCAGCTGCTGCGCATCAACAATTCGCCCATCTTTGATAATGTCGTTCTCCAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAAt                            >  1:147401/2‑235 (MQ=255)                                                                                                                                                                                     cTCCACGGTAACAACGTTTTAGCCAGCTGCTGCGCATCAACAATTCGCCCATCTTTGATAATGTCGTTCTCCAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCt                                                                                                                                                         >  1:109253/1‑108 (MQ=255)                                                                                                                                                                                     cTCCACGGTAACAACGTTTNAGCCAGCGGNNGCNCATCANCNNNTNNNCCNTCTTTGATAATGTCGTTCTNCANNGGCAACNGCCACCAGCGTTGCAAAAAGCATTCt                                                                                                                                                         <  2:109253/108‑1 (MQ=255)                                                                                                                                                                                              aacaacGTTTTAGCCAGCTGCTGCGCATCAACAATTCGCCCATCTTTGATAATGTCGTTCTCCAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACcg  <  1:356742/250‑1 (MQ=255)                                                                                                                                                                                                                                 |                                                                                                                                                                                                                       TCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGCCAGGCCGTTTGCTCCCGCTCACCAAGGGACATCGACGGGCGCGGAAATGACCGCTGTAATGTGCGACTGGCGGGAAACGCCAACATAATGTGATGACGCTGCGGCAGTTCGCGACTCCACGGTAACAACGTTTTAGCCAGCTGCTGCGCATCAACAATCCGCCCATCTTTGATAATGTCGTTCTCCAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCG  >  NC_000913/3522384‑3522822

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