Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,813,656	G→A	38.4%	G507G (GGC→GGT)	alaS ←	alanyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,813,656	0	G	A	38.4%	17.5 / 22.5	39	G507G (GGC→GGT)	alaS	alanyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (11/13);  new base A (13/2);  total (24/15)
Fisher's exact test for biased strand distribution p-value = 1.73e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTATTTCTGCGTATCTTCCACCGCAAAGGAGAAGTTAGCGCCTTTCAGTTCGCCTTTATCGCCAACCTGACCGCCGGATTCCGCATAGAATGGCGTTTGATCCAGCACGACCACAGCTTCCTGGCC  >  minE/1813595‑1813721                                                              |                                                                  cGTATTTCTGCGTATCTTGCACCGCAAAGGAGAAGTTAGCTCGTTTCAGTTCGCCTTTATCGCCAACCt                                                            <  1:635252/69‑1 (MQ=255) cGTATTTCTGCGTATCTTCCACCGCAAAGGAGAAGTTAGCGCCTTTCAGTTCGCCTTTATCGCCAACCt                                                            <  1:2393649/69‑1 (MQ=255)  gTATTTCTGCGTATCTTCCACCGCGAAGGAGAAGTTAGCGCCTTTCAGTTCGCCTTTATCACCAACCTg                                                           <  1:3435713/69‑1 (MQ=255)  gTATTTCTGCGTATCTTCCACCGCGAAGGAGAAGTTAGCGCCTTTCAGTTCGCCTTTATCACCAACCTg                                                           <  1:2484706/69‑1 (MQ=255)  gTATTTCTGCGTATCTTCCACCGCAAAGGAGAAGTTAGCGCCTTTCAGTTCGCCTTTATCGCCAACCt                                                            <  1:1870591/68‑1 (MQ=255)     tttCTGCGTATCTTCCACCGCAAAGGAGAAGTTAGCGCCTTTCAGTTCGCCTTTATCGCCACCCTGAcc                                                        >  1:13416/1‑69 (MQ=255)     tttCTGCGTATCTTCCACCGCAAAGGAGAAGTTAGCGCCTTTCAGTTCGCCTTTATCGCCAACCTGAcc                                                        >  1:3107227/1‑69 (MQ=255)     tttCTGCGTATCTTCCACCGCAAAGGAGAAGTTAGCGCCTTTCAGTTCGCCTTTATCGCCAACCTGAcc                                                        >  1:3313758/1‑69 (MQ=255)          gCGTATCTTCCACCGCGAAGGAGAAGTTAGCGCCTTTCAGTTCGCCTTTATCACCAACCTGAccgccg                                                    >  1:3334548/1‑68 (MQ=255)          gCGTATCTTCCACCGCGAAGGAGAAGTTAGCGCCTTTCAGTTCGCCTTTATCACCAACCTGACCGCCgg                                                   >  1:3420148/1‑69 (MQ=255)          gCGTATCTTCCACCGCGAAGGAGAAGTTAGCGCCTTTCAGTTCGCCTTTATCACCAACCTGACCGCCgg                                                   >  1:2829757/1‑69 (MQ=255)          gCGTATCTTCCACCGCGAAGGAGAAGTTAGCGCCTTTCAGTTCGCCTTTATCACCAACCTGACCGCCgg                                                   >  1:3214685/1‑69 (MQ=255)          gCGTATCTTCCACCGCGAAGGAGAAGTTAGCGCCTTTCAGTTCGCCTTTATCACCAACCTGACCGCCgg                                                   >  1:271304/1‑69 (MQ=255)          gCGTATCTTCCACCGCGAAGGAGAAGTTAGCGCCTTTCAGTTCGCCTTTATCACCAACCTGACCGCCgg                                                   >  1:2487088/1‑69 (MQ=255)          gCGTATCTTCCACCGCGAAGGAGAAGTTAGCGCCTTTCAGTTCGCCTTTATCACCAACCTGACCGCCgg                                                   >  1:2465006/1‑69 (MQ=255)          gCGTATCTTCCACCGCGAAGGAGAAGTTAGCGCCTTTCAGTTCGCCTTTATCACCAACCTGACCGCCgg                                                   >  1:3434467/1‑69 (MQ=255)          gCGTATCTTCCACCGCGAAGGAGAAGTTAGCGCCTTTCAGTTCGCCTTTATCACCAACCTGACCGCCgg                                                   >  1:1895460/1‑69 (MQ=255)          gCGTATCTTCCACCGCGAAGGAGAAGTTAGCGCCTTTCAGTTCGCCTTTATCACCAACCTGACCGCCgg                                                   >  1:1630804/1‑69 (MQ=255)          gCGTATCTTCCACCGCGAAGGAGAAGTTAGCGCCTTTCAGTTCGCCTTTATCACCAACCTGACCGCCgg                                                   >  1:1317083/1‑69 (MQ=255)          gCGTATCTTCCACCGCGAAGGAGAAGTTAGCGCCTTTCAGTTCGCCTTTATCACCAACCTGACCGCCgg                                                   >  1:884503/1‑69 (MQ=255)          gCGTATCTTCCACCGCGAAGGAGAAG‑TAGCGCCTTTCAGTTCGCCTTTATCACCAACCTGACCGCCgg                                                   >  1:685938/1‑68 (MQ=255)                                      aGCGCCTTTCAGTTCGCCTTTATCGCCAACCTGACCGCCGGATTCCGCATAGAATGGCGTTTGATCCAg                       <  1:2356043/69‑1 (MQ=255)                                                    cGCCTTTATCGCCAACCTGACCGCCGGATTCCGCATAg                                        >  1:2061223/1‑38 (MQ=255)                                                    cGCCTTTATCGCCAACCTGACCGCCGGATTCCGCATAg                                        >  1:2639716/1‑38 (MQ=255)                                                    cGCCTTTATCGCCAACCTGACCGCCGGATTCCGCATAg                                        >  1:3138384/1‑38 (MQ=255)                                                    cGCCTTTATCGCCAACCTGACCGCCGGATTCCGCATAg                                        >  1:3083131/1‑38 (MQ=255)                                                    cGCCTTTATCGCCAACCTGACCGCCGGATTCCGCATAGAATGGCGTTTGATCCAGCACGACCACAGCtt         <  1:1211113/69‑1 (MQ=255)                                                    cGCCTTTATCGCCAACCTGACCGCCGGATTCCGCATAGAATGGCGTTTGATCCAGCACGACCACAGCtt         <  1:2593876/69‑1 (MQ=255)                                                    cGCCTTTATCGCCAACCTGACCGCCGGATTCCGCATAGAATGGCGTTTGATCCAGCACGACCACAGCtt         <  1:2384053/69‑1 (MQ=255)                                                      ccTTTATCGCCAACCTGACCGCCGGATTCCGCATAGAATGGCGTTTGATCCAGCACGACCACAGCTTc        <  1:3128714/68‑1 (MQ=255)                                                      ccTTTATCGCCAACCTGACCGCCGGATTCCGCATAGAATGGCGTTTGATCCAGCACGACCACAGCTTc        <  1:2605960/68‑1 (MQ=255)                                                      ccTTTATCGCCAACCTGACCGCCGGATTCCGCATAGAATGGCGTTTGATCCAGCACGACCACAGCTTc        <  1:2000968/68‑1 (MQ=255)                                                      ccTTTATCGCCAACCTGACCGCCGGATTCCGCATAGAATGGCGTTTGATCCAGCACGACCACAGCTTc        <  1:1794647/68‑1 (MQ=255)                                                      ccTTTATCGCCAACCTGACCGCCGGATTCCGCATAGAATGGCGTTTGATCCAGCACGACCACAGCTTc        <  1:1459442/68‑1 (MQ=255)                                                          tATCGCCAACCTGACCGCCGGATTCCGCATAGAATGGCGTTTGATCCAGCACGACCACAGCTTCCTGGc   >  1:1064356/1‑69 (MQ=255)                                                           aTCGCCAACCTGACCGCCGGATTCCGCATAGAATGGCGTTTGATCCAGCACGACCACAGCTTc        <  1:2722366/63‑1 (MQ=255)                                                           aTCGCCAACCTGACCGCCGGATTCCGCATAGAATGGCGTTTGATCCAGCACGACCACAGCTTCCTGGcc  >  1:2938095/1‑69 (MQ=255)                                                           aTCGCCAACCTGACCGCCGGATTCCGCATAGAATGGCGTTTGATCCAGCACGACCACAGCTTCCTGGcc  >  1:633121/1‑69 (MQ=255)                                                           aTCGCCAACCTGACCGCCGGATTCCGCATAGAATGGCGTTTGATCCAGCACGACCACAGCTTCCTGGc   >  1:1949445/1‑68 (MQ=255)                                                              |                                                                  CGTATTTCTGCGTATCTTCCACCGCAAAGGAGAAGTTAGCGCCTTTCAGTTCGCCTTTATCGCCAACCTGACCGCCGGATTCCGCATAGAATGGCGTTTGATCCAGCACGACCACAGCTTCCTGGCC  >  minE/1813595‑1813721

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