Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,871,051	G→A	42.0%	L326L (CTG→CTA) ‡	yqcE →	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,871,051	0	G	A	42.0%	33.0 / 63.1	88	L326L (CTG→CTA) ‡	yqcE	predicted transporter
Reads supporting (aligned to +/- strand):  ref base G (27/24);  new base A (22/15);  total (49/39)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

AACTGCCCTGCTCGTCACGAACTCTAACCCGCAATCGGTCGCGATGGGGATTGGCCTGATTTTACTGCTGGGATTCACCTGTTACGCCTCACGCGGGCTGTACTGGGCCTGCCCTGGCGAAGCGAG  >  minE/1870985‑1871110                                                                   |                                                            aaCTGCCCTGCTCGTCACGAACTCTAACCCGCAATCGGTCGCGATGGGGATTGGCCTGATTTTActgct                                                           >  1:1454303/1‑69 (MQ=255) aaCTGCCCTGCTCGTCACGAACTCTAACCCGCAATCGGTCGCGATGGGGATTGGCCTGATTTTActgct                                                           >  1:288806/1‑69 (MQ=255) aaCTGCCCTGCTCGTCACGAACTCTAACCCGCAATCGGTCGCGATGGGGATTGGCCTGATTTTActgct                                                           >  1:1972760/1‑69 (MQ=255) aaCTGCCCTGCTCGTCACGAACTCTAACCCGCAATCGGTCGCGATGGGGATTGGCCTGATTTTActgct                                                           >  1:2004141/1‑69 (MQ=255) aaCTGCCCTGCTCGTCACGAACTCTAACCCGCAATCGGTCGCGATGGGGATTGGCCTGATTTTActgct                                                           >  1:2156358/1‑69 (MQ=255)        cTGCTCGTCACGAACTCTAACCCGCAATCGGTCGCGATGGGGATTGGCCTGATTTTACTGCTGGGATTc                                                    <  1:3160353/69‑1 (MQ=255)        cTGCTCGTCACGAACTCTAACCCGCAATCGGTCGCGATGGGGATTGGCCTGATTTTACTGCTGGGATTc                                                    <  1:745815/69‑1 (MQ=255)           cTCGTCACGAACTCTAACCCGCAATCGGTCGCGATGGGGATTGGCCTGATTTTACTGCTGGGATTCAcc                                                 >  1:2641776/1‑69 (MQ=255)           cTCGTCACGAACTCTAACCCGCAATCGGTCGCGATGGGGATTGGCCTGATTTTACTGCTGGGATTCAcc                                                 >  1:1664611/1‑69 (MQ=255)           cTCGTCACGAACTCTAACCCGCAATCGGTCGCGATGGGGATTGGCCTGATTTTACTGCTGGGATTCAcc                                                 >  1:1572049/1‑69 (MQ=255)           cTCGTCACGAACTCTAACCCGCAATCGGTCGCGATGGGGATTGGCCTGATTTTACTGCTGGGATTCAcc                                                 >  1:3019552/1‑69 (MQ=255)           cTCGTCACGAACTCTAACCCGCAATCGGTCGCGATGGGGATTGGCCTGATTTTACTGCTGGGATTCAcc                                                 >  1:2946259/1‑69 (MQ=255)           cTCGTCACGAACTCTAACCCGCAATCGGTCGCGATGGGGATTGGCCTGATTTTACTGCTGGGATTCAcc                                                 >  1:1308941/1‑69 (MQ=255)           cTCGTCACGAACTCTAACCCGCAATCGGTCGCGATGGGGATTGGCCTGATTTTACTGCTGGGATTCAcc                                                 >  1:1305017/1‑69 (MQ=255)           cTCGTCACGAACTCTAACCCGCAATCGGTCGCGATGGGGATTGGCCTGATTTTACTGCTGGGATTCAcc                                                 >  1:2364725/1‑69 (MQ=255)           cTCGTCACGAACTCTAACCCGCAATCGGTCGCGATGGGGATTGGCCTGATTTTACTGCTGGGATTCAcc                                                 >  1:496484/1‑69 (MQ=255)           cTCGTCACGAACTCTAACCCGCAATCGGTCGCGATGGGGATTGGCCTGATTTTACTGCTGGGATTCAcc                                                 >  1:1173558/1‑69 (MQ=255)           cTCGTCACGAACTCTAACCAGCAATCGGTCGCGATGGGGATTGGCCTGATTTTACTGCTGGGATTCAcc                                                 >  1:1995421/1‑69 (MQ=255)                 aCGAACTCTAACCCGCAATCGGTCGCGATGGGGATTGGCCTGATTTTACTGCTGGGATTCACCTGTTAc                                           >  1:18466/1‑69 (MQ=255)                 aCGAACTCTAACCCGCAATCGGTCGCGATGGGGATTGGCCTGATTTTACTGCTGGGATTCACCTGTTAc                                           >  1:2937588/1‑69 (MQ=255)                 aCGAACTCTAACCCGCAATCGGTCGCGATGGGGATTGGCCTGATTTTACTGCTGGGATTCACCTGTTAc                                           >  1:1107903/1‑69 (MQ=255)                 aCGAACTCTAACCCGCAATCGGTCGCGATGGGGATTGGCCTGATTTTACTGCTGGGATTAACCTGTTAc                                           >  1:1271498/1‑69 (MQ=255)                     aCTCTAACCCGCAATCGGTCGCGATGGGGATTGGCCTGATTTTACTGCTGTGATTCACCTGTTACGCCt                                       <  1:2412368/69‑1 (MQ=255)                     aCTCTAACCCGCAATCGGTCGCGATGGGGATTGGCCTGATTTTACTGCTGGGATTCACCTGTTACGCCt                                       <  1:2597350/69‑1 (MQ=255)                     aCTCTAACCCGCAATCGGTCGCGATGGGGATTGGCCTGATTTTACTGCTGGGATTCACCTGTTACGCCt                                       <  1:1238705/69‑1 (MQ=255)                     aCTCTAACCCGCAATCGGTCGCGATGGGGATTGGCCTGATTTTACTGCTGGGATTCACCTGTTACGCCt                                       <  1:2103777/69‑1 (MQ=255)                             ccGCAATCGGTCGCGATGGGGATTGGCCTGATTTTATTACTGGGa                                                       >  1:43302/1‑45 (MQ=255)                             ccGCAATCGGTCGCGATGGGGATTGGCCTGATTTTATTACTGGGCTTCACCTGTTa                                            >  1:1339611/1‑56 (MQ=255)                             ccGCAATCGGTCGCGATGGGGATTGGCCTGATTTTATTACTGGGAtt                                                     >  1:396669/1‑47 (MQ=255)                             ccGCAATCGGTCGCGATGGGGATTGGCCTGATTTTATTACTGGGATTCACCTGTTa                                            >  1:476824/1‑56 (MQ=255)                             ccGCAATCGGTCGCGATGGGGATTGGCCTGATTTTATTACTGGGATTCACCTGTTAc                                           >  1:1652904/1‑57 (MQ=255)                             ccGCAATCGGTCGCGATGGGGATTGGCCTGATTTTATTACTGGGATTCACCTGTTAc                                           >  1:2372790/1‑57 (MQ=255)                             ccGCAATCGGTCGCGATGGGGATTGGCCTGATTTTATTACTGGGATTCACCTGTTAc                                           >  1:731439/1‑57 (MQ=255)                             ccGCAATCGGTCGCGATGGGGATTGGCCTGATTTTATTACTGGGATTCACCTGTTAc                                           >  1:2297331/1‑57 (MQ=255)                             ccGCAATCGGTCGCGATGGGGATTGGCCTGATTTTATTACTGGGATTCACCTGTTAc                                           >  1:227586/1‑57 (MQ=255)                             ccGCAATCGGTCGCGATGGGGATTGGCCTGATTTTATTACTGGGATTCACCTGTTAc                                           >  1:1120072/1‑57 (MQ=255)                             ccGCAATCGGTCGCGATGGGGATTGGCCTGATTTTATTACTGGGATTCACCTGTTAc                                           >  1:271510/1‑57 (MQ=255)                             ccGCAATCGGTCGCGATGGGGATTGGCCTGATTTTATTACTGGGATTCACCTGTTAc                                           >  1:456917/1‑57 (MQ=255)                             ccGCAATCGGTCGCGATGGGGATTGGCCTGATTTTATTACTGGGATTCACCTGTTAc                                           >  1:2839407/1‑57 (MQ=255)                             ccGCAATCGGTCGCGATGGGGATTGGCCTGATTTTATTACTGGGATTCACCTGTTAc                                           >  1:440607/1‑57 (MQ=255)                             ccGCAATCGGTCGCGATGGGGATTGGCCTGATTTTATTACTGGGATTCACCTGTTAc                                           >  1:2868344/1‑57 (MQ=255)                             ccGCAATCGGTCGCGATGGGGATTGGCCTGATTTTATTACTGGGATTCACCTGTTAc                                           >  1:1899940/1‑57 (MQ=255)                             ccGCAATCGGTCGCGATGGGGATTGGCCTGATTTTATTACTGGGATTCACCTGTTAc                                           >  1:1302683/1‑57 (MQ=255)                             ccGCAATCGGTCGCGATGGGGATTGGCCTGATTTTATTACTGGGATTCACCTGTTAc                                           >  1:1779454/1‑57 (MQ=255)                             ccGCAATCGGTCGCGATGGGGATTGGCCTGATTTTATTACTGGGATTCACCTGTTAc                                           >  1:3239649/1‑57 (MQ=255)                             ccGCAATCGGTCGCGATGGGGATTGGCCTGATTTTATTACTGGGATTCACCTGTTAc                                           >  1:163136/1‑57 (MQ=255)                             ccGCAATCGGTCGCGATGGGGATTGGCCTGATTTTATTACTGGGATTCACCTGTTAc                                           >  1:3074515/1‑57 (MQ=255)                             ccGCAATCGGTCGCGATGGGGATTGGCCTGATTTTATTACAGGGATTCACCTGTTa                                            >  1:1031613/1‑56 (MQ=255)                                cAATCGGTCGCGATGGGGATTGGCCTGATTTTACTGCTGGGATTCACCTGTTACGCCTCACGCGGGCTg                            <  1:3099580/69‑1 (MQ=255)                                cAATCGGTCGCGATGGGGATTGGCCTGATTTTACTGCTGGGATTCACCTGTTACGCCTCACGCGGGCTg                            <  1:1454582/69‑1 (MQ=255)                                cAATCGGTCGCGATGGGGATTGGCCTGATTTTACTGCTGGGATTCACCTGTTACGCCTCACGCGGGCTg                            <  1:3395576/69‑1 (MQ=255)                                        cgcgATGGGGATTGGCCTGATTTTATTACTGGGATTCACCTGTTACGCCTCACGCTGGCTGTATTGGGc                    <  1:2852027/69‑1 (MQ=255)                                        cgcgATGGGGATTGGCCTGATTTTATTACTGGGATTCACCTGTTACGCCTCACGCGGGCTGTATTGGGc                    <  1:472904/69‑1 (MQ=255)                                        cgcgATGGGGATTGGCCTGATTTTATTACTGGGATTCACCTGTTACGCCTCACGCGGGCTGTATTGGGc                    <  1:3058456/69‑1 (MQ=255)                                        cgcgATGGGGATTGGCCTGATTTTATTACTGGGATTCACCTGTTACGCCTCACGCGGGCTGTATTGGGc                    <  1:3106797/69‑1 (MQ=255)                                        cgcgATGGGGATTGGCCTGATTTTATTACTGGGATTCACCTGTTACGCCTCACGCGGGCTGTATTGGGc                    <  1:130275/69‑1 (MQ=255)                                        cgcgATGGGGATTGGCCTGATTTTATTACTGGGATTCACCTGTTACGCCTCACGCGGGCTGTATTGGGc                    <  1:1873119/69‑1 (MQ=255)                                        cgcgATGGGGATTGGCCTGATTTTATTACTGGGATTCACCTGTTACGCCTCACGCGGGCTGTATTGGGc                    <  1:1184850/69‑1 (MQ=255)                                        cgcgATGGGGATTGGCCTGATTTTATTACTGGGATTCACCTGTTACGCCTCACGCGGGCTGTATTGGGc                    <  1:1111677/69‑1 (MQ=255)                                        cgcgATGGGGATTGGCCTGATTTTATTACTGGGATTCACCTGTTACGCCTCACGCGGGCTGTATTGGGc                    <  1:2695198/69‑1 (MQ=255)                                        cgcgATGGGGATTGGCCTGATTTTATTACTGGGATTCACCTGTTACGCCTCACGCGGGCTGTATTGGGc                    <  1:661013/69‑1 (MQ=255)                                        cgcgATGGGGATTGGCCTGATTTTATTACTGGGATTCACCTGTTACGCCTCACGCGGGCTGTATTGGGc                    <  1:2623742/69‑1 (MQ=255)                                        cgcgATGGGGATTGGCCTGATTTTATTACTGGGATTCACCTGTTACGCCTCACGCGGGCTGTATTGGGc                    <  1:1045437/69‑1 (MQ=255)                                        cgcgATGGGGATTGGCCTGATTTTATTACTGGGATTCACCTGTTACGCCTCACGCGGGCTGTATTGGGc                    <  1:2470630/69‑1 (MQ=255)                                        cgcgATGGGGATTGGCCTGATTTTATTACTGGGATTCACCTGTTACGCCTCACGCGGGCTGTATTCGGc                    <  1:1639031/69‑1 (MQ=255)                                        cgcgATGGGGATTGGCCTGATTTTATTACTGGGATTAACCTGTTACGCCTCACGCGGGCTGTATTGGGc                    <  1:3041189/69‑1 (MQ=255)                                           gATGGGGATTGGCCTGATTTTACTGCTGGGATTCACCTGTTACGCCTCACGCGGGCTGTACTGGGCCTg                 <  1:536516/69‑1 (MQ=255)                                           gATGGGGATTGGCCTGATTTTACTGCTGGGATTCACCTGTTACGCCTCACGCGGGCTGTACTGGGCCTg                 <  1:3127405/69‑1 (MQ=255)                                                       ccTGATTTTACTGCTGGGATTCACCTGTTACGCCTCACGCGGGCtt                            <  1:1028319/46‑2 (MQ=255)                                                       ccTGATTTTACTGCTGGGATTCACCTGTTACGCCTCACGCGGGCTg                            <  1:975465/46‑1 (MQ=255)                                                       ccTGATTTTACTGCTGGGATTCACCTGTTACGCCTCACGCGGGCTg                            <  1:127302/46‑1 (MQ=255)                                                       ccTGATTTTACTGCTGGGATTCACCTGTTACGCCTCACGCGGGCTg                            <  1:2543287/46‑1 (MQ=255)                                                       ccTGATTTTACTGCTGGGATTCACCTGTTACGCCTCACGCGGGCTg                            <  1:2889901/46‑1 (MQ=255)                                                        cTGATTTTACTGCTGGGATTCACCTGTTACGCCTCACGCGGGCTGTACTGGGCCTGCCCTGGCGaa       >  1:1713082/1‑66 (MQ=255)                                                        cTGATTTTACTGCTGGGATTCACCTGTTACGCCTCACGCGGGCTGTACTGGGCCTGCCCTGGCGaa       >  1:2454545/1‑66 (MQ=255)                                                        cTGATTTTACTGCTGGGATTCACCTGTTACGCCTCACGCGGGCTGTACTGGGCCTGCCCTGGCGaa       >  1:2111407/1‑66 (MQ=255)                                                        cTGATTTTACTGCTGGGATTCACCTGTTACGCCTCACGCGGGCTGTACTGGGCCTGCCCTGGCGaa       >  1:3286396/1‑66 (MQ=255)                                                          gATTTTACTGCTGGGATTCACCTGTTACGCCTCACGCGGGCTGTACTGGGCCTGCCCTGGCGAAGCgag  <  1:1640180/69‑1 (MQ=255)                                                          gATTTTACTGCTGGGATTCACCTGTTACGCCTCACGCGGGCTGTACTGGGCCTGCCCTGGCGAAGCgag  <  1:1434974/69‑1 (MQ=255)                                                          gATTTTACTGCTGGGATTCACCTGTTACGCCTCACGCGGGCTGTACTGGGCCTGCCCTGGCGAAGCgag  <  1:748917/69‑1 (MQ=255)                                                          gATTTTACTGCTGGGATTCACCTGTTACGCCTCACGCGGGCTGTACTGGGCCTGCCCTGGCGAAGCgag  <  1:2868854/69‑1 (MQ=255)                                                          gATTTTACTGCTGGGATTCACCTGTTACGCCTCACGCGGGCTGTACTGGGCCTGCCCTGGCGAAGCgag  <  1:1207802/69‑1 (MQ=255)                                                          gATTTTACTGCTGGGATTCACCTGTTACGCCTCACGCGGGCTGTACTGGGCCTGCCCTGGCGAAGCgag  <  1:3003301/69‑1 (MQ=255)                                                          gATTTTACTGCTGGGATTCACCTGTTACGCCTCACGCGGGCTGTACTGGGCCTGCCCTGGCGAAGCgag  <  1:2159689/69‑1 (MQ=255)                                                           aTTTTACTGCTGGGATTCACCTGTTACGCCTCACGCGGGCTGTACTgg                      >  1:1999683/1‑48 (MQ=255)                                                           aTTTTACTGCTGGGATTCACCTGTTACGCCTCACGCGGGCTGTACTgg                      >  1:855361/1‑48 (MQ=255)                                                           aTTTTACTGCTGGGATTCACCTGTTACGCCTCACGCGGGCTGTACTgg                      >  1:2959384/1‑48 (MQ=255)                                                           aTTTTACTGCTGGGATTCACCTGTTACGCCTCACGCGGGCTGTACTgg                      >  1:2824313/1‑48 (MQ=255)                                                           aTTTTACTGCTGGGATTCACCTGTTACGCCTCACGCGGGCTGTACTgg                      >  1:1955977/1‑48 (MQ=255)                                                           aTTTTACTGCTGGGATTCACCTGTTACGCCTCACGCGGGCTGTACTgg                      >  1:1714698/1‑48 (MQ=255)                                                           aTTTTACTGCTGGGATTCACCTGTTACGCCTCACGCGGGCTGTACTgg                      >  1:364222/1‑48 (MQ=255)                                                           aTTTTACTGCTGGGATTCACCTGTTACGCCTCACGCGGGCTGTACTgg                      >  1:1786131/1‑48 (MQ=255)                                                           aTTTTACTGCTGGGACTCACCTGTTACGCCTCACGCGGGCTGTACTGGGCCTGCCCTGGCGAAGCgag  <  1:2184205/68‑1 (MQ=255)                                                                   |                                                            AACTGCCCTGCTCGTCACGAACTCTAACCCGCAATCGGTCGCGATGGGGATTGGCCTGATTTTACTGCTGGGATTCACCTGTTACGCCTCACGCGGGCTGTACTGGGCCTGCCCTGGCGAAGCGAG  >  minE/1870985‑1871110

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