Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,829,664	T→C	47.7%	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	47.7%	5.9 / 129.8	91	N109S (AAC→AGC)	aspA	aspartate ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base T (17/30);  new base C (18/25);  total (36/55)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

AGGTACTGATATTCACCTTTTTGGTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGTTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCGGGAACTGATCCATGCATTTTCCGT  >  minE/2829595‑2829733                                                                      |                                                                      aGGTACTGATATTCACCTTTTTGGTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGgtgt                                                                       <  1:1070832/70‑1 (MQ=255)  ggTACTGATATTCACCTTTTTGGTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGTTc                                                                     <  1:344766/71‑1 (MQ=255)     aCTGATATTCACCTTTTTGGTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATg                                                                  >  1:2550975/1‑71 (MQ=255)     aCTGATATTCACCTTTTTGGTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATg                                                                  >  1:1272463/1‑71 (MQ=255)       tGATATTCACCTTTTTGGTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGTTCATGtt                                                                >  1:2460718/1‑71 (MQ=255)         atatTCACCTTTTTGGTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGttt                                                               <  1:1058332/70‑1 (MQ=255)           atTCACCTTTTTGGTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGTTCATGTTTACg                                                            >  1:1589552/1‑71 (MQ=255)           atTCACCTTTTTGGTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGGTCATGTTTACg                                                            >  1:1693364/1‑71 (MQ=255)            ttCACCTTTTTGGTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGTTc                                                                     <  1:3248316/61‑1 (MQ=255)            ttCACCTTTTTGGTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGTTc                                                                     <  1:1588781/61‑1 (MQ=255)            ttCACCTTTTTGGTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACgg                                                           <  1:1028546/71‑1 (MQ=255)            ttCACCTTTTTGGTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACgg                                                           <  1:3496853/71‑1 (MQ=255)               aCCTTTTTGGTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGTTCATGTTTACGGaa                                                         <  1:83025/70‑1 (MQ=255)               aCCTTTTTGGTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGTTCATGTTTACGGaa                                                         <  1:3145357/70‑1 (MQ=255)                ccTTTTTGGTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACg                                                            >  1:3667395/1‑66 (MQ=255)                ccTTTTTGGTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGt                                                       <  1:951062/71‑1 (MQ=255)                   ttttGGTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGTTCATGTTTACGGAAGTAcc                                                    >  1:2813405/1‑71 (MQ=255)                   ttttGGTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTAcc                                                    >  1:3240276/1‑71 (MQ=255)                   ttttGGTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTAcc                                                    >  1:815707/1‑71 (MQ=255)                   ttttGGTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTAcc                                                    >  1:2065317/1‑71 (MQ=255)                    tttGGTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGTTCATGTTTACGGAAGTAcc                                                    <  1:3653644/70‑1 (MQ=255)                    tttGGTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGTTCATGTTTACGGAAGTAcc                                                    <  1:3372040/70‑1 (MQ=255)                      tGGTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGTTCATGTTTACGGAAGTACCTgc                                                 >  1:1145451/1‑71 (MQ=255)                      tGGTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGaa                                                         >  1:3213944/1‑63 (MQ=255)                      tGGTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTg                                                  <  1:3426224/70‑1 (MQ=255)                       ggTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGTTCATGTTTACGGAAGTACCTgc                                                 <  1:409164/70‑1 (MQ=255)                        gTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTgcg                                                <  1:1069682/70‑1 (MQ=255)                           aCCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGTTCATGTTTACGGAAGTACCTGcgccgc                                            >  1:316201/1‑71 (MQ=255)                           aCCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGTTCATGTTTACGGAAGTACCTGcgccg                                             <  1:674419/70‑1 (MQ=255)                            cccATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACgg                                                           <  1:381302/55‑1 (MQ=255)                            cccATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACgg                                                           <  1:3103431/55‑1 (MQ=255)                            cccATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACgg                                                           <  1:3062248/55‑1 (MQ=255)                            cccATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCgccgcc                                           >  1:366935/1‑71 (MQ=255)                             ccATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGTTCATGTTTACGGAAGTACCTGCGCCGccc                                          >  1:2877738/1‑71 (MQ=255)                               aTCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTg                                        >  1:3159714/1‑71 (MQ=255)                                tCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGTTCATGTTTACg                                                            <  1:3542850/50‑1 (MQ=255)                                      ccAGACCGATATTGGCCAGCACTTCGTTGGTGTTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGAc                                  <  1:2621808/70‑1 (MQ=255)                                      ccAGACCGATATTGGCCAGCACTTCGTTGGTGTTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACg                                 >  1:1984098/1‑71 (MQ=255)                                      ccAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCgccgcc                                           >  1:2949712/1‑61 (MQ=255)                                      ccAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGAc                                  <  1:3575475/70‑1 (MQ=255)                                      ccAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACg                                 >  1:198737/1‑71 (MQ=255)                                      ccAGACCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACg                                 >  1:1499921/1‑71 (MQ=255)                                         gACCGATATTGGCCAGCACTTCGTTGGTGTTCATGTTTACGGa                                                          <  1:3645950/43‑1 (MQ=255)                                         gACCGATATTGGCCAGCACTTCGTTGGTGTTCATGTTTACGGa                                                          <  1:586872/43‑1 (MQ=255)                                         gACCGATATTGGCCAGCACTTCGTTGGTGTTCATGTTTACGGa                                                          <  1:1161602/43‑1 (MQ=255)                                          aCCGATATTGGCCAGCACTTCGTTGGTGTTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTa                             <  1:3605444/71‑1 (MQ=255)                                          aCCGATATTGGCCAGCACTTCGTTGGTGTTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTa                             <  1:1049079/71‑1 (MQ=255)                                          aCCGATATTGGCCAGCACTTCGTTGGTGTTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTa                             <  1:1560081/71‑1 (MQ=255)                                          aCCGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGt                                >  1:166682/1‑68 (MQ=255)                                           ccGATATTGGCCAGCACTTCGTTGGTGTTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTa                             <  1:1584387/70‑1 (MQ=255)                                           ccGATATTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTAc                            <  1:2463048/71‑1 (MQ=255)                                             gATATTGGCCAGCACTTCGTTGGTGTTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCg                          >  1:3717055/1‑71 (MQ=255)                                               tatTGGCCAGCACTTCGTTGGTGTTCATGTTTACGGAAGTACCTgcg                                                >  1:2305547/1‑47 (MQ=255)                                               tatTGGCCAGCACTTCGTTGGTGTTCATGTTTACGGAAGTACCTGCGCCGccc                                          <  1:962205/53‑1 (MQ=255)                                               tatTGGCCAGCACTTCGTTGGTGTTCATGTTTACGGAAGTACCTGCGCCGccc                                          <  1:1373657/53‑1 (MQ=255)                                               tatTGGCCAGCACTTCGTTGGTGTTCATGTTTACGGAAGTACCTGCGCCGccc                                          <  1:2807065/53‑1 (MQ=255)                                               tatTGGCCAGCACTTCGTTGGTGTTCATGTTTACGGAAGTACCTGCGCCGccc                                          <  1:3103582/53‑1 (MQ=255)                                                atTGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCgg                         >  1:3012364/1‑69 (MQ=255)                                                 ttGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGcgccgc                                            >  1:286798/1‑49 (MQ=255)                                                 ttGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCGGGaa                      <  1:2687950/71‑1 (MQ=255)                                                 ttGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCGGGaa                      <  1:1429113/71‑1 (MQ=255)                                                 ttGGCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCGGGaa                      >  1:197903/1‑71 (MQ=255)                                                   ggCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCGGGAACt                    >  1:1463155/1‑71 (MQ=255)                                                    gCCAGCACTTCGTTGGTGTTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCGGGAACTg                   <  1:3718172/71‑1 (MQ=255)                                                    gCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCGGGAACTg                   <  1:1820197/71‑1 (MQ=255)                                                    gCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCGGGAACTg                   <  1:397181/71‑1 (MQ=255)                                                    gCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCGGGAACTg                   <  1:3693172/71‑1 (MQ=255)                                                    gCCAGCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCGGGAACTg                   <  1:789997/71‑1 (MQ=255)                                                      cAGCACTTCGTTGGTGTTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTa                                     <  1:3563048/51‑1 (MQ=255)                                                       aGCACTTCGTTGGTGTTCATGTTTACGGAAGTACCTGcgccgc                                            >  1:173392/1‑43 (MQ=255)                                                       aGCACTTCGTTGGTGTTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGAc                                  <  1:337187/53‑1 (MQ=255)                                                       aGCACTTCGTTGGTGTTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCGGGAACTGATc                >  1:1669852/1‑71 (MQ=255)                                                        gCACTTCGTTGGTGTTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCGGGaa                      >  1:3638999/1‑64 (MQ=255)                                                        gCACTTCGTTGGTGTTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCGGGaa                      >  1:2697194/1‑64 (MQ=255)                                                        gCACTTCGTTGGTGTTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCGGGAACTGATcc               <  1:3350635/71‑1 (MQ=255)                                                        gCACTTCGTTGGTGTTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCGGGAACTGATcc               <  1:2506627/71‑1 (MQ=255)                                                        gCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCGGGAACTGATcc               <  1:1573829/71‑1 (MQ=255)                                                        gCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCGGGAACTGATcc               <  1:1375878/71‑1 (MQ=255)                                                        gCACTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCGGGAACTGATcc               <  1:1237643/71‑1 (MQ=255)                                                         cACTTCGTTGGTGTTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCGGGAACTGATCCa              >  1:3522003/1‑71 (MQ=255)                                                         cACTTCGTTGGTGTTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCGGGAACTGATCCa              >  1:1947890/1‑71 (MQ=255)                                                         cACTTCGTTGGTGTTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCGGGAACTGATCCa              <  1:73860/71‑1 (MQ=255)                                                           cTTCGTTGGTGTTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCGGGAACTGATCCATg            >  1:896535/1‑71 (MQ=255)                                                           cTTCGTTGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACg                                 >  1:3202530/1‑50 (MQ=255)                                                                ttGGTGTTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCGGGAACTGATCCATGCAttt       <  1:1210513/71‑1 (MQ=255)                                                                ttGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACg                                 <  1:2639355/45‑1 (MQ=255)                                                                ttGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACg                                 <  1:1375183/45‑1 (MQ=255)                                                                ttGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACg                                 <  1:2020535/45‑1 (MQ=255)                                                                ttGGTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACg                                 <  1:212277/45‑1 (MQ=255)                                                                 tGGTGTTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCGGGAACTGATCCATGCAtttt      >  1:1004662/1‑71 (MQ=255)                                                                  ggTGTTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCGGGAACTGATCCATGCAtttt      <  1:2885127/70‑1 (MQ=255)                                                                  ggTGCTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCGGGAACTGATCCATGCAtttt      <  1:40016/70‑1 (MQ=255)                                                                      ttCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCGGGAACTGATCCATGCATTTTCCgt  <  1:3282582/70‑1 (MQ=255)                                                                      |                                                                      AGGTACTGATATTCACCTTTTTGGTGACCCATCAGTTCCAGACCGATATTGGCCAGCACTTCGTTGGTGTTCATGTTTACGGAAGTACCTGCGCCGCCCTGGTAGACGTCTACCGGGAACTGATCCATGCATTTTCCGT  >  minE/2829595‑2829733

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