Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	1,352,427	T→C	N101D (AAC→GAC)	thiD ←	bifunctional hydroxy‑methylpyrimidine kinase and hydroxy‑phosphomethylpyrimidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,352,427	0	T	C	100.0%	143.7 / NA	40	N101D (AAC→GAC)	thiD	bifunctional hydroxy‑methylpyrimidine kinase and hydroxy‑phosphomethylpyrimidine kinase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (2/38);  total (2/38)

AACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAG  >  minE/1352362‑1352487                                                                  |                                                             aaCCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGa                                                         <  1:277766/71‑1 (MQ=255) aaCCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGa                                                         <  1:2007355/71‑1 (MQ=255) aaCCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGa                                                         <  1:1600966/71‑1 (MQ=255)         aaGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGATCTGATaa                                                 <  1:38104/71‑1 (MQ=255)         aaGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGATCTGATaa                                                 <  1:3231118/71‑1 (MQ=255)         aaGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGATCTGATaa                                                 <  1:3037675/71‑1 (MQ=255)         aaGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGATCTGATaa                                                 <  1:154918/71‑1 (MQ=255)         aaGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGATCTGATaa                                                 <  1:2884738/71‑1 (MQ=255)         aaGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGATCTGATaa                                                 <  1:157138/71‑1 (MQ=255)         aaGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGATCTGATaa                                                 <  1:853346/71‑1 (MQ=255)         aaGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGATCTGATaa                                                 <  1:1677827/71‑1 (MQ=255)         aaGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGATCTGATaa                                                 <  1:2620317/71‑1 (MQ=255)         aaGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGATCTGATaa                                                 <  1:1681583/71‑1 (MQ=255)         aaGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGATCTGATaa                                                 <  1:2395998/71‑1 (MQ=255)         aaGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGATCTGATaa                                                 <  1:1712273/71‑1 (MQ=255)         aaGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGATCTGATaa                                                 <  1:914043/71‑1 (MQ=255)          aGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGATCTGATaa                                                 <  1:101450/70‑1 (MQ=255)             tGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGATCTGATaa                                                 <  1:2241090/67‑1 (MQ=255)                   cagcGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGATCTGATaa                                                 <  1:1477014/61‑1 (MQ=255)                    agcGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGATCTGATaa                                                 <  1:2312550/60‑1 (MQ=255)                                                  gtgtCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCa                              <  1:1333650/49‑1 (MQ=255)                                                   tgtCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTcc       <  1:276749/71‑1 (MQ=255)                                                   tgtCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTcc       <  1:791820/71‑1 (MQ=255)                                                   tgtCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTcc       <  1:759664/71‑1 (MQ=255)                                                   tgtCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTcc       <  1:1222604/71‑1 (MQ=255)                                                   tgtCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTcc       <  1:3111058/71‑1 (MQ=255)                                                   tgtCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTcc       <  1:3011189/71‑1 (MQ=255)                                                   tgtCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTcc       <  1:1567522/71‑1 (MQ=255)                                                   tgtCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTcc       <  1:1919200/71‑1 (MQ=255)                                                   tgtCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTcc       <  1:1214026/71‑1 (MQ=255)                                                   tgtCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTcc       <  1:2746160/71‑1 (MQ=255)                                                   tgtCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTcc       <  1:1680197/71‑1 (MQ=255)                                                   tgtCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTcc       <  1:2569365/71‑1 (MQ=255)                                                   tgtCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTcc       <  1:2457281/71‑1 (MQ=255)                                                   tgtCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTcc       <  1:1698834/71‑1 (MQ=255)                                                   tgtCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTcc       <  1:1757667/71‑1 (MQ=255)                                                   tgtCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTcc       <  1:1041021/71‑1 (MQ=255)                                                        aGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTcc       <  1:1600317/66‑1 (MQ=255)                                                         gTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAg  >  1:2752969/1‑70 (MQ=255)                                                         gTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAg  >  1:2158193/1‑70 (MQ=255)                                                                  |                                                             AACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAG  >  minE/1352362‑1352487

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