Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,753,942	C→T	55.0%	A20A (GCG→GCA)	mscL ←	mechanosensitive channel

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,753,942	0	C	T	55.0%	14.0 / 75.1	91	A20A (GCG→GCA)	mscL	mechanosensitive channel
Reads supporting (aligned to +/- strand):  ref base C (25/16);  new base T (26/24);  total (51/40)
Fisher's exact test for biased strand distribution p-value = 4.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

GATATCGGCAACCAGTGAAGAGACAATCTTCCCGAATGCCGCACCGATAATGACACCCACCGCCAAATCCACCACGTTCCCGCGCATCGCAAATTCGCGAAATTCTTTAATAATGCTCATGTTA  >  minE/2753882‑2754005                                                             |                                                                gatATCGGCAACCAGTGAAGAGACAATCTTTCCGAATGCCGCACCGATAATGACACCCACCGCCAAATc                                                         >  1:1067149/1‑69 (MQ=255) gatATCGGCAACCAGTGAAGAGACAATCTTCCCGAATGCCGCACCGATAATGACACCCACCGCCAAATc                                                         <  1:1094447/69‑1 (MQ=255)  atatCGGCAACCAGTGAAGAGACAATCTTCCCGAATGCCGCACCGATAATGACACCCACCGCCAAATcc                                                        <  1:204313/69‑1 (MQ=255)  atatCGGCAACCAGTGAAGAGACAATCTTCCCGAATGCCGCACCGATAATGACACCCACCGCCAAATcc                                                        <  1:586605/69‑1 (MQ=255)   tatCGGCAACCAGTGAAGAGACAATCTTCCCGAATGCCGCACCGATAATGACACCCACCGCCAATTcc                                                        <  1:1770280/68‑1 (MQ=255)       ggCAACCAGTGAAGATACAATCTTCCCGAATGCCGCTCCGATAATGACACGCACTGCCAAATCcaccac                                                   <  1:2011598/69‑1 (MQ=255)       ggCAACCAGTGAAGATACAATCTTCCCGAATGCCGCTCCGATAATGACACCCACTGCCAAATCcaccac                                                   <  1:2453868/69‑1 (MQ=255)       ggCAACCAGTGAAGATACAATCTTCCCGAATGCCGCTCCGATAATGACACCCACTGCCAAATCcaccac                                                   <  1:2008578/69‑1 (MQ=255)       ggCAACCAGTGAAGATACAATCTTCCCGAATGCCGCTCCGATAATGACACCCACTGCCAAATCcaccac                                                   <  1:2871292/69‑1 (MQ=255)       ggCAACCAGTGAAGATACAATCTTCCCGAATGCCGCTCCGATAATGACACCCACTGCCAAATCcaccac                                                   <  1:2874712/69‑1 (MQ=255)       ggCAACCAGTGAAGATACAATCTTCCCGAATGCCGCTCCGATAATGACACCCACTGCCAAATCcaccac                                                   <  1:3373518/69‑1 (MQ=255)       ggCAACCAGTGAAGATACAATCTTCCCGAATGCCGCTCCGATAATGACACCCACTGCCAAATCcaccac                                                   <  1:2930337/69‑1 (MQ=255)       ggCAACCAGTGAAGATACAATCTTCCCGAATGCCGCTCCGATAATGACACCCACTGCCAAATCcaccac                                                   <  1:2936950/69‑1 (MQ=255)       ggCAACCAGTGAAGATACAATCTTCCCGAATGCCGCTCCGATAATGACACCCACTGCCAAATCcaccac                                                   <  1:1724681/69‑1 (MQ=255)       ggCAACCAGTGAAGATACAATCTTCCCGAATGCCGCTCCGATAATGACACCCACTGCCAAATCcaccac                                                   <  1:266170/69‑1 (MQ=255)       ggCAACCAGTGAAGATACAATCTTCCCGAATGCCGCTCCGATAATGACACCCACTGCCAAATCcaccac                                                   <  1:661714/69‑1 (MQ=255)       ggCAACCAGTGAAGATACAATCTTCCCGAATGCCGCTCCGATAATGACACCCACTGCCAAATCcaccac                                                   <  1:1527115/69‑1 (MQ=255)       ggCAACCAGTGAAGATACAATCTTCCCGAATGCCGCTCCGATAATGACACCCACTGCCAAATCcaccac                                                   <  1:2624878/69‑1 (MQ=255)       ggCAACCAGTGAAGATACAATCTTCCCGAATGCCGCTCCGATAATGACACCCACTGCCAAATCcaccac                                                   <  1:2994302/69‑1 (MQ=255)       ggCAACCAGTGAAGATACAATCTTCCCGAATGCCGCTCCGATAATGACACCCACTGCCAAATCcaccac                                                   <  1:3346893/69‑1 (MQ=255)       ggCAACCAGTGAAGATACAATCTTCCCGAATGCCGCTCCGATAATGACACCCACTGCCAAATCcaccac                                                   <  1:1361559/69‑1 (MQ=255)       ggCAACCAGTGAAGATACAATCTTCCCGAATGCCGCTCCGATAATGACACCCACTGCCAAATCcaccac                                                   <  1:711651/69‑1 (MQ=255)       ggCAACCAGTGAAGATACAATCTTCCCGAATGCCGCTCCGATAATGACACCCACTGCCAAATCcaccac                                                   <  1:1282564/69‑1 (MQ=255)       ggCAACCAGTGAAGATACAATCTTCCCGAATGCCGCTCCGATAATGACACCCACTGCCAAATCcaacac                                                   <  1:1244958/69‑1 (MQ=255)       ggCAACCAGTGAAGATACAATCTTCACGAATGCCGCTCCGATAATGACACCCACTGCCAAATCcagcac                                                   <  1:1302854/69‑1 (MQ=255)       ggCAACCAGTGAAGAGACAATCTTCCCGAATGCCGCACCGATAATGACACCCACCGCCAAATCcaccac                                                   >  1:2831921/1‑69 (MQ=255)        gCAACCAGTGAAGATACAATCTTCCCGAATGCCGCTCCGATAATGACACCCACTGCCAAATCcaccac                                                   <  1:2678957/68‑1 (MQ=255)        gCAACCAGTGAAGATACAATCTTCCCGAATGCCGCTCCGATAATGACACCCACTGCCAAATCcaccac                                                   <  1:231010/68‑1 (MQ=255)        gCAACCAGTGAAGATACAATCTTCCCGAATGCCGCTCCGATAATGACACCCACTGCCAAATCcaccac                                                   <  1:2821390/68‑1 (MQ=255)        gCAACCAGTGAAGAGACAATCTTCCCGAATGCCGCACCGATAATGCCCCCCACCGCCAAATccacca                                                    >  1:1392513/1‑67 (MQ=255)        gCAACCAGTGAAGAGACAATCTTCCCGAATGCCGCACCGATAATGACACCCACCGCCAAATCCACCACg                                                  >  1:1668151/1‑69 (MQ=255)        gCAACCAGTGAAGAGACAATCTTCCCGAATGCCGCACCGATAATGACACCCACCGCCAAATCCACCACg                                                  >  1:3097453/1‑69 (MQ=255)          aaCCAGTGAAGATACAATCTTCCCGAATGCCGCTCCGATAATGACACCCACTGCCAAATCcaccac                                                   <  1:515834/66‑1 (MQ=255)                tGAAGAGACAATCTTCCCGAATGCCGCACCGATAATGACACCCACCGCCAAATCCACCACGTTCCcgcg                                          >  1:2430123/1‑69 (MQ=255)                tGAAGAGACAATCTTCCCGAATGCCGCACCGATAATGACACCCACCGCCAAATCCACCACGTTCCcgcg                                          <  1:1524001/69‑1 (MQ=255)                 gAAGAGACAATCTTCCCGAATGCCGCACCGATAATGACACCCACCGCCAAATCCACCACGTTCCcgcg                                          >  1:2438050/1‑68 (MQ=255)                 gAAGAGACAATCTTCCCGAATGCCGCACCGATAATGACACCCACCGCCAAATCCACCACGTTCCCgcgc                                         >  1:1846783/1‑69 (MQ=255)                  aaGAGACAATCTTCCCGAATGCCGCACCGATAATGACACCCACCGCCAAATCCACCACGTTCCCGCGCa                                        >  1:960452/1‑69 (MQ=255)                  aaGAGACAATCTTCCCGAATGCCGCACCGATAATGACACCCACCGCCAAATCCACCACGTTCCCGCGCa                                        <  1:1155510/69‑1 (MQ=255)                       aCAATCTTCCCGAATGCCTCACCGATAATGACACCCACCGCCAAATCCACCACGTTCCCGCGCATCGCa                                   <  1:2373270/69‑1 (MQ=255)                       aCAATCTTCCCGAATGCCGCACCGATAATGACACCCACCGCCAAATCCACCACGTTCCCGCGCATCGCa                                   <  1:2603897/69‑1 (MQ=255)                        cAATCTTCCCGAATGCCGCACCGATAATGACACCCACCGCCAAATCCACCACGTTCCCGCGCATCGCa                                   <  1:295662/68‑1 (MQ=255)                            cTTCCCGAATGCCGCACCGTTAATGACACCCACCGCCaaa                                                           >  1:3394850/1‑40 (MQ=255)                            cTTCCCGAATGCCGCACCGATAATGACACCCACCGCCAAATCCACCACGTTCCCGCGCATCGCAAAtt                               >  1:529175/1‑68 (MQ=255)                            cTTCCCGAATGCCGCACCGATAATGACACCCACCGCCAAATCCACCACGTTCCCGCGCATCGCAAATTc                              >  1:2928335/1‑69 (MQ=255)                            cTTCCCGAATGCCGCACCGATAATGACACCCACCGCCAAATCCACCACGTTCCCGCGCATCGCAAATTc                              >  1:3315425/1‑69 (MQ=255)                            cTTCCCGAATGCCGCACCGATAATGACACCCACCGCCAAATCCACCACGTTCCCGCGCATCGCAAATTc                              >  1:3173942/1‑69 (MQ=255)                            cTTCCCGAATGCCGCACCGATAATGACACCCACCGCCAAATCCACCACGTTCCCGCGCATCGCAAATTc                              >  1:2091822/1‑69 (MQ=255)                            cTTCCCGAATGCCGCACCGATAATGACACCCACCGCCAAATCCACCACGTTCCCGCGCATCGCAAATTc                              >  1:1749425/1‑69 (MQ=255)                            cTTCCCGAATGCCGCACCGATAATGACACCCACCGCCAAATCCACCACGTTCCCGCGCATCGCAAATTc                              >  1:624770/1‑69 (MQ=255)                            cTTCCCGAATGCCGCACCGATAATGACACCCACCGCCAAATCCACCACGTTCCCGCGCATCGCAAATTc                              >  1:78572/1‑69 (MQ=255)                            cTTCCCGAATGCCGCACCGATAATGACACCCACCGCCAAATCCACCACGTTCCCGCGCATCGCAAATTc                              >  1:733248/1‑69 (MQ=255)                            cTTCCCGAATGCCGCACCGATAATGACACCCACCGCCAAATCCACCACGTTCCCGCGCATCGCAAATTc                              >  1:606431/1‑69 (MQ=255)                                     tGCCGCACCGATAATGACTCCCACCGCCAAATCCACCACGTTCCCGCGCATCGCAAATTCGCGAAATTc                     >  1:2385604/1‑69 (MQ=255)                                     tGCCGCACCGATAATGACACCCACCGCCAAATCCACCACGTTCCCGCGCATCGCAAATTCGCGAAATTc                     >  1:158524/1‑69 (MQ=255)                                     tGCCGCACCGATAATGACACCCACCGCCAAATCCACCACGTTCCCGCGCATCGCAAATTCGCGAAATTc                     >  1:1115092/1‑69 (MQ=255)                                     tGCCGCACCGATAATGACACCCACCGCCAAATCCACCACGTTCCCGCGCATCGCAAATTCGCGAAATTc                     >  1:2115305/1‑69 (MQ=255)                                           aCCGATAATGACACCCACCGCCAAATCCACCACGTTCCCGCGCATCGCAAATTCGCGAAATTCTTtaat               <  1:2653116/69‑1 (MQ=255)                                           aCCGATAATGACACCCACCGCCAAATCCACCACGTTCCCGCGCATCGCAAATTCGCGAAATTCTTtaat               <  1:1176747/69‑1 (MQ=255)                                           aCCGATAATGACACCCACCGCCAAATCCACCACGTTCCCGCGCATCGCAAATTCGCGAAATTCTTtaat               <  1:2690023/69‑1 (MQ=255)                                           aCCGATAATGACACCCACCGCCAAATCCACCACGTTCCCGCGCATCGCAAATTCGCGAAATTCTTtaat               <  1:2701228/69‑1 (MQ=255)                                           aCCGATAATGACACCCACCGCCAAATCCACCACGTTCCCGCGCATCGCAAATTCGCGAAATTCTTtaat               <  1:2424545/69‑1 (MQ=255)                                           aCCGATAATGACACCCACCGCCAAATCCACCACGTTCCCGCGCATCGCAAATTCGCGAAATTCTTtaat               <  1:2338977/69‑1 (MQ=255)                                            ccGATAATGACACCCACCGCCAAATCCACCACGTTCCCGCGCATCGCAAATTCGCGAAATTCTTtaat               <  1:1768934/68‑1 (MQ=255)                                                tAATGACACCTACTGCCAAATCCACCACGTTCCCGCGCATCGCAAATTCGCGAAATTCTTTAATAATg           >  1:1734239/1‑68 (MQ=255)                                                tAATGACACCCACTGCCAAATCCACCACGTTCCCGCGCATCGCAAATTCGCGAAATTCTTtaataa             >  1:1891060/1‑66 (MQ=255)                                                tAATGACACCCACTGCCAAATCCACCACGTTCCCGCGCATCGCAAATTCGCGAAATTCTTTaataat            >  1:141810/1‑67 (MQ=255)                                                tAATGACACCCACTGCCAAATCCACCACGTTCCCGCGCATCGCAAATTCGCGAAATTCTTTAATAATg           >  1:642546/1‑68 (MQ=255)                                                tAATGACACCCACTGCCAAATCCACCACGTTCCCGCGCATCGCAAATTCGCGAAATTCTTTAATAATg           >  1:1077781/1‑68 (MQ=255)                                                tAATGACACCCACTGCCAAATCCACCACGTTCCCGCGCATCGCAAATTCGCGAAATTCTTTAATAATg           >  1:480220/1‑68 (MQ=255)                                                tAATGACACCCACTGCCAAATCCACCACGTTCCCGCGCATCGCAAATTCGCGAAATTCTTTAATAATg           >  1:673722/1‑68 (MQ=255)                                                tAATGACACCCACTGCCAAATCCACCACGTTCCCGCGCATCGCAAATTCGCGAAATTCTTTAATAATg           >  1:400633/1‑68 (MQ=255)                                                tAATGACACCCACTGCCAAATCCACCACGTTCCCGCGCATCGCAAATTCGCGAAATTCTTTAATAATg           >  1:341233/1‑68 (MQ=255)                                                tAATGACACCCACTGCCAAATCCACCACGTTCCCGCGCATCGCAAATTCGCGAAATTCTTTAATAATg           >  1:801223/1‑68 (MQ=255)                                                tAATGACACCCACTGCCAAATCCACCACGTTCCCGCGCATCGCAAATTCGCGAAATTCTTTAATAATg           >  1:1962737/1‑68 (MQ=255)                                                tAATGACACCCACTGCCAAATCCACCACGTTCCCGCGCATCGCAAATTCGCGAAATTCTTTAATAATg           >  1:1467907/1‑68 (MQ=255)                                                tAATGACACCCACTGCCAAATCCACCACGTTCCCGCGCATCGCAAATTCGCGAAATTCTTTAATAATg           >  1:1578317/1‑68 (MQ=255)                                                tAATGACACCCACTGCCAAATCCACCACGTTCCCGCGCATCGCAAATTCGCGAAATTCTTTAATAATg           >  1:1746211/1‑68 (MQ=255)                                                tAATGACACCCACTGCCAAATCCACCACGTTCCCGCGCATCGCAAATTCGCGAAATTCTTTAATAATg           >  1:135792/1‑68 (MQ=255)                                                tAATGACACCCACTGCCAAATCCACCACGTTCCCGCGCATCGCAAATTCGCGAAATTCTTTAATAATg           >  1:1767212/1‑68 (MQ=255)                                                tAATGACACCCACTGCCAAATCCACCACGTTCCCGCGCATCGCAAATTCGCGAAATTCTTTAATAATg           >  1:1785190/1‑68 (MQ=255)                                                tAATGACACCCACTGCCAAATCCACCACGTTCCCGCGCATCGCAAATTCGCGAAATTCTTTAATAATg           >  1:178873/1‑68 (MQ=255)                                                tAATGACACCCACTGCCAAATCCACCACGTTCCCGCGCATCGCAAATTCGCGAAATTCTTTAATAATg           >  1:1830883/1‑68 (MQ=255)                                                tAATGACACCCACTGCCAAATCCACCACGTTCCCGCGCATCGCAAATTCGCGAAATTCTTTAATAATg           >  1:2057735/1‑68 (MQ=255)                                                tAATGACACCCACTGCCAAATCCACCACGTTCCCGCGCATCGCAAATTCGCGAAATTCTTTAATAATg           >  1:2331686/1‑68 (MQ=255)                                                tAATGACACCCACTGCCAAATCCACCACGTTCCCGCGCATCGCAAATTCGCGAAATTCTTTAATAATg           >  1:2535616/1‑68 (MQ=255)                                                tAATGACACCCACTGCCAAATCCACCACGTTCCCGCGCATCGCAAATTCGCGAAATTCTTTAATAATg           >  1:1338675/1‑68 (MQ=255)                                                tAATGACACCCACTGCCAAATCCACCACGTTCCCGCGCATCGCAAATTCGCGAAATTCTTTAATAATg           >  1:2883291/1‑68 (MQ=255)                                                tAATGACACCCACTGCCAAATCCACCACGTTCCCGCGCATCGCAAATTCGCGAAATTCTTTAATAATg           >  1:1249619/1‑68 (MQ=255)                                                tAATGACACCCACTGCCAAATCCACCACGTTCCCGCGCATCGCAAATTCGCGAAATTCTTTAATAATg           >  1:3132710/1‑68 (MQ=255)                                                  aTGACACCCACCGCCAAATCCACCACGTTCCCGCGCATCGCAAATTCGCGAAATTCTTTAATAATGCTc        >  1:3034002/1‑69 (MQ=255)                                                       acccaccGCCAAAT‑CACCACGTTCCCGCGCATCGCAAATTCGCGAAATTCTTTAATAATGCTCATGTTa  <  1:2871037/69‑1 (MQ=255)                                                        cccaccGCCAAATCCACCACGTTCCCGCGCATCGCAAATTCGCGAAATTCTTTAATAATGCTCATGTTa  <  1:631022/69‑1 (MQ=255)                                                        cccaccGCCAAATCCACCACGTTCCCGCGCATCGCAAATTCGCGAAATTCTTTAATAATGCTCATGTTa  <  1:3352769/69‑1 (MQ=255)                                                        cccaccGCCAAATCCACCACGTTCCCGCGCATCGCAAATTCGCGAAATTCTTTAATAATGCTCATGTTa  <  1:1178452/69‑1 (MQ=255)                                                        cccaccGCCAAATCCACCACGTTCCCGCGCATCGCAAATTCGCGAAATTCTTTAATAATGCTCATGTTa  <  1:1126771/69‑1 (MQ=255)                                                          caccGCCAAATCCACCACGTTCCCGCGCATCGCAAATTCGCGAAATTCTTTAATAATGCTCATGTTa  <  1:964630/67‑1 (MQ=255)                                                             |                                                                GATATCGGCAACCAGTGAAGAGACAATCTTCCCGAATGCCGCACCGATAATGACACCCACCGCCAAATCCACCACGTTCCCGCGCATCGCAAATTCGCGAAATTCTTTAATAATGCTCATGTTA  >  minE/2753882‑2754005

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