Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,590,559	T→G	33.7%	Y54D (TAT→GAT)	yfeK →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,590,559	0	T	G	33.7%	63.6 / 61.5	80	Y54D (TAT→GAT)	yfeK	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (24/29);  new base G (20/7);  total (44/36)
Fisher's exact test for biased strand distribution p-value = 1.81e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

GCTGGAGGGATTAGCACAGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATACCCGTAACCGCATTGACACTGCCG  >  minE/1590495‑1590622                                                                 |                                                                gCTGGAGGGATTAGCACAGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGa                                                             >  1:2724214/1‑69 (MQ=255)  cTGGAGGGATTAGCACAGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGaa                                                            >  1:1343214/1‑69 (MQ=255)          aTTAGCACAGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGttt                                                     >  1:2587559/1‑68 (MQ=255)          aTTAGCACAGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGttt                                                     >  1:1709642/1‑68 (MQ=255)            tAGCACAGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTtctc                                                  >  1:2234349/1‑69 (MQ=255)            tAGCACAGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTtctc                                                  >  1:1543823/1‑69 (MQ=255)             aGCACAGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCa                                                 >  1:275256/1‑69 (MQ=255)             aGCACAGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCa                                                 >  1:1574601/1‑69 (MQ=255)               cacaGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCGATGAAGCGGTTTCTCAt                                                <  1:3336813/68‑1 (MQ=255)               cacaGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCGATGAAGCGGTTTCTCAt                                                <  1:3044786/68‑1 (MQ=255)               cacaGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCGATGAAGCGGTTTCTCAt                                                <  1:640374/68‑1 (MQ=255)               cacaGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCGATGAAGCGGTTTCTCAt                                                <  1:2115598/68‑1 (MQ=255)               cacaGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCGATGAAGCGGTTTCTCAt                                                <  1:1911939/68‑1 (MQ=255)               cacaGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCGATGAAGCGGTTTCTCAt                                                <  1:1754735/68‑1 (MQ=255)               cacaGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCGATGAAGCGGTTTCTCAt                                                <  1:1607572/68‑1 (MQ=255)                      aaaGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTc                                         >  1:347913/1‑68 (MQ=255)                        aGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGa                                      >  1:1145860/1‑69 (MQ=255)                        aGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGa                                      >  1:981686/1‑69 (MQ=255)                          gATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAg                                    <  1:1720497/69‑1 (MQ=255)                          gATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAg                                    <  1:663169/69‑1 (MQ=255)                          gATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAg                                    <  1:878267/69‑1 (MQ=255)                          gATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAg                                    <  1:3375914/69‑1 (MQ=255)                          gATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAg                                    <  1:3023882/69‑1 (MQ=255)                          gATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCACATCTGCGTCTGAAg                                    <  1:2160082/69‑1 (MQ=255)                            tttGATATTTGTGCGCAACGGTGATGAACATACCTGCGATGAAGCg                                                         >  1:2890352/1‑46 (MQ=255)                            tttGATATTTGTGCGCAACGGTGATGAACATACCTGCGATGAAGCg                                                         >  1:3350057/1‑46 (MQ=255)                            tttGATATTTGTGCGCAACGGTGATGAACATACCTGCGATGAAGCg                                                         >  1:3250564/1‑46 (MQ=255)                            tttGATATTTGTGCGCAACGGTGATGAACATACCTGCGATGAAGCg                                                         >  1:1419399/1‑46 (MQ=255)                            tttGATATTTGTGCGCAACGGTGATGAACATACCTGCGATGAAGCg                                                         >  1:3029496/1‑46 (MQ=255)                            tttGATATTTGTGCGCAACGGTGATGAACATACCTGCGATGAAGCg                                                         >  1:2630974/1‑46 (MQ=255)                            tttGATATTTGTGCGCAACGGTGATGAACATACCTGCGATGAAGCg                                                         >  1:346168/1‑46 (MQ=255)                            tttGATATTTGTGCGCAACGGTGATGAACATACCTGCGATGAAGCg                                                         >  1:1252895/1‑46 (MQ=255)                            tttGATATTTGTGCGCAACGGTGATGAACATACCTGCGATGAAGCg                                                         >  1:3310198/1‑46 (MQ=255)                            tttGATATTTGTGCGCAACGGTGATGAACATACCTGCGATGAAGCg                                                         >  1:2532009/1‑46 (MQ=255)                            tttGATATTTGTGCGCAACGGTGATGAACATACCTGCGATGAAGCg                                                         >  1:449158/1‑46 (MQ=255)                            tttGATATTTGTGCGCAACGGTGATGAACATACCTGCGATGAAGCg                                                         >  1:1221394/1‑46 (MQ=255)                            tttGATATTTGTGCGCAACGGTGATGAACATACCTGCGATGAAGCg                                                         >  1:700469/1‑46 (MQ=255)                            tttGATATTTGTGCGCAACGGTGATGAACATACCTGCGATGAAGCg                                                         >  1:2255165/1‑46 (MQ=255)                            tttGATATTTGTGCGCAACGGTGATGAACATACCTGCGATGAAGCg                                                         >  1:1457826/1‑46 (MQ=255)                            tttGATATTTGTGCGCAACGGTGATGAACATACCTGCGATGAAGCg                                                         >  1:2179396/1‑46 (MQ=255)                            tttGATATTTGTGCGCAACGGTGATGAACATACCTGCGATGAAGCg                                                         >  1:811436/1‑46 (MQ=255)                            tttGATATTTGTGCGCAACGGTGATGAACATACCTGCGATGAAGCg                                                         >  1:148087/1‑46 (MQ=255)                            tttGATATTTGTGCGCAACGGTGATGAACATACCTGCGATGAAGCg                                                         >  1:1604563/1‑46 (MQ=255)                            tttGATATTTGTGCGCAACGGTGATGAACATACCTGCGATGAAGCg                                                         >  1:889922/1‑46 (MQ=255)                             ttGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTc                                 >  1:2215818/1‑69 (MQ=255)                             ttGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTc                                 >  1:2886652/1‑69 (MQ=255)                             ttGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTc                                 >  1:2099351/1‑69 (MQ=255)                             ttGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTc                                 >  1:2927424/1‑69 (MQ=255)                             ttGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTATCATCTGCGTCTGAAGCTc                                 >  1:2416755/1‑69 (MQ=255)                                 tatTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGc                              <  1:3425952/68‑1 (MQ=255)                                 tatTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGc                              <  1:3044512/68‑1 (MQ=255)                                  atTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCaa                            <  1:2564225/69‑1 (MQ=255)                                  atTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCaa                            <  1:3312990/69‑1 (MQ=255)                                  atTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCaa                            <  1:3352823/69‑1 (MQ=255)                                     tgtgCGCAACGGTGATGATCATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCATTa                          <  1:1439923/68‑1 (MQ=255)                                     tgtgCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATa                          <  1:1318081/68‑1 (MQ=255)                                     tgtgCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATa                          <  1:741657/68‑1 (MQ=255)                                        gcgcAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATa                          <  1:472146/65‑1 (MQ=255)                                                  gatgaACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGc                              <  1:1472888/51‑1 (MQ=255)                                                     gatcATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATACCCGTAACCGCATTGa          <  1:234342/65‑1 (MQ=255)                                                     gaACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCTGCAATACCCGTAACCGCATTGa          <  1:2632032/68‑1 (MQ=255)                                                     gaACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATACCCGTAACCGCATTGa          <  1:1345415/68‑1 (MQ=255)                                                     gaACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATACCCGTAACCGCATTGa          <  1:2300943/68‑1 (MQ=255)                                                     gaACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATACCCGTAACCGCATTGa          <  1:2698971/68‑1 (MQ=255)                                                     gaACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATACCCGTAACCGCATTGa          <  1:1735154/68‑1 (MQ=255)                                                     gaACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATACCCGTAACCGCATTGa          <  1:3048851/68‑1 (MQ=255)                                                      aaCATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGc                                   >  1:2276962/1‑42 (MQ=255)                                                      aaCATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGc                              >  1:3102237/1‑47 (MQ=255)                                                      aaCATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGc                              >  1:68517/1‑47 (MQ=255)                                                      aaCATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGc                              >  1:504374/1‑47 (MQ=255)                                                      aaCATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGc                              >  1:3399783/1‑47 (MQ=255)                                                      aaCATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATACCCGTAACCGCATTGac         <  1:606272/68‑1 (MQ=255)                                                      aaCATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATACCCGTAACCGCATTGac         <  1:2471705/68‑1 (MQ=255)                                                      aaCATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATACCCGTAACCGCATTGac         <  1:3264892/68‑1 (MQ=255)                                                      aaCATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATACCCGTAACCGCATTGac         <  1:1723193/68‑1 (MQ=255)                                                       aCATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATACCCGTAACCGCATTGac         <  1:1310991/67‑1 (MQ=255)                                                           aCCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATACCCGTAAcc                 >  1:1611981/1‑55 (MQ=255)                                                           aCCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATACCCGTAACCGCATTGACACTg     >  1:1095025/1‑67 (MQ=255)                                                            ccTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATAc                         <  1:3183946/46‑1 (MQ=255)                                                            ccTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATACCCGTAACCGCATTGACACTGCCg  >  1:3183641/1‑69 (MQ=255)                                                                 |                                                                GCTGGAGGGATTAGCACAGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATACCCGTAACCGCATTGACACTGCCG  >  minE/1590495‑1590622

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑