Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,585,538	T→C	N429S (AAC→AGC)	ggt ←	glutathione hydrolase proenzyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,585,538	0	T	C	100.0%	161.8 / NA	41	N429S (AAC→AGC)	ggt	glutathione hydrolase proenzyme
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/25);  total (16/25)

TTACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGGCGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAATGGGTAGTTTGATTACTCTCATA  >  NC_000913/3585403‑3585668                                                                                                                                        |                                                                                                                                   ttACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGGCGGAGAAATCATCCATCTGGctatt                                                                                                                                <  2:319935/140‑1 (MQ=255) ttACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGGCGGAGAAATCATCCATCTGGctatt                                                                                                                                <  2:64375/140‑1 (MQ=255) ttACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGGCGGAGAAATCATCCATCTGGctatt                                                                                                                                <  1:534571/140‑1 (MQ=255)   aCCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGGCGGAGAAATCATCCATCTGGCTATTaa                                                                                                                              >  1:549695/1‑140 (MQ=255)        cTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGGCGGAGAAATCATCCATCTGGctatt                                                                                                                                >  2:422995/1‑133 (MQ=255)        cTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGGCGGAGAAATCATCCATCTGGctatt                                                                                                                                <  1:422995/133‑1 (MQ=255)         tttCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGGCGGAGAAATCATCCATCTGGCTATTAAGCAGaa                                                                                                                        >  1:495953/1‑140 (MQ=255)           tCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGGCGGAGAAATCATCCATCTGGCTATTAAGCAGAATa                                                                                                                      <  1:587117/140‑1 (MQ=255)                 caATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGGCGGAGAAATCATCCATCTGGCTATTAAGCAGAATACCGctc                                                                                                                <  1:186323/140‑1 (MQ=255)                  aaTGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGGCGGAGAAATCATCCATCTGGCTATTAAGCAGAATACCGctct                                                                                                               >  2:424035/1‑140 (MQ=255)                   aTGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGGCGGAGAAATCATCCATCTGGCTATTAAGCAGAATACCGCtctc                                                                                                              >  2:98202/1‑140 (MQ=255)                   aTGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGGCGGAGAAATCATCCATCTGGCTATTAAGCAGAATACCGCtctc                                                                                                              >  1:567827/1‑140 (MQ=255)                    tGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGGCGGAGAAATCATCCATCTGGCTATTAAGCAGAATACCGCTCTCg                                                                                                             >  1:581638/1‑140 (MQ=255)                    tGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGGCGGAGAAATCATCCATCTGGCTATTAAGCAGAATACCGCTCTCg                                                                                                             >  1:65762/1‑140 (MQ=255)                      gTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGGCGGAGAAATCATCCATCTGGCTATTAAGCAGAATACCGCTCTCGcc                                                                                                           <  1:344242/140‑1 (MQ=255)                       tCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGGCGGAGAAATCATCCATCTGGCTATTAAGCAGAATACCGCTCTCGccc                                                                                                          <  2:73169/140‑1 (MQ=255)                        cGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGGCGGAGAAATCATCCATCTGGCTATTAAGCAGAATACCGCTCTCGCCcg                                                                                                         <  2:581638/140‑1 (MQ=255)                         ggCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGGCGGAGAAATCATCCATCTGGCTATTAAGCAGAATACCGCTCTCGCCcgc                                                                                                        <  2:193564/140‑1 (MQ=255)                         ggCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGGCGGAGAAATCATCCATCTGGCTATTAAGCAGAATACCGCTCTCGCCcga                                                                                                        >  1:575838/1‑139 (MQ=255)                                tCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGGCGGAGAAATCATCCATCTGGCTATTAAGCAGAATACCGCTCTCGCCCGCGACAATg                                                                                                 >  1:581759/1‑140 (MQ=255)                                     gacAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGGCGGAGAAATCATCCATCTGGCTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGt                                                                                            <  1:424035/140‑1 (MQ=255)                                          cGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGGCGGAGAAATCATCCATCTGGCTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGa                                                                                       <  1:121998/140‑1 (MQ=255)                                             gCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGGCGGAGAAATCATCCATCTGGCTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAgg                                                                                    <  2:495953/140‑1 (MQ=255)                                                 ttGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGGCGGAGAAATCATCCATCTGGCTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAggtggt                                                                                <  2:65762/140‑1 (MQ=255)                                                             aCGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGGCGGAGAAATCATCCATCTGGCTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGtata                                                                    >  1:500456/1‑140 (MQ=255)                                                                  gTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGGCGGAGAAATCATCCATCTGGCTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTca                                                               <  1:418106/140‑1 (MQ=255)                                                                           aCCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGGCGGAGAAATCATCCATCTGGCTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTcaccgccaccg                                                      <  2:549695/140‑1 (MQ=255)                                                                           aCCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGGCGGAGAAATCATCCATCTGGCTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTcaccgcca                                                         <  2:31250/137‑1 (MQ=255)                                                                           aCCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGGCGGAGAAATCATCCATCTGGCTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTcaccgcca                                                         >  1:31250/1‑137 (MQ=255)                                                                                ccaccaGCCCGTAAACGTTCGGTACGCCCGGTTTGGCGGAGAAATCATCCATCTGGCTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTa                                                 >  2:597369/1‑140 (MQ=255)                                                                                   ccaGCCCGTAAACGTTCGGTACGCCCGGTTTGGCGGAGAAATCATCCATCTGGCTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCa                                              >  1:190816/1‑140 (MQ=255)                                                                                      gCCCGTAAACGTTCGGTACGCCCGGTTTGGCGGAGAAATCATCCATCTGGCTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCt                                           <  2:495318/140‑1 (MQ=255)                                                                                         cGTAAACGTTCGGTACGCCCGGTTTGGCGGAGAAATCATCCATCTGGCTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTa                                        <  2:523890/140‑1 (MQ=255)                                                                                             aaCGTTCGGTACGCCCGGTTTGGTGGAGAAATCATCCATCTGGCTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATcca                                    <  2:537086/140‑1 (MQ=255)                                                                                                 ttCGGTACGCCCGGTTTGGCGGAGAAATCATCCATCTGGCTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCt                                           <  2:177021/129‑1 (MQ=255)                                                                                                 ttCGGTACGCCCGGTTTGGCGGAGAAATCATCCATCTGGCTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCt                                           >  1:177021/1‑129 (MQ=255)                                                                                                  tCGGTACGCCCGGTTTGGCGGAGAAATCATCCATCTGGCTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACt                               <  2:500456/140‑1 (MQ=255)                                                                                                  tCGGTACGCCCGGTTTGGCGGAGAAATCATCCATCTGGCTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACt                               <  2:570511/140‑1 (MQ=255)                                                                                                         gCCCGGTTTGGCGGAGAAATCATCCATCTGGCTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAAt                        <  1:98202/140‑1 (MQ=255)                                                                                                          cccGGTTTGGCGGAGAAATCATCCATCTGGCTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAATg                       >  2:429523/1‑140 (MQ=255)                                                                                                            cGGTTTGGCGGAGAAATCATCCATCTGGCTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAATgtg                     <  2:140477/140‑3 (MQ=255)                                                                                                            cGGTTTGGCGGAGAAATCATCCATCTGGCTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAATggg                     <  1:597369/140‑1 (MQ=255)                                                                                                             ggTTTGGCGGAGAAATCATCCATCTGGCTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAATGGGt                    <  2:190816/140‑1 (MQ=255)                                                                                                              gTTTGGCGGAGAAATCATCCATCTGGCTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAATGGGTa                   >  1:67938/1‑140 (MQ=255)                                                                                                                     ggAGAAATCATCCATCTGGCTATTAAGCAGAATACCGCTCTCGCCCGAGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAATGGGTAGTTTGAt            <  2:575838/140‑1 (MQ=255)                                                                                                                               tccatcTGGCTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAATGGGTAGTTTGATTACTCTCATa  <  2:67938/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                   TTACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGGCGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAATGGGTAGTTTGATTACTCTCATA  >  NC_000913/3585403‑3585668

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: ‑