Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	2,829,664	T→C	N109S (AAC→AGC)	aspA ←	aspartate ammonia‑lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,829,664	0	T	C	100.0%	308.3 / NA	83	N109S (AAC→AGC)	aspA	aspartate ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (36/47);  total (36/47)

CTGATATTCACCTTTTTGGTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGTTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCGGGAACTGATCCATGCATTTTCCG  >  minE/2829600‑2829732                                                                 |                                                                     cTGATATTCACCTTTTTGGTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGt                                                                >  1:3068481/1‑71 (MQ=255)    atatTCACGTTTTTGGTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTa                                                             <  1:2012913/71‑1 (MQ=255)    atatTCACCTTTTTGGTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTa                                                             <  1:1513327/71‑1 (MQ=255)    atatTCACCTTTTTGGTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTa                                                             <  1:340889/71‑1 (MQ=255)      atTCACCTTTTTGGTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACg                                                           >  1:3141045/1‑71 (MQ=255)       ttCACCTTTTTGGTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACgg                                                          <  1:1217741/71‑1 (MQ=255)       ttCACCTTTTTGGTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACgg                                                          <  1:319947/71‑1 (MQ=255)        tCACCTTTTTGGTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACgg                                                          <  1:1572475/70‑1 (MQ=255)         cACCTTTTTGGTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACgg                                                          >  1:3157793/1‑69 (MQ=255)           ccTTTTTGGTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAg                                                       >  1:1322490/1‑70 (MQ=255)            cTTTTTGGTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGttt                                                              <  1:1974686/62‑1 (MQ=255)            cTTTTTGGTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGttt                                                              <  1:2251754/62‑1 (MQ=255)            cTTTTTGGTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGa                                                         <  1:1872385/67‑1 (MQ=255)            cTTTTTGGTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGa                                                         <  1:1400599/67‑1 (MQ=255)            cTTTTTGGTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGt                                                      <  1:2701745/70‑1 (MQ=255)                ttGGTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTc                                                                    >  1:2911359/1‑52 (MQ=255)                ttGGTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCt                                                  <  1:1996440/70‑1 (MQ=255)                ttGGTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCt                                                  <  1:172087/70‑1 (MQ=255)                ttGGTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTg                                                 >  1:33248/1‑71 (MQ=255)                     gACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGcgccg                                            >  1:7677/1‑71 (MQ=255)                      aCCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGcgccgc                                           >  1:3123335/1‑71 (MQ=255)                      aCCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGcgccgc                                           >  1:2127777/1‑71 (MQ=255)                      aCCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGcgccgc                                           >  1:3119587/1‑71 (MQ=255)                      aCCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGcgccgc                                           >  1:412791/1‑71 (MQ=255)                      aCCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGcgccg                                            <  1:3140268/70‑1 (MQ=255)                       cccATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCgccgcc                                          >  1:1653887/1‑71 (MQ=255)                        ccATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGccc                                         >  1:2838615/1‑71 (MQ=255)                        ccATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGccc                                         <  1:791732/71‑1 (MQ=255)                          aTCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTg                                       >  1:1669071/1‑71 (MQ=255)                          aTCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTg                                       >  1:1016121/1‑71 (MQ=255)                          aTCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTg                                       >  1:1076234/1‑71 (MQ=255)                          aTCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTg                                       >  1:2156435/1‑71 (MQ=255)                           tCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTgg                                      <  1:996951/71‑1 (MQ=255)                            cAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGcgcc                                             <  1:1653134/63‑1 (MQ=255)                             aGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGt                                     <  1:2099709/70‑1 (MQ=255)                              gTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTa                                    <  1:2499022/70‑1 (MQ=255)                               ttCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGa                                  <  1:2267080/71‑1 (MQ=255)                               ttCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGa                                  <  1:1234166/71‑1 (MQ=255)                                 ccAGACCGATATTGGCCAGCACTTCGTTTGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGAc                                 <  1:1666006/70‑1 (MQ=255)                                 ccAGACCGATATTGGCCAGCACTTCGTTGGTGCTCa                                                                   >  1:120781/1‑36 (MQ=255)                                 ccAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGaa                                                        >  1:1341961/1‑47 (MQ=255)                                 ccAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGaa                                                        >  1:1721104/1‑47 (MQ=255)                                 ccAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGAc                                 <  1:1329646/70‑1 (MQ=255)                                 ccAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACg                                >  1:2189162/1‑71 (MQ=255)                                 ccAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACg                                >  1:696271/1‑71 (MQ=255)                                  cAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTa                                                             <  1:2502048/41‑1 (MQ=255)                                  cAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTa                                                             <  1:804377/41‑1 (MQ=255)                                  cAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTa                                                             <  1:1431034/41‑1 (MQ=255)                                  cAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACg                                >  1:819822/1‑70 (MQ=255)                                    gACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCt                             <  1:2459397/71‑1 (MQ=255)                                     aCCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGccc                                         <  1:3116861/58‑1 (MQ=255)                                       cGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTgg                                      <  1:2616877/59‑1 (MQ=255)                                        gATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCg                         <  1:1815136/71‑1 (MQ=255)                                         atatTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTg                                       <  1:2459401/56‑1 (MQ=255)                                         atatTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCgg                        >  1:264858/1‑71 (MQ=255)                                           atTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGCCGTCTACCggg                       >  1:585482/1‑70 (MQ=255)                                            ttGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAg                                                       <  1:1691491/37‑1 (MQ=255)                                            ttGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAg                                                       <  1:513341/37‑1 (MQ=255)                                            ttGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCGGGaa                     >  1:2269744/1‑71 (MQ=255)                                               gCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCGGGAACTg                  <  1:2527172/71‑1 (MQ=255)                                                ccAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCGGGAACTGa                 <  1:1309705/71‑1 (MQ=255)                                                 cAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGa                                  <  1:705523/53‑1 (MQ=255)                                                  aGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGt                               <  1:3087045/55‑1 (MQ=255)                                                  aGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCGGGAACt                   >  1:1642122/1‑67 (MQ=255)                                                  aGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCGGGAACTGATc               >  1:387862/1‑71 (MQ=255)                                                     aCTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCGGGa                      <  1:2326936/61‑1 (MQ=255)                                                     aCTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCGGGa                      <  1:2369994/61‑1 (MQ=255)                                                     aCTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCGGGAACTGATCCAt            >  1:3179421/1‑71 (MQ=255)                                                     aCTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCGGGAACTGATCCAt            >  1:1479203/1‑71 (MQ=255)                                                       ttCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCg                         <  1:27838/56‑1 (MQ=255)                                                       ttCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCg                         <  1:2342774/56‑1 (MQ=255)                                                       ttCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCg                         <  1:2286134/56‑1 (MQ=255)                                                       ttCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCg                         <  1:1769161/56‑1 (MQ=255)                                                       ttCGTTGGTGCTCATGTTTACGCAAGTACCTgcgc                                              <  1:364417/35‑1 (MQ=37)                                                         cGTTGGTGCTCATGTTTACGGAAGTACCTGcgccg                                            >  1:2002060/1‑35 (MQ=255)                                                         cGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCGGGAACTGATCCATGCAt        >  1:210973/1‑71 (MQ=255)                                                         cGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCGGGAACTGATCCATGCAt        >  1:149483/1‑71 (MQ=255)                                                         cGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCGGGAACTGATCCATGCAt        >  1:627144/1‑71 (MQ=255)                                                           ttGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCGGGAACTGATCCATGCAttt      <  1:2490596/71‑1 (MQ=255)                                                            tGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCGGGAACTGATCCATGCAtttt     >  1:1426530/1‑71 (MQ=255)                                                            tGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCGGGAACTGATCCATGCAttt      >  1:124380/1‑70 (MQ=255)                                                             ggTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCGGGAACTGATCCATGCAtttt     <  1:2608549/70‑1 (MQ=255)                                                             ggTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCGGGAACTGATCCATGCAtttt     <  1:102959/70‑1 (MQ=255)                                                               tgtTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCGGGAACTGATCCATGCATTTTCCg  <  1:286631/71‑1 (MQ=255)                                                                 |                                                                     CTGATATTCACCTTTTTGGTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGTTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCGGGAACTGATCCATGCATTTTCCG  >  minE/2829600‑2829732

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