Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	2,078,255	C→T	D88N (GAT→AAT)	uxaA ←	altronate hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,078,255	0	C	T	98.1%	164.8 / ‑3.9	52	D88N (GAT→AAT)	uxaA	altronate hydrolase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (28/23);  total (28/24)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

GATCTGCACTTCACGATCTGCCGCTTGCGCAGGCAGATCCTGAAAATCAGGTTGATAGCGATACTGATCCAGATCGCTCAGATTCGTGCGCGTATTGTGGGCGTGAACGTGAACACCTGCCGCAATATCCGCCA  >  minE/2078187‑2078320                                                                     |                                                                  gATCTGCACTTCACGATCTGCCGCTTGCGCAGGCAGATCCTGAAAATCAGGTTGATAGCGATACTGATTCa                                                                 <  1:813701/71‑1 (MQ=255)       cACTTCACGATCTGCCGCTTGCGCAGGCAGATCCTGAAAATCAGGTTGATAGCGATACTGATTCAg                                                                <  1:1509767/66‑1 (MQ=255)       cACTTCACGATCTGCCGCTTGCGCAGGCAGATCCTGAAAATCAGGTTGATAGCGATACTGATTCAGATCg                                                            <  1:1301338/70‑1 (MQ=255)        aCTTCACGATCTGCCGCTTGCGCAGGCAGATCCTGAAAATCAGGTTGATAGCGATACTGATTCAGATCGCt                                                          >  1:69956/1‑71 (MQ=255)        aCTTCACGATCTGCCGCTTGCGCAGGCAGATCCTGAAAATCAGGTTGATAGCGATACTGATTCAGATCGCt                                                          >  1:325990/1‑71 (MQ=255)        aCTTCACGATCTGCCGCTTGCGCAGGCAGATCCTGAAAATCAGGTTGATAGCGATACTGATTCAGATCGCt                                                          >  1:455331/1‑71 (MQ=255)          ttCACGATCTGCCGCTTGCGCAGGCAGATCCTGAAAATCAGGTTGATAGCGATACTGATTCAGATCg                                                            >  1:846169/1‑67 (MQ=255)          ttCACGATCTGCCGCTTGCGCAGGCAGATCCTGAAAATCAGGTTGATAGCGATACTGATTCAGATCGCt                                                          >  1:2116594/1‑69 (MQ=255)               gATCTGCCGCTTGCGCAGGCAGATCCTGAAAATCAGGTTGATAGCGATACTGATTGAGATCg                                                            <  1:316271/62‑1 (MQ=255)                     ccGCTTGCGCAGGCAGATCCTGAAAATCAGGTTGATAGCGATACTGATTCAGATCGCTCAGATTc                                                   <  1:1689883/65‑1 (MQ=255)                     ccGCTTGCGCAGGCAGATCCTGAAAATCAGGTTGATAGCGATACTGATTCAGATCGCTCAGATTc                                                   <  1:123962/65‑1 (MQ=255)                     ccGCTTGCGCAGGCAGATCCTGAAAATCAGGTTGATAGCGATACTGATTCAGATCGCTCAGATTc                                                   <  1:1458919/65‑1 (MQ=255)                     ccGCTTGCGCAGGCAGATCCTGAAAATCAGGTTGATAGCGATAATGATTCAGATCGCTCAGATTCGt                                                 <  1:716447/67‑1 (MQ=255)                         ttGCGCAGGCAGATCCTGAAAATCAGGTTGATAGCGATACTGATTCAGATCGCTCa                                                        >  1:1052197/1‑56 (MQ=255)                         ttGCGCAGGCAGATCCTGAAAATCAGGTTGATAGCGATACTGATTCAGATCGCTCa                                                        >  1:151151/1‑56 (MQ=255)                         ttGCGCAGGCAGATCCTGAAAATCAGGTTGATAGCGATACTGATTCAGATCGCTCAGATTCGTg                                                >  1:478593/1‑64 (MQ=255)                         ttGCGCAGGCAGATCCTGAAAATCAGGTTGATAGCGATACTGATTCAGATCGCTCAGATTCGTg                                                >  1:29289/1‑64 (MQ=255)                         ttGCGCAGGCAGATCCTGAAAATCAGGTTGATAGCGATACTGATTCAGATCGCTCAGATTCGTGCGCGTAt                                         >  1:13879/1‑71 (MQ=255)                          tGCGCAGGCAGATCCTGAAAATCAGGTTGATAGCGATACTGATTCAGATCGCTCAGATTCGTGCGCGTAtt                                        >  1:1323486/1‑71 (MQ=255)                           gcgcAGGCAGATCCTGAAAATCAGGTTGATAGCGATACTGATTCAGATCGCTCAGATTCGTGCGCGTAt                                         <  1:947444/69‑1 (MQ=255)                            cgcAGGCAGATCCTGAAAATCAGGTTGATAGCGATACTGATTCAGATCGCTCAGATTCGTGCGCGTATtg                                       <  1:271252/70‑1 (MQ=255)                            cgcAGGCAGATCCTGAAAATCAGGTTGATAGCGATACTGATTCAGATCGCTCAGATTCGTGCGCGTATtg                                       <  1:2401164/70‑1 (MQ=255)                            cgcAGGCAGATCCTGAAAATCAGGTTGATAGCGATACTGATTCAGATCGCTCAGATTCGTGCGCGTATtg                                       <  1:2057810/70‑1 (MQ=255)                              caggcagATCCTGAAAATCAGGTTGATAGCGATACTGATTCAGATCGCTCAGATTc                                                   <  1:919851/56‑1 (MQ=255)                               aggcagATCCTGAAAATCAGGTTGATAGCGATACTGATTCAGATCGCTCAGATTCGTGCGCGTATTGTggg                                   >  1:333266/1‑71 (MQ=255)                                ggcagATCCTGAAAATCAGGTTGATAGCGATACTGATTCAGATCGCTCAGATTCGTGCGCGTATTGTGGGc                                  <  1:1668732/71‑1 (MQ=255)                                 gcagATCCTGAAAATCAGGTTGATAGCGATACTGATTCAGATCGCTCAGATTCGTGCGCGTATTGTGGGcg                                 >  1:1643533/1‑71 (MQ=255)                                   agATCCTGAAAATCAGGTTGATAGCGATACTGATTCAGATCGCTCAGATTCGTGCGCGTATTGTGGGcgt                                >  1:1003941/1‑70 (MQ=255)                                   agATCCTGAAAATCAGGTTGATAGCGATACTGATTCAGATCGCTCAGATTCGTGCGCGTATTGTGGGcgt                                >  1:1705371/1‑70 (MQ=255)                                    gATCCTGAAAATCAGGTTGATAGCGATACTGATTCAGATCGCTCAGATTCGTGCGCGTATTGTGGGcgtga                              <  1:1228538/71‑1 (MQ=255)                                    gATCCTGAAAATCAGGTTGATAGCGATACTGATTCAGATCGCTCAGATTCGTGCGCGTATTGTGGGcgtga                              <  1:1370396/71‑1 (MQ=255)                                      tCCTGAAAATCAGGTTGATAGCGATACTGATTCAGATCGCTCAGATTCGTGCGCGTATTGTGGGcgtgaac                            >  1:1933482/1‑71 (MQ=255)                                       ccTGAAATTCAGGTTGATAGCGATACTGATTCAGATCGc                                                           >  1:1441840/1‑39 (MQ=37)                                        cTGAAAATCAGGTTGATAGCGATACTGATTCAGATCGCTCAGATTCGTGCGCGTATTGTGGGcgtgaacg                           <  1:834132/70‑1 (MQ=255)                                         tGAAAATCAGGTTGATAGCGATACTGATTCAGATCGCTCAGATTCGTGCGCGTa                                          >  1:1634067/1‑54 (MQ=255)                                         tGAAAATCAGGTTGATAGCGATACTGATTCAGATCGCTCAGATTCGTGCGCGTa                                          >  1:184939/1‑54 (MQ=255)                                         tGAAAATCAGGTTGATAGCGATACTGATTCAGATCGCTCAGATTCGTGCGCGTa                                          >  1:378083/1‑54 (MQ=255)                                         tGAAAATCAGGTTGATAGCGATACTGATTCAGATCGCTCAGATTCGTGCGCGTa                                          >  1:55187/1‑54 (MQ=255)                                          gAAAATCAGGTTGATAGCGATACTGATTCAGATCGCTCAGATTCGTGCGCGTATTGTGGGcgtgaacgtg                         >  1:2114146/1‑70 (MQ=255)                                          gAAAATCAGGTTGATAGCGATACTGATCCAGATCGCTCAGAt                                                     <  1:556978/42‑1 (MQ=255)                                           aaaaTCAGGTTGATAGCGATACTGATTCAGATCGCTCAGATTCGTGCGCGTATTGTGGGcgtgaacgtgaa                       <  1:1301840/71‑1 (MQ=255)                                           aaaaTCAGGTTGATAGCGATACTGATTCAGATCGCTCAGATTCGTGCGCGTATTGTGGGcgtgaacgtgaa                       <  1:2054051/71‑1 (MQ=255)                                           aaaaTCAGGTTGATAGCGATACTGATTCAGATCGCTCAGATTCGTGCGCGTATTGTGGGcgtgaacgtg                         <  1:693568/69‑1 (MQ=255)                                               tCAGGTTGATAGCGATACTGATTCAGATCGCTCAGATTCGTGCGCGTAtt                                        >  1:610313/1‑50 (MQ=255)                                                 aGGTTGATAGCGATACTGATTCAGATCGCTCAGATTCGTGCGCGTATTGTGGGCGTGAACGTGAACACCt                  <  1:475619/70‑1 (MQ=255)                                                    ttGATAGCGATACTGATTCAGATCGCTCAGATTCGTGCGCGTATTGTGGGCGTGAACGTGGACACCTGCCg              >  1:2117373/1‑71 (MQ=255)                                                     tGATAGCGATACTGATTCAGATCGCTCAGATTCGTGCGCGTATTGTGGGCGTGAACGTGAACACCTGc                >  1:1477130/1‑68 (MQ=255)                                                        tAGCGATACTGATTCAGATCGCTCAGATTCGTGCGCGTATTGTGGGCGTGAACGTGAACACCTGCCGCa            >  1:1965823/1‑69 (MQ=255)                                                               aCTGATTCAGATCGCTCAGATTCGTGCGCGTATTGTGGGCGTGAACGTGAACGCCTGCCGCAATATCCGc    >  1:1603956/1‑70 (MQ=255)                                                                cTGATTCAGATCGCTCAGATTCGTGCGCGTATTGTGGGCGTGAACGTGAACACCTGCCGCAATATCCGcc   >  1:171605/1‑70 (MQ=255)                                                                cTGATTCAGATCGCTCAGATTCGTGCGCGTATTGTGGGCGTGAACGTGAACACCTGCCGCAATATCCGCCa  <  1:1641754/71‑1 (MQ=255)                                                                cTGATTCAGATCGCTCAGATTCGTGCGCGTATTGTGGGCGTGAACGTGAACACCTGCCGCAATATCCGCCa  <  1:1436260/71‑1 (MQ=255)                                                                     |                                                                  GATCTGCACTTCACGATCTGCCGCTTGCGCAGGCAGATCCTGAAAATCAGGTTGATAGCGATACTGATCCAGATCGCTCAGATTCGTGCGCGTATTGTGGGCGTGAACGTGAACACCTGCCGCAATATCCGCCA  >  minE/2078187‑2078320

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