Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	986,858	G→A	5.0%	A367T (GCA→ACA)	ynfE →	oxidoreductase subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	986,858	0	G	A	5.0%	268.2 / 3.1	79	A367T (GCA→ACA)	ynfE	oxidoreductase subunit
Reads supporting (aligned to +/- strand):  ref base G (30/45);  new base A (2/2);  total (32/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.23e-01

GGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCG  >  minE/986794‑986927                                                                 |                                                                      ggCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAAc                                                                 >  1:4532775/1‑71 (MQ=255)       tcAAATTCCCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAAc                                                                 <  1:4747949/64‑1 (MQ=255)       gCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCt                                                            >  1:4503400/1‑70 (MQ=255)       gCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCt                                                            >  1:3640360/1‑70 (MQ=255)        cAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTta                                                          >  1:82667/1‑71 (MQ=255)        cAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTta                                                          >  1:3865580/1‑71 (MQ=255)          aaTTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTtata                                                        <  1:794874/71‑1 (MQ=255)            ttACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCg                                                              <  1:2205220/63‑1 (MQ=255)            ttACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCg                                                              <  1:4723767/63‑1 (MQ=255)            ttACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCg                                                              <  1:1826402/63‑1 (MQ=255)              aCCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCt                                                            <  1:118993/63‑1 (MQ=255)              aCCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCt                                                            <  1:3042946/63‑1 (MQ=255)              aCCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCt                                                            <  1:2832605/63‑1 (MQ=255)               ccGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTtat                                                         <  1:719319/65‑1 (MQ=255)               ccGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGc                                                   >  1:3120262/1‑71 (MQ=255)                cGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAACAAAACCCGCTTATATCTGc                                                   >  1:2599372/1‑70 (MQ=255)                  gTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAg                                                <  1:1544393/71‑1 (MQ=255)                  gTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAg                                                <  1:839579/71‑1 (MQ=255)                   tATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAgg                                               <  1:1216201/71‑1 (MQ=255)                       ccGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTg                                           <  1:1594160/71‑1 (MQ=255)                       ccGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTg                                           >  1:1786436/1‑71 (MQ=255)                       ccGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTg                                           >  1:67009/1‑71 (MQ=255)                       ccGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTg                                           <  1:4018452/71‑1 (MQ=255)                       ccGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTg                                           <  1:1122201/71‑1 (MQ=255)                            ggaCCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCt                                                            <  1:2109419/49‑1 (MQ=255)                               ccGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCcac                                   <  1:1361094/71‑1 (MQ=255)                                cGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCcaca                                  <  1:4062813/71‑1 (MQ=255)                                 gTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACaa                                 >  1:3938914/1‑71 (MQ=255)                                   atcatcAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTtat                                                         <  1:3443428/45‑1 (MQ=255)                                   atcatcAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTtat                                                         <  1:3462886/45‑1 (MQ=255)                                   atcatcAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTtat                                                         <  1:3627229/45‑1 (MQ=255)                                   atcatcAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTtat                                                         <  1:1121038/45‑1 (MQ=255)                                   atcatcAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTtat                                                         <  1:1351497/45‑1 (MQ=255)                                   atcatcAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACg                               >  1:1201584/1‑71 (MQ=255)                                   atcatcAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACg                               >  1:1884112/1‑71 (MQ=255)                                    tcatcaAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCcac                                   <  1:4657888/66‑1 (MQ=255)                                    tcatcaAACTGGCGCGTGAAATTGGCACAACAAAACCCGCTTATATCTGc                                                   <  1:3555306/50‑1 (MQ=255)                                    tcatcaAACTGGCGCGTGAAATTGGCACAACAAAACCCGCTTATATCTGc                                                   <  1:220675/50‑1 (MQ=255)                                     catcaAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTgg                                          >  1:941787/1‑58 (MQ=255)                                     catcaAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTgg                                          >  1:2067299/1‑58 (MQ=255)                                     catcaAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACaa                                 >  1:358061/1‑67 (MQ=255)                                     catcaAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACaa                                 >  1:4021346/1‑67 (MQ=255)                                      atcaAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCcac                                   <  1:2523512/64‑1 (MQ=255)                                      atcaAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCcac                                   <  1:3914432/64‑1 (MQ=255)                                      atcaAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCcac                                   <  1:3805545/64‑1 (MQ=255)                                      atcaAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCcac                                   <  1:3729432/64‑1 (MQ=255)                                      atcaAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGc                              >  1:3549555/1‑69 (MQ=255)                                      atcaAACTGGCGCGTGAAATTGGCACAACAAAACCCGCTTATATCTGCCa                                                 >  1:1365393/1‑50 (MQ=255)                                        caAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTTCCAgg                                               <  1:1660594/50‑1 (MQ=255)                                        caAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAgg                                               <  1:789167/50‑1 (MQ=255)                                        caAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAgg                                               <  1:1538556/50‑1 (MQ=255)                                        caAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAgg                                               <  1:655810/50‑1 (MQ=255)                                         aaaCTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACaa                                 >  1:1626486/1‑63 (MQ=255)                                         aaaCTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACaa                                 >  1:2661052/1‑63 (MQ=255)                                         aaaCTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAgg                          >  1:2958060/1‑70 (MQ=255)                                         aaaCTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAgg                          >  1:2923355/1‑70 (MQ=255)                                          aaCTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTgg                                          >  1:2825743/1‑53 (MQ=255)                                          aaCTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTgg                                          >  1:3963246/1‑53 (MQ=255)                                          aaCTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGc                         >  1:241574/1‑70 (MQ=255)                                          aaCTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCa                        <  1:1598855/71‑1 (MQ=255)                                          aaCTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCa                        <  1:2414732/71‑1 (MQ=255)                                            cTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCaaa                      <  1:3851813/71‑1 (MQ=255)                                            cTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCaaa                      >  1:1470352/1‑71 (MQ=255)                                             tGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAAc                     >  1:1968850/1‑71 (MQ=255)                                               gcgcGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACgg                   >  1:99538/1‑71 (MQ=255)                                               gcgcGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACgg                   >  1:4056986/1‑71 (MQ=255)                                                  cgTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGa                <  1:4459972/71‑1 (MQ=255)                                                      aaaTTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAactg            <  1:2651802/71‑1 (MQ=255)                                                      aaaTTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAact             <  1:2724016/70‑1 (MQ=255)                                                      aaaTTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAact             <  1:1520697/70‑1 (MQ=255)                                                      aaaTTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAact             >  1:3073236/1‑70 (MQ=255)                                                       aaTTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAactga           >  1:3032638/1‑71 (MQ=255)                                                       aaTTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAactga           <  1:4329950/71‑1 (MQ=255)                                                       aaTTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAactga           <  1:2150548/71‑1 (MQ=255)                                                         ttGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAactgact         <  1:1088857/71‑1 (MQ=255)                                                           ggCACAGCAAAACCCGCTTATATCTGCCAGGGCTGCGGGCCACAACGCCAGGCAAACGGCGAactgactg        >  1:2684770/1‑70 (MQ=255)                                                              acaGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCcac                                   <  1:3221374/40‑1 (MQ=255)                                                               caGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGc                              <  1:4522497/44‑1 (MQ=255)                                                                aGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCAcgcg  <  1:2066670/71‑1 (MQ=255)                                                                 gCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCcac                                   >  1:1231129/1‑37 (MQ=255)                                                                 gCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCa                            >  1:3333450/1‑44 (MQ=255)                                                                 gCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCa                            >  1:2635201/1‑44 (MQ=255)                                                                 |                                                                      GGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCG  >  minE/986794‑986927

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