Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,963,226	G→T	A335D (GCT→GAT)	gpp ←	guanosine‑5'‑triphosphate,3'‑diphosphate phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,963,226	0	G	T	100.0%	44.8 / NA	14	A335D (GCT→GAT)	gpp	guanosine‑5'‑triphosphate,3'‑diphosphate phosphatase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (5/9);  total (5/9)

TAAGCAGCGTGTTGCGGCGCTTGTTTGAAGTCAACGCTCAGGCCGATTTCATGAAGCTGACAGGCGCTGATGAGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGC  >  NC_000913/3963094‑3963362                                                                                                                                     |                                                                                                                                         tAAGCAGCGTGTTGCGGCGCTTGTTTGAAGTCAACGCTCAGGCCGATTTCATGAAGCTGACAGGCGCTGATGAGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTATCCGCAAc                                                                                                                                   <  2:292549/140‑1 (MQ=255)         gtgtTGCGGCGCTTGTTTGAAGTCAACGCTCAGGCCGATTTCATGAAGCTGACAGGCGCTGATGAGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTATCCGCAACTTTGGCTa                                                                                                                           >  2:190433/1‑140 (MQ=255)                   gcTTGTTTGAAGTCAACGCTCAGGCCGATTTCATGAAGCTGACAGGCGCTGATGAGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTATCCGCAACTTTGGCTACGCGCTGTGc                                                                                                                 <  1:165560/140‑1 (MQ=255)                                 aaCGCTCAGGCCGATTTCATGAAGCTGACAGGCGCTGATGAGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTATCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGa                                                                                                   >  2:165383/1‑140 (MQ=255)                                 aaCGCTCAGGCCGATTTCATGAAGCTGACAGGCGCTGATGAGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTATCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGa                                                                                                   <  2:225505/140‑1 (MQ=255)                                                ttCATGAAGCTGACAGGCGCTGATGAGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTATCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGtt                                                                                    <  1:15609/140‑1 (MQ=255)                                                                 cgcTGATGAGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTATCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTa                                                                   <  2:259656/140‑1 (MQ=255)                                                                                       tATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTATCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGa                                             <  2:17809/140‑1 (MQ=255)                                                                                         ttGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTATCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACg                                           <  1:190433/140‑1 (MQ=255)                                                                                               cAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTATCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGa                                     >  2:32999/1‑140 (MQ=255)                                                                                                  gATGCCATTCGTTTTCCACCTGATCGAAGAAGTTATCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATACCCTGCTCGACGGCCAGATGc                                  >  2:78786/1‑140 (MQ=255)                                                                                                                  cACCTGATCGAAGAAGTTATCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCa                  <  1:32999/140‑1 (MQ=255)                                                                                                                            aagaagTTATCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACg        <  1:165383/140‑1 (MQ=255)                                                                                                                                  ttATCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAgcgc  >  1:347859/1‑140 (MQ=255)                                                                                                                                     |                                                                                                                                         TAAGCAGCGTGTTGCGGCGCTTGTTTGAAGTCAACGCTCAGGCCGATTTCATGAAGCTGACAGGCGCTGATGAGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGC  >  NC_000913/3963094‑3963362

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