Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,323,532	A→C	49.4%	D274E (GAT→GAG)	pfkA ←	6‑phosphofructokinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,323,532	0	A	C	49.4%	3.2 / 103.3	81	D274E (GAT→GAG)	pfkA	6‑phosphofructokinase I
Reads supporting (aligned to +/- strand):  ref base A (24/17);  new base C (17/23);  total (41/40)
Fisher's exact test for biased strand distribution p-value = 1.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGAACCAGCTGTTCGTTCTGGATACCTACACAACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCATACGGGAAGCCAGAATACGGTCGTAAGGCACCGGAGAACCACCGCGCTG  >  minE/2323471‑2323600                                                              |                                                                     ggtgAACCAGCTGTTCGTTCTGGATACCTACACAACGACCGCCGTAACCTGCCAGCAGCAGATCGATAg                                                               <  1:1858990/69‑1 (MQ=255)    gAACCAGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGt                                                            <  1:1708824/69‑1 (MQ=255)    gAACCAGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGt                                                            <  1:2254338/69‑1 (MQ=255)    gAACCAGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGt                                                            <  1:2104254/69‑1 (MQ=255)    gAACCAGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGt                                                            <  1:1992493/69‑1 (MQ=255)    gAACCAGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGt                                                            <  1:2615969/69‑1 (MQ=255)    gAACCAGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGt                                                            <  1:1827933/69‑1 (MQ=255)    gAACCAGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGt                                                            <  1:1807854/69‑1 (MQ=255)    gAACCAGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGt                                                            <  1:38808/69‑1 (MQ=255)    gAACCAGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGt                                                            <  1:648882/69‑1 (MQ=255)    gAACCAGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGt                                                            <  1:1241346/69‑1 (MQ=255)    gAACCAGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGt                                                            <  1:145281/69‑1 (MQ=255)    gAACCAGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGt                                                            <  1:1400149/69‑1 (MQ=255)    gAACCAGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGt                                                            <  1:1470022/69‑1 (MQ=255)    gAACCAGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCAAGCAGCAGCTCGATAGCGt                                                            <  1:1526756/69‑1 (MQ=255)     aaCCAGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGt                                                            <  1:2480164/68‑1 (MQ=255)         aGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAgcg                                                       >  1:184856/1‑69 (MQ=255)         aGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAgcg                                                       >  1:93357/1‑69 (MQ=255)         aGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAgcg                                                       >  1:2065230/1‑69 (MQ=255)         aGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAgcg                                                       >  1:2268084/1‑69 (MQ=255)         aGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAgcg                                                       >  1:2521671/1‑69 (MQ=255)         aGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAgcg                                                       >  1:1354856/1‑69 (MQ=255)         aGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAgcg                                                       >  1:1359702/1‑69 (MQ=255)         aGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAgcg                                                       >  1:3200115/1‑69 (MQ=255)         aGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAgcg                                                       >  1:1894148/1‑69 (MQ=255)         aGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAgcg                                                       >  1:1848141/1‑69 (MQ=255)         aGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAgcg                                                       >  1:66379/1‑69 (MQ=255)         aGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAgcg                                                       >  1:362979/1‑69 (MQ=255)         aGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAgcg                                                       >  1:1543635/1‑69 (MQ=255)         aGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAgcg                                                       >  1:1494352/1‑69 (MQ=255)         aGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAgcg                                                       >  1:601687/1‑69 (MQ=255)         aGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAg                                                         >  1:879082/1‑67 (MQ=255)         aGCTGTTCGTTCTGGATACCTAC‑GAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAgcg                                                       >  1:416705/1‑68 (MQ=255)             gttcgttcTGGATACCTACACAACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCa                                                   >  1:1245800/1‑69 (MQ=255)             gttcgttcTGGATACCTACACAACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCa                                                   >  1:2036442/1‑69 (MQ=255)                cgttcTGGATACCTACACAACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCATa                                                 <  1:2187042/68‑1 (MQ=255)                cgttcTGGATACCTACACAACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCATa                                                 <  1:3055639/68‑1 (MQ=255)                cgttcTGGATACCTACACAACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCATa                                                 <  1:3045183/68‑1 (MQ=255)                cgttcTGGATACCTACACAACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCATa                                                 <  1:2018933/68‑1 (MQ=255)                cgttcTGGATACCTACACAACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCATa                                                 <  1:2843954/68‑1 (MQ=255)                cgttcTGGATACCTACACAACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCATAc                                                >  1:837423/1‑69 (MQ=255)                cgttcTGGATACCTACACAACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCATAc                                                >  1:1725362/1‑69 (MQ=255)                cgttcTGGATACCTACACAACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCATAc                                                >  1:2225426/1‑69 (MQ=255)                cgttcTGGATACCTACACAACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCATAc                                                >  1:2315724/1‑69 (MQ=255)                cgttcTGGATACCTACACAACGACCGCAGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCATAc                                                >  1:2110762/1‑69 (MQ=255)                    cTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGt                                                            <  1:2291226/53‑1 (MQ=38)                             tACACAACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCATACGGGAAGCCAGAAt                                   >  1:1897715/1‑69 (MQ=255)                              acacAACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCa                                                   <  1:984062/52‑1 (MQ=255)                              acacAACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCa                                                   <  1:589294/52‑1 (MQ=255)                              acacAACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCa                                                   <  1:3297257/52‑1 (MQ=255)                              acacAACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCa                                                   <  1:2560595/52‑1 (MQ=255)                              acacAACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCa                                                   <  1:2420389/52‑1 (MQ=255)                               cacaACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCg                                                             >  1:726326/1‑41 (MQ=255)                               cacaACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTa                                                           >  1:624244/1‑43 (MQ=255)                               cacaACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTa                                                           >  1:219305/1‑43 (MQ=255)                               cacaACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTa                                                           >  1:313893/1‑43 (MQ=255)                               cacaACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTa                                                           >  1:3101123/1‑43 (MQ=255)                                 caACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCATACGGGAAGCCAGAATACg                                <  1:433744/68‑1 (MQ=255)                                 caACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCATACGGGAAGCCAGAATACg                                <  1:3258550/68‑1 (MQ=255)                                 caACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCATACGGGAAGCCAGAATACg                                <  1:1482090/68‑1 (MQ=255)                                 caACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCATACGGGAAGCCAGAATACg                                <  1:837087/68‑1 (MQ=255)                                 caACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCATACGGGAAGCCAGAATACg                                <  1:3016590/68‑1 (MQ=255)                                 caACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCATACGGGAAGCCAGAATACg                                <  1:267176/68‑1 (MQ=255)                                   aCGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAGCGCCCATACggg                                             <  1:818720/53‑1 (MQ=39)                                   aCGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAGCGCCCATACggg                                             <  1:219928/53‑1 (MQ=39)                                   aCGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAGCGCCCATACggg                                             <  1:555699/53‑1 (MQ=39)                                   aCGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAGCGCCCATACggg                                             <  1:2181920/53‑1 (MQ=39)                                   aCGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAGCGCCCATACggg                                             <  1:2047111/53‑1 (MQ=39)                                   aCGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAGCGCCCATACggg                                             <  1:3135311/53‑1 (MQ=39)                                   aCGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAGCGCCCATACggg                                             <  1:3080509/53‑1 (MQ=39)                                                    ccAGCAGCAGATCGATAGCGTAAGCGCCCATACGGGAAGCCAGAATACGGTCGTAAGGCACCGGAGAac            >  1:1704023/1‑69 (MQ=255)                                                    ccAGCAGCAGATCGATAGCGTAAGCGCCCATACGGGAAGCCAGAATACGGTCGTAAGGCACCGGAGAac            >  1:22110/1‑69 (MQ=255)                                                    ccAGCAGCAGATCGATAGCGTAAGCGCCCATACGGGAAGCCAGAATACGGTCGTAAGGCACCGGAGAac            >  1:1405836/1‑69 (MQ=255)                                                    ccAGCAGCAGATCGATAGCGTAAGCGCCCATACGGGAAGCCAGAATACGGTCGTAAGGCACCGGAGAac            >  1:714193/1‑69 (MQ=255)                                                    ccAGCAGCAGATCGATAGCGTAAGCGCCCATACGGGAAGCCAGAATACGGTCGTAAGGCACCGGAGAac            >  1:3065929/1‑69 (MQ=255)                                                    ccAGCAGCAGATCGATAGCGTAAGCGCCCATACGGGAAGCCAGAATACGGTCGTAAGGCACCGGAGAac            >  1:2426335/1‑69 (MQ=255)                                                    ccAGCAGCAGATCGATAGCGTAAGCGCCCATACGGGAAGCCAGAATACGGTCGTAAGGCACCGGAGAac            >  1:2481221/1‑69 (MQ=255)                                                    ccAGCAGCAGATCGATAGCGTAAGCGCCCATACGGGAAGCCAGAATACGGTCGTAAGGCACCGGAGAac            >  1:1395800/1‑69 (MQ=255)                                                    ccAGCAGCAGATCGATAGCGTAAGCGCCCATACGGGAAGCCAGAATACGGTCGTAAGGCACCGGAGAac            >  1:1962896/1‑69 (MQ=255)                                                    ccAGCAGCAGATCGATAGCGTAAGCGCCCATACGGGAAACCAGaa                                    >  1:12813/1‑45 (MQ=255)                                                    ccAGCAGCAGATCGATAGCGTAAGCGCCCAAACGGGAAGCCAGAATACGGTCGTAAGGCACCGGAGAac            >  1:163406/1‑69 (MQ=255)                                                              aTCGATAGCGTAAGCGCCCATACGGGAAGCCAGAATACGGTCGTAAGGCACCGGAGAACCACCGCGCtt  <  1:149352/69‑2 (MQ=255)                                                              |                                                                     GGTGAACCAGCTGTTCGTTCTGGATACCTACACAACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCATACGGGAAGCCAGAATACGGTCGTAAGGCACCGGAGAACCACCGCGCTG  >  minE/2323471‑2323600

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