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%	172.5 / NA	50	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 (25/25);  total (25/25)

AGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCG  >  minE/1352359‑1352494                                                                     |                                                                    aGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGtctt                                                                   >  1:2642700/1‑71 (MQ=255)   cAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTg                                                                  >  1:2088795/1‑70 (MQ=255)    aaCCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGa                                                                <  1:2146560/71‑1 (MQ=255)     aCCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGAt                                                               >  1:915319/1‑71 (MQ=255)       cgccgAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGATCt                                                             <  1:1622140/71‑1 (MQ=255)       cgccgAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGATCt                                                             <  1:2264467/71‑1 (MQ=255)          cgAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGATCTGa                                                           >  1:2593400/1‑70 (MQ=255)            aaGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGATCTGATaa                                                        <  1:1719862/71‑1 (MQ=255)             aGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGATCTGATAAc                                                       <  1:978715/71‑1 (MQ=255)              ggTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGATCTGATAACg                                                      <  1:546499/71‑1 (MQ=255)              ggTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGATCTGATAACg                                                      <  1:2301160/71‑1 (MQ=255)               gTGAAAGCAGCGGGTCGCCGATTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGATCTGATAACg                                                      <  1:608610/70‑1 (MQ=255)                tGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTgg                                                                 <  1:2711479/58‑1 (MQ=255)                tGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTgg                                                                 <  1:713594/58‑1 (MQ=255)                  aaaGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGAt                                                               >  1:2206034/1‑58 (MQ=255)                   aaGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGATCTGATAACGTTGc                                                  >  1:1287515/1‑70 (MQ=255)                       agcGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGAt                                                               >  1:568519/1‑53 (MQ=255)                            gTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGtt                                           <  1:1325007/68‑1 (MQ=255)                              cgccgcTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTg                                        >  1:635688/1‑69 (MQ=255)                              cgccgcTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGc                                       <  1:1422508/70‑1 (MQ=255)                               gccgcTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTc                                          <  1:2836119/66‑1 (MQ=255)                                 cgcTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTgg                                                                 <  1:2370449/41‑1 (MQ=255)                                  gcTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGATCTGATaa                                                        >  1:2866316/1‑49 (MQ=255)                                  gcTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCt                                         >  1:2005480/1‑64 (MQ=255)                                  gcTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCt                                         >  1:981706/1‑64 (MQ=255)                                   cTTTTTGCCAGCATAACGGTGTCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGc                                 >  1:984149/1‑71 (MQ=255)                                      tttGCCAGCATAACGGTGTCGAGTACCACGTCTTGGATCTGATAACg                                                      >  1:1522359/1‑47 (MQ=255)                                       ttGCCAGCATAACGGTGTCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTc                              <  1:2220650/70‑1 (MQ=255)                                           cAGCATAACGGTGTCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGt                                            <  1:11605/52‑1 (MQ=255)                                            aGCATAACGGTGTCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTc                              >  1:2896642/1‑65 (MQ=255)                                             gCATAACGGTGTCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAata                       <  1:1374516/71‑1 (MQ=255)                                              cATAACGGTGTCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAatat                      >  1:706159/1‑71 (MQ=255)                                               aTAACGGTGTCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATc                     <  1:1084974/71‑1 (MQ=255)                                               aTAACGGTGTCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATc                     >  1:738518/1‑71 (MQ=255)                                                 aaCGGTGTCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCgg                   >  1:1751412/1‑71 (MQ=255)                                                 aaCGGTGTCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCgg                   >  1:1134124/1‑71 (MQ=255)                                                 aaCGGTGTCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCgg                   >  1:875671/1‑71 (MQ=255)                                                  aCGGTGTCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCg                                  >  1:291342/1‑55 (MQ=255)                                                   cGGTGTCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCAc                                    >  1:221523/1‑52 (MQ=255)                                                   cGGTGTCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCAc                                    >  1:250542/1‑52 (MQ=255)                                                    ggTGTCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGttt                <  1:3040928/71‑1 (MQ=255)                                                     gtgtCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAatat                      <  1:1438109/64‑1 (MQ=255)                                                       gtCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTc                                          >  1:1249111/1‑42 (MQ=255)                                                        tCGAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGc            >  1:1411745/1‑71 (MQ=255)                                                          gAGTACCACGTCTTGGATCTGATAACGTTGCAACCGtt                                           >  1:2571168/1‑38 (MQ=255)                                                          gAGTACCACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGcc           <  1:90317/70‑1 (MQ=255)                                                               ccACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCAt      <  1:1280467/70‑1 (MQ=255)                                                               ccACGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCAt      <  1:1081295/70‑1 (MQ=255)                                                                cACGTCTTGGATCTGATAACGTTGCAACCGGTCTGCCACCGCTTCAACAata                       <  1:1077620/52‑1 (MQ=255)                                                                  cGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATAcc   >  1:148547/1‑70 (MQ=255)                                                                  cGTCTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCg  <  1:2917974/71‑1 (MQ=255)                                                                     |                                                                    AGCAACCGCCGAAGGTGAAAGCAGCGGGTCGCCGCTTTTTGCCAGCATAACGGTGTCGAGTACCACGTTTTGGATCTGATAACGTTGCAACCGTTCTGCCACCGCTTCAACAATATCGGTTTCCGCCAGCATACCG  >  minE/1352359‑1352494

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