Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,393,281	C→T	100%	E157E (GAG→GAA)	pldA ←	outer membrane phospholipase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,393,281	0	C	T	83.5%	44.5 / 2.0	30	E157E (GAG→GAA)	pldA	outer membrane phospholipase A
Reads supporting (aligned to +/- strand):  ref base C (4/1);  new base T (8/17);  total (12/18)
Fisher's exact test for biased strand distribution p-value = 1.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AAGGCGGTTCCAGCTGCGGGAGGTCGGGTCGGAACGCCCGTTAGAGTCGTGGTTATACCCCATCTCCACATCGCGCAGCGTCCAACCTGCAAAACGGTAATCGGTGGCAAAACCGAGGA  >  minE/2393218‑2393336                                                                |                                                        aaGGCGGTTCCAGCTGCGGGAGGTCGGGTCGGAACGCCCGTTAGAGTCGTGGTTATACCCCATCTCca                                                     <  1:1640940/68‑1 (MQ=255)            aGCTGCGGGAAGTCGGGTCTGAACGCCCGTTTGAGTCGTGGTTGTAGCCCATTTcc                                                      >  1:686892/1‑56 (MQ=255)            aGCTGCGGGAAGTCGGGTCTGAACGCCCGTTTGAGTCGTGGTTGTAGCCCATTTCCACATCgcgc                                             >  1:564188/1‑65 (MQ=255)            aGCTGCGGGAAGTCGGGTCTGAACGCCCGTTTGAGTCGTGGTTGTAGCCCATTTCCACATCGCGCAGCg                                         >  1:2607814/1‑69 (MQ=255)            aGCTGCGGGAAGTCGGGTCTGAACGCCCGTTTGAGTCGTGGTTGTAGCCCATTTCCACATCGCGCAGCg                                         >  1:975535/1‑69 (MQ=255)            aGCTGCGGGAAGTCGGGTCTGAACGCCCGTTTGAGTCGTGGTTGTAGCCCATTTCCACATCGCGCAGCg                                         >  1:954261/1‑69 (MQ=255)            aGCTGCGGGAAGTCGGGTCTGAACGCCCGTTTGAGTCGTGGTTGTAGCCCATTTCCACATCGCGCAGCg                                         >  1:1042924/1‑69 (MQ=255)            aGCTGCGGGAAGTCGGGTCTGAACGCCCGTTTGAGTCGTGGTTGTAGCCCATTTCCACATCGCGCAGCg                                         >  1:1933345/1‑69 (MQ=255)            aGCTGCGGGAAGTCGGGTCTGAACGCCCGTTTGAGTCGTGGTTGTAGCCCATTTCCACATCGCGCAGCg                                         >  1:1733921/1‑69 (MQ=255)                                           tgagTCGTGGTTGTAGCCCATTTCCACATCGCGCAGCGTCCAGCCAGCAAAGTTGTAATCCGTGGCaaa          <  1:1989989/68‑1 (MQ=255)                                           tgagTCGTGGTTGTAGCCCATTTCCACATCGCGCAGCGTCCAGCCAGCAAAGTTGTAATCCGTGGCaaa          <  1:1237532/68‑1 (MQ=255)                                           tgagTCGTGGTTGTAGCCCATTTCCACATCGCGCAGCGTCCAGCCAGCAAAGTTGTAATCCGTGGCaaa          <  1:1457671/68‑1 (MQ=255)                                           tgagTCGTGGTTGTAGCCCATTTCCACATCGCGCAGCGTCCAGCCAGCAAAGTTGTAATCCGTGGCaaa          <  1:709688/68‑1 (MQ=255)                                           tgagTCGTGGTTGTAGCCCATTTCCACATCGCGCAGCGTCCAGCCAGCAAAGTTGTAATCCGTGGCaaa          <  1:693219/68‑1 (MQ=255)                                           tgagTCGTGGTTGTAGCCCATTTCCACATCGCGCAGCGTCCAGCCAGCAAAGTTGTAATCCGTGGCaaa          <  1:3345300/68‑1 (MQ=255)                                           tgagTCGTGGTTGTAGCCCATTTCCACATCGCGCAGCGTCCAGCCAGCAAAGTTGTAATCCGTGGCaaa          <  1:2865636/68‑1 (MQ=255)                                           tgagTCGTGGTTGTAGCCCATTTCCACATCGCGCAGCGTCCAGCCAGCAAAGTTGTAATCCGTGGCaaa          <  1:2603235/68‑1 (MQ=255)                                           tgagTCGTGGTTGTAGCCCATTTCCACATCGCGCAGCGTCCAGCCAGCAAAGTTGTAATCCGTGGCaaa          <  1:2521861/68‑1 (MQ=255)                                            gagTCGTGGTTGTAGCCCATTTCCACATCGCGCAGCGTCCAGCCAGCAAAGTTGTAATCCGTGGCaaa          <  1:859031/68‑1 (MQ=255)                                               tCGTGGTTGTAGCCCATTTCCACATCGCGCAGCGTCCAGCCAGCAAAGTTGTAATCCGTGGCAAAAcc       <  1:2627825/68‑1 (MQ=255)                                               tCGTGGTTGTAGCCCATTTCCACATCGCGCAGCGTCCAGCCAGCAAAGTTGTAATCCGTGGCAAAAcc       <  1:1732805/68‑1 (MQ=255)                                               tCGTGGTTGTAGCCCATTTCCACATCGCGCAGCGTCCAGCCAGCAAAGTTGTAATCCGTGGCAAAAcc       <  1:3096281/68‑1 (MQ=255)                                               tCGTGGTTGTAGCCCATTTCCACATCGCGCAGCGTCCAGCCAGCAAAGTTGTAATCCGTGGCAAAAcc       <  1:3100300/68‑1 (MQ=255)                                               tCGTGGTTGTAGCCCATTTCCACATCGCGCAGCGTCCAGCCAGCAAAGTTGTAATCCGTGGCAAAAcc       <  1:1714971/68‑1 (MQ=255)                                               tCGTGGTTGTAGCCCATTTCCACATCGCGCAGCGTCCAGCCAGCAAAGTTGTAATCCGTGGCAAAAcc       <  1:1628669/68‑1 (MQ=255)                                               tCGTGGTTGTAGCCCATTTCCACATCGCGCAGCGTCCAGCCAGCAAAGTTGTAATCCGTGGCAAAAcc       <  1:2598970/68‑1 (MQ=255)                                                    gTTATACCCCATCTCCACATCGCGCAGCGTCCAACCTGCAAAACGGTAATCGGTGGCAAAACCGAGGa  >  1:3276127/1‑68 (MQ=255)                                                    gTTATACCCCATCTCCACATCGCGCAGCGTCCAACCTGCAAAACGGTAATCGGTGGCAAAACCGAGGa  >  1:3387111/1‑68 (MQ=255)                                                    gTTATACCCCATCTCCACATCGCGCAGCGTCCAACCTGCAAAACGGTAATCGGTGGCAAAACCGAGGa  >  1:2254978/1‑68 (MQ=255)                                                    gTTATACCCCATCTCCACATCGCGCAGCGTCCAACCTGCAAAACGGTAATCGGTGGCAAAACCGAGGa  >  1:1993437/1‑68 (MQ=255)                                                                |                                                        AAGGCGGTTCCAGCTGCGGGAGGTCGGGTCGGAACGCCCGTTAGAGTCGTGGTTATACCCCATCTCCACATCGCGCAGCGTCCAACCTGCAAAACGGTAATCGGTGGCAAAACCGAGGA  >  minE/2393218‑2393336

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