Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	7,456	G→T	28.8%	R168R (CGC→CGA)	yaaJ ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	7,456	0	G	T	28.8%	70.0 / 26.8	52	R168R (CGC→CGA)	yaaJ	predicted transporter
Reads supporting (aligned to +/- strand):  ref base G (19/18);  new base T (8/7);  total (27/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

AGTAAAGACAGCGAGAATAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGAC  >  minE/7393‑7518                                                                |                                                               agTAAAGACAGCGAGAATAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGc                                                           <  1:2055533/69‑1 (MQ=255) agTAAAGACAGCGAGAATAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGc                                                           <  1:1745951/69‑1 (MQ=255)   tAAAGACAGCGAGAATAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCaa                                                         >  1:1685297/1‑69 (MQ=255)   tAAAGACAGCGAGAATAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCaa                                                         >  1:2097720/1‑69 (MQ=255)   tAAAGACAGCGAGAATAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCaa                                                         >  1:3256733/1‑69 (MQ=255)   tAAAGACAGCGAGAATAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCaa                                                         >  1:3254242/1‑69 (MQ=255)   tAAAGACAGCGAGAATAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCaa                                                         >  1:2375650/1‑69 (MQ=255)   tAAAGACAGCGAGAATAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCaa                                                         >  1:2466525/1‑69 (MQ=255)   tAAAGACAGCGAGAATAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCaa                                                         >  1:1162779/1‑69 (MQ=255)   tAAAGACAGCGAGAATAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCaa                                                         >  1:1341071/1‑69 (MQ=255)   tAAAGACAGCGAGAATAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCaa                                                         >  1:1291513/1‑69 (MQ=255)   tAAAGACAGCGAGAATAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCaa                                                         >  1:957541/1‑69 (MQ=255)   tAAAGACAGCGAGAATAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCa                                                          >  1:3141709/1‑68 (MQ=255)               gaATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGcc                                              >  1:2809046/1‑67 (MQ=255)               gaATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGcc                                              >  1:3056560/1‑67 (MQ=255)               gaATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGcc                                              >  1:2746716/1‑67 (MQ=255)               gaATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGcc                                              >  1:364056/1‑67 (MQ=255)               gaATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGcc                                              >  1:390109/1‑67 (MQ=255)               gaATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGcc                                              >  1:1872903/1‑67 (MQ=255)               gaATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGcc                                              >  1:172994/1‑67 (MQ=255)               gaATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGcc                                              >  1:1622251/1‑67 (MQ=255)                aataatGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTg                                            <  1:456271/69‑1 (MQ=255)                aataatGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTg                                            <  1:95376/69‑1 (MQ=255)                    atGCCTGTCCCCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAAc                                         >  1:1503811/1‑68 (MQ=255)                    atGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACt                                        >  1:34314/1‑69 (MQ=255)                    atGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACt                                        >  1:197792/1‑69 (MQ=255)                      gCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTcc                                      >  1:2151620/1‑69 (MQ=255)                      gCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTcc                                      >  1:641618/1‑69 (MQ=255)                            cACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGaa                                >  1:2714375/1‑69 (MQ=255)                            cACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGaa                                >  1:1853294/1‑69 (MQ=255)                             aCCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTg                                            <  1:58078/56‑1 (MQ=255)                             aCCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTg                                            <  1:483692/56‑1 (MQ=255)                             aCCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTg                                            <  1:3119713/56‑1 (MQ=255)                             aCCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTg                                            <  1:2415785/56‑1 (MQ=255)                             aCCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGaa                                >  1:3105866/1‑68 (MQ=255)                               cAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAAt                             <  1:1223678/69‑1 (MQ=255)                               cAGCGGGGGAAAATCAAAAGAAAAACTCAGGACGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAAt                             <  1:2612660/69‑1 (MQ=255)                                aGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAAt                             <  1:80123/68‑1 (MQ=255)                                    ggggAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGc                        <  1:2000316/69‑1 (MQ=255)                                             cAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTg                                            <  1:3099417/40‑1 (MQ=255)                                                     aaaCTCAGGGCTCGCGCAACGGCGTTCGCCTGAACTCCGCTGAAAATTATGGCGTAGGCGATGAGCaaa       <  1:2231443/69‑1 (MQ=255)                                                     aaaCTCAGGGCTCGCGCAACGGCGTTCGCCTGAACTCCGCTGAAAATTATGGCGTAGGCGATGAGCaaa       <  1:2144471/69‑1 (MQ=255)                                                     aaaCTCAGGGCTCGCGCAACGGCGTTCGCCTGAACTCCGCTGAAAATTATGGCGTAGGCGATGAGCaaa       <  1:3074767/69‑1 (MQ=255)                                                     aaaCTCAGGGCTCGCGCAACGGCGTTCGCCTGAACTCCGCTGAAAATTATGGCGTAGGCGATGAGCaaa       <  1:679474/69‑1 (MQ=255)                                                      aaCTCAGGGCTCGCGCAACGGCGTTCGCCTGAACTCCGCTGAAAATTATGGCGTAGGCGATGAGCaaa       <  1:2216594/68‑1 (MQ=255)                                                      aaCTCAGGGCTCGCGCAACGGCGTTCGCCTGAACTCCGCTGAAAATTATGGCGTAGGCGATGAGCaaa       <  1:901386/68‑1 (MQ=255)                                                      aaCTCAGGGCTCGCGCAACGGCGTTCGCCTGAACTCCGCTGAAAATTATGGCGTAGGCGATGAGCaaa       <  1:1511347/68‑1 (MQ=255)                                                        cTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAg    <  1:2900651/69‑1 (MQ=255)                                                          cAGGGCGCTGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGAc  <  1:740452/69‑1 (MQ=255)                                                          cAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGAc  <  1:44685/69‑1 (MQ=255)                                                          cAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGAc  <  1:2372362/69‑1 (MQ=255)                                                          cAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGAc  <  1:164911/69‑1 (MQ=255)                                                                |                                                               AGTAAAGACAGCGAGAATAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGAC  >  minE/7393‑7518

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