Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,876,341	A→G	66.4%	D361D (GAT→GAC)	eno ←	enolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,876,341	0	A	G	66.4%	68.4 / 45.7	86	D361D (GAT→GAC)	eno	enolase
Reads supporting (aligned to +/- strand):  ref base A (12/17);  new base G (12/45);  total (24/62)
Fisher's exact test for biased strand distribution p-value = 7.38e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGTCAGCGATGGTAGCGTCTTCAGTTTCGCCAGAACGGTGAGAGATAACTGCAGTGTAGCCAGCATCTTTCGCCATCTTGATTGCAGCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGA  >  minE/1876275‑1876408                                                                   |                                                                    cAGGTCAGCGATGGTAGCGTCTTCAGTTTCGCCAGAACGGTGAGAGATAACTGCAGTGTAGCCAGCATc                                                                   >  1:939305/1‑69 (MQ=255) cAGGTCAGCGATGGTAGCGTCTTCAGTTTCGCCAGAACGGTGAGAGATAACTGCAGTGTAGCCAGCATc                                                                   >  1:2830133/1‑69 (MQ=255) cAGGTCAGCGATGGTAGCGTCTTCAGTTTCGCCAGAACGGTGAGAGATAACTGCAGTGTAGCCAGCATc                                                                   >  1:2122638/1‑69 (MQ=255) cAGGTCAGCGATGGTAGCGTCTTCAGTTTCGCCAGAACGGTGAGAGATAACTGCAGTGTAGCCAGCATc                                                                   >  1:1572560/1‑69 (MQ=255)            tGGTATCATCTTCAGTTTCACCAGAACGGTGAGAGATAACTGCAGTGTAGCCAGCGTCTTTCGCCATCt                                                        <  1:1774560/69‑1 (MQ=255)            tGGTAGCATCTTCAGTTTCACCAGAACGGTGAGAGATAACTTCAGTGTAGCCAGCGTCTTTCGCCATCt                                                        <  1:2405895/69‑1 (MQ=255)            tGGTAGCATCTTCAGTTTCACCAGAACGGTGAGAGATAACTGCAGTGTAGCCAGCGTCTTTCGCCATCt                                                        <  1:837357/69‑1 (MQ=255)            tGGTAGCATCTTCAGTTTCACCAGAACGGTGAGAGATAACTGCAGTGTAGCCAGCGTCTTTCGCCATCt                                                        <  1:1961099/69‑1 (MQ=255)            tGGTAGCATCTTCAGTTTCACCAGAACGGTGAGAGATAACTGCAGTGTAGCCAGCGTCTTTCGCCATCt                                                        <  1:1998035/69‑1 (MQ=255)            tGGTAGCATCTTCAGTTTCACCAGAACGGTGAGAGATAACTGCAGTGTAGCCAGCGTCTTTCGCCATCt                                                        <  1:3310568/69‑1 (MQ=255)            tGGTAGCATCTTCAGTTTCACCAGAACGGTGAGAGATAACTGCAGTGTAGCCAGCGTCTTTCGCCATCt                                                        <  1:2027808/69‑1 (MQ=255)            tGGTAGCATCTTCAGTTTCACCAGAACGGTGAGAGATAACTGCAGTGTAGCCAGCGTCTTTCGCCATCt                                                        <  1:3354647/69‑1 (MQ=255)            tGGTAGCATCTTCAGTTTCACCAGAACGGTGAGAGATAACTGCAGTGTAGCCAGCGTCTTTCGCCATCt                                                        <  1:283390/69‑1 (MQ=255)            tGGTAGCATCTTCAGTTTCACCAGAACGGTGAGAGATAACTGCAGTGTAGCCAGCGTCTTTCGCCATCt                                                        <  1:1804023/69‑1 (MQ=255)            tGGTAGCATCTTCAGTTTCACCAGAACGGTGAGAGATAACTGCAGTGTAGCCAGCGTCTTTCGCCATCt                                                        <  1:392426/69‑1 (MQ=255)            tGGTAGCATCTTCAGTTTCACCAGAACGGTGAGAGATAACTGCAGTGTAGCCAGCGTCTTTCGCCATCt                                                        <  1:2652011/69‑1 (MQ=255)            tGGTAGCATCTTCAGTTTCACCAGAACGGTGAGAGATAACTGCAGTGTAGCCAGCGTCTTTCGCCATCt                                                        <  1:1229240/69‑1 (MQ=255)            tGGTAGCATCTTCAGTTTCACCAGAACGGTGAGAGATAACTGCAGTGTAGCCAGCGTCTTTCGCCATCt                                                        <  1:634469/69‑1 (MQ=255)            tGGTAGCATCTTCAGTTTCACCAGAACGGTGAGAGATAACTGCAGTGTAGCCAGCGTCTTTCGCCATCt                                                        <  1:2440246/69‑1 (MQ=255)            tGGTAGCATCTTCAGTTTCACCAGAACGGTGAGAGATAACTGCAGTGTAGCCAGCGTCTTTCGCCATCt                                                        <  1:1301244/69‑1 (MQ=255)            tGGTAGCATCTTCAGTTTCACAAGAACGGTGAGAGATAACTGCAGTGTAGCCAGCGTCTTTCGGCATCt                                                        <  1:678369/69‑1 (MQ=255)             ggTAGCGTCTTCAGTTTCTCCAGAACGGTGAGAGATAACTGCAGTGTAGCCAGCATCTTTCGCCATCtt                                                       <  1:2455140/69‑1 (MQ=255)             ggTAGCGTCTTCAGTTTCGCCAGAACGGTGAGAGATAAGTGCAGTGTAGCCAGCATCTTTCGCCATCtt                                                       <  1:1555690/69‑1 (MQ=255)             ggTAGCGTCTTCAGTTTCGCCAGAACGGTGAGAGATAACTGCAGTGTAGCCAGCATCTTTCGCCATCtt                                                       <  1:1951864/69‑1 (MQ=255)             ggTAGCGTCTTCAGTTTCGCCAGAACGGTGAGAGATAACTGCAGTGTAGCCAGCATCTTTCGCCATCtt                                                       <  1:3347626/69‑1 (MQ=255)             ggTAGCATCTTCAGTTTCACCAGAACGGTGAGAGATAACTGCAGTGTAGCCAGCGTCTTTCGCCATCt                                                        <  1:777378/68‑1 (MQ=255)             ggTAGCATCTTCAGTTTCACCAGAACGGTGAGAGATAACTGCAGTGTAGCCAGCGTCTTTCGCCATCt                                                        <  1:433223/68‑1 (MQ=255)             ggTAGCATCTTCAGTTTCACCAGAACGGTGAGAGATAACTGCAGTGTAGCCAGCGTCTTTCGCCATCt                                                        <  1:2623542/68‑1 (MQ=255)                            ttCGCCAGAACGGTGAGAGATAACTGCAGTGTAGCCAGCATCTTTCGCCATCTTGATTGCAGCCAGAGt                                        >  1:1445991/1‑69 (MQ=255)                                    aaCGGTGAGAGATAACTGCAGTGTAGCCAGCGTCTTTCGCCATCt                                                        <  1:2694423/45‑1 (MQ=255)                                           agagaTAACTGCAGTTTAGCCAGCGTCTTTCGCCATCTTAATAGCAGCCAGGGTTTCGGt                                  <  1:2400707/60‑1 (MQ=255)                                           agagaTAACTGCAGTGTAGCCAGCGTCTTTGGCCATCTTAATAGCAGCCAGGGTTTCGGt                                  <  1:905711/60‑1 (MQ=255)                                           agagaTAACTGCAGTGTAGCCAGCGTCTTTCGCCATCTTAATAGCAGCCAGGGTTTCGGt                                  <  1:3144923/60‑1 (MQ=255)                                           agagaTAACTGCAGTGTAGCCAGCGTCTTTCGCCATCTTAATAGCAGCCAGGGTTTCGGt                                  <  1:2905728/60‑1 (MQ=255)                                           agagaTAACTGCAGTGTAGCCAGCGTCTTTCGCCATCTTAATAGCAGCCAGGGTTTCGGt                                  <  1:2893072/60‑1 (MQ=255)                                           agagaTAACTGCAGTGTAGCCAGCGTCTTTCGCCATCTTAATAGCAGCCAGGGTTTCGGt                                  <  1:372551/60‑1 (MQ=255)                                           agagaTAACTGCAGTGTAGCCAGCGTCTTTCGCCATCTTAATAGCAGCCAGGGTTTCGGt                                  <  1:826317/60‑1 (MQ=255)                                           agagaTAACTGCAGTGTAGCCAGCGTCTTTCGCCATCTTAATAGCAGCCAGGGTTTCGGt                                  <  1:2740610/60‑1 (MQ=255)                                           agagaTAACTGCAGTGTAGCCAGCGTCTTTCGCCATCTTAATAGCAGCCAGGGTTTCGGt                                  <  1:766802/60‑1 (MQ=255)                                           agagaTAACTGCAGTGTAGCCAGCGTCTTTCGCCATCTTAATAGCAGCCAGGGTTTCGGt                                  <  1:2583812/60‑1 (MQ=255)                                           agagaTAACTGCAGTGTAGCCAGCGTCTTTCGCCATCTTAATAGCAGCCAGGGTTTCGGt                                  <  1:212237/60‑1 (MQ=255)                                           agagaTAACTGCAGTGTAGCCAGCGTCTTTCGCCATCTTAATAGCAGCCAGGGTTTCGGt                                  <  1:1878094/60‑1 (MQ=255)                                           agagaTAACTGCAGTGTAGCCAGCGTCTTTCGCCATCTTAATAGCAGCCAGGGTTTCGGt                                  <  1:1876623/60‑1 (MQ=255)                                           agagaTAACTGCAGTGTAGCCAGCGTCTTTCGCCATCTTAATAGCAGCCAGGGTTTCGGt                                  <  1:1801120/60‑1 (MQ=255)                                           agagaTAACTGCAGTGTAGCCAGCGTCTTTCGCCATCTTAATAGCAGCCAGGGTTTCGGt                                  <  1:1734726/60‑1 (MQ=255)                                           agagaTAACTGCAGTGTAGCCAGCGTCTTTCGCCATCTTAATAGCAGCCAGGGTTTCGGt                                  <  1:1344561/60‑1 (MQ=255)                                           agagaTAACTGCAGTGTAGCCAGCGTCTTTCGCCATCTTAATAGCAGCCAGGGTTTCGGt                                  <  1:1331990/60‑1 (MQ=255)                                           agagaTAACTGCAGTGTAGCCAGCGTCTTTCGCCATCTTAATAGCAGCCAGGGTTTCGGt                                  <  1:1257693/60‑1 (MQ=255)                                           agagaTAACTGCAGTGTAGCCAGCGTCTTTCGCCATCTTAATAGCAGCCAGGGTTTCGGt                                  <  1:2273455/60‑1 (MQ=255)                                           agagaTAACTGCAGTGTAGCCAGCGTCTTTCGCCATCTTAATAGCAGCCAGGGTTTCGGt                                  <  1:2301444/60‑1 (MQ=255)                                           agagaTAACTGCAGTGTAGCCAGCGTCTTTCGCCATCTTAATAGCAGCCAGGGTTTCGGt                                  <  1:2395280/60‑1 (MQ=255)                                           agagaTAACTGCAGTGTAGCCAGCGTCTTTCGCCATCTTAATAGCAGCCAGGGCTTCGGt                                  <  1:438216/60‑1 (MQ=255)                                           agagaTAACTCCAGTGTAGCCAGCGTCTTTCGCCATCTTAATAGCAGCCAGGGTTTCGGt                                  <  1:2499520/60‑1 (MQ=255)                                            gagaTAACTGCAGTGTAGCCAGCGTCTTTCGCCATCTTAATAGCAGCCAGGGTTTCGGt                                  <  1:590673/59‑1 (MQ=255)                                             agaTAACTGCAGTGTAGCCAGCATCTTTCGCCATCTTGATTGCAGCCAGAGTTTCGGTCAGAGAACCg                        <  1:804237/68‑1 (MQ=255)                                             agaTAACTGCAGTGTAGCCAGCATCTTTCGCCATCTTGATTGCAGCCAGAGTTTCGGTCAGAGAACCg                        <  1:1904490/68‑1 (MQ=255)                                             agaTAACTGCAGTGTAGCCAGCATCTTTCGCCATCTTGATTGCAGCCAGAGTTTCGGTCAGAGAACCg                        <  1:2049076/68‑1 (MQ=255)                                             agaTAACTGCAGTGTAGCCAGCATCTTTCGCCATCTTGATTGCAGCCAGAGTTTCGGTCAGAGAACCg                        <  1:597209/68‑1 (MQ=255)                                                 aaCTGCAGTGTAGCCAGCATCTTTCGCCATCTTGATTGCAGCCAGAGTTTCGGTCAGAGAACCg                        >  1:538399/1‑64 (MQ=255)                                                 aaCTGCAGTGTAGCCAGCATCTTTCGCCATCTTGATTGCAGCCAGAGTTTCGGTCAGAGAACCGATCTg                   >  1:1183545/1‑69 (MQ=255)                                                 aaCTGCAGTGTAGCCAGCATCTTTCGCCATCTTGATTGCAGCCAGAGTTTCGGTCAGAGAACCGATCTg                   >  1:848196/1‑69 (MQ=255)                                                 aaCTGCAGTGTAGCCAGCATCTTTCGCCATCTTGATTGCAGCCAGAGTTTCGGTCAGAGAACCGATCTg                   >  1:3166068/1‑69 (MQ=255)                                                 aaCTGCAGTGTAGCCAGCATCTTTCGCCATCTTGATTGCAGCCAGAGTTTCGGTCAGAGAACCGATCTg                   >  1:1414452/1‑69 (MQ=255)                                                         tgtAGCCAGCGTCTTTCGCCATCTTAATAGCAGCCAGGGTTTCGGTCAGAGAACCGATCTGGTTGAAt            >  1:2405456/1‑68 (MQ=255)                                                         tgtAGCCAGCGTCTTTCGCCATCTTAATAGCAGCCAGGGTTTCGGTCAGAGAACCGATCTGGTTGAAt            >  1:928691/1‑68 (MQ=255)                                                         tgtAGCCAGCGTCTTTCGCCATCTTAATAGCAGCCAGGGTTTCGGTCAGAGAACCGATCTGGTTGAAt            >  1:101648/1‑68 (MQ=255)                                                         tgtAGCCAGCGTCTTTCGCCATCTTAATAGCAGCCAGGGTTTCGGTCAGAGAACCGATCTGGTTGAAt            >  1:1007080/1‑68 (MQ=255)                                                         tgtAGCCAGCGTCTTTCGCCATCTTAATAGCAGCCAGGGTTTCGGTCAGAGAACCGATCTGGTTGAAt            >  1:1095071/1‑68 (MQ=255)                                                         tgtAGCCAGCGTCTTTCGCCATCTTAATAGCAGCCAGGGTTTCGGTCAGAGAACCGATCTGGTTGAAt            >  1:2436807/1‑68 (MQ=255)                                                         tgtAGCCAGCGTCTTTCGCCATCTTAATAGCAGCCAGGGTTTCGGTCAGAGAACCGATCTGGTTGAAt            >  1:3007449/1‑68 (MQ=255)                                                         tgtAGCCAGCGTCTTTCGCCATCTTAATAGCAGCCAGGGTTTCGGTCAGAGAACCGATCTGGTTGAAt            >  1:2670739/1‑68 (MQ=255)                                                         tgtAGCCAGCGTCTTTCGCCATCTTAATAGCAGCCAGGGTTTCGGTCAGAGAACCGATCTGGTTGAAt            >  1:2199138/1‑68 (MQ=255)                                                         tgtAGCCAGCGTCTTTCGCCATCTTAATAGCAGCCAGGGTTTCGGTCAGAGAACCGATCTGGTTGAAt            >  1:2703956/1‑68 (MQ=255)                                                         tgtAGCCAGCGTCTTTCGCCATCTTAATAGCAGCCAGGGTTTCGGTCAGAGAACCGATCTGGTTGAAt            >  1:2133811/1‑68 (MQ=255)                                                         tgtAGCCAGCGTCTTTCGCCATCTTAATAGCAGCCAGGGTTTCGGTCAGAGAACCGATCTGGTTGAAt            >  1:2154485/1‑68 (MQ=255)                                                          gtAGCCAGCATCTTTCGCCATCTTGATTGCAGCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAAttt          <  1:296716/69‑1 (MQ=255)                                                          gtAGCCAGCATCTTTCGCCATCTTGATTGCAGCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAAttt          <  1:627555/69‑1 (MQ=255)                                                           tAGCCAGCATCTTTCGCCATCTTGATTGCAGCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTg         <  1:1001729/69‑1 (MQ=255)                                                           tAGCCAGCATCTTTCGCCATCTTGATTGCAGCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTg         >  1:3335871/1‑69 (MQ=255)                                                           tAGCCAGCATCTTTCGCCATCTTGATTGCAGCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTg         >  1:2853494/1‑69 (MQ=255)                                                                aGCATCTTTCGCCATCTTGATTGCAGCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAg    <  1:2011525/69‑1 (MQ=255)                                                                aGCATCTTTCGCCATCTTGATTGCAGCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAg    <  1:23038/69‑1 (MQ=255)                                                                aGCATCTTTCGCCATCTTGATTGCAGCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAg    <  1:883792/69‑1 (MQ=255)                                                                aGCATCTTTCGCCATCTTGATTGCAGCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAg    <  1:868922/69‑1 (MQ=255)                                                                aGCATCTTTCGCCATCTTGATTGCAGCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAg    <  1:2947611/69‑1 (MQ=255)                                                                 gCATCTTTCGCCATCTTGATTGCAGCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAg    <  1:428262/68‑1 (MQ=255)                                                                   aTCTTTCGCCATCTTGATTGCAGCCAGAGTTTCGGTCAgaga                            >  1:524215/1‑42 (MQ=255)                                                                   aTCTTTCGCCATCTTGATTGCAGCCAGAGTTTCGGTCAgaga                            >  1:3068952/1‑42 (MQ=255)                                                                   aTCTTTCGCCATCTTGATTGCAGCCAGAGTTTCGGTCAgaga                            >  1:1926606/1‑42 (MQ=255)                                                                   aTCTTTCGCCATCTTGATTGCAGCCAGAGTTTCGGTCAgaga                            >  1:2570030/1‑42 (MQ=255)                                                                   aTCTTTCGCCATCTTGATTGCAGCCAGAGTTTCGGTCAgaga                            >  1:3149122/1‑42 (MQ=255)                                                                   aTCTTTCGCCATCTTGATTGCAGCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGa  >  1:2952894/1‑68 (MQ=255)                                                                   aTCTTTCGACATCTTGATTGCAGCCAGAGTTTCGGTCAgaga                            >  1:1810599/1‑42 (MQ=255)                                                                   |                                                                    CAGGTCAGCGATGGTAGCGTCTTCAGTTTCGCCAGAACGGTGAGAGATAACTGCAGTGTAGCCAGCATCTTTCGCCATCTTGATTGCAGCCAGAGTTTCGGTCAGAGAACCGATCTGGTTGAATTTGATCAGGA  >  minE/1876275‑1876408

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