Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,829,174	A→G	26.0%	S272S (TCT→TCC)	aspA ←	aspartate ammonia‑lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,829,174	0	A	G	26.0%	71.3 / 23.1	65	S272S (TCT→TCC)	aspA	aspartate ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base A (22/26);  new base G (6/11);  total (28/37)
Fisher's exact test for biased strand distribution p-value = 5.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCTTCACAGCCAGGCGTTTCAGCGCGCCGTGAACCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTT  >  minE/2829118‑2829239                                                         |                                                                  aTCTTCACAGCCAGGCGTTTCAGCGCGCCGTGAACCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCg                                                       <  1:2465595/69‑1 (MQ=255)    ttCACAGCCAGGCGTTTCAGCGCGCCGTGAACCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATc                                                    >  1:1071732/1‑69 (MQ=255)    ttCACAGCCAGGCGTTTCAGCGCGCCGTGAACCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATc                                                    >  1:1352611/1‑69 (MQ=255)        caGCCAGGCGTTTCAGCGCGCCGTGAACCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATCAg                                                  <  1:2206857/67‑1 (MQ=255)        caGCCAGGCGTTTCAGCGCGCCGTGAACCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATCAg                                                  <  1:2619329/67‑1 (MQ=255)        caGCCAGGCGTTTCAGCGCGCCGTGAACCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATCAg                                                  <  1:729790/67‑1 (MQ=255)        caGCCAGGCGTTTCAGCGCGCCGTGAACCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATCAg                                                  <  1:1792776/67‑1 (MQ=255)        caGCCAGGCGTTTCAGCGCGCCGTGAACCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATCAg                                                  <  1:734068/67‑1 (MQ=255)        caGCCAGGCGTTTCAGCGCGCCGTGAACCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATCAg                                                  <  1:432917/67‑1 (MQ=255)        caGCCAGGCGTTTCAGCGCGCCGTGAACCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATCAg                                                  <  1:1634710/67‑1 (MQ=255)        caGCCAGGCGTTTCAGCGCGCCGTGAACCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATCAg                                                  <  1:1577404/67‑1 (MQ=255)        caGCCAGGCGTTTCAGCGCGCCGTGAACCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATCAg                                                  <  1:1275995/67‑1 (MQ=255)        caGCCAGGCGTTTCAGCGCGCCGTGAACCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATCAg                                                  <  1:908146/67‑1 (MQ=255)        caGCCAGGCGTTTCAGCGCGCCGTGAACCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATCAg                                                  <  1:3409890/67‑1 (MQ=255)          gCCAGGCGTTTCAGCGCGCCGTGATCCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATCAGGt                                                >  1:3195288/1‑67 (MQ=255)          gCCAGGCGTTTCAGCGCGCCGTGAACCATAACATAATCGCCGCAGTCAGAGGTCGCTTCGATCAGGTc                                               >  1:642771/1‑68 (MQ=255)          gCCAGGCGTTTCAGCGCGCCGTGAACCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTc                                               >  1:1562224/1‑68 (MQ=255)          gCCAGGCGTTTCAGCGCGCCGTGAACCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTc                                               >  1:470350/1‑68 (MQ=255)          gCCAGGCGTTTCAGCGCGCCGTGAACCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTc                                               >  1:2968817/1‑68 (MQ=255)          gCCAGGCGTTTCAGCGCGCCGTGAACCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTc                                               >  1:1805382/1‑68 (MQ=255)          gCCAGGCGTTTCAGCGCGCCGTGAACCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTc                                               >  1:2383052/1‑68 (MQ=255)          gCCAGGCGTTTCAGCGCGCCGTGAACCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTc                                               >  1:2621677/1‑68 (MQ=255)          gCCAGGCGTTTCAGCGCGCCGTGAACCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTc                                               >  1:1585371/1‑68 (MQ=255)          gCCAGGCGTTTCAGCGCGCCGTGAACCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTc                                               >  1:333093/1‑68 (MQ=255)                cGTTTCAGCGCGCCGTGAACCATAACATAAGCGCCGCAGTCagag                                                                <  1:1333940/45‑1 (MQ=255)                    tCAGCGCGCCGTGAACCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTa                                    >  1:2274965/1‑69 (MQ=255)                    tCAGCGCGCCGTGAACCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTa                                    >  1:235263/1‑69 (MQ=255)                    tCAGCGCGCCGTGAACCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTa                                    >  1:2700578/1‑69 (MQ=255)                    tCAGCGCGCCGTGAACCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTa                                    >  1:445536/1‑69 (MQ=255)                        cgcgcCGTGAACCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCa                                <  1:2305165/69‑1 (MQ=255)                        cgcgcCGTGAACCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCa                                <  1:502690/69‑1 (MQ=255)                        cgcgcCGTGAACCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCa                                <  1:1498514/69‑1 (MQ=255)                        cgcgcCGTGAACCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCa                                <  1:3346791/69‑1 (MQ=255)                                 aaCCATAACATAAGCGCCGCAGTCGGAGGTCGCTTCGATCATGTCTTCAGCCGGTACGCATGGGAAGcc                       <  1:2752185/69‑1 (MQ=255)                                 aaCCATAACATAAGCGCCGCAGTCGGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGcc                       <  1:2802112/69‑1 (MQ=255)                                 aaCCATAACATAAGCGCCGCAGTCGGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGcc                       >  1:2808085/1‑69 (MQ=255)                                 aaCCATAACATAAGCGCCGCAGTCGGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGcc                       <  1:3103106/69‑1 (MQ=255)                                 aaCCATAACATAAGCGCCGCAGTCGGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGcc                       >  1:3429093/1‑69 (MQ=255)                                 aaCCATAACATAAGCGCCGCAGTCGGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGcc                       <  1:2158693/69‑1 (MQ=255)                                 aaCCATAACATAAGCGCCGCAGTCGGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGcc                       >  1:2305680/1‑69 (MQ=255)                                 aaCCATAACATAAGCGCCGCAGTCGGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGcc                       <  1:2049209/69‑1 (MQ=255)                                 aaCCATAACATAAGCGCCGCAGTCGGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGcc                       >  1:2022593/1‑69 (MQ=255)                                 aaCCATAACATAAGCGCCGCAGTCGGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGcc                       <  1:1745416/69‑1 (MQ=255)                                 aaCCATAACATAAGCGCCGCAGTCGGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGcc                       <  1:496784/69‑1 (MQ=255)                                 aaCCATAACATAAGCGCCGCAGTCGGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGcc                       <  1:1668290/69‑1 (MQ=255)                                 aaCCATAACATAAGCGCCGCAGTCGGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGcc                       <  1:1572312/69‑1 (MQ=255)                                 aaCCATAACATAAGCGCCGCAGTCGGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGcc                       <  1:131258/69‑1 (MQ=255)                                 aaCCATAACATAAGCGCCGCAGTCGGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGcc                       >  1:1280766/1‑69 (MQ=255)                                 aaCCATAACATAAGCGCCGCAGTCGGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGcc                       >  1:652725/1‑69 (MQ=255)                                 aaCCATAACATAAGCGCCGCAGTCGGAGGGCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGcc                       <  1:2869463/69‑1 (MQ=255)                                    cataacataaGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGt                    <  1:975692/69‑1 (MQ=255)                                    cataacataaGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGt                    <  1:624948/69‑1 (MQ=255)                                    cataacataaGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGt                    <  1:2885273/69‑1 (MQ=255)                                    cataacataaGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGt                    <  1:789279/69‑1 (MQ=255)                                      taacataaGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTCTTc                                            >  1:301927/1‑43 (MQ=255)                                      taacataaGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTaa                  >  1:318196/1‑69 (MQ=255)                                      taacataaGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTaa                  >  1:1710500/1‑69 (MQ=255)                                      taacataaGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTaa                  >  1:13/1‑69 (MQ=255)                                      taacataaGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTaa                  >  1:707573/1‑69 (MQ=255)                                        acataaGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATggg                            >  1:1613046/1‑57 (MQ=255)                                                      gTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGtttttt  <  1:2450402/69‑1 (MQ=255)                                                      gTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGtttttt  <  1:1979563/69‑1 (MQ=255)                                                      gTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGtttttt  <  1:1738069/69‑1 (MQ=255)                                                      gTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGtttttt  <  1:1711864/69‑1 (MQ=255)                                                      gTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGtttttt  <  1:1207608/69‑1 (MQ=255)                                                       tCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGtttttt  <  1:1385825/68‑1 (MQ=255)                                                         |                                                                  ATCTTCACAGCCAGGCGTTTCAGCGCGCCGTGAACCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTT  >  minE/2829118‑2829239

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