Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,522,408	T→C	D115G (GAC→GGC)	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,408	0	T	C	100.0%	50.8 / NA	16	D115G (GAC→GGC)	hofM	DNA utilization protein HofM
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (5/11);  total (5/11)

GTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGCCAGGCCGTTTGCTCCCGCTCACCAAGGGACATCGACGGGCGCGGAAATGACCGCTGTAATGTGCGACTGGCGGGAAACGC  >  NC_000913/3522279‑3522529                                                                                                                                  |                                                                                                                          gTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGCCCg                                                                                                                        <  2:538877/133‑1 (MQ=255) gTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGCCCg                                                                                                                        <  2:393425/133‑1 (MQ=255)   gATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGCCCGGa                                                                                                                      >  1:555902/1‑133 (MQ=255)         gcTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGCCCGGCTCCata                                                                                                                >  1:610761/1‑133 (MQ=255)                       gCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGTCTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGCCCGGGTCCATATTCACCTCCCGGGc                                                                                                  >  2:18467/1‑133 (MQ=255)                       gCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGCCCGGATCCATATCCAGCTCGCGGGc                                                                                                  >  2:9716/1‑133 (MQ=255)                       gCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGCCCGGATCCATATCCAGCTCGCGGGc                                                                                                  >  1:545964/1‑133 (MQ=255)                         aaCCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGCCCGGATCCATATCCAGCTCGCGGGCCa                                                                                                <  1:599602/133‑1 (MQ=255)                             tttCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGCCCGGATCCATATCCAGCTCGCGGGCCATCGt                                                                                            <  2:555902/133‑1 (MQ=255)                                           agcagcGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGCCCGGATCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGCCAgg                                                                              <  1:9716/133‑1 (MQ=255)                                                 gTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGGCTTCGCTGTAATCGAAGCGCAGGGAGCCCGGCTCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGGCAGGGCGttt                                                                        >  2:9614/1‑133 (MQ=255)                                                          tctTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGCCCGGATCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGCCAGGCCGTTTGCTCCCGCt                                                               <  2:545964/133‑1 (MQ=255)                                                            tttgctttgcGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGCCCGGATCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGCCAGGCCGTTTGCTCCCGCTCa                                                             <  2:414608/133‑1 (MQ=255)                                                                                       ggCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGCCCGGATCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGCCAGGCCGTTTGCTCCCGCTCACCAAGGGACATCGACGGGCGCGGAAAt                                  <  1:9614/133‑1 (MQ=255)                                                                                       ggCGGGACTGAGTGAATCATCGTTGTAATCGAAGCGCAGGGAGCCCGGATCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGCCAGGCCGTTTGCTCCCGCTCACCAAGGGACATCGACGGGCGCGGAAAt                                  <  1:259201/133‑1 (MQ=255)                                                                                                           cGCTGTAATCGAAGCGCAGGGAGCCCGGATCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGCCAGGCCGTTTGCTCCCGCTCACCAAGGGACATCGACGGGCGCGGAAATGACCACTGTAATGTGCGACt              >  1:217200/1‑133 (MQ=255)                                                                                                                 aaTCGAAGCGCAGGGAGCCCGGATCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGCCAGGCCGTTTGCTCCCGCTCACCAAGGGACATCGACGGGCGCGGAAATGACCGCTGTAATGTGCGACTGGCggg        >  1:313271/1‑133 (MQ=255)                                                                                                                       aGCGCAGGGAGCCCGGATCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGCCAGGCCGTTTGCTCCCGCTCACCAAGGGACATCGACGGGCGCGGAAATGACCGCTGTAATGTGCGACTGGCGGGAAACGc  <  1:357446/133‑1 (MQ=255)                                                                                                                          gcAGGGAGCCCGGATCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGCCAGGCCGTTTGCTCCCGCTCACCAAGGGACATCGACGGGCGCGGAAATGACCGCTGTAATGTGCGACTGGCGGGaaa     >  2:35392/1‑127 (MQ=255)                                                                                                                          gcAGGGAGCCCGGATCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGCCAGGCCGTTTGCTCCCGCTCACCAAGGGACATCGACGGGCGCGGAAATGACCGCTGTAATGTGCGACTGGCGGGaaa     <  1:35392/127‑1 (MQ=255)                                                                                                                                  |                                                                                                                          GTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGCCAGGCCGTTTGCTCCCGCTCACCAAGGGACATCGACGGGCGCGGAAATGACCGCTGTAATGTGCGACTGGCGGGAAACGC  >  NC_000913/3522279‑3522529

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