Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,253,716	C→T	34.1%	A241A (GCG→GCA)	aspS ←	aspartyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,253,716	0	C	T	34.1%	49.2 / 38.4	44	A241A (GCG→GCA)	aspS	aspartyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base C (9/20);  new base T (8/7);  total (17/27)
Fisher's exact test for biased strand distribution p-value = 1.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCACACCCTTCACTTCCAGCCACAGATGACGCACCAGCGCTTCCATCACTTCACGCACTTGCGGCGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTAAATTCAGGCTGACGGTCAGCACGCAGGTCT  >  minE/1253651‑1253782                                                                  |                                                                   tCCACACCCTTCACTTCCAGCCACAGATGACGCACCAGCGCTTCCATCACTTCACGCACTTGCGGCGCGGt                                                               >  1:1033670/1‑71 (MQ=255)  ccACACCCTTCACTTCCAGCCACAGATGACGCACCAGCGCTTCCATCACTTCACGCACTTgcggtgcggt                                                               >  1:957286/1‑70 (MQ=255)  ccACACCCTTCACTTCCAGCCACAGATGACGCACCAGCGCTTCCATCACTTCACGCACTTgcggtgcggt                                                               >  1:1777133/1‑70 (MQ=255)      acCCTTCACTTCCAGCCACAGATGACGCACCAGCGCTTCCATCACTTCACGCACTTGCGGTGCGGTCATGa                                                          <  1:1690905/71‑1 (MQ=255)      acCCTTCACTTCCAGCCACAGATGACGCACCAGCGCTTCCATCACTTCACGCACTTGCGGCGCGGTCATGa                                                          <  1:1769098/71‑1 (MQ=255)      acCCTTCACTTCCAGCCACAGATGACGCACCAGCGCTTCCATCACTTCACGCACTTGCGGCGCGGTCATGa                                                          <  1:1223927/71‑1 (MQ=255)         cTTCACTTCCAGCCACAGATGACGCACCAGCGCTTCCATCACTTCACGCACTTGCGGTGCGGTCATGAaag                                                       <  1:1435026/71‑1 (MQ=255)           tCACTTCCAGCCACAGATGACGCACCAGCGCTTCCATCACTTCACGCACTTGCGGTGCGGTCATGaa                                                         >  1:629705/1‑67 (MQ=255)           tCACTTCCAGCCACAGATGACGCACCAGCGCTTCCATCACTTCACGCACTTGCGGCGCGGTCATGAaag                                                       >  1:1091582/1‑69 (MQ=255)              cTTCCAGCCACAGATGACGCACCAGCGCTTCCATCACTTCACGCACTTGCGGTGCGGTCATGAAAGAAGtt                                                  <  1:1985570/71‑1 (MQ=255)               ttCCAGCCACAGATGACGCACCAGCGCTTCCATCACTTCACGCACTTGCGGCGCGGTCATGAaagaa                                                     >  1:1566383/1‑67 (MQ=255)                       acaGATGACGCACCAGCGCTTCCATCACTTCACGCACTTGCGGCGCGGTCATGAAAGAAGTTTCCACATCg                                         <  1:1141743/71‑1 (MQ=255)                        caGATGACGCACCAGCGCTTCCATCACTTCACGCACTTGCGGCGCGGTCATGAAAGAAGTTTCCACATCGa                                        <  1:2095687/71‑1 (MQ=255)                         aGATGACGCACCAGCGCTTCCATCACTTCACGCACTTGCGGTGCGGTCATGAAAGAAGTTTCCACATCGAt                                       >  1:2147974/1‑71 (MQ=255)                           aTGACGCACCAGCGCTTCCATCACTTCACGCACTTGCGGCGCGGTCATGAaagaag                                                    <  1:1319850/56‑1 (MQ=255)                                gCACCAGCGCTTCCATCACTTCACGCACTTGCGGCGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTa                                <  1:1089275/71‑1 (MQ=255)                                gCACCAGCGCTTCCATCACTTCACGCACTTGCGGCGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTa                                <  1:1826546/71‑1 (MQ=255)                                  aCCAGCGCTTCCATCACTTCACGCACTTGCGGCGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTaaa                              <  1:1116951/71‑1 (MQ=255)                                  aCCAGCGCTTCCATCACTTCACGCACTTGCGGCGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTaaa                              >  1:2084883/1‑71 (MQ=255)                                      gcgcTTCCATCACTTCACGCACTTGCGGTGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTAAATTCa                          <  1:1200770/71‑1 (MQ=255)                                      gcgcTTCCATCACTTCACGCACTTGCGGTGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTAAATTCa                          <  1:1460070/71‑1 (MQ=255)                                      gcgcTTCCATCACTTCACGCACTTGCGGTGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTAAATTCa                          <  1:937098/71‑1 (MQ=255)                                      gcgcTTCCATCACTTCACGAACTTGCGGTGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTAAAttta                          <  1:2094978/71‑3 (MQ=255)                                        gcTTCCATCACTTCACGCACTTGCGGCGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTAAATTCAgg                        >  1:188148/1‑71 (MQ=255)                                            ccATCACTTCACGCACTTGCGGCGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTAAATTCAGGCTGa                    >  1:663105/1‑71 (MQ=255)                                            ccATCACTTCACGCACTTGCGGCGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTAAATTCAGGCTGa                    >  1:350532/1‑71 (MQ=255)                                            ccATCACTTCACGCACTTGCGGCGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTAAATTCAGGCTGa                    >  1:1195897/1‑71 (MQ=255)                                                 aCTTCACGCACTTGCGGCGCGGTCATGAAAGAAGTTTCCACATCg                                         <  1:1808050/45‑1 (MQ=255)                                                 aCTTCACGCACTTGCGGCGCGGTCATGAAAGAAGTTTCCACATCg                                         <  1:1139482/45‑1 (MQ=255)                                                 aCTTCACGCACTTGCGGCGCGGTCATGAAAGAAGTTTCCACATCg                                         <  1:1986370/45‑1 (MQ=255)                                                 aCTTCACGCACTTGCGGCGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTAAATTCAGGCTGACGGt                >  1:2050477/1‑70 (MQ=255)                                                 aCTTCACGCACTTGCGGCGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTAAATTCAGGCTGACGGt                <  1:505805/70‑1 (MQ=255)                                                 aCTTCACGCACTTGCGGCGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTAAATTCAGGCTGACGGt                <  1:715687/70‑1 (MQ=255)                                                   ttCACGCACTTGCGGCGCGGTCATGAAAGAAGTTTCcac                                             <  1:1160249/39‑1 (MQ=255)                                                   ttCACGCACTTGCGGCGCGGTCATGAAAGAAGTTTCcac                                             <  1:2180078/39‑1 (MQ=255)                                                   ttCACGCACTTGCGGCGCGGTCATGAAAGAAGTTTCcac                                             <  1:633780/39‑1 (MQ=255)                                                    tCACGCACTTGCGGTGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTAAATTCAGGCTGACGGTCAg             >  1:845593/1‑70 (MQ=255)                                                    tCACGCACTTGCGGTGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTAAATTCAGGCTGACGGTCAg             >  1:848079/1‑70 (MQ=255)                                                     cACGCACTTGCGGTGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTAAATTCAGGCTGACGGt                >  1:48957/1‑66 (MQ=255)                                                     cACGCACTTGCGGTGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTAAATTCAGGCTGACGGTCAGCa           >  1:1699898/1‑71 (MQ=255)                                                     cACGCACTTGCGGTGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTAAATTCAGGCTGACGGTCAGCa           >  1:938408/1‑71 (MQ=255)                                                     aaCGCACTTGCGGTGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTAAATTCAGGCTGACGGt                >  1:1887692/2‑66 (MQ=255)                                                        gCACTTGCGGCGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTAAATTCAGGCTGACGGTCAGCACGc        <  1:2011798/71‑1 (MQ=255)                                                        gCACTTGCGGCGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTAAATTCAGGCTGACGGTCAGCACGc        <  1:107449/71‑1 (MQ=255)                                                           cTTGCGGCGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTAAATTCa                          <  1:1041886/50‑1 (MQ=255)                                                            ttGCGGCGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTaa                               <  1:1439085/44‑1 (MQ=255)                                                               cggcgCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTAAATTCAGGCTGACGGTCAGCACGCAGGTCt  <  1:2179843/70‑1 (MQ=255)                                                                  |                                                                   TCCACACCCTTCACTTCCAGCCACAGATGACGCACCAGCGCTTCCATCACTTCACGCACTTGCGGCGCGGTCATGAAAGAAGTTTCCACATCGATCTGAGTAAATTCAGGCTGACGGTCAGCACGCAGGTCT  >  minE/1253651‑1253782

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