Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	7,459	G→C	28.8%	A167A (GCC→GCG)	yaaJ ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	7,459	0	G	C	28.8%	66.0 / 30.7	52	A167A (GCC→GCG)	yaaJ	predicted transporter
Reads supporting (aligned to +/- strand):  ref base G (19/18);  new base C (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.36e-01

AGTAAAGACAGCGAGAATAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAG  >  minE/7393‑7527                                                                   |                                                                     agTAAAGACAGCGAGAATAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGc                                                                    <  1:2055533/69‑1 (MQ=255) agTAAAGACAGCGAGAATAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGc                                                                    <  1:1745951/69‑1 (MQ=255)   tAAAGACAGCGAGAATAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCaa                                                                  >  1:3254242/1‑69 (MQ=255)   tAAAGACAGCGAGAATAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCaa                                                                  >  1:3256733/1‑69 (MQ=255)   tAAAGACAGCGAGAATAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCaa                                                                  >  1:2466525/1‑69 (MQ=255)   tAAAGACAGCGAGAATAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCaa                                                                  >  1:2375650/1‑69 (MQ=255)   tAAAGACAGCGAGAATAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCaa                                                                  >  1:2097720/1‑69 (MQ=255)   tAAAGACAGCGAGAATAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCaa                                                                  >  1:1685297/1‑69 (MQ=255)   tAAAGACAGCGAGAATAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCaa                                                                  >  1:957541/1‑69 (MQ=255)   tAAAGACAGCGAGAATAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCaa                                                                  >  1:1291513/1‑69 (MQ=255)   tAAAGACAGCGAGAATAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCaa                                                                  >  1:1341071/1‑69 (MQ=255)   tAAAGACAGCGAGAATAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCaa                                                                  >  1:1162779/1‑69 (MQ=255)   tAAAGACAGCGAGAATAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCa                                                                   >  1:3141709/1‑68 (MQ=255)               gaATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGcc                                                       >  1:2809046/1‑67 (MQ=255)               gaATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGcc                                                       >  1:1872903/1‑67 (MQ=255)               gaATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGcc                                                       >  1:2746716/1‑67 (MQ=255)               gaATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGcc                                                       >  1:390109/1‑67 (MQ=255)               gaATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGcc                                                       >  1:364056/1‑67 (MQ=255)               gaATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGcc                                                       >  1:172994/1‑67 (MQ=255)               gaATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGcc                                                       >  1:3056560/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:197792/1‑69 (MQ=255)                    atGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACt                                                 >  1:34314/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:3119713/56‑1 (MQ=255)                             aCCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTg                                                     <  1:483692/56‑1 (MQ=255)                             aCCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTg                                                     <  1:2415785/56‑1 (MQ=255)                             aCCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTg                                                     <  1:58078/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:3074767/69‑1 (MQ=255)                                                     aaaCTCAGGGCTCGCGCAACGGCGTTCGCCTGAACTCCGCTGAAAATTATGGCGTAGGCGATGAGCaaa                <  1:679474/69‑1 (MQ=255)                                                     aaaCTCAGGGCTCGCGCAACGGCGTTCGCCTGAACTCCGCTGAAAATTATGGCGTAGGCGATGAGCaaa                <  1:2144471/69‑1 (MQ=255)                                                     aaaCTCAGGGCTCGCGCAACGGCGTTCGCCTGAACTCCGCTGAAAATTATGGCGTAGGCGATGAGCaaa                <  1:2231443/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)                                                                   ggCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAg  <  1:1409568/69‑1 (MQ=255)                                                                   ggCAACGGCGTTCGCTTGAACTCCCCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAg  <  1:2131063/69‑1 (MQ=255)                                                                   |                                                                     AGTAAAGACAGCGAGAATAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAG  >  minE/7393‑7527

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