Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,187,797	G→C	27.2%	R20R (CGC→CGG)	nanT ←	sialic acid transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,187,797	0	G	C	27.2%	94.6 / 37.2	66	R20R (CGC→CGG)	nanT	sialic acid transporter
Reads supporting (aligned to +/- strand):  ref base G (22/26);  new base C (8/10);  total (30/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.98e-01

CCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACTCAT  >  minE/2187733‑2187856                                                                 |                                                            ccAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGcc                                                         >  1:2164307/1‑69 (MQ=255) ccAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGcc                                                         >  1:2070469/1‑69 (MQ=255)   aGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCAt                                                       >  1:3146591/1‑69 (MQ=255)   aGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCAt                                                       >  1:431576/1‑69 (MQ=255)   aGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCAt                                                       >  1:3415450/1‑69 (MQ=255)   aGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCAt                                                       >  1:2397954/1‑69 (MQ=255)   aGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCAt                                                       >  1:2758262/1‑69 (MQ=255)   aGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCAt                                                       >  1:1987841/1‑69 (MQ=255)   aGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCAt                                                       >  1:2786052/1‑69 (MQ=255)   aGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCAt                                                       >  1:1760634/1‑69 (MQ=255)   aGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCAt                                                       >  1:1242759/1‑69 (MQ=255)   aGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCAt                                                       >  1:3250256/1‑69 (MQ=255)   aGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCAt                                                       >  1:1409159/1‑69 (MQ=255)   aGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCAt                                                       >  1:171103/1‑69 (MQ=255)      gCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATtg                                                     >  1:1724221/1‑68 (MQ=255)       cGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATtgtg                                                   <  1:2234036/69‑1 (MQ=255)            aaaaCGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACgg                                              >  1:1407769/1‑69 (MQ=255)            aaaaCGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACgg                                              >  1:319546/1‑69 (MQ=255)            aaaaCGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACgg                                              >  1:524967/1‑69 (MQ=255)            aaaaCGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACgg                                              >  1:700019/1‑69 (MQ=255)            aaaaCGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACgg                                              >  1:755855/1‑69 (MQ=255)                cGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTga                                          <  1:3342317/69‑1 (MQ=255)                cGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTga                                          <  1:194453/69‑1 (MQ=255)                cGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTga                                          <  1:2790611/69‑1 (MQ=255)                cGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTga                                          <  1:3049398/69‑1 (MQ=255)                cGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTga                                          <  1:2675798/69‑1 (MQ=255)                cGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTga                                          <  1:2346047/69‑1 (MQ=255)                cGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTga                                          <  1:185497/69‑1 (MQ=255)                cGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTga                                          <  1:238780/69‑1 (MQ=255)                cGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTga                                          <  1:825191/69‑1 (MQ=255)                cGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTga                                          <  1:2665236/69‑1 (MQ=255)                cGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTga                                          <  1:1611631/69‑1 (MQ=255)                      cAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGc                                    <  1:1147196/69‑1 (MQ=255)                         aaCCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGCGa                                  <  1:3375760/68‑1 (MQ=255)                         aaCCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGCGa                                  <  1:1882090/68‑1 (MQ=255)                                     gATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGCGATACCACGGGatat                     <  1:2626913/69‑1 (MQ=255)                                     gATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGCGATACCACGGGatat                     <  1:2646313/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:844700/69‑1 (MQ=255)                                                     aGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACt     <  1:3025876/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:1362843/69‑1 (MQ=255)                                                     aGCGGAAAATGCCCGCCATTGTGCGCGGTTAAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACt     <  1:2524875/69‑1 (MQ=255)                                                     aGCGGAAAATGCCCGCCATTGTGCGCGGTTAAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACt     <  1:2513549/69‑1 (MQ=255)                                                     aGCGGAAAATGCCCGCCATTGTGCGCGGTTAAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACt     <  1:1423826/69‑1 (MQ=255)                                                     aGCGGAAAATGCCCGCCATTGTGCGCGGTTAAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACt     <  1:2434710/69‑1 (MQ=255)                                                     aGCGGAAAATGCCCGCCATTGTGCGCGGTTAAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACt     <  1:1791290/69‑1 (MQ=255)                                                     aGCGGAAAATGCCCGCCATTGTGCGCGGTTAAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACt     <  1:2036585/69‑1 (MQ=255)                                                      gCGGAAAATGCCCGCCATTGTGCGCGGTTAAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACt     <  1:573486/68‑1 (MQ=255)                                                        ggAAAATGCCCGCCATTGTGCGCGGTTAAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACTCa   >  1:972614/1‑68 (MQ=255)                                                        ggAAAATGCCCGCCATTGTGCGCGGTTAAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACTCAt  >  1:1127655/1‑69 (MQ=255)                                                        ggAAAATGCCCGCCATTGTGCGCGGTTAAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACTCAt  >  1:3130235/1‑69 (MQ=255)                                                        ggAAAATGCCCGCCATTGTGCGCGGTTAAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACTCAt  >  1:2010899/1‑69 (MQ=255)                                                        ggAAAATGCCCGCCATTGTGCGCGGTTAAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACTCAt  >  1:2342137/1‑69 (MQ=255)                                                        ggAAAATGCCCGCCATTGTGCGCGGTTAAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACTCAt  >  1:2442623/1‑69 (MQ=255)                                                        ggAAAATGCCCGCCATTGTGCGCGGTTAAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACTCAt  >  1:2624744/1‑69 (MQ=255)                                                        ggAAAATGCCCGCCATTGTGCGCGGTTAAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACTCAt  >  1:294650/1‑69 (MQ=255)                                                            aaTGCGCGCCATTGTGCACGGTTGAGATGGCGATACCACGGGa                        <  1:2386160/43‑1 (MQ=255)                                                            aaTGCGCGCCATTGTGCACGGTTGAGATGGCGATACCACGGGa                        <  1:277017/43‑1 (MQ=255)                                                                 |                                                            CCAGGGCGATTAAAACGAAATCAAAACCGTCAAGCAGATATCCCAACCAGGCAGCGGAAAATGCGCGCCATTGTGCACGGTTGAGATGGCGATACCACGGGATATTCTGGGTTGTAGTACTCAT  >  minE/2187733‑2187856

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