Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,971,608	C→A	H4Q (CAC→CAA)	tas →	putative NADP(H)‑dependent aldo‑keto reductase Tas

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,971,608	0	C	A	100.0%	108.5 / NA	32	H4Q (CAC→CAA)	tas	putative NADP(H)‑dependent aldo‑keto reductase Tas
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (13/19);  total (13/19)

CAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCACCGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCTTGGCACGATGACGTTTGGTGAACAGAACAGCGAAGCCGACGCCCACGCACAACTCGACTATGCCGTCGCTCA  >  NC_000913/2971472‑2971721                                                                                                                                         |                                                                                                                  cAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCAACGt                                                                                                                <  2:43820/140‑1 (MQ=255)  aaaTTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCAACGta                                                                                                               <  2:543031/140‑1 (MQ=255)                tAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCAACGTATACCCCACAGTTCg                                                                                                 >  1:52845/1‑140 (MQ=255)                tAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCAACGTATACCCCACAGTTCg                                                                                                 >  1:762324/1‑140 (MQ=255)                     aaTAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCAACGTATACCCCACAGTTCGCTGGa                                                                                            >  2:284450/1‑140 (MQ=255)                          ggTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCAACGTATACCCCACAGtt                                                                                                   <  1:275399/128‑1 (MQ=255)                          ggTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCAACGTATACCCCACAGtt                                                                                                   >  2:275399/1‑128 (MQ=255)                                cATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCAACGTATACCCCACAGTTCGCTGGAAGTCa                                                                                       >  2:518756/1‑134 (MQ=255)                                cATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCAACGTATACCCCACAGTTCGCTGGAAGTCa                                                                                       <  1:518756/134‑1 (MQ=255)                                cATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCAACGTATACCCCACAGTTCGCTGGAAGTCAGCACGc                                                                                 >  1:139370/1‑140 (MQ=255)                                     ggCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCAACGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTgggg                                                                            >  1:700865/1‑140 (MQ=255)                                       cTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCAACGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCt                                                                          <  2:700865/140‑1 (MQ=255)                                           ccTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCAACGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCTTGGc                                                                      <  1:376750/140‑1 (MQ=255)                                           ccTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCAACGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCTTGGc                                                                      >  1:775414/1‑140 (MQ=255)                                                aGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCAACGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCTTGGCACGAt                                                                 <  1:158048/140‑1 (MQ=255)                                                    ttttCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCAACGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCTTGGCACGATGACg                                                             >  2:641086/1‑140 (MQ=255)                                                             tcccttttcccttCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAGGGAAACGGCTATGCAATATCAACGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCTTGGCACGATGACGTTTGGTgaa                                                    <  1:745418/140‑1 (MQ=255)                                                                   ttcccttcccTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCAACGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCTTGGCACGATGACGTTTGGTgaacagaac                                              >  2:537670/1‑140 (MQ=255)                                                                              ctGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCAACGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCTTGGCACGATGACGTTTGGTGAACAGAACAGCGAAGCCGa                                   <  1:184347/140‑1 (MQ=255)                                                                                      tttATATTCGTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCAACGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCTTGGCACGATGACGTTTGGTGAACAGAACAGCGAAGCCGACGCCCACg                           <  1:677404/140‑1 (MQ=255)                                                                                      tttATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCAACGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCTTGGCACGATGACGTTTGGTGAACAGAACAGCGAAGCCGACGCCCACg                           <  2:139053/140‑1 (MQ=255)                                                                                      tttATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCAACGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCTTGGCACGATGACGTTTGGTGAACAGAACAGCGAAGCCGACGCCCACg                           <  2:701969/140‑1 (MQ=255)                                                                                      tttATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCAACGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCTTGGCACGATGACGTTTGGTGAACAGAACAGCGAAGCCGACGCCCACg                           <  2:374135/140‑1 (MQ=255)                                                                                      tttATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCAACGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCTTGACACGATGACGTTTGGTGAACAGAACAGCGAAGCCGACGCCCACg                           <  2:166824/140‑1 (MQ=255)                                                                                           atTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCAACGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCTTGGCACGATGACGTTTGGTGAACAGAACAGCGAAGCCGACGCCCACGCACaa                      <  2:610497/140‑1 (MQ=255)                                                                                                 tATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCAACGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCTTGGCACGATGACGTTTGGTGAACAGAACAGCGAAGCCGACGCCCACGCACAACTCGAc                <  1:641086/140‑1 (MQ=255)                                                                                                       gTGATTATAAAAAGGAAACGGCTATGCAATATCAACGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCTTGGCACGATGACGTTTGGTGAACAGAACAGCGAAGCCGACGCCCACGCACAACTCGACTATGcc          >  2:13364/1‑140 (MQ=255)                                                                                                        tGATTATAAAAAGGAAACGGCTATGCAATATCAACGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCTTGGCACGATGACGTTTGGTGAACAGAACAGCGAAGCCGACGCCCACGCACAACTCGACTATGCCg         >  1:82813/1‑140 (MQ=255)                                                                                                        tGATTATAAAAAGGAAACGGCTATGCAATATCAACGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCTTGGCACGATGACGTTTGGTGAACAGAACAGCGAAGCCGACGCCCACGCACAACTCGACTATGCCg         <  2:762324/140‑1 (MQ=255)                                                                                                          aTTATAAAAAGGAAACGGCTATGCAATATCAACGTATACCCCACAGTTCGCTGGAAGTCAGCACGGTGGGGCTTGGCACGATGACGTTTGGTGAACAGAACAGCGAAGCCGACGCCCACGCACAACTCGACTATGCCGTc       >  2:351056/1‑140 (MQ=255)                                                                                                               aaaaaGGAAACGGCTATGCAATATCAACGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCTTGGCACGATGACGTTTGGTGAACAGAACAGCGAAGCCGACGCCCACGCACAACTCGACTATGCCGTCGCTCa  <  2:775414/140‑1 (MQ=255)                                                                                                               aaaaaGGAAACGGCTATGCAATATCAACGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCTTGGCACGATGACGTTTGGTGAACAGAACAGCGAAGCCGACGCCCACGCACAACTCGACTATGCCGTCGCTCa  <  2:82813/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                  CAAATTATTGAACGTTAACAAATAAGGTCAGCATCCGGCTGGCCTTAAGATTTTTCTCTTTCCCTTTTCCCTTCCCTCTGCCATTTTTATATTCCTTATGTCGTGATTATAAAAAGGAAACGGCTATGCAATATCACCGTATACCCCACAGTTCGCTGGAAGTCAGCACGCTGGGGCTTGGCACGATGACGTTTGGTGAACAGAACAGCGAAGCCGACGCCCACGCACAACTCGACTATGCCGTCGCTCA  >  NC_000913/2971472‑2971721

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