Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,586,164	G→A	D220D (GAC→GAT)	ggt ←	glutathione hydrolase proenzyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,586,164	0	G	A	100.0%	86.7 / NA	26	D220D (GAC→GAT)	ggt	glutathione hydrolase proenzyme
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (14/12);  total (14/12)

TTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTGTCGCCCTTTTTCAGCGGCTCGCCCTCTTTCCAGAAGATAGCTTTACTGTTTTCGTGATTCGGCAACACTTCGCTACCGTAGGTTTTGAGATCGTCAGCCAGCGCGTCGTTAACGATA  >  NC_000913/3586032‑3586281                                                                                                                                     |                                                                                                                      tttAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTATCGCCCt                                                                                                                >  1:378932/1‑140 (MQ=255)        aTCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTATCGCCCTTTTTCAg                                                                                                         >  2:211702/1‑140 (MQ=255)         tCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTATCGCCCTTTTTCAGc                                                                                                        <  1:211702/140‑1 (MQ=255)             gCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTATCGCCCTTTTTCAGCGGCt                                                                                                    >  1:187405/1‑140 (MQ=255)             gCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTATCGCCCTTTTTCAGCGGCt                                                                                                    >  2:135357/1‑140 (MQ=255)              ccACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTATCGCCCTTTTTCAGCGGCTc                                                                                                   <  1:135357/140‑1 (MQ=255)                 ccGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTATCGCCCTTTTTCAGCGGCTCGcc                                                                                                <  2:378932/140‑1 (MQ=255)                  cGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTATCGCCCTTTTTCAGCGGCTCGccc                                                                                               >  1:203185/1‑140 (MQ=255)                   gTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTATCGCCCTTTTTCAGCGGCTCGCCct                                                                                              <  2:203185/140‑1 (MQ=255)                       ttCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTATCGCCCTTTTTCAGc                                                                                                        >  1:245786/1‑126 (MQ=255)                       ttCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTATCGCCCTTTTTCAGc                                                                                                        <  2:245786/126‑1 (MQ=255)                       ttCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTATCGCCCTTTTTCAGCGGCTCGCCCTCttt                                                                                          >  1:351419/1‑140 (MQ=255)                                          tCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTATCGCCCTTTTTCAGCGGCTCGCCCTCTTTCCAGAAGATAGCTTTACTg                                                                       >  2:385251/1‑140 (MQ=255)                                                        gTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTATCGCCCTTTTTCAGCGGCTCGCCCTCTTTCCAGAAGATAGCTTTACTGTTTTCGTGATTCgg                                                         >  2:223709/1‑140 (MQ=255)                                                                ataGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTATCGCCCTTTTTCAGCGGCTCGCCCTCTTTCCAGAAGATAGCTTTACTGTTTTCGTGATTCGGCAACACtt                                                 <  2:318027/140‑1 (MQ=255)                                                                  aGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTATCGCCCTTTTTCAGCGGCTCGCCCTCTTTCCAGAAGATAGCTTTACTGTTTTCGTGATTCGGCAACACTTCg                                               <  2:187405/140‑1 (MQ=255)                                                                      ttCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTATCGCCCTTTTTCAGCGGCTCGCCCTCTTTCCAGAAGATAGCTTTACTGTTTTCGTGATTCGGCAACACTTCGCTAc                                           >  2:111368/1‑140 (MQ=255)                                                                       tCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTATCGCCCTTTTTCAGCGGCTCGCCCTCTTTCCAGAAGATAGCTTTACTGTTTTCGTGATTCGGCAACACTTCGCTAcc                                          >  2:245986/1‑140 (MQ=255)                                                                        cGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTATCGCCCTTTTTCAGCGGCTCGCCCTCTTTCCAGAAGATAGCTTTACTGTTTTCGTGATTCGGCAACACTTCGCTACCg                                         <  1:385251/140‑1 (MQ=255)                                                                        cGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTATCGCCCTTTTTCAGCGGCTCGCCCTCTTTCCAGAAGATAGCTTTACTGTTTTCGTGATTCGGCAACACTTCGCTACCg                                         <  1:111368/140‑1 (MQ=255)                                                                         gTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTATCGCCCTTTTTCAGCGGCTCGCCCTCTTTCCAGAAGATAGCTTTACTGTTTTCGTGATTCGGCAACACTTCGCTACCGt                                        <  2:351419/140‑1 (MQ=255)                                                                                                 ctcCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTATCGCCCTTTTTCAGCGGCTCGCCCTCTTTCCAGAAGATAGCTTTACTGTTTTCGTGATTCGGCAACACTTCGCTACCGTAGGTTTTGAGATCGTCAGCCAgcg                <  1:11479/140‑1 (MQ=255)                                                                                                    cAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTATCGCCCTTTTTCAGCGGCTCGCCCTCTTTCCAGAAGATAGCTTTACTGTTTTCGTGATTCGGCAACACTTCGCTACCGTAGGTTTTGAGATCGTCAGCCAGCGcgt             >  1:232823/1‑140 (MQ=255)                                                                                                    cAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTATCGCCCTTTTTCAGCGGCTCGCCCTCTTTCCAGAAGATAGCTTTACTGTTTTCGTGATTCGGCAACACTTCGCTACCGTAGGTTTTGAGATCGTCAGCCAGCGcgt             >  1:118409/1‑140 (MQ=255)                                                                                                                   gTTCGCCTGCACCAGCGTATCGCCCTTTTTCAGCGGCTCGCCCTCTTTCCAGAAGATAGCTTTACTGTTTTCGTGATTCGGCAACACTTCGCTACCGTAGGTTTTGAGATCGTCAGCCAGCGCGTCGTTAACGATa  <  1:68189/136‑1 (MQ=255)                                                                                                                   gTTCGCCTGCACCAGCGTATCGCCCTTTTTCAGCGGCTCGCCCTCTTTCCAGAAGATAGCTTTACTGTTTTCGTGATTCGGCAACACTTCGCTACCGTAGGTTTTGAGATCGTCAGCCAGCGCGTCGTTAACGATa  >  2:68189/1‑136 (MQ=255)                                                                                                                                     |                                                                                                                      TTTAGTGATCAAGCCACCGTTTTTCTGCATCTCCTGGGCGATCTGTTCCGCAATCGTGCCTTTATAGAATTCGTCCGGGCCGTTTTCAGCAATCATCTCCAGGCTCTTTGCCAGGTTCGCCTGCACCAGCGTGTCGCCCTTTTTCAGCGGCTCGCCCTCTTTCCAGAAGATAGCTTTACTGTTTTCGTGATTCGGCAACACTTCGCTACCGTAGGTTTTGAGATCGTCAGCCAGCGCGTCGTTAACGATA  >  NC_000913/3586032‑3586281

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