Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,187,815	G→A	51.4%	L14L (CTC→CTT)	nanT ←	sialic acid transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,187,815	0	G	A	51.4%	‑5.6 / 29.4	35	L14L (CTC→CTT)	nanT	sialic acid transporter
Reads supporting (aligned to +/- strand):  ref base G (2/15);  new base A (8/10);  total (10/25)
Fisher's exact test for biased strand distribution p-value = 5.99e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACTCATTGTGAGTCTCCCGCGGTGGGCGA  >  minE/2187748‑2187879                                                                    |                                                                 cGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTga                                                                 <  1:2665236/69‑1 (MQ=255) cGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTga                                                                 <  1:2675798/69‑1 (MQ=255) cGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTga                                                                 <  1:2790611/69‑1 (MQ=255) cGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTga                                                                 <  1:1611631/69‑1 (MQ=255) cGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTga                                                                 <  1:825191/69‑1 (MQ=255) cGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTga                                                                 <  1:185497/69‑1 (MQ=255) cGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTga                                                                 <  1:194453/69‑1 (MQ=255) cGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTga                                                                 <  1:3049398/69‑1 (MQ=255) cGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTga                                                                 <  1:3342317/69‑1 (MQ=255) cGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTga                                                                 <  1:238780/69‑1 (MQ=255) cGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTga                                                                 <  1:2346047/69‑1 (MQ=255)       cAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGc                                                           <  1:1147196/69‑1 (MQ=255)          aaCCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGCGa                                                         <  1:1882090/68‑1 (MQ=255)          aaCCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGCGa                                                         <  1:3375760/68‑1 (MQ=255)                      gATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGCGATACCACGGGatat                                            <  1:2646313/69‑1 (MQ=255)                      gATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGCGATACCACGGGatat                                            <  1:2626913/69‑1 (MQ=255)                         atCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGCGATACCACGGGATATTCt                                         >  1:469169/1‑69 (MQ=255)                         atCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGCGATACCACGGGATATTCt                                         >  1:130087/1‑69 (MQ=255)                               aCCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGCGATACCACGGGATATTCTGGGTTg                                   <  1:2038854/69‑1 (MQ=255)                                ccAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGCGATAGCACGGGATATTCTGGGTTgt                                  <  1:588187/69‑1 (MQ=255)                                      aGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACt                            <  1:468024/69‑1 (MQ=255)                                      aGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACt                            <  1:2286436/69‑1 (MQ=255)                                      aGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACt                            <  1:3025876/69‑1 (MQ=255)                                      aGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACt                            <  1:844700/69‑1 (MQ=255)                                      aGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACt                            <  1:1151458/69‑1 (MQ=255)                                      aGCGGAAAATGCCCGCCATTGTGCGCGGTTAAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACt                            <  1:825250/69‑1 (MQ=255)                                      aGCGGAAAATGCCCGCCATTGTGCGCGGTTAAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACt                            <  1:97576/69‑1 (MQ=255)                                      aGCGGAAAATGCCCGCCATTGTGCGCGGTTAAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACt                            <  1:2434710/69‑1 (MQ=255)                                      aGCGGAAAATGCCCGCCATTGTGCGCGGTTAAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACt                            <  1:1362843/69‑1 (MQ=255)                                      aGCGGAAAATGCCCGCCATTGTGCGCGGTTAAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACt                            <  1:1423826/69‑1 (MQ=255)                                      aGCGGAAAATGCCCGCCATTGTGCGCGGTTAAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACt                            <  1:1791290/69‑1 (MQ=255)                                      aGCGGAAAATGCCCGCCATTGTGCGCGGTTAAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACt                            <  1:2036585/69‑1 (MQ=255)                                      aGCGGAAAATGCCCGCCATTGTGCGCGGTTAAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACt                            <  1:2513549/69‑1 (MQ=255)                                      aGCGGAAAATGCCCGCCATTGTGCGCGGTTAAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACt                            <  1:2524875/69‑1 (MQ=255)                                       gCGGAAAATGCCCGCCATTGTGCGCGGTTAAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACt                            <  1:573486/68‑1 (MQ=255)                                         ggAAAATGCCCGCCATTGTGCGCGGTTAAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACTCa                          >  1:972614/1‑68 (MQ=255)                                         ggAAAATGCCCGCCATTGTGCGCGGTTAAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACTCAt                         >  1:2442623/1‑69 (MQ=255)                                         ggAAAATGCCCGCCATTGTGCGCGGTTAAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACTCAt                         >  1:3130235/1‑69 (MQ=255)                                         ggAAAATGCCCGCCATTGTGCGCGGTTAAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACTCAt                         >  1:2624744/1‑69 (MQ=255)                                         ggAAAATGCCCGCCATTGTGCGCGGTTAAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACTCAt                         >  1:2342137/1‑69 (MQ=255)                                         ggAAAATGCCCGCCATTGTGCGCGGTTAAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACTCAt                         >  1:294650/1‑69 (MQ=255)                                         ggAAAATGCCCGCCATTGTGCGCGGTTAAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACTCAt                         >  1:2010899/1‑69 (MQ=255)                                         ggAAAATGCCCGCCATTGTGCGCGGTTAAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACTCAt                         >  1:1127655/1‑69 (MQ=255)                                             aaTGCGCGCCATTGTGCACGGTTGAGATGGCGATACCACGGGa                                               <  1:2386160/43‑1 (MQ=255)                                             aaTGCGCGCCATTGTGCACGGTTGAGATGGCGATACCACGGGa                                               <  1:277017/43‑1 (MQ=255)                                                                ggTTGAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACTCATTGTGAGTCTCCCGCGGTGGGCGa  <  1:2587859/69‑1 (MQ=255)                                                                    |                                                                 CGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACTCATTGTGAGTCTCCCGCGGTGGGCGA  >  minE/2187748‑2187879

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