Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,253,551	A→G	70.7%	A296A (GCT→GCC)	aspS ←	aspartyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,253,551	0	A	G	70.7%	49.8 / 19.5	51	A296A (GCT→GCC)	aspS	aspartyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (10/5);  new base G (19/17);  total (29/22)
Fisher's exact test for biased strand distribution p-value = 5.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

TACAGCAAACTCAACAGATTTCAGCAGATCAGCAACGTCAGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCGTTCTGCTT  >  minE/1253521‑1253618                               |                                                                    tACAGCAAACTCAACAGATTTCAGCAGATCAGCAACGTCAGTCAGTTCCATCGGGTTACGCAGATCCgg                               <  1:163092/69‑1 (MQ=255) tACAGCAAACTCAACAGATTTCAGCAGATCAGCAACGTCAGTCAGTTCCATCGGGTTACGCAGATCCgg                               <  1:1676405/69‑1 (MQ=255)  aCAGCAAACTCAACAGATTTCAGCAGATCAGCAACGTCAGTCAGTTCCATCGGGTTACGCAGATCCgg                               >  1:973646/1‑68 (MQ=255)  aCAGCAAACTCAACAGATTTCAGCAGATCAG‑AACGTCAGTCAGTTCCATCGGGTTACGCAGATCCGGtt                             >  1:1007453/1‑69 (MQ=255)        aaCTCAACAGATTTCAGCAGATCAGCAACGTCAGTCAGTTcc                                                   >  1:1380568/1‑42 (MQ=255)        aaCTCAACAGATTTCAGCAGATCAGCAACGTCAGTCAGTTCCATCGGGTTACGCAGATCCGGTTTAta                         >  1:1422909/1‑67 (MQ=255)        aaCTCAACAGATTTCAGCAGATCAGCAACGTCAGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATc                         >  1:945861/1‑68 (MQ=255)        aaCTCAACAGATTTCAGCAGATCAGCAACGTCAGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATc                         >  1:418632/1‑68 (MQ=255)         aCTCAACAGATTTCAGCAGATCAGCAACGTCAGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAg                       <  1:792072/69‑1 (MQ=255)         aCTCAACAGATTTCAGCAGATCAGCAACGTCAGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAg                       <  1:726609/69‑1 (MQ=255)         aCTCAACAGATTTCAGCAGATCAGCAACGTCAGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAg                       <  1:1917024/69‑1 (MQ=255)                       agcagATCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg          <  1:2316198/68‑1 (MQ=255)                       agcagATCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg          <  1:2375964/68‑1 (MQ=255)                       agcagATCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg          <  1:2828635/68‑1 (MQ=255)                       agcagATCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg          <  1:2954045/68‑1 (MQ=255)                       agcagATCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg          <  1:297186/68‑1 (MQ=255)                       agcagATCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg          <  1:3182686/68‑1 (MQ=255)                       agcagATCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg          <  1:3414039/68‑1 (MQ=255)                       agcagATCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg          <  1:460090/68‑1 (MQ=255)                       agcagATCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg          <  1:1575565/68‑1 (MQ=255)                       agcagATCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg          <  1:1025237/68‑1 (MQ=255)                       agcagATCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg          <  1:1176772/68‑1 (MQ=255)                       agcagATCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg          <  1:1851042/68‑1 (MQ=255)                        gcagATCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg          <  1:1957261/67‑1 (MQ=255)                        gcagATCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg          <  1:788737/67‑1 (MQ=255)                        gcagATCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg          <  1:3311873/67‑1 (MQ=255)                        gcagATCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAAACATAACGGCg          <  1:692211/67‑1 (MQ=255)                        gcagATCGGCAACGTCGGTCAGTTCCATCAGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg          <  1:243877/67‑1 (MQ=255)                            aTCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGc           >  1:1662891/1‑62 (MQ=255)                            aTCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg          >  1:2101798/1‑63 (MQ=255)                            aTCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg          >  1:2232308/1‑63 (MQ=255)                            aTCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg          >  1:1081628/1‑63 (MQ=255)                            aTCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg          >  1:1165223/1‑63 (MQ=255)                            aTCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg          >  1:2172140/1‑63 (MQ=255)                            aTCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg          >  1:625734/1‑63 (MQ=255)                            aTCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg          >  1:571514/1‑63 (MQ=255)                            aTCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg          >  1:129976/1‑63 (MQ=255)                            aTCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg          >  1:1366377/1‑63 (MQ=255)                            aTCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg          >  1:1461619/1‑63 (MQ=255)                            aTCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg          >  1:3396170/1‑63 (MQ=255)                            aTCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg          >  1:1510916/1‑63 (MQ=255)                            aTCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg          >  1:3144496/1‑63 (MQ=255)                            aTCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg          >  1:1698047/1‑63 (MQ=255)                            aTCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg          >  1:2792467/1‑63 (MQ=255)                            aTCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg          >  1:25941/1‑63 (MQ=255)                            aTCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg          >  1:256278/1‑63 (MQ=255)                            aTCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg          >  1:2525754/1‑63 (MQ=255)                              cAGCAACGTCAGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCGTTCTGCtt  >  1:2003440/1‑69 (MQ=255)                              cAGCAACGTCAGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCGTTCTGCtt  >  1:2648270/1‑69 (MQ=255)                              cAGCAACGTCAGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCGTTCTGCtt  >  1:2432146/1‑69 (MQ=255)                              cAGCAACGTCAGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCGTTCTGCtt  >  1:2176278/1‑69 (MQ=255)                               |                                                                    TACAGCAAACTCAACAGATTTCAGCAGATCAGCAACGTCAGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCGTTCTGCTT  >  minE/1253521‑1253618

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