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%	165.7 / NA	45	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 (18/27);  total (18/27)

AACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACC  >  minE/1352362‑1352493                                                                  |                                                                   aaCCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTg                                                                 >  1:3706205/1‑69 (MQ=255) aaCCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGa                                                               >  1:724727/1‑71 (MQ=255)  aCCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGAt                                                              <  1:1017461/71‑1 (MQ=255)  aCCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGAt                                                              <  1:2934513/71‑1 (MQ=255)  aCCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGAt                                                              <  1:1642577/71‑1 (MQ=255)  aCCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGAt                                                              <  1:1975273/71‑1 (MQ=255)  aCCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGAt                                                              <  1:2867626/71‑1 (MQ=255)  aCCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGAt                                                              <  1:434155/71‑1 (MQ=255)  aCCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGAt                                                              <  1:2867554/71‑1 (MQ=255)  aCCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGAt                                                              <  1:2804186/71‑1 (MQ=255)            gTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATCACGGTGTCGAGTACCACGTCTTGGATCTGATAACGt                                                    >  1:62307/1‑71 (MQ=255)            gTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATATCGGTGTCGAGTACCACGTCTTGGATCTGATAACGt                                                    >  1:2421299/1‑71 (MQ=255)            gTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGATCTGATAACGt                                                    >  1:3719327/1‑71 (MQ=255)            gTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGATCTGATAACGt                                                    >  1:4011812/1‑71 (MQ=255)            gTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGATCTGATAACGt                                                    >  1:1433973/1‑71 (MQ=255)            gTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGATCTGATAACGt                                                    >  1:364126/1‑71 (MQ=255)            gTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGATCTGATAACGt                                                    >  1:4101569/1‑71 (MQ=255)            gTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGATCTGATAACGt                                                    >  1:52982/1‑71 (MQ=255)            gTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGATCTGATAACGt                                                    >  1:3411519/1‑71 (MQ=255)            gTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGATCTGATAACGt                                                    >  1:920225/1‑71 (MQ=255)            gTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGATCTGATAACGt                                                    >  1:3270539/1‑71 (MQ=255)                                        cAGCATAACGGTGTCGAGTACCACGTCTTGGAGCTGATAACGTTGCAACCGTTCTGCCACCGCTTcaacaa                        <  1:2898299/71‑1 (MQ=255)                                            aTAACGGTGTCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATc                    <  1:3491330/71‑1 (MQ=255)                                            aTAACGGTGTCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATc                    <  1:370647/71‑1 (MQ=255)                                            aTAACGGTGTCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATc                    <  1:2719415/71‑1 (MQ=255)                                            aTAACGGTGTCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATc                    <  1:3779709/71‑1 (MQ=255)                                            aTAACGGTGTCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATc                    <  1:2289413/71‑1 (MQ=255)                                            aTAACGGTGTCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATc                    <  1:2443918/71‑1 (MQ=255)                                            aTAACGGTGTCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATc                    <  1:4279277/71‑1 (MQ=255)                                                  gtgtCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTc              <  1:3297100/71‑1 (MQ=255)                                                  gtgtCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTc              <  1:895419/71‑1 (MQ=255)                                                  gtgtCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTc              <  1:876364/71‑1 (MQ=255)                                                  gtgtCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTc              <  1:1799846/71‑1 (MQ=255)                                                  gtgtCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTc              <  1:2163143/71‑1 (MQ=255)                                                  gtgtCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTc              <  1:440330/71‑1 (MQ=255)                                                  gtgtCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTc              <  1:3990000/71‑1 (MQ=255)                                                  gtgtCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTc              <  1:375265/71‑1 (MQ=255)                                                  gtgtCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTc              <  1:3663406/71‑1 (MQ=255)                                                  gtgtCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTc              <  1:2825438/71‑1 (MQ=255)                                                  gtgtCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTc              <  1:3063167/71‑1 (MQ=255)                                                              aCGTCTTGGATCTGATAACGTTGCAACCGTTCTGccc                                    >  1:2723030/1‑36 (MQ=38)                                                              aCGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATa    >  1:2362660/1‑69 (MQ=255)                                                              aCGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATAcc  >  1:4361849/1‑71 (MQ=255)                                                              aCGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATAc   >  1:2219708/1‑70 (MQ=255)                                                              aCGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATAc   >  1:1493291/1‑70 (MQ=255)                                                                  |                                                                   AACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACC  >  minE/1352362‑1352493

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