Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,997,980	C→T	43.8%	V384V (GTG→GTA)	pitB ←	phosphate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,997,980	0	C	T	43.8%	21.0 / 54.8	73	V384V (GTG→GTA)	pitB	phosphate transporter
Reads supporting (aligned to +/- strand):  ref base C (30/11);  new base T (3/29);  total (33/40)
Fisher's exact test for biased strand distribution p-value = 2.72e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCGCTACACGACGCCAGCCAATCATGGTGCCAATGCCGAGCGCCAGTGCTACCGCCATGATGATCCACACCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTT  >  minE/1997912‑1998048                                                                     |                                                                     aTCGCTACACGACGCCAGCCAATCTTGGTGCCAATGCCGAGCGCCAGTGCTACCGCCATGATGATCcac                                                                      >  1:2560939/1‑69 (MQ=255) aTCGCTACACGACGCCAGCCAATCATGGTGCCAATGCCGAGCGCCAGTGCTACCGCCATGATGATCcac                                                                      >  1:1411693/1‑69 (MQ=255) aTCGCTACACGACGCCAGCCAATCATGGTGCCAATGCCGAGCGCCAGTGCTACCGCCATGATGATCcac                                                                      >  1:2168748/1‑69 (MQ=255) aTCGCTACACGACGCCAGCCAATCATGGTGCCAATGCCGAGCGCCAGTGCTACCGCCATGATGATCcac                                                                      >  1:2049857/1‑69 (MQ=255) aTCGCTACACGACGCCAGCCAATCATGGTGCCAATGCCGAGCGCCAGTGCTACCGCCATGATGATCcac                                                                      >  1:2502335/1‑69 (MQ=255) aTCGCTACACGACGCCAGCCAATCATGGTGCCAATGCCGAGCGCCAGTGCTACCGCCATGATGATCcac                                                                      >  1:3110986/1‑69 (MQ=255) aTCGCTACACGACGCCAGCCAATCATGGTGCCAATGCCGAGCGCCAGTGCTACCGCCATGATGATCcac                                                                      >  1:3316264/1‑69 (MQ=255) aTCGCTACACGACGCCAGCCAATCATGGTGCCAATGCCGAGCGCCAGTGCTACCGCCATGATGATCcac                                                                      >  1:335720/1‑69 (MQ=255) aTCGCTACACGACGCCAGCCAATCATGGTGCCAATGCCGAGCGCCAGTGCTACCGCCATGATGATCcac                                                                      >  1:1321722/1‑69 (MQ=255) aTCGCTACACGACGCCAGCCAATCATGGTGCCAATGCCGAGCGCCAGTGCTACCGCCATGATGATCcac                                                                      >  1:1198974/1‑69 (MQ=255) aTCGCTACACGACGCCAGCCAATCATGGTGCCAATGCCGAGCGCCAGTGCTACCGCCATGATGATCcac                                                                      >  1:585932/1‑69 (MQ=255) aTCGCTACACGACGCCAGCCAATCATGGTGCCAATGCCGAGCGCCAGTGCTACCGCCATGATGATCcac                                                                      >  1:498916/1‑69 (MQ=255) aTCGCTACACGACGCCAGACAATCATGGTGCCAATGCCGAGCGCCAGTGCTACCGCCATGATGATCcac                                                                      >  1:1535346/1‑69 (MQ=255)        cacGACGCCAGCCAATCATGGTGCCGATGCCGAGCGCCAGTGCTACCGCCATGATGATCCATACCggtg                                                               >  1:307187/1‑67 (MQ=18)        cacGACGCCAGCCAATCATGGTGCCGATGCCGAGCGCCAGTGCTACCGCCATGATGATCCATACCggtg                                                               >  1:2306379/1‑67 (MQ=18)        cacGACGCCAGCCAATCATGGTGCCGATGCCGAGCGCCAGTGCTACCGCCATGATGATCCATACCggtg                                                               >  1:1659311/1‑67 (MQ=18)                    caATCATGGTGCCAATGCCGAGCGCCAGTGCTACCGCCATGATGATCCACACCGGAGCGTACTCAATgg                                                   >  1:670036/1‑69 (MQ=255)                    caATCATGGTGCCAATGCCGAGCGCCAGTGCTACCGCCATGATGATCCACACCGGAGCGTACTCAATgg                                                   >  1:2554296/1‑69 (MQ=255)                    caATCATGGTGCCAATGCCGAGCGCCAGTGCTACCGCCATGATGATCCACACCGGAGCGTACTCAATgg                                                   >  1:410478/1‑69 (MQ=255)                    caATCATGGTGCCAATGCCGAGCGCCAGTGCTACCGCCATGATGATCCACACCGGAGCGTACTCAATgg                                                   >  1:2621226/1‑69 (MQ=255)                    caATCATGGTGCCAATGCCGAGCGCCAGTGCTACCGCCATGATGATCCACACCGGAGCGTACTCAATgg                                                   >  1:2717658/1‑69 (MQ=255)                    caATCATGGTGCCAATGCCGAGCGCCAGTGCTACCGCCATGATGATCCACACCGGAGCGTACTCAATgg                                                   >  1:1404779/1‑69 (MQ=255)                    caATCATGGTGCCAATGCCGAGCGCCAGTGCTACCGCCATGATGATCCACACCGGAGCGTACTCAATg                                                    >  1:1559262/1‑68 (MQ=255)                    caATCATGGTGCCAATGCCGAGCGCCAGTGCTACCGCCATGATGATCCACACCGGAGCGTACTCAATg                                                    >  1:1684335/1‑68 (MQ=255)                    caATCATGGTGCCAATGCCGAGCGCCAGTGCTACCGCCATGATGATCCACACCGGAGCGTACTCAATg                                                    >  1:2295725/1‑68 (MQ=255)                        cATGGTGCCAATGCCGAGCGCCAGTGCTACCGCCATGATGATCCACACCGGAGCGTACTCAATggg                                                  <  1:1402843/66‑2 (MQ=255)                        cATGGTGCCAATGCCGAGCGCCAGTGCTACCGCCATGATGATCCACACCGGAGCGTACTCAATGGt                                                  <  1:1360351/66‑1 (MQ=255)                        cATGGTGCCAATGCCGAGCGCCAGTGCTACCGCCATGATGATCCACACCGGAGCGTACTCAATGGt                                                  <  1:69174/66‑1 (MQ=255)                        cATGGTGCCAATGCCGAGCGCCAGTGCTACCGCCATGATGATCCACACCGGAGCGTACTCAATGGt                                                  <  1:2380815/66‑1 (MQ=255)                        cATGGTGCCAATGCCGAGCGCCAGTGCTACCGCCATGATGATCCACACCGGAGCGTACTCAATGGt                                                  <  1:2062778/66‑1 (MQ=255)                           ggTGCCAATGCCGAGCGCCAGTGCTACCGCCATGATGATCCACACCGGAGCGTACTCAATGGTGCTTaa                                            <  1:3035188/69‑1 (MQ=255)                           ggTGCCAATGCCGAGCGCCAGTGCTACCGCCATGATGATCCACACCGGAGCGTACTCAATGGTGCTTaa                                            <  1:969082/69‑1 (MQ=255)                           ggTGCCAATGCCGAGCGCCAGTGCTACCGCCATGATGATCCACACCGGAGCGTACTCAATGGTGCTTaa                                            <  1:748698/69‑1 (MQ=255)                           ggTGCCAATGCCGAGCGCCAGTGCTACCGCCATGATGATCCACACCGGAGCGTACTCAATGGTGCTTaa                                            <  1:785870/69‑1 (MQ=255)                           ggTGCCAATGCCGAGCGCCAGTGCTACCGCCATGATGATCCACACCGGAGCGTACTCAATGGTGCTTaa                                            <  1:1638468/69‑1 (MQ=255)                               ccAATGCCGAGCGCCAGTGCTACCGCCATGATGATCCACACCGGAGCGTACTCAATGGTGCTTAACata                                        >  1:1750199/1‑69 (MQ=255)                                   tGCCGAGCGCCAGTGCTACCGCCATGATGATCCACACCGGAGCGTACTCAATGGTGCTTAACATATCGc                                    >  1:2619516/1‑69 (MQ=255)                                   tGCCGAGCGCCAGTGCTACCGCCATGATGATCCACACCGGAGCGTACTCAATGGTGCTTAACATATCGc                                    >  1:2571098/1‑69 (MQ=255)                                   tGCCGAGCGCCAGTGCTACCGCCATGATGATCCACACCGGAGCGTACTCAATGGTGCTTAACATATCGc                                    >  1:752295/1‑69 (MQ=255)                                    gCCGAGCGCCAGTGCTACCGCCATGATGATCCATACCGGTGCGTACTCAATGGTGCTTAACATATCGCt                                   <  1:1780472/69‑1 (MQ=255)                                    gCCGAGCGCCAGTGCTACCGCCATGATGATCCATACCGGTGCGTACTCAATGGTGCTTAACATATCGCt                                   <  1:3002454/69‑1 (MQ=255)                                    gCCGAGCGCCAGTGCTACCGCCATGATGATCCATACCGGTGCGTACTCAATGGTGCTTAACATATCGCt                                   <  1:3332648/69‑1 (MQ=255)                                    gCCGAGCGCCAGTGCTACCGCCATGATGATCCATACCGGTGCGTACTCAATGGTGCTTAACATATCGCt                                   <  1:2420383/69‑1 (MQ=255)                                    gCCGAGCGCCAGTGCTACCGCCATGATGATCCATACCGGTGCGTACTCAATGGTGCTTAACATATCGCt                                   <  1:1081243/69‑1 (MQ=255)                                    gCCGAGCGCCAGTGCTACCGCCATGATGATCCATACCGGTGCGTACTCAATGGTGCTTAACATATCGCt                                   <  1:2267903/69‑1 (MQ=255)                                    gCCGAGCGCCAGTGCTACCGCCATGATGATCCATACCGGTGCGTACTCAATGGTGCTTAACATATCGCt                                   <  1:2171386/69‑1 (MQ=255)                                                tGCTACCGCCATGATGATCCACACCGGAGCGTACTCAATGGTGCTTaa                                            <  1:3315516/48‑1 (MQ=255)                                                  cTACCTCCATGATGATCCACACCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTc                      >  1:1140235/1‑68 (MQ=255)                                                  cTACCGCCATGATGATCCACACCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTc                      >  1:1370624/1‑68 (MQ=255)                                                  cTACCGCCATGATGATCCACACCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTc                      >  1:655719/1‑68 (MQ=255)                                                  cTACCGCCATGATGATCCACACCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTc                      >  1:1368330/1‑68 (MQ=255)                                                  cTACCGCCATGATGATCCACACCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTc                      >  1:1014671/1‑68 (MQ=255)                                                  cTACCGCCATGATGATCCACACCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTc                      >  1:894333/1‑68 (MQ=255)                                                  cTACCGCCATGATGATCCACACCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTc                      >  1:2132977/1‑68 (MQ=255)                                                  cTACCGCCATGATGATCCACACCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTc                      >  1:2532430/1‑68 (MQ=255)                                                  cTACCGCCATGATGATCCACACCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTc                      >  1:2565646/1‑68 (MQ=255)                                                  cTACCGCCATGATGATCCACACCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTc                      >  1:1477793/1‑68 (MQ=255)                                                            gatgttCCATACCGGTGCGTACTCAATGGTGCTTAACATATCGCTGCGCAGTTTTTTCAGCAGGTTCTg           <  1:2014350/69‑1 (MQ=255)                                                            gatgatCCATACCGGTGCGTACTCAATGGTGCTTAACATATCGCTGCTCAGTTTTTTCAGCATGTTCTg           <  1:442564/69‑1 (MQ=255)                                                            gatgatCCATACCGGTGCGTACTCAATGGTGCTTAACATATCGCTGCGCAGTTTTTTCAGCAGGTTCTg           <  1:1025319/69‑1 (MQ=255)                                                            gatgatCCATACCGGTGCGTACTCAATGGTGCTTAACATATCGCTGCGCAGTTTTTTCAGCAGGTTCTg           <  1:353295/69‑1 (MQ=255)                                                            gatgatCCATACCGGTGCGTACTCAATGGTGCTTAACATATCGCTGCGCAGTTTTTTCAGCAGGTTCTg           <  1:48508/69‑1 (MQ=255)                                                            gatgatCCATACCGGTGCGTACTCAATGGTGCTTAACATATCGCTGCGCAGTTTTTTCAGCAGGTTCTg           <  1:2462889/69‑1 (MQ=255)                                                            gatgatCCATACCGGTGCGTACTCAATGGTGCTTAACATATCGCTGCGCAGTTTTTTCAGCAGGTTCTg           <  1:1114810/69‑1 (MQ=255)                                                            gatgatCCATACCGGTGCGTACTCAATGGTGCTTAACATATCGCTGCGCAGTTTTTTCAGCAGGTTCTg           <  1:1799254/69‑1 (MQ=255)                                                            gatgatCCATACCGGTGCGTACTCAATGGTGCTTAACATATCGCTGCGCAGTTTTTTCAGCAGGTTCTg           <  1:3408187/69‑1 (MQ=255)                                                            gatgatCCATACCGGTGCGTACTCAATGGTGCTTAACATATCGCTGCGCAGTTTTTTCAGCAGGTTCTg           <  1:2317958/69‑1 (MQ=255)                                                            gatgatCCATACCGGTGCGTACTCAATGGTGCTTAACATATCGCTGCGCAGTTTTTTCAGCAGGTTCTg           <  1:1861730/69‑1 (MQ=255)                                                            gatgatCCATACCGGTGCGTACTCAATGGTGCTTAACATATCGCTGCGCAGTTTTTTCAGCAGGTTCTg           <  1:1007555/69‑1 (MQ=255)                                                            gatgatCCATACCGGTGCGTACTCAATGGTGCTTAACATATCGCTGCGCAGTTTTTTCAGCAGGTTCTg           <  1:88830/69‑1 (MQ=255)                                                            gatgatCCATACCGGTGCGTACTCAATGGTGCTTAACATATCGCTGCGCAGTTTTTTCAGCAGGTTCTg           <  1:1005709/69‑1 (MQ=255)                                                            gatgatCCATACCGGTGCGTACTCAATGGTGCTTAACATATCGCTGCGCAGTTTTTTCAGCAGGTTCTg           <  1:3309791/69‑1 (MQ=255)                                                            gatgatCCATACCGGTGCGTACTCAATGGTGCTTAACATATCGCTGCGCAGTTTTTTCAGCAGGTTCTg           <  1:3268093/69‑1 (MQ=255)                                                            gatgatCCATACCGGTGCGTACTCAATGGTGCTTAACATATCGCTGCGCAGTTTTTTCAGCAGGTTCTg           <  1:2921292/69‑1 (MQ=255)                                                            gatgatCCATACCGGTGCGTACTCAATGGTGCTTAACATATCGCTGCGCAGTTTTTTCAGCAGGTTCTg           <  1:1388411/69‑1 (MQ=255)                                                            gatgatCCATACCGGTGCGTACTCAATGGTGCTTAACATATCGCTGCGCAGTTTTTTCAGCAGGTTCTg           <  1:1438247/69‑1 (MQ=255)                                                            gatgatCCATACCGGTGCGTACTCAATGGTGCTTAACATATCGCTGCGCAGTTTTATCAGCAGGTTCTg           <  1:3100970/69‑1 (MQ=255)                                                            gatgatCCATACCGGTGCATACTCAATGGTGCTTAACATATCGCTGCGCAGTTTTTTCAGCAGGTTCTg           <  1:1338695/69‑1 (MQ=25)                                                             atgatCCATACCGGTGCGTACTCAATGGTGCTTAACATATCGCTGCGCAGTTTTTTCAGCAGGTTCTg           <  1:2542767/68‑1 (MQ=255)                                                             atgatCCACACCGGAGCGTACTCAATGGTGCTTAACATATCGCt                                   >  1:3391361/1‑44 (MQ=255)                                                             atgatCCACACCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTgg          >  1:3191397/1‑69 (MQ=255)                                                             atgatCCACACCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTgg          >  1:929166/1‑69 (MQ=255)                                                             atgatCCACACCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTgg          >  1:1477179/1‑69 (MQ=255)                                                             atgatCCACACCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTgg          >  1:86337/1‑69 (MQ=255)                                                             atgatCCACACCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTgg          >  1:2804410/1‑69 (MQ=255)                                                             atgatCCACACCGGAGAGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTgg          >  1:1219306/1‑69 (MQ=255)                                                                     cacCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTcttctt   <  1:2366210/68‑1 (MQ=255)                                                                     cacCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTcttctt   <  1:3252081/68‑1 (MQ=255)                                                                     cacCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCttt  <  1:2869489/69‑1 (MQ=255)                                                                     |                                                                     ATCGCTACACGACGCCAGCCAATCATGGTGCCAATGCCGAGCGCCAGTGCTACCGCCATGATGATCCACACCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTT  >  minE/1997912‑1998048

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