Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,828,639	C→A	12.4%	A451S (GCG→TCG)	aspA ←	aspartate ammonia‑lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,828,639	0	C	A	12.4%	286.7 / 14.8	88	A451S (GCG→TCG)	aspA	aspartate ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base C (53/24);  new base A (7/4);  total (60/28)
Fisher's exact test for biased strand distribution p-value = 7.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01

TCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACACTCTTACCGGTTTCGGCACAGATTTT  >  minE/2828569‑2828705                                                                       |                                                                   tcaGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGc                                                                    >  1:788381/1‑71 (MQ=255)     tataGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCtt                                                                  >  1:4067481/1‑69 (MQ=255)     tataGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCtt                                                                  >  1:4071816/1‑69 (MQ=255)       taGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAg                                                               >  1:4093842/1‑70 (MQ=255)       taGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGt                                                              >  1:876167/1‑71 (MQ=255)       taGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGATTCAGt                                                              >  1:1224864/1‑71 (MQ=255)       taGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGATTCAGt                                                              >  1:2450392/1‑71 (MQ=255)       taGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGATTCAGt                                                              >  1:392779/1‑71 (MQ=255)       taGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGATTCAGt                                                              >  1:2201739/1‑71 (MQ=255)       taGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGATTCAGt                                                              >  1:2329799/1‑71 (MQ=255)         gCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTc                                                             <  1:2782726/70‑1 (MQ=255)         gCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTc                                                             <  1:2029875/70‑1 (MQ=255)         gCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTc                                                             <  1:128498/70‑1 (MQ=255)         gCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTc                                                             <  1:106747/70‑1 (MQ=255)         gCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTATGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTc                                                             >  1:354067/1‑70 (MQ=255)              ttGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGa                                                       >  1:1971114/1‑71 (MQ=255)               tGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGAc                                                      >  1:476091/1‑71 (MQ=255)                gCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGAcc                                                     >  1:2684321/1‑71 (MQ=255)                     gTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTc                                                             <  1:506571/58‑1 (MQ=255)                      tAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTc                                               >  1:1110161/1‑71 (MQ=255)                         gCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAg                                            <  1:258491/71‑1 (MQ=255)                         gCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAg                                            <  1:3348895/71‑1 (MQ=255)                         gCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGATTCAGTCAACAGACCGCGTTCCAg                                            <  1:3517286/71‑1 (MQ=255)                         gCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGATTCAGTCAACAGACCGCGTTCCAg                                            <  1:3953028/71‑1 (MQ=255)                         gCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGATTCAGTCAACAGACCGCGTTCCAg                                            <  1:3652243/71‑1 (MQ=255)                         gCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGATTCAGTCAACAGACCGCGTTCCAg                                            <  1:3798280/71‑1 (MQ=255)                               tGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGAc                                      >  1:1801524/1‑71 (MQ=255)                               tGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGAc                                      <  1:2302278/71‑1 (MQ=255)                               tGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGAc                                      >  1:47006/1‑71 (MQ=255)                               tGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGAc                                      >  1:51911/1‑71 (MQ=255)                                   tCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACa                                                         <  1:4030347/48‑1 (MQ=255)                                   tCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGa                                                       <  1:3232073/50‑1 (MQ=255)                                   tCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGa                                                       <  1:4044661/50‑1 (MQ=255)                                   tCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGa                                                       <  1:1445259/50‑1 (MQ=255)                                   tCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAAc                                         <  1:1542210/64‑1 (MQ=255)                                   tCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCa                                  <  1:2799756/71‑1 (MQ=255)                                     aGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGCCCGCGTTCCAg                                            >  1:1875853/1‑59 (MQ=255)                                     aGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAg                                            >  1:1002886/1‑59 (MQ=255)                                     aGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAg                                            >  1:3742369/1‑59 (MQ=255)                                     aGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAg                                            >  1:2124157/1‑59 (MQ=255)                                     aGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAg                                            >  1:709787/1‑59 (MQ=255)                                     aGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACg                                >  1:3610859/1‑71 (MQ=255)                                     aGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACg                                >  1:899215/1‑71 (MQ=255)                                     aGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACg                                >  1:897797/1‑71 (MQ=255)                                     aGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACg                                >  1:1827888/1‑71 (MQ=255)                                     aGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACg                                >  1:1780501/1‑71 (MQ=255)                                     aGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACg                                >  1:2338500/1‑71 (MQ=255)                                     aGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACg                                >  1:242642/1‑71 (MQ=255)                                     aGATTCTGTACGGAGAAAATATCGTCAAGTTCAGCTTCAGTCAACAGACCGCGTTCCAg                                            >  1:2805761/1‑59 (MQ=255)                                     aGAATCTGTACCGAGAAAATATCGACAAGTTCCGc                                                                    >  1:1880690/1‑35 (MQ=25)                                      gATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACtt                                    <  1:2786521/66‑1 (MQ=255)                                         tCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTaca                            >  1:2177643/1‑71 (MQ=255)                                          cTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTaca                            <  1:202308/70‑1 (MQ=255)                                          cTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTaca                            <  1:4000702/70‑1 (MQ=255)                                          cTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTaca                            <  1:3808860/70‑1 (MQ=255)                                           tGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGa                                       >  1:1528754/1‑58 (MQ=255)                                           tGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGa                                       >  1:3651323/1‑58 (MQ=255)                                           tGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGa                                       >  1:2458049/1‑58 (MQ=255)                                           tGTACGGAGAAAATATCGTCAAGTTCCGATTCAGTCAACAGACCGCGTTCCAGAACGa                                       >  1:2474559/1‑58 (MQ=255)                                             tACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTacac                           <  1:4020784/68‑1 (MQ=255)                                             tACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACAct                          <  1:628718/69‑1 (MQ=255)                                             tACGGAGAAAATATCGTCAAGTTCCGCTTCAGACAACAGACCGCGTTCCAGAACGACTTCACGTACActct                        <  1:3756387/71‑1 (MQ=255)                                               cGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGAcc                                                     >  1:2748655/1‑40 (MQ=255)                                               cGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGAcc                                                     >  1:3648749/1‑40 (MQ=255)                                               cGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGAcc                                                     >  1:628876/1‑40 (MQ=255)                                               cGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGAcc                                                     >  1:1262195/1‑40 (MQ=255)                                               cGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGAcc                                                     >  1:2848307/1‑40 (MQ=255)                                               cGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGAcc                                                     >  1:2176605/1‑40 (MQ=255)                                               cGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGAcc                                                     >  1:298537/1‑40 (MQ=255)                                               cGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGAcc                                                     >  1:112193/1‑40 (MQ=255)                                                 gagaAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACACTCTTAcc                    >  1:933052/1‑71 (MQ=255)                                                 gagaAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACACTCTTAcc                    >  1:3254212/1‑71 (MQ=255)                                                     aaaTATCGTCAAGTTCCGCTTCAGTCAACAGACCGCTTTCCAGAACGACTTCACGTACACTCTTACCGGt                 <  1:3221886/70‑1 (MQ=255)                                                       atatCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAAc                                         >  1:800001/1‑44 (MQ=255)                                                       atatCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACACTCTTACCGGTTTc              >  1:423027/1‑71 (MQ=255)                                                       atatCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACACTCTTACCGGTTTc              >  1:3524863/1‑71 (MQ=255)                                                       atatCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACACTCTTACCGGTTTc              >  1:3394186/1‑71 (MQ=255)                                                         atCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACACTCTTACCGGTTTCgg            >  1:3821042/1‑71 (MQ=255)                                                         atCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACACTCTTACCGGTTTCgg            >  1:143827/1‑71 (MQ=255)                                                         atCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACACTCTTACCGGTTTCg             >  1:2848002/1‑70 (MQ=255)                                                          tCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACAct                          >  1:454802/1‑56 (MQ=255)                                                          tCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACAct                          >  1:627403/1‑56 (MQ=255)                                                          tCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACAct                          >  1:3156644/1‑56 (MQ=255)                                                          tCGTCAAGTTCCGCTTCAGTCAACAGACCGCCTTCCAGAACGACTTCACGTACAct                          >  1:2247068/1‑56 (MQ=255)                                                                  ttCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACACTCTTACCGGTTTCGGCACAGAttt   <  1:3421249/71‑1 (MQ=255)                                                                   tCCGATTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACACTCTTACCGGTTTCGGCACAGAttt   >  1:1056291/1‑70 (MQ=255)                                                                    ccGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACACTCTTACCGGttt               >  1:3146079/1‑57 (MQ=255)                                                                    ccGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACACTCTTACCGGTTTc              >  1:1464095/1‑58 (MQ=255)                                                                    ccGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACACTCTTACCGGTTTCGGCACAGAtttt  <  1:3194117/70‑1 (MQ=255)                                                                    ccGCTTCAGTCAACAGACCGCGTTCCAGAACAACTTCACGTACACTCTTACCGGTTTc              >  1:3888077/1‑58 (MQ=255)                                                                       |                                                                   TCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAACGACTTCACGTACACTCTTACCGGTTTCGGCACAGATTTT  >  minE/2828569‑2828705

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