Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,971,955	G→A	R120H (CGC→CAC)	tas →	putative NADP(H)‑dependent aldo‑keto reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,971,955	0	G	A	100.0%	121.4 / NA	35	R120H (CGC→CAC)	tas	putative NADP(H)‑dependent aldo‑keto reductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/20);  total (15/20)

GAGCGGACCGTCGCGCAATAATGACAAGGGCATCCGCCCGGATCAGGCGCTGGATCGGAAGAATATCCGCGAAGCGCTGCATGACAGCCTCAAGCGCCTACAGACTGATTACCTCGATCTTTATCAGGTGCACTGGCCGCAGCGCCCGACCAACTGCTTCGGCAAACTCGGTTATAGCTGGACGGATT  >  NC_000913/2971860‑2972047                                                                                                |                                                                                             gAGCGGACCGTCGCGCAATAATGACAAGGGCATCCGCCCGGATCAGGCGCTGGATCGGAAGAATATCCGCGAAGCGCTGCATGACAGCCTCAAGCACCTa                                                                                          <  2:1078981/100‑1 (MQ=255)  aGCGGACCGTCGCGCAATAATGACAAGGGCATCCGCCCGGATCAGGCGCTGGATCGGAAGAATATCCGCGAAGCGCTGCATGACAGCCTCAAGCACCTAc                                                                                         <  1:222882/100‑1 (MQ=255)   gCGGACCGTCGCGCAATAATGACAAGGGCATCCGCCCGGATCAGGCGCTGGATCGGAAGAATATCCGCGAAGCGCTGCATGACAGCCTCAAGCACCTACa                                                                                        <  1:384730/100‑1 (MQ=255)      gACCGTCGCGCAATAATGACAAGGGCATCCGCCCGGATCAGGCGCTGGATCGGAAGAATATCCGCGAAGCGCTGCATGACAGCCTCAAGCACCTACAGAc                                                                                     <  2:229566/100‑1 (MQ=255)        ccGTCGCGCAATAATGACAAGGGCATCCGCCCGGATCAGGCGCTGGATCGGAAGAATATCCGCGAAGCGCTGCATGACAGCCTCAAGCACCTACAGACTg                                                                                   >  2:606410/1‑100 (MQ=255)         cGTCGCGCAATAATGACAAGGGCATCCGCCCGGATCAGGCGCTGGATCGGAAGAATATCCGCGAAGCGCTGCATGACAGCCTCAAGCACCTACAGACTGa                                                                                  <  1:72390/100‑1 (MQ=255)          gTCGCGCAATAATGACAAGGGCATCCGCCCGGATCAGGCGCTGGATCGGAAGAGTATCCGCGAAGCGCTGCATGACAGCCTCAAGCACCTACAGACTGAt                                                                                 >  1:785122/1‑100 (MQ=255)             gcgcAATAATGACAAGGGCATCCGCCCGGATCAGGCGCTGGATCGGAAGAATATCCGCGAAGCGCTGCATGACAGCCTCAAGCACCTACAGACTGATTAc                                                                              <  1:405097/100‑1 (MQ=255)                cAATAATGACAAGGGCATCCGCCCGGATCAGGCGCTGGATCGGAAGAATATCCGCGAAGCGCTGCATGACAGCCTCAAGCACCTACAGACTGATTACCTc                                                                           >  2:438545/1‑100 (MQ=255)                 aataatGACAAGGGCATCCGCCCGGATCAGGCGCTGGATCGGAAGAATATCCGCGAAGCGCTGCATGACAGCCTCAAGCACCTACAGACTGATTACCTCg                                                                          <  2:1071037/100‑1 (MQ=255)                      tGACAAGGGCATCCGCCCGGATCAGGCGCTGGATCGGAAGAATATCCGCGAAGCGCTGCATGACAGCCTCAAGCACCTACAGACTGATTACCTCGATCtt                                                                     <  2:947092/100‑1 (MQ=255)                        aCAAGGGCATCCGCCCGGATCAGGCGCTGGATCGGAAGAATATCCGCGAAGCGCTGCATGACAGCCTCAAGCACCTACAGACTGATTACCTCGATCTTTa                                                                   <  1:1143442/100‑1 (MQ=255)                        aCAAGGGCATCCGCCCGGATCAGGCGCTGGATCGGAAGAATATCCGCGAAGCGCTGCATGACAGCCTCAAGCACCTACAGACTGATTACCTCGATCTTTa                                                                   <  2:408512/100‑1 (MQ=255)                          aaGGGCATCCGCCCGGATCAGGCGCTGGATCGGAAGAATATCCGCGAAGCGCTGCATGACAGCCTCAAGCACCTACAGACTGATTACCTCGATCTTTATc                                                                 <  1:431543/100‑1 (MQ=255)                            gggCATCCGCCCGGATCAGGCGCTGGATCGGAAGAATATCCGCGAAGCGCTGCATGACAGCCTCAAGCACCTACAGACTGATTACCTCGATCTTTATCAg                                                               >  1:315764/1‑100 (MQ=255)                            gggCATCCGCCCGGATCAGGCGCTGGATCGGAAGAATATCCGCGAAGCGCTGCATGACAGCCTCAAGCACCTACAGACTGATTACCTCGATCTTTATCAg                                                               <  2:444628/100‑1 (MQ=255)                            gggCATCCGCCCGGATCAGGCGCTGGATCGGAAGAATATCCGCGAAGCGCTGCATGACAGCCTCAAGCACCTACAGACTGATTACCTCGATCTTTATCAg                                                               >  1:679065/1‑100 (MQ=255)                                aTCCGCCCGGATCAGGCGCTGGATCGGAAGAATATCCGCGAAGCGCTGCATGACAGCCTCAAGCACCTACAGACTGATTACCTCGATCTTTATCAGGTGc                                                           <  2:790284/100‑1 (MQ=255)                                   cGCCCGGATCAGGCGCTGGATCGGAAGAATATCCGCGAAGCGCTGCATGACAGCCTCAAGCACCTACAGACTGATTACCTCGATCTTTATCAGGTGCACt                                                        >  1:60281/1‑100 (MQ=255)                                    gCCCGGATCAGGCGCTGGATCGGAAGAATATCCGCGAAGCGCTGCATGACAGCCTCAAGCACCTACAGACTGATTACCTCGATCTTTATCAGGTGCACTg                                                       >  1:520142/1‑100 (MQ=255)                                                 gcTGGATCGGAAGAATATCCGCGAAGCGCTGCATGACAGCCTCAAGCACCTACAGACTGATTACCTCGATCTTTATCAGGTGCACTGGCCGCAGCGCCCg                                          >  1:780621/1‑100 (MQ=255)                                                       tCGGAAGAATATCCGCGAAGCGCTGCATGACAGCCTCAAGCACCTACAGACTGATTACCTCGATCTTTATCAGGTGCACTGGCCGCAGCGCCCGACCAAc                                    <  2:659413/100‑1 (MQ=255)                                                       tCGGAAGAATATCCGCGAAGCGCTGCATGACAGCCTCAAGCACCTACAGACTGATTACCTCGATCTTTATCAGGTGCACTGGCCGCAGCGCCCGACCAAc                                    <  2:524774/100‑1 (MQ=255)                                                               atatCCGCGAAGCGCTGCATGACAGCCTCAAGCACCTACAGACTGATTACCTCGATCTTTATCAGGTGCACTGGCCGCAGCGCCCGACCAACTGCTTCgg                            <  2:315761/100‑1 (MQ=255)                                                                 atCCGCGAAGCGCTGCATGACAGCCTCAAGCACCTACAGACTGATTACCTCGATCTTTATCAGGTGCACTGGCCGCAGCGCCCGACCAACTGCTTCGGCa                          <  2:641575/100‑1 (MQ=255)                                                                 atCCGCGAAGCGCTGCATGACAGCCTCAAGCACCTACAGACTGATGACCTCGATCTTTATCAGGTGCACTGGCCGCAGCGCCCGACCAACTGCTTCGGCa                          <  1:606414/100‑1 (MQ=255)                                                                   ccGCGAAGCGCTGCATGACAGCCTCAAGCACCTACAGACTGATTACCTCGATCTTTATCAGGTGCACTGGCCGCAGCGCCCGACCAACTGCTTCGGCaaa                        >  2:956250/1‑100 (MQ=255)                                                                    cgcgAAGCGCTGCATGACAGCCTCAAGCACCTACAGACTGATTACCTCGATCTTTATCAGGTGCACTGGCCGCAGCGCCCGACCAACTGCTTCGGCAAAc                       >  2:974952/1‑100 (MQ=255)                                                                    cgcgAAGCGCTGCATGACAGCCTCAAGCACCTACAGACTGATTACCTCGATCTTTATCAGGTGCACTGGCCGCAGCGCCCGACCAACTGCTTCGGCAAAc                       >  2:275292/1‑100 (MQ=255)                                                                      cgAAGCGCTGCATGACAGCCTCAAGCACCTACAGACTGATTACCTCGATCTTTATCAGGTGCACTGGCCGCAGCGCCTGACCAACTGCTTCGGCAAACTc                     <  2:756752/100‑1 (MQ=255)                                                                          gcgcTGCATGACAGCCTCAAGCACCTACAGACTGATTACCTCGATCTTTATCAGGTGCACTGGCCGCAGCGCCCGACCAACTGCTTCGGCAAACTCGGtt                 >  1:407263/1‑100 (MQ=255)                                                                                cATGACAGCCTCAAGCACCTACAGACTGATTACCTCGATCTTTATCAGGTGCACTGGCCGCAGCGCCCGACCAACTGCTTCGGCAAACTCGGTTATAGCt           <  1:655087/100‑1 (MQ=255)                                                                                     cAGCCTCAAGCACCTACAGACTGATTACCTCGATCTTTATCAGGTGCACTGGCCGCAGCGCCCGACCAACTGCTTCGGCAAACTCGGTTATAGCTGGACg      >  2:590364/1‑100 (MQ=255)                                                                                       gCCTCAAGCACCTACAGACTGATTACCTCGATCTTTATCAGGTGCACTGGCCGCAGCGCCCGACCAACTGCTTCGGCAAACTCGGTTATAGCTGGACGGa    >  1:1270277/1‑100 (MQ=255)                                                                                         cTCAAGCACCTACAGACTGATTACCTCGATCTTTATCAGGTGCACTGGCCGCAGCGCCCGACCAACTGCTTCGGCAAACTCGGTTATAGCTGGACGGAtt  >  1:386876/1‑100 (MQ=255)                                                                                                |                                                                                             GAGCGGACCGTCGCGCAATAATGACAAGGGCATCCGCCCGGATCAGGCGCTGGATCGGAAGAATATCCGCGAAGCGCTGCATGACAGCCTCAAGCGCCTACAGACTGATTACCTCGATCTTTATCAGGTGCACTGGCCGCAGCGCCCGACCAACTGCTTCGGCAAACTCGGTTATAGCTGGACGGATT  >  NC_000913/2971860‑2972047

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