Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,876,422	A→G	41.1%	G334G (GGT→GGC)	eno ←	enolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,876,422	0	A	G	41.1%	24.2 / 52.3	73	G334G (GGT→GGC)	eno	enolase
Reads supporting (aligned to +/- strand):  ref base A (17/26);  new base G (15/15);  total (32/41)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGATTGCAGCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAGGTCGTCACCAACCAG  >  minE/1876354‑1876484                                                                     |                                                               ttGATTGCAGCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATa                                                                >  1:1133409/1‑69 (MQ=255) ttGATTGCAGCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATa                                                                >  1:1133747/1‑69 (MQ=255) ttGATTGCAGCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATa                                                                >  1:781148/1‑69 (MQ=255) ttGATTGCAGCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATa                                                                >  1:35409/1‑69 (MQ=255) ttGATTGCAGCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATa                                                                >  1:1847757/1‑69 (MQ=255) ttGATTGCAGCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATa                                                                >  1:2595503/1‑69 (MQ=255) ttGATTGCAGCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATa                                                                >  1:3046917/1‑69 (MQ=255) ttGATTGCAGCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTAGATCAGGATGGAGTTAGCGATa                                                                >  1:159495/1‑69 (MQ=255)          gCCAGAGTTTCGGTCAGAGAACCTATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTc                                                        <  1:1464272/68‑1 (MQ=255)          gCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTc                                                        <  1:2858731/68‑1 (MQ=255)          gCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTc                                                        <  1:2976125/68‑1 (MQ=255)          gCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTc                                                        <  1:2063430/68‑1 (MQ=255)          gCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTc                                                        <  1:3051903/68‑1 (MQ=255)          gCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTc                                                        <  1:2778924/68‑1 (MQ=255)          gCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTc                                                        <  1:2144334/68‑1 (MQ=255)          gCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTc                                                        <  1:2370504/68‑1 (MQ=255)          gCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTc                                                        <  1:1761117/68‑1 (MQ=255)          gCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTc                                                        <  1:2264451/68‑1 (MQ=255)          gCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTc                                                        <  1:3367726/68‑1 (MQ=255)          gCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTc                                                        <  1:345821/68‑1 (MQ=255)          gCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTc                                                        <  1:383187/68‑1 (MQ=255)          gCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTc                                                        <  1:2396592/68‑1 (MQ=255)          gCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTc                                                        <  1:1493667/68‑1 (MQ=255)          gCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTc                                                        <  1:2470151/68‑1 (MQ=255)          gCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTc                                                        <  1:1253278/68‑1 (MQ=255)          gCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTc                                                        <  1:650924/68‑1 (MQ=255)          gCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTc                                                        <  1:1315255/68‑1 (MQ=255)               agTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTCGATAcc                                                  <  1:1445294/69‑1 (MQ=255)               agTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTCGATAcc                                                  <  1:2893355/69‑1 (MQ=255)                         agagAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTc                                                        <  1:1662391/53‑1 (MQ=255)                                            aTTTGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACg                      <  1:2114592/68‑1 (MQ=255)                                            aTTTGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACg                      <  1:2033288/68‑1 (MQ=255)                                            aTTTGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACg                      <  1:1970374/68‑1 (MQ=255)                                             tttGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACg                      <  1:671784/67‑1 (MQ=255)                                               tGATCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTAc                       >  1:2226962/1‑64 (MQ=255)                                               tGATCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAAc                   >  1:2858264/1‑68 (MQ=255)                                               tGATCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAAc                   >  1:2653551/1‑68 (MQ=255)                                               tGATCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACa                  >  1:819559/1‑69 (MQ=255)                                               tGATCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACa                  >  1:2494292/1‑69 (MQ=255)                                               tGATCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACa                  >  1:2711504/1‑69 (MQ=255)                                               tGATCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACa                  >  1:3299307/1‑69 (MQ=255)                                               tGATCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACa                  >  1:386242/1‑69 (MQ=255)                                               tGATCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACa                  >  1:425656/1‑69 (MQ=255)                                               tGATCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACa                  >  1:653834/1‑69 (MQ=255)                                               tGATCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACa                  >  1:224368/1‑69 (MQ=255)                                               tGATCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACa                  >  1:1325567/1‑69 (MQ=255)                                               tGATCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACa                  >  1:151448/1‑69 (MQ=255)                                               tGATCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACa                  >  1:1871569/1‑69 (MQ=255)                                               tGATCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACa                  >  1:1966956/1‑69 (MQ=255)                                                  tCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAgg                <  1:2960419/68‑1 (MQ=255)                                                  tCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAgg                <  1:1565517/68‑1 (MQ=255)                                                  tCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAgg                <  1:1852878/68‑1 (MQ=255)                                                  tCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAgg                <  1:3256479/68‑1 (MQ=255)                                                  tCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAgg                <  1:3255996/68‑1 (MQ=255)                                                  tCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAgg                <  1:1524891/68‑1 (MQ=255)                                                  tCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAgg                <  1:3150808/68‑1 (MQ=255)                                                  tCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAgg                <  1:2936710/68‑1 (MQ=255)                                                  tCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAgg                <  1:1446038/68‑1 (MQ=255)                                                  tCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAgg                <  1:473068/68‑1 (MQ=255)                                                  tCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAgg                <  1:2853160/68‑1 (MQ=255)                                                  tCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAgg                <  1:1321661/68‑1 (MQ=255)                                                  tCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAgg                <  1:217687/68‑1 (MQ=255)                                                  tCAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAgg                <  1:712265/68‑1 (MQ=255)                                                   cAGGATGGAGTTAGCGATGCCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAgg                <  1:1907698/67‑1 (MQ=255)                                                               aGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGt                          >  1:2555928/1‑45 (MQ=255)                                                               aGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGt                          >  1:720576/1‑45 (MQ=255)                                                               aGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGt                          >  1:370081/1‑45 (MQ=255)                                                               aGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGt                          >  1:124057/1‑45 (MQ=255)                                                               aGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGt                          >  1:1549297/1‑45 (MQ=255)                                                               aGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGt                          >  1:1572817/1‑45 (MQ=255)                                                               aGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGt                          >  1:3159068/1‑45 (MQ=255)                                                               aGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGt                          >  1:3097041/1‑45 (MQ=255)                                                               aGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGt                          >  1:1988072/1‑45 (MQ=255)                                                               aGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGt                          >  1:2001478/1‑45 (MQ=255)                                                               aGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGt                          >  1:2064883/1‑45 (MQ=255)                                                               aGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGt                          >  1:2474173/1‑45 (MQ=255)                                                               aGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAGGTCGTCACCAACCAg  >  1:2141385/1‑69 (MQ=255)                                                               aGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAGGTCGTCACCAACCAg  >  1:1273558/1‑69 (MQ=255)                                                               aGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAGGTCGTCACCAACCAg  >  1:2600063/1‑69 (MQ=255)                                                               aGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAGGTCGTCACCAACCAg  >  1:1107602/1‑69 (MQ=255)                                                               aGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAGCTCGTCACCAACCAg  >  1:238009/1‑69 (MQ=255)                                                                     |                                                               TTGATTGCAGCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGATGGAGTTAGCGATACCTTTTTCGATACCTTCTTTCAGGATCTTGGTGTTGGTTACGAACAGGTCGTCACCAACCAG  >  minE/1876354‑1876484

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑